Basis for using thioredoxin as an electron donor by Schizosaccharomyces pombe Gpx1 and Tpx1

Glutathione (GSH) peroxidases (GPxs or GSHPx) and thioredoxin (Trx) peroxidases (TPxs) are two classes of peroxidases that catalyze the reduction of peroxides. GPxs and TPxs generally use GSH or Trx, respectively, to recycle the oxidized cysteine (Cys) residue in the protein. However, it is unclear why unlike human GPxs, the Schizosaccharomyces pombe Gpx1 (spGpx1) prefers Trx over GSH for recycling of the active-site peroxidatic Cys residue. Here, we compared spGpx1 and S. pombe Tpx1 (spTpx1) protein sequences with those of their respective homologs in Saccharomyces cerevisiae and humans. Our analysis revealed that like spTpx1, spGpx1 contains a pair of conserved Cys residues (Cys36 and Cys82). These two conserved Cys residues are named peroxidatic and resolving Cys residues, respectively, and are found only in GPxs and TPxs that prefer Trx as an electron donor. Our analysis suggested that Cys36 and Cys82 in spGpx1 are most likely to form a disulfide bond upon oxidation of Cys36. Molecular modelling predicted that a conformational change might be required for the formation of this disulfide bond. Evolutionary analysis suggested that fungal GPxs and TPxs are related by divergent evolution from a common ancestor. Our analyses support a prediction that while spGpx1 and spTpx1 are phylogenetically and functionally different, they evolved from a common ancestor and use a similar mechanism for recycling of the active-site peroxidatic Cys residue. Supplementary Information The online version contains supplementary material available at 10.1186/s13568-022-01381-2.

Introduction maintenance of physiological functions (Schieber andChandel 2014, Sies andJones 2020). However, the toxic levels of ROS cause damage to lipids, proteins and DNA, and are implicated in aging and numerous diseases, including neurodegenerative diseases and cancer.
To deal with ROS toxicity, cells have evolved mechanisms such as production of enzymatic and nonenzymatic antioxidant defenses (Do et al. 2019). Antioxidant enzymes include superoxide dismutase, which catalyzes the dismutation of superoxide anions into oxygen and H 2 O 2 , catalase that degrades the conversion of H 2 O 2 into water and oxygen, and thiol-dependent peroxidase that uses redox-active cysteines (Cys) or selenocysteine (Sec) to reduce H 2 O 2 or organic hydroperoxides (ROOR) to H 2 O or alcohols (Nordberg andArnér 2001, Tairum et al. 2012). Nonenzymatic antioxidant defenses include glutathione (GSH), ascorbic acid, glutathione-S-transferase (GST) and metallothionein.
Saccharomyces cerevisiae has 3 GPxs [scGpx1, scGpx2 and scGpx3 (also known as Hyr1 or scOrp1)]. scGpx2 and scGpx3 prefer Trx over GHS as an electron donor. Although scGpx1 uses GSH and Trx almost equally as the electron donor in vitro, Trx plays a crucial role in the reduction of spGpx1 in vivo (Ohdate et al. 2010, Ukai et al. 2011. The three proteins are localized in the peroxisome, mitochondrial, and cytosol, respectively (Ohdate and Inoue 2012). Unlike scGpx1 and scGpx2, scGpx3 also functions as the major H 2 O 2 sensor in the cytosol. After oxidized by H 2 O 2 , the Cys36 of scGpx3 forms an intermolecular disulfide bond with the Cys598 of scYap1, the central regulator of oxidative stress response in S. cerevisiae. This disulfide bond is rearranged to form an intramolecular disulfide bond in scYap1, leading to the activation of scYap1 (Delaunay et al. 2002). In addition, an isoform of scGpx3 is localized in mitochondria and facilitates oxidative protein folding (Kritsiligkou et al. 2017). In contrast to most other organisms that contain several GPxs, Schizosaccharomyces pombe has only a single GPx (spGpx1). Like S. cerevisiae GPxs, spGpx1 prefers Trx over GSH for electron donor (Lee et al. 2008, Paulo et al. 2014. Consistent with the view that the human genome shows a tendency to expand the number of ROS-related proteins compared to fungi, humans have 8 GPxs including 5 selenoproteins (hsGPx1-hsGPx4, and hsGPx6) and 3 non-selenoproteins (Matoušková et al. 2018).
The reduction of peroxide by Se-GPxs involves the formation of several intermediate modifications to the Sec active site (Lubos et al. 2011). After reacting with peroxide, the selenol group (Se-H) at the active site is oxidized to selenenic acid (Se-OH), and one molecule of GSH then reduces Se-OH and forms glutathiolated selenol (Se-GS). Subsequently, another molecule of GSH reduces Se-GS, and simultaneously generates the oxidized form of GSH disulfide (GSSG), which is reduced to regenerate GSH by GSH reductase using NADPH as an electron donor. The active sites of non-Se-GPxs contain a Cys residue instead of a Sec residue. In contrast, the reduction of peroxide by non-Se-GPxs involves the formation of an intramolecular disulfide bond between two Cys residues. The first redox-active Cys residue, designated the peroxidatic Cys is oxidized by peroxide, yielding sulfenic acid (Cys-SOH). This oxidized Cys residue then forms a disulfide bond with another Cys residue, designated the resolving Cys. This disulfide bond is reduced by the Trx system, which is composed of NADPH, Trx reductase (TrxR), and Trx.
Peroxiredoxins (Prxs or Prdxs) are a family of thiol peroxidases that catalyze the reduction of H 2 O 2 , organic hydroperoxides, and peroxynitrite (Wood et al. 2003a, b, Wood et al. 2003a, b, Perkins et al. 2015, Rhee 2016, Bolduc et al. 2021. Based on the number of active cysteine residues, Prxs can be categorized into 1-Cys and 2-Cys Prxs. All 2-Cys Prxs share the similar catalytic mechanism. First, a redox-active Cys residue (called the peroxidatic Cys) is oxidized to a sulfenic acid (Cys-SOH) by either H 2 O 2 or organic hydroperoxides. Second, the sulfenic acid forms a disulfide bond with another redoxactive Cys (called the resolving Cys) of the other subunit of the Prx homodimer (typical 2-Cys Prxs), or with the C-terminal resolving Cys of the same monomer (atypical 2-Cys Prxs). Last, this resulting disulfide bond is normally reduced by thioredoxin. Thus 2-Cys Prxs are also named Trx peroxidases (TPxs). For 1-Cys Prxs, the sulfenic acid is directly reduced to thiol due to no nearby Cys to form a disulfide bond.
Prxs are found in both eukaryotic and prokaryotic organisms. S. cerevisiae has five Prxs (scTsa1, scTsa2, scAhp1, scPrx1, and scDot5) (Kim et al. 2010). Mammalian cells have 6 Prxs: 4 typical 2-Cys Prx (huPrx1-huPrx4), 1 atypical 2-Cys Prx (huPrx5) and 1 1-Cys Prx (huPrx6) (Rhee and Kil 2017). S. pombe contains 3 Prxs (spTpx1, spPmp20, and spBcp1). spTpx1 is a 2-Cys Prx that functions as the main H 2 O 2 scavenger during aerobic growth in S. pombe (Jara et al. 2007, Paulo et al. 2014). spTpx1 also functions as the H 2 O 2 sensor that transduces the redox signal to the transcription factor spPap1, which a homolog of scYap1 (Paulo et al. 2014). Upon exposure to lower levels of H 2 O 2 , Cyc48 in spTpx1 is oxidized to sulfenic acid, which can form a disulfide bond with Cyc. The reduction of this labile disulfide bond is mediated by the Trx system (Jara et al. 2007, Lu andHolmgren 2014). The sulfenylated spTpx1 can activate spPap1 by inducing the formation of an intramolecular disulfide bond in spPap1, which then translocates from cytoplasm to nucleus, where it activates genes required for adaptation to oxidative stress (Quinn et al. 2002, Calvo et al. 2013. However, at high levels of oxidants, the sulfenylated spTpx1 can be further oxidized to the sulfinic acid (SO 2 H) or sulfonic acid (SO 3 H) forms, resulting in loss of its peroxidase activity and prevention of Pap1 activation (Bozonet et al. 2005, Vivancos et al. 2005. The hyperoxidized spTpx1 can be reduced to restore the peroxidase activity by sulfiredoxin spSrx1. The functions of spPmp20 and spBcp1 are unclear.
In this study, we conducted a bioinformatics analysis to understand the basis for why both spGpx1 and spTpx1 use Trx as an electron donor. Our analyses revealed that spGpx1 harbors two conserved Cys residues that are found in other GPxs and TPxs that use Trx as the electron donor. Our analysis suggests that like other Trxdependent GPxs and TPxs, spGpx1 may also form a disulfide bond between Cys36 and Cys82 upon oxidation. Our phylogenetic study suggests that fungal GPxs and TPxs descended from a common ancestor.

Proteins data set and sequence analyses
Full-length sequences of the spGpx1 and spTpx1 were obtained from PomBase. We used the spGpx1 and spTpx1 protein sequences as queries in BLAST (Basic Local Alignment Search Tool) searches against human genome databases and fungal genome databases, including the National Center for Biotechnology Information (NCBI), using default parameters, to identify candidate GPxs and TPxs. All protein sequences were retrieved by using the cut-off E-value of 0.05. The UniProt Knowledgebase (UniProtKB) was used to check the presence of the conserved active site in the candidate proteins. The candidate proteins were then examined for the existence of the conserved domains using the Simple Modular Architecture Research Tool (SMART) (Letunic and Bork 2018). Multiple EM for Motif Elicitation (MEME Suite) Version 5.4.1 was used to identified the conserved sequence motifs using default parameters (Peng et al. 2018). The web-based prediction program WoLF PSORT was used to predict protein subcellular localization (Horton et al. 2007).

Homology prediction of the three-dimensional structure of spGpx1 by I-TASSER
The structure of spGpx1 was predicted by the I-TASSER server (IterativeThreadingAssemblyRefinement-https:// zhang lab. ccmb. med. umich. edu/I-TASSER/). The sequence was uploaded in FASTA format and output result was downloaded in PBD format. C-score was measured for the spGpx1 model and the lower rank model, which has a high C-score was selected as a better quality structure and C-score is confidence score for the model (Zhang 2008). UCSF Chimera was used for image visualization and structural analysis.

Phylogenetic analysis
To construct a phylogenetic tree, a multiple amino acid sequence alignment was generated with the MUSCLE algorithm included in the Molecular Evolutionary Genetics Analysis (MEGA X) software (version 10.1.8) (Kumar et al. 2018). The phylogenetic tree was constructed using the maximum-likelihood method and the Jones-Taylor-Thornton (JTT) model in MEGA X, with 1000 bootstrap replications.

Bayesian evolution analysis
To estimate the evolutionary rates of candidate fungal GPxs and TPxs, the Bayesian evolutionary tree was constructed using Beast v1.10.4 (Zhang and Drummond 2020). BEAST XML files were generated for GPxs and TPxs by using BEAUti software which is available in the BEAST package. The substitution model was based on BLOSUM62, and the 4-category gamma model was used for the site-heterogeneity model. A strict clock and constant-size coalescent tree prior were employed. The Markov Chain Monte Carlo (MCMC) sampling algorithm was used to estimate statistical confidence. BEAST analysis was run, and the tree was viewed by using the software FigTree v.1.4,4 (http:// tree. bio. ed. ac. uk/ softw are/ figtr ee/). Tracer v1.7.2 was used to calculate substitution per amino acid per year was calculated as well (Rambaut et al. 2018).

spGpx1 harbours two conserved Cys residues that are likely to form an intramolecular disulfide bond
We conducted BLAST homology analyses using the spGpx1 protein sequence against the human genome database and fungal genome databases as part of our attempts to investigate the structural and functional diversity of GPxs. For further analysis, we chose 12 GPxs: 1 from S. pombe (spGpx1), 3 from S. cerevisiae (scGpx1, scGpx2 and scGpx3) and 8 from humans (hsGPx1-hsGPx8). spGpx1 is more similar to S. cerevisiae GPxs than to human GPxs (Table 1). It shares the highest similarity with scGpx3 (71% identity). To examine the evolutionary relationships among these GPxs, the phylogenetic tree of GPxs was constructed based on the similarity, using the maximum-likelihood method (Additional file 1: Fig. S1). Our analysis suggested that spGpx1 is most closely related to scGpx3.
To understand the structural and functional diversity of GPxs, we performed multiple sequence alignment of S. pombe, S. cervisiae and human GPxs. Multiple sequence alignment revealed the existence of a conserved Cys/ Sec residue in all selected GPxs (Cys in spGpx1, scGpx1, scGpx2, scGpx3, hsGPx5, hsGPx7 and hsGPx8, and Sec in hsGPx1-hsGPx4 and hsGPx6). This Cys/Sec residue in GPxs has been demonstrated to be oxidized by peroxide (Delaunay et al. 2002, Tanaka et al. 2005, Ohdate et al. 2010). Our analysis also revealed an exclusively conserved Cys residue (Cys82 in spGpx1) in S. pombe and S. cerevisiae GPxs (Fig. 1). It has been shown an intramolecular disulfide bond formed by these two Cys residues is required for the peroxidase activity of the S. cerevisiae GPxs (Delaunay et al. 2002, Tanaka et al. 2005, Ohdate et al. 2010. Using the SMART database, we found that S. pombe, S. cerevisiae and human GPxs all contain the GSHPx domain that is involved in catalyzing the reduction of peroxide (Fig. 1). Using the MEME Suite, we identified 3 highly conserved motifs in these GPxs (Fig. 2). These three motifs are important for catalysis as they contain residues that participate in catalysis (Fig. 2).

Homology modelling of spGpx1
We carried out homology modelling of spGpx1 to determine whether the two conserved cysteines (Cys36 and Cys82) are in close proximity. Because spGpx1 shares the highest sequence similarity with scGpx3, whose structure has been determined (Zhang et al. 2008), we selected the structure of scGpx3 to build a three-dimensional model of spGpx1 using the online tool I-TASSER. The predicted spGpx1 structure is shown in Fig. 3. In the structure of spGpx1, Cys36 and Cys82 are in a distant of 15.89 Ǻ and is thus not in the close proximity to form a disulfide bond. Similarly, Cys36 and Cys82 are in a distant of 13.24 Ǻ in the scGpx3 structure. This analysis suggests that a conformational change is required to form a disulfide bond between these two Cys residues upon oxidation.

spTpx1 contains N-terminal AhpC/TSA domain and C-terminal 1-cysPrx C domain.
A BLAST search using the spTpx1 protein sequence as a query revealed that spTpx1 homologs are present in S. cerevisiae and humans. Three spTpx1 homologs (scTsa1, scTsa2, and scPrx1) were found in S. cerevisiae, and five (hsPrdx1-hsPrdx4 and hsPrdx6) were found in humans. spTpx1 shares the highest sequence similarity (66% identity) with hsPrdx2 (Table 2). Phylogenetic analysis by maximum likelihood revealed that hsPrxd2 and scTsa1 were the most similar proteins to spTpx1 (Additional file 1: Fig. S2). SMART database searches showed that spTpx1 and representative Tpx1 homologs have an N-terminal AhpC/TSA domain and the C-terminal 1-cys Prx C domain, which are characteristic of Prxs (Fig. 4).
Multiple sequence alignment of spTpx1 and its homologs in S. cerevisiae and humans (Fig. 4) revealed that except for scPrx1 and hsPrdx6, which has oneconserved Cys residue, all Prxs examined contain two characteristic Cys residues; an absolutely conserved peroxidatic Cys (Cys48 in scTsa1) and one highly conserved resolving Cys (Cys146 in scTsa1). Notably, the second Cys residues from spTpx1 and hsPrdx1-hsPrdx4 are not aligned with those in scTsa1 and scTsa2. The active sites of Prxs contain a catalytic triad composed of a peroxidatic Cys, a Thr (or a Ser) and an Arg (Tairum et al. 2016). This catalytic triad (Thr45-Cys48-Arg124 in ScTsa1) are found in spTpx1 and all other Prxs examined (Fig. 4). Glu51 and Arg147 in scTsa1 that are important for the scTsa1-Trx interaction are absolutely conserved in all Prxs examined (Tairum, de Oliveira et al. 2012).
Like scTsa1 and scTsa2, spTpx1 is a typical 2-Cys Prx that uses Trx to reduce the disulfide bond in Prxs. The reductant responsible for the reduction of the disulfide bond in 1-Cys Prx in vivo has not been definitively resolved. Candidate reductants include ascorbate, Trx, Trx reductase, glutaredoxin and GSH (Monteiro et al. 2007, Shadel Gerald and Horvath Tamas, 2015, Pedrajas et al. 2016. spGpx1 and spTpx1 are probably related via divergent evolution from a common ancestor A Bayesian evolutionary analysis was performed using a data set composed of 28 GPxs and 28 TPxs belonging to taxonomically diverse fungi to determine the evolutionary relationships among fungal GPxs and TPxs (Additional file 1: Table S1 and S2). The evolutionary tree suggested that both fungal GPxs and TPxs share a common ancestor (Fig. 5)  The motifs were detected by the MEME Suite (Version 5.4.1). Motifs 1-3 are labelled in green, red and blue, respectively. B Sequence logos for Motifs 1-3. The sequence logos of the motifs derived from S. pombe, S. cerevisiae and human GPxs were generated using the MEME Suite. The height of each amino acid indicates the level of conservation at that position. Amino acid residues are highlighted according to the biochemical properties of the amino acid residues. The catalytic residues Cys, Gln, Trp and Asn (Tosatto et al. 2008) are indicated by stars substitution/aa/year, respectively. Next, 28 GPXs were used to estimate the mutation rate, the evolutionary rate for the GPxs with the strict clock model calculated was 35.9718 substitutions/aa/year, whereas with the relaxed clock model, it was 36.9963 substitutions/aa/year. Furthermore, the mean evolutionary rate for the 28 TPxs was 12.6302 substitutions/aa/year with the strict clock model, and was 12.5194 substitutions/aa/year with the relaxed clock model. Our results showed that the TPxs have a relatively low substitution rate compared to GPxs (35.9718 substitutions/aa/year). However, the mean substitution rate for GPxs was relatively higher (35.9718 substitutions/aa/year) compared with GPXs/TPXs (28.73076 substitutions/aa/year) (Table 3).

Discussion
The GSH and Trx systems are two major thiol-dependent antioxidant systems, which are also involved in regulating the redox state of GPxs and TPxs, respectively. Unlike human GPxs, which use GSH as an electron donor, many fungal GPxs prefer Trx over GSH as an electron donor in vivo. In this study, we carried out a bioinformatics analysis to understand the basis for the selection of Trx by both spGpx1 and spTpx1. Sequence analysis revealed that unlike human GPxs, which contain a single conserved Cys/ Sec residue (called peroxidatic Cys/ Sec), S. pombe and S. cerevisiae GPxs all contain two conserved Cys residues called the peroxidatic and resolving Cys residues, respectively. It has been shown that all three S. cerevisiae GPxs function as atypical 2-Cys Prxs which form an intramolecular disulfide bond between the two conserved Cys residues and use Trx as an electron donor for recycling of the oxidized peroxidatic Cys residue in vivo (Delaunay et al. 2002, Tanaka et al. 2005, Ohdate et al. 2010). Our analysis predicts that spGpx1 is likely to function as a 2-Cys Prx using Trx as an electron donor to reduce the disulfide bond between Cys36 and Cys82 residues. Molecular modelling of spGpx1 suggests that a conformational change upon Cys36 oxidation may be required to bring these two Cys residues together. This is consistent with the idea that upon oxidation of the peroxidatic Cys residue, a conformational change in the active sites of Prxs is required for the disulfide bond formation between the peroxidatic and resolving Cys residues (Hall et al. 2011, Zhou et al. 2016). However, we cannot predict whether Cys36 and Cys82 residues in spGpx1 form disulfide bond intramolecularly or intermolecularly. spGpx1 is phylogenetically and functionally distinct from spTpx1. No sequence similarity between spGpx1 and spTpx1 was detected. While both proteins are H 2 O 2 scavengers, they play different roles in cells. spGpx1 is localized in both the cytosol and mitochondria, and is the predominant peroxidase in stationary phase (Lee et al. 2008). In contrast, spTpx1 is a ROS sensor that transduces the redox signal to the transcription factor Pap1 (Calvo et al. 2013). Unlike spGpx1, spTpx1 is required for the normal aerobic growth of S. pombe (Jara et al. 2007). Despite these differences between spGpx1 and spTpx1, our phylogenetic analysis suggests that the two proteins may share a common ancestor and use a similar mechanism for the regeneration of the active-site peroxidatic Cys residue. associated taxa clustered together is shown next to the branches. Initial tree(s) for the heuristic search were obtained automatically by applying Neighbor-Join and BioNJ algorithms to a matrix of pairwise distances estimated using the JTT model, and then selecting the topology with superior log likelihood value. There was a total of 305 positions in the final dataset. Table S1. Comparison of spGpx1 with selected candidate fungal GPxs used in the evolutionary analysis. Table S2. Comparison of spTpx1 with selected candidiate fungal TPxs used in the evolutionary analysis