Hyperoxidation of Peroxiredoxin 6 Induces Alteration from Dimeric to Oligomeric State

Peroxiredoxins(Prdx), the family of non-selenium glutathione peroxidases, are important antioxidant enzymes that defend our system from the toxic reactive oxygen species (ROS). They are thiol-based peroxidases that utilize self-oxidation of their peroxidatic cysteine (Cp) group to reduce peroxides and peroxidized biomolecules. However, because of its high affinity for hydrogen peroxide this peroxidatic cysteine moiety is extremely susceptible to hyperoxidation, forming peroxidase inactive sulfinic acid (Cys-SO2H) and sulfonic acid (Cys-SO3H) derivatives. With the exception of peroxiredoxin 6 (Prdx6), hyperoxidized sulfinic forms of Prdx can be reversed to restore peroxidase activity by the ATP-dependent enzyme sulfiredoxin. Interestingly, hyperoxidized Prdx6 protein seems to have physiological significance as hyperoxidation has been reported to dramatically upregulate its calcium independent phospholipase A2 activity. Using biochemical studies and molecular dynamic (MD) simulation, we investigated the roles of thermodynamic, structural and internal flexibility of Prdx6 to comprehend the structural alteration of the protein in the oxidized state. We observed the loosening of the hydrophobic core of the enzyme in its secondary and tertiary structures. These changes do not affect the internal dynamics of the protein (as indicated by root-mean-square deviation, RMSD and root mean square fluctuation, RMSF plots). Native-PAGE and dynamic light scattering experiments revealed the formation of higher oligomers of Prdx6 under hyperoxidation. Our study demonstrates that post translational modification (like hyperoxidation) in Prdx6 can result in major alterations of its multimeric status.


Introduction
Reactive oxygen species (ROS) are generated as a result of normal cellular metabolism such as aerobic respiration, photosynthesis, and exposure to environmental stimuli [1]. Depending on their concentration in biological systems, ROS can either be harmful or beneficial to living systems [2]. Prolonged presence of these imbalanced reactive oxygen intermediates further initiates a chain reaction where any macromolecule coming into contact with them is oxygenated, thereby exponentially enhancing the ROS concentration and cellular damage. In certain cases, it is very difficult to distinguish whether oxidative stress is the cause or an effect of a disease. functioning of both activities without hampering any other biological process. As Prdx6 has affinity for oxidized substrate both at acidic and neutral pH, it shows aiPLA2 activity at both pH to play an important role in the repair of peroxidized cell membranes [16][17][18].
Like 2-cys Prdx, the conserved peroxidatic catalytic residue Cys-47 of Prdx6 is also prone to hyperoxidation forming hyperoxidized forms, Cys-SO 2 H and Cys-SO 3 H [19]. It should be noted that hyperoxidized Prdx6 is not a substrate for sulfiredoxin and therefore cannot be regenerated like other Prdxs, even from the Cys-SO 2 H form, making hyperoxidation irreversible and functional inactivation for Prdx6. Interestingly, it is already reported that H 2 O 2 -mediated hyperoxidation of Prdx6 has been shown to induce cell cycle arrest at the G2/M transition through up-regulation of its aiPLA 2 activity though the hyperoxidized forms do not have peroxidase activity [19]. This effect of hyperoxidation on PLA2 activity is confirmed in a report by Zhou et.al. Our aim is to comprehend the structural alterations that incur due to hyperoxidation of Prdx6's peroxidatic cysteine residue. These structural changes could be the mechanism that allows Prdx6 to up-regulate its aiPLA2 activity by accepting reduced phospholipid as its substrate within the cytosol. Until now, to our best knowledge, only a limited number of studies have been conducted to comparatively analyze the hyperoxidized and reduced forms of Prdx6. Here, we used molecular dynamic (MD) simulations along with various spectroscopic methods to investigate the molecular, conformational and thermodynamic studies of the different oxidation states of rat Prdx6. We found that hyperoxidation of Cys47 of Prdx6 induces changes at secondary, and tertiary as well as quaternary levels in the structure of enzyme, which might be responsible for the upregulation of its total cellular aiPLA 2 activity. For the first time, we are reporting the existence of an oligomeric form of Prdx6 enzyme.

Materials
Analytical grade Tris-HCl, Dithiothreitol (DTT), and IPTG samples were purchased from the Sigma-Aldrich Corp. (St. Louis, MO, USA). H 2 O 2 was obtained from MP Biomedicals (Santa Ana, CA, USA). NaCl was from Merck (Kenilwort, NJ, USA). This and other chemicals were analytical grade reagents and were used without further purification. 8-Anilinonaphthalene-1-sulfonate (ANS) was also purchased from the Sigma-Aldrich Corp. (St. Louis, MO, USA). Chitin resin beads were purchased from New England Biolabs (NEB) (Ipswich, MA, USA).

Preparation of Recombinant Prdx6
The cloning of the target gene of rat Prdx6 into the pTyB1 vector was done as described previously [20,21]. The construct contained Nde1 (at 5 end) and sap1 (at 3 end) as restriction sites. The restriction site of Sap1 is not generated after recombinant cloning. Thus, the nucleotide sequence of rPrdx6 is immediately followed by that of an intein tag derived from pTyB1 vector ( Figure 1A,B). The recombinant plasmid was constructed to over-express the protein conjugated with an intein tag. During the purification, the protein was eluted by cleaving the intein tag, the cleavage occurred between the last residue of rPrdx and the Cys residue of the intein tag. The purified protein consists of 224 amino acids, as in the cDNA, without any vector derived residues. The recombinant vector was transformed into the E. coli BL21 (DE3) expression strain. A single colony of transformant was selected and inoculated in luria bertani (LB) medium containing ampicillin (50 µg/mL). Cells were grown at 37 • C in a shaker incubator overnight. When the optical density (at 600 nm) of the growing cells reached 0.6-0.8, isopropyl-β-D-thiogalactopyranoside (0.6 mM) was added for its induction. Then, the cells were grown at 20 • C overnight. After harvesting by centrifugation, the cell pellet was suspended in lysis buffer (20 mM Tris-HCl pH 7.5, 500 mM NaCl and 1 mM Ethylenediaminetetraacetic acid (EDTA)) and sonicated on ice (Bandelin sonicator, Bandelin, Germany) with five pulses of 10 s at an interval of 5 min. The lysate after sonication was centrifuged at 10,000 rpm for 20 min. The supernatant obtained after centrifugation was loaded onto a chitin affinity column equilibrated with 20 mM Tris-HCl pH 7.5, 500 mM NaCl and 1 mM EDTA. Induction of the on-column cleavage was performed by quickly flushing the column with cleavage buffer containing 20 mM Tris-HCl pH 8.5, 500 mM NaCl in presence of 80 mM DTT. After 48 h incubation of the inducted column at room temperature, the target protein was eluted with column buffer (20 mM Tris-HCl pH 7.5, 500 mM NaCl and 1 mM EDTA). The purified protein was analyzed by running 10% sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) and found to be more than 90% pure as shown in Figure 1C. Induction of the on-column cleavage was performed by quickly flushing the column with cleavage buffer containing 20 mM Tris-HCl pH 8.5, 500 mM NaCl in presence of 80 mM DTT. After 48 h incubation of the inducted column at room temperature, the target protein was eluted with column buffer (20 mM Tris-HCl pH 7.5, 500 mM NaCl and 1 mM EDTA). The purified protein was analyzed by running 10% sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) and found to be more than 90% pure as shown in Figure 1C. Details of nucleotide sequence of plasmid expressing recombinant rat peroxiredoxin 6 (rPrdx6). The 5′ terminal underlined and italics indicate the derived primer that was used. The 3′ terminal sequence is vector derived nucleotides after recombinant cloning. Inset view shows the different nucleotides of recombinant plasmid to cDNA sequence of rPrdx6 at 3′ terminal (A). The amino acid sequence of the prepared construct contains an intein tag for purification. Underlined sequence encodes for the vector derived amino acids. The cleavage site of the protein from purification tag is shown with an arrow (B). Sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) of purified Prdx6 using chitin affinity chromatography. The purity is more than 90% (C).

Preparation of Reduced and Hyperoxidized Prdx6
The purified protein was used for preparation of reduced Prdx6 by adding 1 mM DTT in standard buffer i.e., 50 mM Tris-HCl, 100 mM NaCl pH 7.4. We incubated the protein with 500 µM H2O2 to get the hyperoxidized Prdx6 (C47 sulfinic acid or sulfonic acid) at room temperature for 30 min to 1 h in the standard buffer. At this concentration of H2O2 the antibody against the oxidized protein detects both the sulfinic (-SO2H) and sulfonic (-SO3H) states of the protein [15,19].

Circular Dichroism (CD) Measurements
CD measurements of reduced Prdx6 and hyperoxidized Prdx6 were recorded in 50 mM Tris-HCl, 100 mM NaCl, pH 7.4 using a Jasco-750 CD Spectropolarimeter. All the measurements were done at 25 °C in a thermoelectric cell holder, with the temperature being maintained using a Peltier element. Spectra were recorded in the far-ultraviolet region (190-260 nm), with a bandwidth of 1.0 nm, a step size of 1 nm, an integration time of 30 s, and with three repeats. A fused quartz cell with a pathlength of 0.1 cm was used. The protein concentration was 10 µM. The results of all the CD measurements are expressed as mean residue ellipticity ([θ]λ) in deg cm 2 dmol −1 at a given wavelength λ (nm) using the relation: [θ]λ = θλMo/10cl, where θλ is the observed ellipticity in Figure 1. Details of nucleotide sequence of plasmid expressing recombinant rat peroxiredoxin 6 (rPrdx6). The 5 terminal underlined and italics indicate the derived primer that was used. The 3 terminal sequence is vector derived nucleotides after recombinant cloning. Inset view shows the different nucleotides of recombinant plasmid to cDNA sequence of rPrdx6 at 3 terminal (A). The amino acid sequence of the prepared construct contains an intein tag for purification. Underlined sequence encodes for the vector derived amino acids. The cleavage site of the protein from purification tag is shown with an arrow (B). Sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) of purified Prdx6 using chitin affinity chromatography. The purity is more than 90% (C).

Preparation of Reduced and Hyperoxidized Prdx6
The purified protein was used for preparation of reduced Prdx6 by adding 1 mM DTT in standard buffer i.e., 50 mM Tris-HCl, 100 mM NaCl pH 7.4. We incubated the protein with 500 µM H 2 O 2 to get the hyperoxidized Prdx6 (C47 sulfinic acid or sulfonic acid) at room temperature for 30 min to 1 h in the standard buffer. At this concentration of H 2 O 2 the antibody against the oxidized protein detects both the sulfinic (-SO 2 H) and sulfonic (-SO 3 H) states of the protein [15,19].

Circular Dichroism (CD) Measurements
CD measurements of reduced Prdx6 and hyperoxidized Prdx6 were recorded in 50 mM Tris-HCl, 100 mM NaCl, pH 7.4 using a Jasco-750 CD Spectropolarimeter. All the measurements were done at 25 • C in a thermoelectric cell holder, with the temperature being maintained using a Peltier element. Spectra were recorded in the far-ultraviolet region (190-260 nm), with a bandwidth of 1.0 nm, a step size of 1 nm, an integration time of 30 s, and with three repeats. A fused quartz cell with a pathlength of 0.1 cm was used. The protein concentration was 10 µM. The results of all the CD measurements are expressed as mean residue ellipticity ([θ] λ ) in deg cm 2 dmol −1 at a given wavelength λ (nm) using the relation: [θ] λ = θ λ M o /10cl, where θ λ is the observed ellipticity in millidegrees at wavelength λ, Mo is the mean residual weight of the protein, c is the protein concentration (mg/cm 3 ), and l is the path length (cm).

Fluorescence Measurements
Fluorescence spectroscopy was performed with a PTI spectrofluorometer (Photon Technology International, Inc., Lawrenceville, NJ, USA) equipped with a single photon counting system for fluorescence intensity detection, dual fluorescence and absorbance channels using excitation and emission slits of 1 nm each. All measurements were performed at 25 • C in 50 mM Tris-HCl, 100 mM NaCl, pH 7.4 buffer. For tryptophan fluorescence measurements, the emission spectra were recorded in the wavelength range 310-450 nm after excitation at 295 nm to avoid tyrosine fluorescence.
8-Anilinonaphthalene-1-sulfonate (ANS) fluorescence spectra of Prdx6 were collected from 400 to 600 nm after excitation at 360 nm. The protein concentration used for Trp and ANS fluorescence measurements was 1 µM. The ANS concentration in ANS fluorescence was 16 µM. All the measurements were done in micro quartz fluorescence cuvettes with pathlength of 0.3 cm (Sterna).

Thermal-Induced Denaturation
Heat-induced denaturation of Prdx6 was carried out in a Jasco-750 CD Spectropolarimeter equipped with a Peltier-type temperature controller with a heating rate of 1 • C/min, a scan rate providing adequate time for equilibration. Thermal denaturation was recorded following changes in [θ] 220 from 20 • C to 80 • C at a rate of 1 • C/min in Jasco-750 CD Spectropolarimeter. A fused quartz cell with a pathlength of 1 cm was used. The protein concentration was 1 µM. The CD measured mean residual ellipticity at [θ] 220 as a function of temperature is normalized at the start of the temperature scan. As shown earlier [22], each thermal denaturation curve was fitted to a two-state unfolding model after assuming a linear dependence of pre-and post-transition baselines on temperature [23,24]: y(T) is the observed mean residue ellipticity at a given temperature, m N and m D are slopes of the native and denatured baselines, while y N and y D are the intercept of the native and denatured baselines, respectively. T is the temperature in kelvin, T m is the melting temperature in kelvin and ∆H m is the enthalpy change of denaturation at melting temperature and R is the universal gas constant. Curve fitting of the data points was done in Sigma Plot 13.0 software (Systat Software, Inc., San Jose, CA, USA) using Equation (1) in the dynamic fit.

Dynamic Light Scattering Measurements
The size distribution of the particles present in the protein sample was obtained using a Zetasizer Micro V/ZMV 2000 (Malvern, UK). The reduced and hyperoxidized samples were in the standard buffer, i.e., 50 mM Tris-HCl, 100 mM NaCl pH 7.4. The protein concentration was 2 mg/mL. All measurements were performed at 25 • C. Measurements were made at a fixed angle of 90 • using an incident laser beam of 689 nm. Fifteen measurements were made with an acquisition time of 30 s for each sample at sensitivity of 10%. The data was analyzed using Zetasizer software provided by the manufacturer to get hydrodynamic diameters of the particles and the volume fraction of particles associated with that particular diameter. The hydrodynamic data obtained from this measurement was used to calculate the apparent molecular weight of reduced and hyperoxidized Prdx6.

Native PAGE
Reduced and hyperoxidized Prdx6 in the standard buffer, i.e., 50 mM Tris-HCl, 100 mM NaCl pH 7.4 were loaded in a continuous 7.5% polyacrylamide gel under non-reducing and non-denaturing condition. Proteins were electrophoresed using running buffer of 30 mM beta-alanine and 20 mM lactic acid (85-90%), pH 3.8 at a constant voltage (180 V) for 2 h with reverse polarity. The gel was run in the presence of two markers, bovine albumin (Mr 66.5 kDa, pI 4.7-4.9) and high molecular weight horse ferritin (Mr 450 kDa, pI 4.4) under non-reducing and non-denaturing conditions. Coomassie brilliant blue G-250 was used for the detection of the proteins on the gel.

Homology Modelling
The complete sequence of rat Prdx6 (consisting of 224 amino acids) was retrieved from UniProtKB database (accession number O35244). The best three templates, based on sequence identity (91.52%, 91.52%, 91.52%) and coverage (100%, 100% and 100%), of models with dimer option were generated from the selected templates (PDB codes: 1PRX, 5B6M and 5B6N), respectively. The reduced Cys47 was converted to hyperoxidized (Cys47-SH to Cys47-SO 3 H) form in Maestro [25]. To remove any steric clashes, energy minimization was performed using the GROMOS96 force field implemented in Swiss-Pdb Viewer [26].
PROCHECK [27] was used to verify the backbone conformation of the modeled structures, and the quality of the models was evaluated by qualitative model energy analysis (QMEAN) [28]. The QMEAN score gives the global score of the model and it ranges from 0 to 1, the higher the scores, the more reliable are the models. Further, the QMEAN Z-score provides an estimation of the "degrees of nativeness" of the structural features observed in a model and indicates model quality compared to experimental structures. The model having the best validation results was selected for further studies.

Molecular Dynamics Simulation
Desmond [29] was used to perform MD simulations by implementing OPLS-AA 2005 force field [30,31] for reduced and hyperoxidized rat Prdx6 in its dimeric state. Protein Preparation wizard was used to prepare the protein structure in Schrodinger. TIP3P water model was selected to solvate the protein, the orthorhombic periodic boundary box was chosen to construct the required systems for subsequent MD simulations. Appropriate numbers of counter ions were added to the system to maintain charge neutrality. Further, the distance between the box and wall was set to be greater than 10 Å in order to avoid direct interaction with its own periodic image of the protein complex. Steepest descent method was selected to minimize the potential energies of the systems by applying a maximum of 5000 steps with a gradient threshold of 25 kcal mol −1 Å −1 , followed by L-BFGS (low-memory Broyden-Fletcher-Goldfarb-Shanno quasi-Newtonian) minimizer until a convergence criterion of 1 kcal mol −1 Å −1 was achieved.
Default parameters were selected to equilibrate the pressure and temperature. Thereafter, the final MD simulation was done on the equilibrated systems for 50 ns at a constant temperature of 300 K, and constant pressure of 1 atm with a time step of 2 fs. To compute long-range electrostatic interactions, the particle-mesh Ewald method (PME) [32] was applied with a grid spacing of 0.8 Å. The van der Waals and short-range electrostatic interactions were smoothly truncated at 9.0 Å.
Root mean square deviation (RMSD), root mean square fluctuation (RMSF), radius of gyration, solvent accessible surface area (SASA), internal hydrogen bond contacts of the reduced and hyperoxidized Prdx6 and PCA were calculated using VMD tcl scripts.

Results and Discussion
It is universally known that the ultimate folded native state conformation, not the primary structure of a protein determines the functional capacity and efficiency of a protein. Therefore, for multifunctional proteins (such as Prdx6, aconitase, activating transcription factor 2, ERK2, MAP kinase, glutamate racemase, etc.) usage of native state conformational change as a mechanism to switch between two different activities seems very common [22,33]. To comprehend structural alterations as a possibility for upregulation of aiPLA2 activity due to the hyperoxidation of peroxidatic cysteine residue of Prdx6 [19], we first did the in vitro comparison of the secondary and tertiary structure of hyperoxidized and reduced Prdx6 via various spectroscopic methods (Figure 2). Tryptophan fluorescence emission is a very sensitive tool for measuring the polarity of the local environment of tryptophan residues [34], and hence, mostly used as a diagnostic probe for analyzing conformational state of proteins [35]. Prdx6 has three tryptophan (Trp) residues, Trp33, Trp82 and Trp181. The solvent exposure of these Trp residues under different redox conditions is assessed by measuring tryptophan fluorescence of reduced and hyperoxidized Prdx6. As shown in Figure 2A, a decrease in fluorescence intensity and red shift from 328 nm to 334 nm of the hyperoxidized Prdx6 in comparison to that of reduced Prdx6 clearly reveals that Trp residues of hyperoxidized Prdx6 (relative to reduced protein) have been shifted to a more polar solvent environment. We were further interested to see the impact of oxidation mediated conformational change on the thermodynamic stability of Prdx6. To investigate this, we measured the thermodynamic stability of reduced and hyperoxidized Prdx6 via an optical method, i.e., by monitoring the changes in molar ellipticity at 220 nm over the temperature range of 20-80 °C ( Figure 2D). The transition curve obtained from the thermal-induced denaturation was analyzed using a two-state unfolding model, given above in equation (1). The analysis gave the change in melting temperature (Tm) and change in enthalpy (ΔHm) at the melting temperature of Prdx6, as shown in Table 1  To confirm the likelihood of greater exposition of hydrophobic amino acids in hyperoxidized Prdx6 (as suggested by the Trp fluorescence results), we investigated the ANS binding behaviour of reduced and hyperoxidized Prdx6. The fluorescent hydrophobic probe 8-Anilinonaphthalene-1-sulfonate (ANS) is usually used to detect the accessible hydrophobic surfaces on proteins [36,37]. In aqueous solvent, the quantum yield of the ANS probe is very weak, which on binding to a hydrophobic surface increases several fold [38]. As seen in Figure 2B, ANS fluorescence measurement in aqueous solvent shows an emission maximum at~520 nm, and at 519 nm and 515 nm with reduced and hyperoxidized Prdx6, respectively. And there is increase in fluorescence intensity of the hyperoxidized Prdx6 in comparison to that of the aqueous solvent and reduced Prdx6. The blue shift in emission maxima and increase in the fluorescence intensity of ANS fluorescence with hyperoxidized Prdx6 as compared to the aqueous solvent is a clear signature for ANS binding, and indicates the exposure of hydrophobic amino acids to the solvent in the hyperoxidized Prdx6. Such exposure of hydrophobic amino acids in oxidized Prdx6, suggests conformational alterations with opening of the hydrophobic core.
We further analyzed whether the disparity in the hydrophobic core's packing within reduced and hyperoxidized Prdx6 is limited to only tertiary alterations or also has its root in the changes at secondary structure level. The far-UV region of the circular dichroism (CD) spectrum monitors the changes in the secondary structure of protein [39,40]. Here we evaluated the effect of redox state on the secondary structure of Prdx6 by performing CD measurements with reduced and hyperoxidized Prdx6 in the far-UV region. The mean residual ellipticity at θ 222 , the measure of secondary structure content of Prdx6, is shown in Table 1. Figure 2C and Table 1 show the far-UV CD spectra of both proteins, and the decrease in secondary structure content of hyperoxidized Prdx6 (in comparison with reduced Prdx6) indicates the influence of the redox state of Prdx6 for the functional switching of the protein by altering its secondary as well as tertiary structure. Taken together, the results indicate that there is a conformational difference between the reduced and hyperoxidized protein. However, we do not rule out the possibility that oligomer formation may also contribute to the observed spectral properties of the hyperoxidized protein.
We were further interested to see the impact of oxidation mediated conformational change on the thermodynamic stability of Prdx6. To investigate this, we measured the thermodynamic stability of reduced and hyperoxidized Prdx6 via an optical method, i.e., by monitoring the changes in molar ellipticity at 220 nm over the temperature range of 20-80 • C ( Figure 2D). The transition curve obtained from the thermal-induced denaturation was analyzed using a two-state unfolding model, given above in Equation (1). The analysis gave the change in melting temperature (T m ) and change in enthalpy (∆H m ) at the melting temperature of Prdx6, as shown in Table 1. As seen in Figure 2D and Table 1, the hyperoxidized Prdx6 shows a decrease of 4 • C in its T m as compared to that of reduced Prdx6. This decreased stability shows the conformational change of the native state of Prdx6 upon hyperoxidation, which may lead to changes in its activity under physiological condition. To investigate the actual molecular level differences in the tertiary interaction, we performed molecular dynamic (MD) simulations for both the proteins (see Figure 3). The initial 3D structure models used while initiating simulation were run at Prdx6 homology models with their peroxidatic cysteine either reduced (Cys47-SH) or oxidized (Cys47-SO 3 H). The structure validation for the modeled rPrdx6 protein was done via PROCHECK (Ramachandran plot generator, Figure S1) and QMEAN ( Figure S2) in-silico tools demonstrated the predicted model had good correlation with the experimental structures and were acceptable for further use. In MD simulation, root mean square deviation (RMSD) changes in Cα atoms of the reduced Prdx6 showed a slow increment during the simulation and finally converged at the end of the simulation around 2.5 Å. At the same time, hyperoxidized Prdx6 showed fluctuations at the start, reached 3 Å by 200 ns, and ultimately converging at 2.5 Å, indicating more fluctuation in oxidized Prdx6 compared to reduced Prdx6 ( Figure 3A). The root mean square fluctuation (RMSF) analysis is an important criterion to measure the stability of the protein under study. As seen in Figure 3B, the hyperoxidized Prdx6 displayed more fluctuations in 4 regions (in comparison to reduced Prdx6) localized in the conserved thioredoxin fold. It must be noted that apart from these 4 peaks, the fluctuations showed by side chains of other amino acids were almost similar, and highly flexible regions seems to mostly lie in the loops. In terms of conformational analysis, the radius of gyration is described as the moment of inertia of the group of atoms from their center of mass. As shown in Figure 3C, we observed, that the hyperoxidized Prdx6 showed a higher gyration radius than that of the reduced Prdx6 throughout the simulation run. Further details on the structural packing of the proteins were assessed by measuring SASA. SASA provides insight about the compactness of the hydrophobic core, which is an important factor for determining protein stability. As illustrated in Figure 3D, the SASA profile of the hyperoxidized Prdx6 simulation run exhibited a higher SASA value than that of reduced Prdx6, which could explain the differences in polarity of the local environment of tryptophan residues in both forms of the protein. However, it seems that the hyperoxidation of Prdx6 did not affect the overall dynamics or stability of the protein, indeed, it somehow disrupted the packing of the hydrophobic core, leading to solvent exposition of hydrophobic amino acids as evidenced by the higher gyration radii, SASA profile, and ANS fluorescence of hyperoxidized Prdx6 as compared to reduced Prdx6.  Secondary structure elements, such as alpha-helices and beta-sheets were also monitored throughout the simulation period (see Figure 4). There are 23.11% helix and 19.17% beta-sheets in reduced Prdx6, while the hyperoxidized Prdx6 has 21.74% helix and 18.09% beta-strands, thereby, the total percentage difference in secondary structure elements in both proteins is only 2%. To assess the contribution of non-covalent interaction on the flexibility of Prdx6, hydrogen bond formation during the simulation period was monitored. For each trajectory frame, total hydrogen bonds present in the protein were calculated. The hydrogen bond analysis revealed a decrease in the total number of hydrogen bonds in hyperoxidized Prdx6 ( Figure 5) as compared to reduced Prdx6 during the simulation period. The average hydrogen bonds throughout the simulation for hyperoxidized Prdx6 was calculated to be 123.515 while for the reduced protein it was 126.823. However, since this loss does not seem to affect the internal dynamics of the protein (as indicated by RMSD and RMSF plots), the reason for up-regulation of aiPLA2 activity in hyperoxidized Prdx6 remains ambiguous. Secondary structure elements, such as alpha-helices and beta-sheets were also monitored throughout the simulation period (see Figure 4). There are 23.11% helix and 19.17% beta-sheets in reduced Prdx6, while the hyperoxidized Prdx6 has 21.74% helix and 18.09% beta-strands, thereby, the total percentage difference in secondary structure elements in both proteins is only 2%. To assess the contribution of non-covalent interaction on the flexibility of Prdx6, hydrogen bond formation during the simulation period was monitored. For each trajectory frame, total hydrogen bonds present in the protein were calculated. The hydrogen bond analysis revealed a decrease in the total number of hydrogen bonds in hyperoxidized Prdx6 ( Figure 5) as compared to reduced Prdx6 during the simulation period. The average hydrogen bonds throughout the simulation for hyperoxidized Prdx6 was calculated to be 123.515 while for the reduced protein it was 126.823. However, since this loss does not seem to affect the internal dynamics of the protein (as indicated by RMSD and RMSF plots), the reason for up-regulation of aiPLA 2 activity in hyperoxidized Prdx6 remains ambiguous.
assess the contribution of non-covalent interaction on the flexibility of Prdx6, hydrogen bond formation during the simulation period was monitored. For each trajectory frame, total hydrogen bonds present in the protein were calculated. The hydrogen bond analysis revealed a decrease in the total number of hydrogen bonds in hyperoxidized Prdx6 ( Figure 5) as compared to reduced Prdx6 during the simulation period. The average hydrogen bonds throughout the simulation for hyperoxidized Prdx6 was calculated to be 123.515 while for the reduced protein it was 126.823. However, since this loss does not seem to affect the internal dynamics of the protein (as indicated by RMSD and RMSF plots), the reason for up-regulation of aiPLA2 activity in hyperoxidized Prdx6 remains ambiguous.   Prdxs are inherently oligomeric proteins whose multimeric status is known to be influenced by the oxidation states of the peroxidatic cysteine. For instance, hyperoxidation of peroxidatic cysteine to a Cys-SO2H and Cys-SO3H under oxidative stress has been reported to alter the quaternary structure of a typical 2-Cys Prdx from a dimer into higher oligomers, thus, forming molecular complexes that led to inactivation of its peroxidase activity [41,42]. In fact, in a homolog of Prdx, this structural change from low molecular weight multimer to a high molecular weight oligomeric complex, has been observed to cause a functional change from peroxidase to molecular chaperone [43]. In addition, subunit associations are often preceded by conformational transitions that cause the exposition of hydrophobic amino acids on the surface for stronger interactions at the oligomeric interface. Taking all this into consideration, it is highly probable that hyperoxidation of Prdx6 causes a similar effect (as that of other Prdxs) on the quaternary structure of the enzyme. To analyze the native state oligomeric composition of Prdx6, the native-PAGE of reduced and hyperoxidized Prdx6 was performed (see Figure 6). It must be noted here that, the crystallographic studies as well as equilibrium sedimentation and DuoLink analyses have shown Prdx6 in reduced (Cys47-SH) and sulfinic forms (Cys47-SO2H) to exist as dimers [44][45][46]. In native PAGE gel we also observed the dimeric nature of reduced Prdx6, as indicated by the protein band (monomer molecular weight, Mr: 25 KDa), which is close to that of serum albumin (Mr: 64kDa). However, the hyperoxidized Prdx6 showed a thick band similar to that of horse ferritin, having Mr of 450 kDa. Prdxs are inherently oligomeric proteins whose multimeric status is known to be influenced by the oxidation states of the peroxidatic cysteine. For instance, hyperoxidation of peroxidatic cysteine to a Cys-SO 2 H and Cys-SO 3 H under oxidative stress has been reported to alter the quaternary structure of a typical 2-Cys Prdx from a dimer into higher oligomers, thus, forming molecular complexes that led to inactivation of its peroxidase activity [41,42]. In fact, in a homolog of Prdx, this structural change from low molecular weight multimer to a high molecular weight oligomeric complex, has been observed to cause a functional change from peroxidase to molecular chaperone [43]. In addition, subunit associations are often preceded by conformational transitions that cause the exposition of hydrophobic amino acids on the surface for stronger interactions at the oligomeric interface. Taking all this into consideration, it is highly probable that hyperoxidation of Prdx6 causes a similar effect (as that of other Prdxs) on the quaternary structure of the enzyme. To analyze the native state oligomeric composition of Prdx6, the native-PAGE of reduced and hyperoxidized Prdx6 was performed (see Figure 6). It must be noted here that, the crystallographic studies as well as equilibrium sedimentation and DuoLink analyses have shown Prdx6 in reduced (Cys47-SH) and sulfinic forms (Cys47-SO 2 H) to exist as dimers [44][45][46]. In native PAGE gel we also observed the dimeric nature of reduced Prdx6, as indicated by the protein band (monomer molecular weight, M r : 25 KDa), which is close to that of serum albumin (M r : 64kDa). However, the hyperoxidized Prdx6 showed a thick band similar to that of horse ferritin, having M r of 450 kDa. This varied mobility of both proteins on a native PAGE gel, indicates that hyperoxidized Prdx6 exists in a higher oligomer state. These observations were confirmed by the dynamic light scattering (DLS) measurements at physiological pH 7.4 (See Table 1). The reduced Prdx6 was found to have a hydrodynamic radius of 5.36 nm (corresponding to that of a protein with Mr 50 KDa) while the hyperoxidized Prdx6 had a hydrodynamic radius of 36.19 nm.; the apparent molecular weight calculated from the hydrodynamic radii of reduced and hyperoxidized Prdx6 are shown in Table 1. Thus, we concluded that the quaternary state of reduced and hyperoxidized Prdx6 is different. This change in the quaternary structure of the proteins is preceded by the conformational change shown by the above biophysical studies (CD and fluorescence) for its subunit association during oligomerization [43,47]. The exposure of hydrophobic residues in hyperoxidized Prdx6 (detected by ANS fluorescence) is necessary for non-covalent interaction among subunits in high oligomers. However, the lesser change in the internal dynamics of the protein (as indicated by RMSD and RMSF plots) might be because MD simulation of hyperoxidized Prdx6 was performed in its dimeric state.
During the redox cycle, the hydrophobic residues (Leu 145, Ile 147, Leu 148 and Tyr 149) between β-strands (βb7 from each monomer) in Prdx6 play a significant role on the stabilization of Prdx6-Prdx6 homodimerization and Prdx6-πGST heterodimerization (an important step on its peroxidase activity) [44]. Our finding that the formation of higher oligomer of Prdx6 on the hyperoxidation may not allow the formation of Prdx6-πGST hetero-dimer to recycle the hyperoxidized Prdx6, leading to loss of its peroxidase activity but it might overexpose the catalytic triad for aiPLA2 activity located at the protein surface, inducing the upregulation of its phosphopholipase activity, which has previously been shown in [19]. Therefore, the results lead us to propose that the change in the oligomeric state of Prdx6 due to hyperoxidation might induce the switching of its peroxidase activity to aiPLA2 activity.

Conclusions
In summary, our study demonstrates that post-translational modifications (like hyperoxidation) cause the alteration of reduced Prdx6 dimer to a multimeric status. The formation of a high order oligomer is consequently due to a decrease in protein stability as a result of the opening of some of the hydrophobic groups to the polar solvent. In future, it will be important to further identify the properties of this oligomeric protein and the physiological significance thereof. This varied mobility of both proteins on a native PAGE gel, indicates that hyperoxidized Prdx6 exists in a higher oligomer state. These observations were confirmed by the dynamic light scattering (DLS) measurements at physiological pH 7.4 (See Table 1). The reduced Prdx6 was found to have a hydrodynamic radius of 5.36 nm (corresponding to that of a protein with M r 50 KDa) while the hyperoxidized Prdx6 had a hydrodynamic radius of 36.19 nm.; the apparent molecular weight calculated from the hydrodynamic radii of reduced and hyperoxidized Prdx6 are shown in Table 1. Thus, we concluded that the quaternary state of reduced and hyperoxidized Prdx6 is different. This change in the quaternary structure of the proteins is preceded by the conformational change shown by the above biophysical studies (CD and fluorescence) for its subunit association during oligomerization [43,47]. The exposure of hydrophobic residues in hyperoxidized Prdx6 (detected by ANS fluorescence) is necessary for non-covalent interaction among subunits in high oligomers. However, the lesser change in the internal dynamics of the protein (as indicated by RMSD and RMSF plots) might be because MD simulation of hyperoxidized Prdx6 was performed in its dimeric state.
During the redox cycle, the hydrophobic residues (Leu 145, Ile 147, Leu 148 and Tyr 149) between β-strands (βb7 from each monomer) in Prdx6 play a significant role on the stabilization of Prdx6-Prdx6 homodimerization and Prdx6-πGST heterodimerization (an important step on its peroxidase activity) [44]. Our finding that the formation of higher oligomer of Prdx6 on the hyperoxidation may not allow the formation of Prdx6-πGST hetero-dimer to recycle the hyperoxidized Prdx6, leading to loss of its peroxidase activity but it might overexpose the catalytic triad for aiPLA2 activity located at the protein surface, inducing the upregulation of its phosphopholipase activity, which has previously been shown in [19]. Therefore, the results lead us to propose that the change in the oligomeric state of Prdx6 due to hyperoxidation might induce the switching of its peroxidase activity to aiPLA2 activity.

Conclusions
In summary, our study demonstrates that post-translational modifications (like hyperoxidation) cause the alteration of reduced Prdx6 dimer to a multimeric status. The formation of a high order oligomer is consequently due to a decrease in protein stability as a result of the opening of some of the hydrophobic groups to the polar solvent. In future, it will be important to further identify the properties of this oligomeric protein and the physiological significance thereof.