IMR Press / FBL / Volume 24 / Issue 6 / DOI: 10.2741/4770
Interactions between human hemoglobin subunits and peroxiredoxin 2
Show Less
1 Laboratory of Hemoglobin, Baotou Medical College, 014060, Baotou, China
2 Key Laboratory of Carcinogenesis and Translational Research (Ministry of Education/Beijing), Central Laboratory, Peking University Cancer Hospital & Institute, 100142, Beijing, China
3 State Key Laboratory of Heavy Oil Processing and Department of Materials Science and Engineering, China University of Petroleum, 102249, Beijing, China
*Correspondence: (Yan Su)
Front. Biosci. (Landmark Ed) 2019, 24(6), 1085–1096;
Published: 1 March 2019
(This article belongs to the Special Issue Leader sequences of coronavirus are altered during infection)

Red blood cells (RBCs) are exposed to exogenous reactive oxygen species in the circulatory system. To this end, the interactions between the different hemoglobin (Hb) subunits and peroxiredoxin 2, which is a ubiquitous member of the antioxidant enzymes that also controls the cytokine-induced peroxide levels, were assessed. We predicted by the increment of diversity with quadratic discriminant analysis (IDQD) that peroxiredoxin2 (Prx2) could interact with the hemoglobin alpha, beta and gamma subunits but not with the delta subunit. Coimmunoprecipitation (co-IP), electrospray ionization quadrupole time of flight (ESI-Q-TOF) mass spectrometry, Western blotting and X-ray absorption fine structure (XAFS) spectroscopy were performed to verify these predictions. The results showed that Prx2 was a member of the beta-globin immunoprecipitating complex that existed in hemoglobin A, hemolysate-hemoglobin A, hemoglobin A-hemoglobin A2, hemolysate-hemoglobin A-hemoglobin A2 and hemoglobin A2 but not in hemolysate-hemoglobin A2. Adding Prx2 to hemoglobin A altered the second shell of iron embedded in hemoglobin A. Therefore, Prx2 interacts with hemoglobin A (Alpha2Beta2) and hemoglobin F (Alpha2Gamma2) but not with hemoglobin A2 (Alpha2Delta2).

Peroxiredoxin 2
Increment Of Diversity With Quadratic Discriminant Analysis
X-Ray Absorption Fine Structure

In living cells, biological functions are generally not carried out by individual proteins but by networks of intricate interactions between numerous proteins. Therefore, a deeper knowledge of the functional organization of such networks has crucial significance in understanding the biological activities of proteins (1-3). At present, many experimental techniques have been developed to detect protein-protein interactions (PPIs), such as the yeast two-hybrid-based method, coimmunoprecipitation (co-IP) and far-Western blot analysis (4-6). As most of these experiments are time-consuming and expensive, reliable computational methods need to be developed to predict the PPIs before experimentation. The prediction can improve the experimental efficiency and reduce the upfront costs (7). Several algorithms, such as neural networks (NNs) (8), support vector machines (SVMs) (9) and random forests (RFs) (10), have been reported. Park et al. attempted to identify the interacting partners of a protein by viewing the interactions of evolutionarily related domains that belong to the same structural families of the investigated protein (11). Sprinzak and Margalit put forward another indirect interaction prediction method based on digging out the signature feature related to interactions rather than the domain interaction information from the protein sequences via protein classification (12). The algorithm of the increment of diversity with quadratic discriminant analysis (IDQD), focusing on the recognition of sequence patterns, was proposed by Zhang and Luo in 2003 (13). This algorithm was successfully applied in the recognition of the transcriptional start site (14), intron splice site (13), DNase I hypersensitive site (15), protein classification (16, 17), nucleosome positioning prediction (18-20) and recombination hotspots (21). However, successful cases of protein interaction prediction by IDQD methods have not yet been reported.

Hemoglobin is the ideal protein model for structure and function research. Each hemoglobin molecule contains one pair of alpha globin family chains (including alpha and zeta globin chains) and one pair of beta-globin family chains (including beta, gamma, delta and epsilon globin chains). There are three types of Hb in normal adult RBCs, which are HbA (alpha2beta2, 96%), HbA2 (Alpha2Delta2, 2-3%) and HbF (Alpha2Gamma2, <1%) (22). RBCs are exposed to more oxidative stress than any other cell type due to the abundance of heme iron and oxygen (23). The autoxidation of Hb produces a large number of reactive oxygen species (ROS), which can damage the lipids, proteins and cytoskeleton of the membrane and can decrease the deformability of RBCs if the ROS cannot be cleared in time (24, 25). As a nonnucleated cell, a RBC cannot synthesize proteins by itself; therefore, a series of antioxidant enzymes, such as superoxide dismutase (SOD), catalase, glutathione peroxidase (GSH-Px), and peroxiredoxins (Prxs), have evolved to protect RBCs from oxidative damage and maintain their physiological functions during erythropoiesis (26). Peroxiredoxins (Prxs) are a newly discovered antioxidant enzyme family. Prx2 is the main type of Prx and plays an important role in preventing the damage mediated by ROS in RBCs (27-29). Our previous study (30) and a recent report (31) suggested that Prx2 might interact with Hb; however, the mechanism of the interaction between Prx2 and Hb has not been thoroughly elucidated. In this study, the algorithm of IDQD was applied to predict the interaction between different human Hb subunits. Prx2, co-IP and Western blotting were subsequently used to confirm the prediction results, and X-ray absorption fine structure (XAFS) spectroscopy was used to study the interaction mechanism.

3.1. IDQD prediction method

Protein sequences were obtained from the following sources: Hemoglobin subunit alpha (Homo sapiens) (NCBI Reference Sequence: NP_000549.1), hemoglobin subunit beta (Homo sapiens) (NCBI Reference Sequence: NP_000509.1), hemoglobin subunit gamma (Homo sapiens) (Swiss-Prot: P68068.2), hemoglobin subunit delta (Homo sapiens) (GenBank: ABK79072.1) and Peroxiredoxin2 (Homo sapiens) (GenBank: AAH00452.1). The original PPI data of humans used in this work was downloaded from the Database of Interacting Proteins (DIP) ( (32). The corresponding protein sequence information was extracted from the Universal Protein Resource (UniProt) (33). There were 1,962 protein pairs (involving 1,703 proteins) in both the positive and negative sets. The prediction of PPI was based on protein sequences. A descriptor using a conjoint quintuple (5-mer), which regarded any five continuous amino acids as a unit, was first proposed by us (34). The classification of amino acids and IDQD analysis was introduced in our previous work. The sensitivity (Sn), specificity (Sp), false positive rate (FPR), total accuracy (TA), correlation coefficient (CC), receiver operator characteristics (ROC) and the area under ROC curve (auROC) were used to access the accuracy of the PPI predictions. A three-fold cross-check is combined with the above four indicators, as well as the ROC curve and auROC, to evaluate the predictive power of the IDQD algorithm. The specific method is as follows: The data set is randomly divided into 3 parts, each of which contains an equal amount of positive and negative set data. Two of them were used as training sets, and the remaining one was used as a test set. The rotation was carried out three times in turn, and the average of the various indicators was used to evaluate the model performance.

3.2. Preparation of RBC suspension and hemolysate

All the experimental protocols of this study were approved by the Ethics Committee of Baotou Medical College. Fresh, normal anti-coagulated blood was collected from volunteers. Before the blood samples were collected, written consent was obtained from all of the volunteers. Then, the RBC suspensions and hemolysates were prepared in accordance with the approved guidelines (35).

3.3. Purification of different Hb components from starch-agarose mixed gel

A two percent starch-agarose mixed gel (starch: agarose=4:1) was prepared with 1× TEB buffer (pH 8.6) (36). Eight microliters of RBC suspension and hemolysate were loaded into the gel, and electrophoresis was performed at 5 V/cm for approximately 2 hours. After the electrophoresis, the red bands of RBC-HbA (RA), hemolysate-HbA (HA), RBC-HbA-A2 (RA-A2), hemolysate-HbA-A2 (HA-A2), RBC-HbA2 (RA2) and hemolysate-HbA2 (HA2) were cut out separately and frozen at -80°C for at least 30 min. Before use, the gel was taken out and thawed at room temperature. After centrifuging at 10,000 g for 10 min, the supernatants were pipetted into new Eppendorf tubes and concentrated with Sephadex G25.

3.4. Coimmunoprecipitation (co-IP)

Seven microliters of PMSF and 20 μL of agarose conjugate hemoglobin Beta (37-8) (Santa Cruz, sc-21757AC) were added into 700 μL of HA and RA samples, respectively. A negative control was prepared by adding 7 μL of PMSF and 20 μL of agarose conjugated mouse IgG into 700 μL HA and RA solutions, respectively. After incubating the samples overnight at 4°C on a flat platform, they were centrifuged at 10,000 g for 1 min at 4°C. Next, the beads were collected, and the supernatants were removed carefully. The beads were washed with 700 μL 1×PBS (pH 7.4) buffer six times and 0.1×PBS (pH 7.4) buffer one time. Forty microliters of 1× loading buffer was added into each sample, which were then heated at 100°C for 5 min. After centrifuging at 10,000 g for 5 min, the supernatant was ready for SDS-PAGE.


The samples of supernatant (20 μL) prepared by the co-IP complex procedure described above were separated by 5-12% SDS-PAGE (36). After electrophoresis, the gel was stained overnight with 0.08% Coomassie Brilliant Blue G250 and was destained in Milli-Q water until the background staining of the gel was low.

3.6. In-gel digestion

The desired gel bands were excised into 1 mm3 pieces and further destained and digested according to a routine protocol (36). The digested peptide samples were extracted with 50-100 mL of 5% trifluoroacetic acid (TFA) solution at 40°C for 1 hour and later extracted with the same volume of 50% CH3CN/2.5% TFA solution at 30°C for 1 hour. Finally, the samples were extracted ultrasonically with 50 mL of CH3CN solution. The extracted solutions were pooled and dried in a SpeedVac vacuum dryer and resuspended in 3-5 mL of 0.1% formic acid for MS analysis.

3.7. Identification of proteins by liquid chromatography tandem mass spectrometry

The peptide mixture was analyzed on a nanoliquid chromatography tandem mass spectrometry (LC-MS/MS) system, which consists of an Ultimate HPLC system (Dionex) and a quadrupole time of flight (Q-TOF) mass spectrometer (micrOTOF-Q II mass spectrometer, Bruker) equipped with a nano-ESI source. The MS/MS data were processed using Data Analysis 4.0 and then searched against the Swissprot protein sequence databases by MASCOT ( The TOF mass analyzer of the instrument was calibrated by using 10 mM Na TFA (pos). For the MS/MS ion search, a peptide charge state of +2, +3 or +4 and peptide/fragment mass tolerance of ±0.1 Da was used. Probability-based MASCOT scores were estimated by comparing the search results against the estimated random match population and reported as −10 × log (P), where P is the absolute probability. The significance threshold was set at P < 0.05.

3.8. Western blotting

Twenty microliters of HA, RA, HA2 and RA2 were first separated by 5-12% SDS-PAGE, and then the protein was transferred to the PVDF membrane. The PVDF membrane was subsequently blocked in 5% skim milk for 1 hour at room temperature and then incubated with rabbit polyclonal anti-Prx2 antibody (Abcam, ab86295) at a 1:1000 dilution overnight at 4°C. Next, the membranes were hybridized with anti-rabbit IgG-conjugated horseradish peroxidase (Santa Cruz, sc-2004) at a 1:2000 dilution for 1 hour at room temperature. Finally, the protein signals were detected by the Pierce ECL Western Blotting Substrate and visualized by a Tanon 4200 (Tanon Science & Technology Co., Ltd.).

3.9. XAFS experiment

HA samples were collected from the starch-agarose gel by the freeze-thaw method after electrophoresis of the hemolysate (30). Peroxiredoxin 2 was purchased from Abcam (Abcam, ab8). XAFS experiments were carried out at the National Synchrotron Radiation Laboratory (USRL), University of Science and Technology of China. The XAFS data of the iron K-edge was collected at the U-7C end station in transmission mode. The electron energy in the storage ring was 0.8 GeV, and the ring current was approximately 160 mA. Higher harmonics were reduced sufficiently by detuning the double-crystal Si(III) monochromator at the beam line. The raw data were analyzed by using the NSRLXAFS 3.0 program. The basic EXAFS data, χ (k), was fitted to the following equation:

χ j ( k ) = S 0 2 N j k R j 2 F j ( π , k ) exp ( 2 k 2 Δ σ j 2 ) exp ( 2 R j / λ j ) sin ( 2 k R j + ϕ j + Σ j ) (1)


Σ j = 4 Δ σ j 2 R j K 4 3 σ j 3 k 3 (2)

where k is the photoelectron wavenumber, S0 is the amplitude reduction factor, Nj is the coordination number of the jth shell, Rj is the distance between the coordination atoms and central atoms in the shell, Fj (π,k), λj and фj are the jth atom backscattering factors, Δσj2 is the Debye-Waller factor, and σj3 is the third cumulant (37).

4.1. Discrimination of the interactions between human hemoglobin subunits and Prx2 based on the IDQD model

In our approach, the PPI prediction model was constructed based on the training set. Three training sets and three testing sets were prepared by the sampling method described previously (34). Each training set consisted of 2,616 protein pairs, with half of the protein pairs randomly selected from the data of positive PPI pairs and the other half randomly selected from the negative protein pairs. Each testing set was constructed with another 1,308 protein pairs. Thus, three prediction models were generated for the three training sets of data. The best prediction results are listed in Tab 1, and the ROC curves are shown in Figure 1. The performance of the IDQD modules for the prediction of the PPIs in human samples was trained on different types of tests. The values showed that the performance was at the best threshold (ξ0). The ξ0 of the self-consistency test and the 3-fold cross-validation are –0.45 and –0.08, respectively. The PPI prediction results between each two Hb subunits (Alpha, Beta, Gamma, Delta) and between each Hb subunit and Prx2 (P) were analyzed with an IDQD model. As shown in Figure 2, there were interactions between Beta and Beta, Gamma and Gamma, Beta and Alpha, Alpha and Alpha, P and Beta, Gamma and Alpha, P and Gamma, P and P, P and Alpha, Gamma and Beta, and Delta and Alpha subunits; however, there were no interactions between Beta and Delta, Delta and Delta, Delta and Gamma, or Delta and P subunits.

Table 1.IDQD prediction results of protein-protein interactions
Tests Performance Measures
Sn (percent) Sp (percent) TA (percent) CC
Self-consistency test 76.04 75.74 75.89 0.52
3-fold cross-validation 64.22 64.68 64.45 0.29
Figure 1.

Evaluation of the performance of IDQD models. ROC curves for IDQD model were plotted for the self-consistency test and the 3-fold cross-validation test. A mean auROC of 0.82 was obtained in the self-consistency test, while a mean auROC of 0.67 was obtained in the 3-fold cross-validation experiment.

Figure 2.

IDQD prediction results of PPI between Hb subunits (Alpha, Beta, Gamma, Delta) and Prx2 subunit. The scores were calculated by deducting -0.45, which is the threshold value (ξ0), from the performance values between every two subunits. A positive score represents a positive interaction, while a negative score represents a negative interaction.

4.2. Comparison of IDQD-predicted interactions between different Hb subunits with the experimental results

Hb subunits have been reported to interact with each other (38-47). To assess the predictive accuracy of the IDQD model, the predicted results of the Hb subunits were compared with the reported experimental results (Tab 2). The total predictive accuracy is 80% (the interactions of BetaBeta, GammaGamma, AlphaBeta, AlphaAlpha, AlphaGamma, AlphaDelta, BetaDelta and DeltaGamma are correct, but those of BetaGamma and DeltaDelta are incorrect), and this result demonstrates that the accuracy of the algorithm for the PPI prediction is good. The interactions between the Hb subunits (Alpha, Beta, Gamma, Delta) and Prx2 were then analyzed by this algorithm. The results showed that Prx2 had interactions with the Alpha, Beta and Gamma subunits of Hb but had no interaction with the Delta subunit (Figure 2). These results indicate that Prx2 should interact primarily with HbA, which consists of Alpha and Beta subunits.

Table 2.Comparison of IDQD-predicted interactions between different Hb subunits with the reported experimental results
Classification Predicted interactions Proven interactions (38-47)
Positive interactions BetaBeta, GammaGamma, AlphaBeta, AlphaAlpha, AlphaGamma, BetaGamma, AlphaDelta BetaBeta, GammaGamma, AlphaBeta, AlphaAlpha, AlphaGamma, AlphaDelta, DeltaDelta
Negative interactions BetaDelta, DeltaDelta, DeltaGamma BetaDelta, BetaGamma, DeltaGamma
4.3. Interaction between HbA and Prx2

The interaction between HbA and Prx2 predicted by IDQD was confirmed by co-IP. Co-IP was performed with anti- beta globin antibody to pull down the interacting complex of HbA from HA and RA samples. HA and RA represent the HbA band of the hemolysate and red blood cell samples during starch-agarose mixed gel electrophoresis, respectively (Figure 3). The immunoprecipitating complexes were separated by 5%-12% sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) (Figure 4). Five protein bands were pulled down by the beta-globin antibody from both the HA and RA samples. These five bands were subsequently cut out, digested, and detected by ESI-Q-TOF. The MS detection results showed that the ~52 kD, ~22 kD and ~16 kD bands were the immunoglobulin heavy chain, Prx2 and beta subunits of hemoglobin, respectively (Tab 3). Thus, Prx2 was confirmed to be a composite component of the beta-globin immunoprecipitating complex. However, the detection results of the ~20 kD and ~30 kD bands were not satisfactory, so they could not be identified in this study.

Figure 3.

Starch-agarose mixed gel electrophoresis of human Hb. H represents a hemolysate sample, and R represents a RBC sample. HA represents hemolysate-HbA, RA represents RBC-HbA, HA2 represents hemolysate-HbA2, and RA2 represents RBC-HbA2.

Figure 4.

Separating Beta-globin immunoprecipitating complexes by SDS-PAGE. Negative control was prepared by immunoprecipitating RA and HA samples with mouse IgG; HA and RA represent hemolysate-HbA and RBC-HbA; M represents a protein marker.

Table 3.MS results of the target band
NO. Name Mass Score Matches Peptide
1 Immunoglobulin heavy chain 51254 2393 52 (39)
2 Undetected - - - -
3 PRDX2 22049 1578 123 (83) K.DVLTITLTPK.V
4 Undetected - - - -
5 Hemoglobin subunit beta 16102 3689 95 (66) K.SAVTALWGK.V

To determine whether Prx2 can interact with HbA2, the freeze-thawed protein extracted from the starch-agarose gel of hemolysate-HbA, RBC-HbA, hemolysate-HbA2, RBC-HbA2, hemolysate-HbA-HbA2 and RBC-HbA-HbA2 were immunoblotted with anti-Prx2 and anti-Alpha-globin (internal control) antibodies, respectively. The results showed that Prx2 was present in hemolysate-HbA, RBC-HbA, RBC-HbA2, hemolysate-HbA-HbA2 and RBC-HbA-HbA2 but not in hemolysate-HbA2 (Figure 5).

Figure 5.

Detecting Prx2 in different Hb components by Western blot. HbA, HbA2 and HbA-A2 represent the protein samples extracted from the starch-agarose gels of HbA, HbA2 and Hb located between HbA and HbA2; H represents hemolysate; and R represents RBCs.

4.4. XAFS to clarify the interaction mechanism between HbA and Prx2

Figure 6A shows the iron K-edge absorption curve of HbA, while Figure 6B shows the iron K-edge absorption curve of the HbA-Prx2 complex. Both curves show the leading-edge absorption peaks in the X-ray absorption near-edge structure (XANES) spectrum of each sample. The formation of a prepeak feature is due to a 1s→3d transition, which is forbidden in the octahedral coordination but is allowed in coordination geometries (distorted octahedral and tetrahedral), which is lacking a center of inversion (48). Behind the main peak, there is a flattened peak that represents the simple sinusoidal behavior due to the octahedral symmetry of the first coordination sphere (49). Our results showed that heme-Fe had the same local geometry, regardless of whether HbA binds with Prx2. Figures 6C and 6D are the Fe K-edge Fourier transform spectra of HbA and the HbA-Prx2 complex. The results showed that the distance of the first shell (the first shell containing a planar ligand formed by four Fe-N ligands and one axial ligand composed of oxygen and nitrogen (48, 50)) had not been changed. However, after the addition of Prx2, the distance of the second and other extended shells increased significantly.

Figure 6.

Detecting the change in the iron K-edge of HbA by XAFS. (A) The iron K-edge absorption curve of the HbA group. (B) The iron K-edge absorption curve of the HbA-Prx2 group. (C) The iron K-edge Fourier transform spectra of the HbA group. (D) The iron K-edge Fourier transform spectra of the HbA-Prx2 group.


The functional analysis of proteins is one of the most extensive topics in biology (51), and the biological functions of proteins are mainly dependent on PPIs. However, experimental methods for studying PPIs, such as the yeast two-hybrid method and co-IP, are time-consuming and expensive. In contrast, computational algorithms can provide high-throughput and low-cost methods for PPI prediction. In this study, the IDQD algorithm was successfully used to predict the interactions between different Hb subunits and between Hb subunits and Prx2 through their protein sequences. At present, there are three main types of testing methods for predictive models: the self-consistent test, jack-knife test and k-fold cross-test. Although the results of the self-consistent test and jack-knife test are better than that of the k-fold cross-test, the k-fold cross-test is more realistic and is considered to be a more rigorous and objective test method in statistics. In this study, the sensitivity, specificity and correlation coefficient of the algorithm are 76.04%, 75.74% and 0.52, respectively, which are better than the previously reported domain-based combination information prediction performance results of the support vector machine (SVM) or other algorithms (52). The total prediction accuracy is 80%, which is obtained by predicting the interactions between different Hb subunits. In addition, Prx2 is thought to interact with itself to form polymers, a prediction that has been confirmed by other studies (27, 53).

The predicted results indicated that Prx2 might interact primarily with HbA but not with HbA2, and these predictions were confirmed by co-IP and Western blotting. Co-IP was performed with a beta-globin antibody to pull down the protein complex from the HA and RA samples, and Prx2 and beta-globin were identified by ESI-Q-TOF in both samples. For the other pulldown bands, the ~52 kD band was identified as an immunoglobulin heavy chain from the antibodies used in this experiment and appeared in the IgG co-IP group. The results of the ~20 kD and ~30 kD bands are not satisfactory. Due to the high abundance of immunoglobulins, even after several rounds of detection, we could not find proteins that match their molecular weight. In addition, the ~16 kD band appeared not only in the beta-globin co-IP group but also in the IgG co-IP group and markers. This phenomenon might be caused by the electrophoresis position of bromophenol blue (electrophoretic indicator) very close to beta-globin. The Western blotting result showed that Prx2 existed in hemolysate-HbA, RBC-HbA, RBC-HbA2, hemolysate-HbA-HbA2 and RBC-HbA-HbA2 but not in hemolysate-HbA2. Our previous study demonstrated that RA2 was a complex of HbA and HbA2 (30); therefore, Prx2 existed in RBC HbA2, and it was presumed to bind to HbA. The interaction mechanism between HbA and Prx2 was studied by XAFS, which was a powerful tool for probing local atomic structures of protein in solution (54, 55) that had been used to detect the structural changes of iron in heme. The Fourier transform spectra of HbA showed that the distance of the first shell of iron had not been changed, but the second and other extension shells of iron were increased significantly after adding Prx2. These results indicate that Prx2 affected the second shell and other extended shells but did not affect the first iron shell in hemoglobin (the first shell contains four Fe-Ns in the planar ligand, and the axial ligand contains oxygen and nitrogen) (48, 50)). This finding indicates that Prx2 interacts with HbA by inserting itself into the hydrophobic gap of heme.

Our experimental results also demonstrated that the IDQD prediction algorithm had significant prediction accuracy. In the future, the accuracy and scope of the IDQD algorithm were expected to be further improved, and the method could be developed into a reliable and accurate computer algorithm for protein function research and the discovery of new drug targets.

There are reversible changes of Prx2 between the oxidized (dimer or decamer) and reductive states (monomer) (56). Generally, the reductive state of Prx2 can interact with Hb, and the oxidized state will bind to cell membranes through band3 and other membrane proteins (57). The change in the redox status of Prx2 is significantly correlated with the RBC oxidative stress caused by different types of anemia (58-60), hereditary spherocytosis (61) and blood storage (62-63). Further studies of the interaction between Prx-2 and Hb will provide a new and more effective method for evaluating the oxidative stress level during some systemic diseases, RBC-related diseases and RBC storage.


Qiang Ma, Liang An, Huifang Tian and Jia Liu equally contributed to this paper contribution. This work was supported by grants from the National Natural Science Foundation of China (8160214, 81860029), Natural Science Foundation of Inner Mongolia (2016MS0801), Scientific Research Foundation of Health and Family Planning Committee of Inner Mongolia (201701090), and Natural Science Foundation of Beijing (7143172). We also especially acknowledge the volunteers who donated blood samples for our research.

Abbreviations: Reactive oxygen species (ROS), hemoglobin (Hb), peroxiredoxin2 (Prx2), increment of diversity with quadratic discriminant analysis (IDQD), co-immunoprecipitation (co-IP), electrospray ionization quadrupole-time of flight (ESI-Q-TOF), X-ray absorption fine structure (XAFS), protein-protein interactions (PPIs), support vector machines (SVMs), random forests (RFs), red blood cells (RBCs), reactive oxygen species (ROSs), glutathione peroxidase (GSH-Px), peroxiredoxins (Prxs), X-ray absorption near-edge structure (XANES), sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE), trifluoroacetic acid (TFA).

I. Sendina-Nadal , Y. Ofran , J. A. Almendral , J. M. Buldu , I. Leyva , D. Q. Li , S. Havlin and S. Boccaletti : Unveiling protein functions through the dynamics of the interaction network. PloS one 6 ( 3 ), e17679 ( 2011) DOI: 10.1371/journal.pone.0017679
D. Bu , Y. Zhao , L. Cai , H. Xue , X. Zhu , H. Lu , J. Zhang , S. Sun , L. Ling , N. Zhang , G. Li and R. Chen : Topological structure analysis of the protein-protein interaction network in budding yeast. Nucleic Acids Res 31 ( 9 ), 2443 2450 ( 2003) DOI: 10.1093/nar/gkg340
T. Hase , H. Tanaka , Y. Suzuki , S. Nakagawa and H. Kitano : Structure of protein interaction networks and their implications on drug design. Plos Comput Biol 5 ( 10 ), e1000550 ( 2009) DOI: 10.1371/journal.pcbi.1000550
T. Schutze , A. K. Ulrich , L. Apelt , C. L. Will , N. Bartlick , M. Seeger , G. Weber , R. Lührmann , U. Stelzl , and M. C. Wahl : Multiple protein-protein interactions converging on the Prp38 protein during activation of the human spliceosome. RNA (New York, NY) 22 ( 2 ), 265 277 ( 2015) DOI: 10.1261/rna.054296.115
J. Song , J. Li , H. D. Liu , W. Liu , Y. Feng , X. T. Zhou and J. D. Li : Snapin interacts with G-protein coupled receptor PKR2. Biochemi Bioph Res Co 469 , 501 506 ( 2016) DOI: 10.1016/j.bbrc.2015.12.023
R. Hegde , S. M. Srinivasula , Z. Zhang , R. Wassell , R. Mukattash , L. Cilenti , G. DuBois , Y. Lazebnik , A. S. Zervos , T. Fernandes-Alnemri and E. S. Alnemri : Identification of Omi/HtrA2 as a mitochondrial apoptotic serine protease that disrupts inhibitor of apoptosis protein-caspase interaction. J Biol Chem 277 ( 1 ), 432 438 . ( 2002). DOI: 10.1074/jbc.M109721200
G. H. Liu , H. B. Shen and D. J. Yu : Prediction of protein-protein interaction sites with machine-learning-based data-cleaning and post-filtering procedures. J Membr Biol 249 , 141 153 ( 2016) DOI: 10.1007/s00232-015-9856-z
P. Fariselli , F. Pazos , A. Valencia and R. Casadio : Prediction of protein-protein interaction sites in heterocomplexes with neural networks. European journal of biochemistry / FEBS 269 ( 5 ), 1356 1361 ( 2002) DOI: 10.1046/j.1432-1033.2002.02767.x
B. Wang , P. Chen , D. S. Huang , J. J. Li , T. M. Lok and M. R. Lyu : Predicting protein interaction sites from residue spatial sequence profile and evolution rate. FEBS letters 580 ( 2 ), 380 384 ( 2006) DOI: 10.1016/j.febslet.2005.11.081
J. Jia , X. Xiao and B. Liu : Prediction of Protein-Protein Interactions with Physicochemical Descriptors and Wavelet Transform via Random Forests. J Lab Autom 21 ( 3 ), 368 377 ( 2016) DOI: 10.1177/2211068215581487
S. H. Park , J. A. Reyes , D. R. Gilbert , J. W. Kim and S. Kim : Prediction of protein-protein interaction types using association rule based classification. BMC bioinformatics 10 , 36 ( 2009) DOI: 10.1186/1471-2105-10-36
E. Sprinzak and H. Margalit : Correlated sequence-signatures as markers of protein-protein interaction. J Mol Biol 11 ( 4 ), 681 692 ( 2001). DOI: 10.1006/jmbi.2001.4920
L. Zhang and L. Luo : Splice site prediction with quadratic discriminant analysis using diversity measure. Nucleic Acids Res 31 ( 21 ), 6214 6220 ( 2003) DOI: 10.1093/nar/gkg805
J. Lu , L. Luo : Prediction for human transcription start site using diversity measure with quadratic discriminant. Bioinformation 2 ( 7 ), 316 321 ( 2008) DOI: 10.6026/97320630002316
W. Chen , Luo , L. Zhang and H. Lin : Recognition of DNase I hypersensitive sites in multiple cell lines. Int J Logist-Res App 5 ( 4 ), 378 384 ( 2009) DOI: 10.1504/IJBRA.2009.027508
Y. Feng , L. Luo : Use of tetrapeptide signals for protein secondary-structure prediction. Amino Acids 35 ( 3 ), 607 614 ( 2008) DOI: 10.1007/s00726-008-0089-7
W. Chen , L. Luo : Classification of antimicrobial peptide using diversity measure with quadratic discriminant analysis. J Microbiol Methods 78 ( 1 ), 94 96 ( 2009) DOI: 10.1016/j.mimet.2009.03.013
W. Chen , L. Luo and L. Zhang : The organization of nucleosomes around splice sites. Nucleic Acids Res 38 ( 9 ), 2788 2798 ( 2010) DOI: 10.1093/nar/gkq007
Y. Xing , X. Zhao and L. Cai : Prediction of nucleosome occupancy in Saccharomyces cerevisiae using position-correlation scoring function. Genomics 98 ( 5 ), 359 366 ( 2011) DOI: 10.1016/j.ygeno.2011.07.008
X. Zhao , Z. Pei , J. Liu , S. Qin and L. Cai : Prediction of nucleosome DNA formation potential and nucleosome positioning using increment of diversity combined with quadratic discriminant analysis. Chromosome Res 18 ( 7 ), 777 85 ( 2010) DOI: 10.1007/s10577-010-9160-9
G. Liu , J. Liu , X. Cui and L. Cai : Sequence-dependent prediction of recombination hotspots in Saccharomyces cerevisiae. J theoretical biol 293 , 49 54 ( 2012) DOI: 10.1016/j.jtbi.2011.10.004
Y. Su , G. Shao , L. Gao , L. Zhou , L. Qin and W. Qin : RBC electrophoresis with discontinuous power supply a newly established hemoglobin release test. Electrophoresis 30 ( 17 ), 3041 43 ( 2009) DOI: 10.1002/elps.200900176
T. H. Lee , S. U. Kim , S. L. Yu , S. H. Kim , D. S. Park , H. B. Moon , S. H. Dho , K. S. Kwon , H. J. Kwon , Y. H. Han , S. Jeong , S. W. Kang , H. S. Shin , K. K. Lee , S. G. Rhee and D. Y. Yu : Peroxiredoxin II is essential for sustaining life span of erythrocytes in mice. Blood 101 ( 12 ), 5033 38 ( 2003) DOI: 10.1182/blood-2002-08-2548
D. G. Kakhniashvili , Jr. Bulla L. A. and S. R. Goodman : The human erythrocyte proteome: analysis by ion trap mass spectrometry. Mol Cell Proteomics 3 ( 5 ), 501 509 ( 2004) DOI: 10.1074/mcp.M300132-MCP200
R. M. Johnson , Jr. Goyette G. , Y. Ravindranath and Y. S. Ho : Hemoglobin autoxidation and regulation of endogenous H 2O 2 levels in erythrocytes. Free Radic Biol Med 39 ( 11 ), 1407 1417 ( 2005) DOI: 10.1016/j.freeradbiomed.2005.07.002
K. M. Stuhlmeier , J. J. Kao , P. Wallbrandt , M. Lindberg , B. Hammarstrom , H. Broell and B. Paigen : Antioxidant protein 2 prevents methemoglobin formation in erythrocyte hemolysates. Eur J Biochem / FEBS 270 ( 2 ), 334 341 ( 2003) DOI: 10.1046/j.1432-1033.2003.03393.x
F. M. Low , M. B. Hampton and C. C. Winterbourn : Peroxiredoxin 2 and peroxide metabolism in the erythrocyte. Antioxid Redox Signal 10 ( 9 ), 1621 1630 ( 2008) DOI: 10.1089/ars.2008.2081
M. Brizuela , H. M. Huang , C. Smith , G. Burgio , S. J. Foote and B. J. McMorran : Treatment of erythrocytes with the 2-cys peroxiredoxin inhibitor, Conoidin A, prevents the growth of Plasmodium falciparum and enhances parasite sensitivity to chloroquine. PloS one 9 ( 4 ), e92411 ( 2014) DOI: 10.1371/journal.pone.0092411
S. B. Bayer , G. Maghzal , R. Stocker , M. B. Hampton and C. C. Winterbourn : Neutrophil-mediated oxidation of erythrocyte peroxiredoxin 2 as a potential marker of oxidative stress in inflammation. FASEB J 27 ( 8 ), 3315 3322 ( 2013) DOI: 10.1096/fj.13-227298
Y. Su , L. Gao , Q. Ma , L. Zhou , L. Qin , L. Han and W. B. Qin : Interactions of hemoglobin in live red blood cells measured by the electrophoresis release test. Electrophoresis 31 ( 17 ), 2913 20 ( 2010) DOI: 10.1002/elps.201000034
A. Basu and A. Chakrabarti : Hemoglobin interacting proteins and implications of spectrin hemoglobin interaction. J Proteomics 128 : 469 75 ( 2015) DOI: 10.1016/j.jprot.2015.06.014
L. Salwinski , C. S. Miller , A. J. Smith , F. K. Pettit , J. U. Bowie and D. Eisenberg : The Database of Interacting Proteins: 2004 update. Nucleic Acids Res 32 ( Database issue ), D449 451 ( 2004) DOI: 10.1093/nar/gkh086
B. E. Suzek , H. Huang , P. McGarvey , R. Mazumder and C. H. WuUniRef : comprehensive and non-redundant UniProt reference clusters. Bioinformatics 23 ( 10 ), 1282 1288 ( 2007) DOI: 10.1093/bioinformatics/btm098
L. Cai , J. Liu and X. Zhao : Protein-Protein Interaction Prediction Using Increment of Diversity Combined with Quadratic Discriminant Analysis. The 2010 International Congress on Computer Applications and Computational Science; Singapore; 613 616 ( 2010)
Y. Su , L. Gao and W. Qin : Interactions of hemoglobin in live red blood cells measured by the electrophoresis release test. Methods Mol Biol 869 , 393 402 ( 2012) DOI: 10.1007/978-1-61779-821-4_32
Y. Su , J. Shen , L. Gao , H. F. Tian , Z. H Tian and W. B. Qin : Molecular interactions of re-released proteins in electrophoresis of human erythrocytes. Electrophoresis, 33 ( 9-10 ), 1402 05 ( 2012) DOI: 10.1002/elps.201100644
J. M. Tranquada and R. Ingalls : Extended x-ray\char22{}absorption fine-structure study of anharmonicity in CuBr. Physical Review B, 28 ( 6 ), 3520 28 ( 1983) DOI: 10.1103/PhysRevB.28.3520
L. R. Manning , J. E. Russell , J. C. Padovan , B. T. Chait , A. Popowicz , R. S. Manning and J. M. Manning : Human embryonic, fetal, and adult hemoglobins have different subunit interface strengths. Correlation with lifespan in the red cell. Protein Sci 16 ( 8 ), 1641 1658 ( 2007) DOI: 10.1110/ps.072891007
L. R. Manning , J. E. Russell , A. M. Popowicz , R. S. Manning , J. C. Padovan and J. M. Manning : Energetic differences at the subunit interfaces of normal human hemoglobins correlate with their developmental profile. Biochemistry, 48 ( 32 ), 7568 74 ( 2009) DOI: 10.1021/bi900857r
MJ. McDonald : Assembly of human adult and sickle hemoglobins from their oxygenated subunits. Differential rates of beta chain tetramer dissociation. J Biol Chem 256 ( 12 ), 6487 6490 ( 1981)
E. R. Huehns and E. M. Shooter : Human Haemoglobins. J Med Genet 2 ( 1 ), 48 90 ( 1965) DOI: 10.1136/jmg.2.1.48
E. R. Huehns , G. H. Beaven and B. L. Stevens : Reaction of haemoglobin alpha-A with haemoglobins beta-A4, gamma-F4 and delta-A2. Biochem J 92 ( 2 ): 444 448 ( 1964) DOI: 10.1042/bj0920444
R. M. Nalbandian , R. L. Henry , L. F. J. R. Camp , P. L. Wolf and T. N. Evans : Embryonic, fetal, and neonatal hemoglobin synthesis: relationship to abortion and thalassemia. Obstet Gynecol Surv 26 ( 2 ), 184 191 ( 1971) DOI: 10.1097/00006254-197102000-00027
U. Sen , J. Dasgupta , D. Choudhury , P. Datta , A. Chakrabarti , S. B. Chakrabarty , A. Chakrabarty and J. K Dattagupta : Crystal structures of HbA2 and HbE and modeling of hemoglobin delta 4: interpretation of the thermal stability and the antisickling effect of HbA2 and identification of the ferrocyanide binding site in Hb. Biochemistry, 43 ( 39 ), 12477 12488 ( 2004) DOI: 10.1021/bi048903i
L. R. Manning , A. M. Popowicz , J. Padovan , B. T. Chait , J. E. Russell and J. M. Manning : Developmental expression of human hemoglobins mediated by maturation of their subunit interfaces. Protein Sci 19 ( 8 ), 1595 1599 ( 2010) DOI: 10.1002/pro.441
E. P. Vichinsky : Clinical manifestations of alpha-thalassemia. Cold Spring Harb Perspect Med 3 ( 5 ), a011742 ( 2013) DOI: 10.1101/cshperspect.a011742
D. R. Higgs : The molecular basis of alpha-thalassemia. Cold Spring Harb Perspect Med 3 ( 1 ), a011718 ( 2013) DOI: 10.1101/cshperspect.a011718
C. Jin , Y. Li , Y. L. Li , Y. Zou , G. L. Zhang , M. Normura and G. Y. Zhu : Blood lead: Its effect on trace element levels and iron structure in hemoglobin. Nucl Instrum Meth B 266 ( 16 ), 3607 3613 ( 2008) DOI: 10.1016/j.nimb.2008.05.087
P. D’Angelo and M. Benfatto : Effect of Multielectronic ConFigurations on the XAFS Analysis at the Fe K Edge. J Phys Chem A 108 ( 20 ), 4505 4514 ( 2004) DOI: 10.1021/jp0499732
P. Eisenberger , R. G. Shulman , B. M. Kincaid , G. S. Brown and S. Ogawa : Extended X-ray absorption fine structure determination of iron nitrogen distances in haemoglobin. Nature, 274 ( 5666 ): 30 34 ( 1978) DOI: 10.1038/274030a0
L. Laraia , G. McKenzie , D. R. Spring , A. R. Venkitaraman and D. J. Huggins : Overcoming Chemical, Biological, and Computational Challenges in the Development of Inhibitors Targeting Protein-Protein Interactions, Chem Biol 22 ( 6 ), 689 703 ( 2015) DOI: 10.1016/j.chembiol.2015.04.019
S. Martin , D. Roe and J. L. Faulon : Predicting protein-protein interactions using signature products. Bioinformatics 21 ( 2 ), 218 226 ( 2005) DOI: 10.1093/bioinformatics/bth483
E. Schroder , J. A. Littlechild , A. A. Lebedev , N. Errington , A. A. Vagin and M. N : Isupov Crystal structure of decameric 2-Cys peroxiredoxin from human erythrocytes at 1.7 A resolution. Structure 8 ( 6 ), 605 15 ( 2000) DOI: 10.1016/S0969-2126(00)00147-7
T. L. Sorensen , K. E. McAuley , R. Flaig and E. M. Duke : New light for science synchrotron radiation in structural medicine. Trends in biotechnology 24 ( 11 ), 500 08 ( 2006) DOI: 10.1016/j.tibtech.2006.09.006
Z. Wu , R. E. Benfield , Y. Wang , L. Guo , M. Tan , H. Zhang , Y. Ge and D. Grandiean : EXAFS study on the local atomic structures around iron in glycosylated haemoglobin. Phys Med Biol 46 ( 3 ), N71 7 7 ( 2001) DOI: 10.1088/0031-9155/46/3/403
E. Nagababu , J. G. Mohanty , J. S. Friedman and J. M. Rifkind : Role of peroxiredoxin-2 in protecting RBCs from hydrogen peroxide-induced oxidative stress. Free Radic Res 47 ( 3 ), 164 71 ( 2013). DOI: 10.3109/10715762.2012.756138
A. Matte , M. Bertoldi , N. Mohandas , X. L. An , A. Bugatti , A. M. Brunati , M. Rusnati , E. Tibaldi , A. Siciliano , F. Turrini , S. Perrotta and L. D. Franceschi : Membrane association of peroxiredoxin-2 in red cells is mediated by the N-terminal cytoplasmic domain of band 3. Free Radic Biol Med 55 , 27 35 ( 2013). DOI: 10.1016/j.freeradbiomed.2012.10.543
Y. H. Han , S. U. Kim , T. H. Kwon , D. S. Lee , H. L. Ha , D. S. Park , E. J. Woo , S. H. Lee , J. M. Kim , H. B. Chae , S. Y. Lee , B. Y. Kim , D. Y. Yoon , S. G. Rhee , E. Fibach and D. Y. Yu : Peroxiredoxin II is essential for preventing hemolytic anemia from oxidative stress through maintaining hemoglobin stability. Biochem Biophys Res Commun 426 ( 3 ), 427 32 (2012). DOI: 10.1016/j.bbrc.2012.08.113
F. C. Cheah , A. V. Peskin , F. L. Wong , A. Ithnin , A. Othman and C. C. Winterbourn : Increased basal oxidation of peroxiredoxin 2 and limited peroxiredoxin recycling in glucose-6-phosphate dehydrogenase-deficient erythrocytes from newborn infants. FASEB J, 28 ( 7 ): 3205 10 ( 2014). DOI: 10.1096/fj.14-250050
A. Mattea , P. S. Lowb , F. Turrinic , M. Bertoldid , M. E. Campanellab , D. Spanoe , A. Pantaleoa,f , A. Sicilianoa and L. D. Franceschia : Peroxiredoxin-2 expression is increased in beta-thalassemic mouse red cells but is displaced from the membrane as a marker of oxidative stress. Free Radic Biol Med, 49 ( 3 ): 457 66 ( 2010). DOI: 10.1016/j.freeradbiomed.2010.05.003
S. Rocha , R. M. P. Vitorino , F. M. Lemos-Amado , E. B. Castro , P. Rocha-Pereira , J. Barbot , E. Cleto , F. Ferreira , A. Quintanilha , L. Belo , A. Santos-Silva : Presence of cytosolic peroxiredoxin 2 in the erythrocyte membrane of patients with hereditary spherocytosis. Blood Cells Mol Dis, 41 ( 1 ): 5 9 ( 2008). DOI: 10.1016/j.bcmd.2008.02.008
S. Rinalducci , G. M. D’Amici , B. Blasi , S. Vaglio , G. Grazzini and L. Zolla : Peroxiredoxin-2 as a candidate biomarker to test oxidative stress levels of stored red blood cells under blood bank conditions. Transfusion, 51 ( 7 ): 1439 49 ( 2011). DOI: 10.1111/j.1537-2995.2010.03032.x
J. Y. Oh , V. Harper , C. W. Sun , R. Stapley , L. Wilson , M. B. Marques , S. Barnes , T. Townes and R. P. PaTel : Peroxiredoxin-2 recycling is inhibited during erythrocyte storage. Antioxid Redox Signal, 22 ( 4 ): 294 307 ( 2015). DOI: 10.1089/ars.2014.5950
Back to top