key: cord-0012607-bbb5myc3 authors: Protto, Virginia; Tramutola, Antonella; Fabiani, Marco; Marcocci, Maria Elena; Napoletani, Giorgia; Iavarone, Federica; Vincenzoni, Federica; Castagnola, Massimo; Perluigi, Marzia; Di Domenico, Fabio; De Chiara, Giovanna; Palamara, Anna Teresa title: Multiple Herpes Simplex Virus-1 (HSV-1) Reactivations Induce Protein Oxidative Damage in Mouse Brain: Novel Mechanisms for Alzheimer’s Disease Progression date: 2020-06-29 journal: Microorganisms DOI: 10.3390/microorganisms8070972 sha: cbf2a8ce743bc19fbf8b89bccaec836b25897cfe doc_id: 12607 cord_uid: bbb5myc3 Compelling evidence supports the role of oxidative stress in Alzheimer’s disease (AD) pathophysiology. Interestingly, Herpes simplex virus-1 (HSV-1), a neurotropic virus that establishes a lifelong latent infection in the trigeminal ganglion followed by periodic reactivations, has been reportedly linked both to AD and to oxidative stress conditions. Herein, we analyzed, through biochemical and redox proteomic approaches, the mouse model of recurrent HSV-1 infection we previously set up, to investigate whether multiple virus reactivations induced oxidative stress in the mouse brain and affected protein function and related intracellular pathways. Following multiple HSV-1 reactivations, we found in mouse brains increased levels of oxidative stress hallmarks, including 4-hydroxynonenal (HNE), and 13 HNE-modified proteins whose levels were found significantly altered in the cortex of HSV-1-infected mice compared to controls. We focused on two proteins previously linked to AD pathogenesis, i.e., glucose-regulated protein 78 (GRP78) and collapsin response-mediated protein 2 (CRMP2), which are involved in the unfolded protein response (UPR) and in microtubule stabilization, respectively. We found that recurrent HSV-1 infection disables GRP78 function and activates the UPR, whereas it prevents CRMP2 function in mouse brains. Overall, these data suggest that repeated HSV-1 reactivation into the brain may contribute to neurodegeneration also through oxidative damage. protein 78 (GRP78), a protein involved in the unfolded protein response (UPR) [30] , and collapsin response-mediated protein 2 (CRMP2), a protein able to bind and stabilize microtubules [31] . We found that recurrent HSV-1 infection disables GRP78 function and activates UPR, whereas it prevents CRMP2 function in mouse brains. The authors certify that all experimental protocols used in the present study are in agreement with the Italian and European legislation on animal experimentation (Decreto L.vo 26/ 2014; Direttiva n. 63/2010/EU). Experimental protocols were authorized by the Italian Ministry of Health (protocol numbers 801/2016-B, 745/2016-PR). We used brain tissue from the mouse model of recurrent HSV-1 infection we recently set up and characterized [25] . Briefly, 2 µl of HSV-1 suspension (equivalent to 10 6 particle-forming units, pfu) or mock suspension were inoculated in 6-8-week-old female BALB/c mice (purchased from Harlan Laboratories) by lip scarification, following anesthesia with intraperitoneally injection of ketamine (80mg/Kg) + xylazine (5mg/Kg). HSV-1-infected mice (HSV1-M) and mock-infected ones (CTRL-M) were marked and kept in separate cages. When latent infection was established (6 weeks postinfection, p.i.) mice were individually subjected to hyperthermia (40-42 °C for 15 min: thermal stress = TS) to reactivate the latent virus, according to protocols from Sawtell et al [32] . TS was repeated at 6-8-week intervals for 7 times during mouse life. Following the 7th TS, a set of animals per group was sacrificed and analyzed as described below (see experimental schedule in Figure 1 ). In this study, we used brain tissues from the mice assessed for AD hallmarks in De Chiara et al. 2019 [25] . For further details about the mouse model, see De Chiara et al. (2019) [25] . Six weeks after the infection, mice underwent several cycles of thermal stress (TS) that were spaced 6-8 weeks from each other. After the 7th TS, the mice were sacrificed, and brain tissues were collected for biochemical analyses (redox proteomics, Western blot, WB, and immunoprecipitation, IP) or perfused with paraformaldehyde (PFA) for confocal immunofluorescence analysis (IF). Brain tissues (frontal cortex) from CTRL-M and HSV1-M (n = 6 per group) were homogenized by sonication on ice with RIPA buffer (20 mM Tris, 150 mM NaCl, 1% Triton X-100, 1% sodium deoxycholate, 0.1% SDS) supplied with protease and phosphatase inhibitor cocktails (Sigma-Aldrich, St Louis, MO, USA) and centrifuged at 15,000 g at 4 °C for 20 min to remove cellular debris. The supernatant was extracted to determine the total protein concentration by the the BCA method Figure 1 . Experimental design. We infected two-month-old BALB/c mice with Herpes simplex virus-1 (HSV-1, HSV1-M) or a mock solution (CTRL-M). Six weeks after the infection, mice underwent several cycles of thermal stress (TS) that were spaced 6-8 weeks from each other. After the 7th TS, the mice were sacrificed, and brain tissues were collected for biochemical analyses (redox proteomics, Western blot, WB, and immunoprecipitation, IP) or perfused with paraformaldehyde (PFA) for confocal immunofluorescence analysis (IF). Brain tissues (frontal cortex) from CTRL-M and HSV1-M (n = 6 per group) were homogenized by sonication on ice with RIPA buffer (20 mM Tris, 150 mM NaCl, 1% Triton X-100, 1% sodium deoxycholate, 0.1% SDS) supplied with protease and phosphatase inhibitor cocktails (Sigma-Aldrich, St Louis, MO, USA) and centrifuged at 15,000 g at 4 • C for 20 min to remove cellular debris. The supernatant was extracted to determine the total protein concentration by the the BCA method (Pierce, Rockford, IL, USA). Equal amounts of proteins (20 µg) were separated on SDS-PAGE using the TGX precast gels technology to obtain a Total Load gel used for the normalization of specific protein signals or using homemade 10% polyacrylamide gels. Subsequently, the gels were blotted onto nitrocellulose membranes (Amersham Protran 0.45 µm, Hybond ECL, GE Healthcare 0.22 µm, Chicago, IL, USA). The membranes were then blocked for 1-3 h in 10% non-fat milk or 3% bovine serum albumin (BSA, SERVA Electrophoresis GmbH, Heidelberg, Germany), as requested by the specific primary antibodies, in 0.1% tween-100 TBS (T-TBS) and incubated overnight at 4 • C with suitable dilutions of specific primary antibodies (in 5% non-fat milk or 3% BSA Secondary antibodies horseradish peroxidase-conjugated anti-mouse, anti-rabbit, or anti-goat IgG (1:10,000, Bio-Rad, Hercules, CA, USA) were incubated for 1 h at room temperature (RT). The membranes were developed with Clarity enhanced chemiluminescence (ECL) substrate (Bio-Rad, Hercules, CA, USA, #1705061), acquired using Chemi-Doc MP System (Bio-Rad, Hercules, CA, USA), and analyzed with Image Lab software (Bio-Rad, Hercules, CA, USA). For carbonylated protein detection, after blocking, the membranes were equilibrated in 20% methanol, incubated in 2 N HCl, and finally derivatized with 0.5 mM 2,4-dinitrophenylhydrazone (DNP). Finally, DNP-protein adducts were detected on nitrocellulose using a primary DNP-specific antibody (MerckMillipore, Merck KGaA, Darmstadt, Germany). Following 3 washes in T-TBS, the membranes were incubated with horseradish peroxidase-conjugated antibodies (Jackson ImmunoResearch Laboratories, West Grove, PA, USA, 1:5000 in T-TBS). Brain samples proteins (frontal cortex, 150 µg) from CTRL-M and HSV1-M (n = 6 per group) were isolated by trichloroacetic acid (TCA, 15% final volume) in ice for 10 min. Subsequently, the samples were centrifuged for 5 min at 10,000 g and the resulting pellets were rinsed three times in ice-cold ethanol/ethyl acetate 1:1 solution. The pellets were then resuspended in 200 µl of rehydration buffer (8 M urea, 20 mM dithiothreitol (DTT), 2.0% (w/v) Chaps, 0.2% Bio-Lyte, 2M thiourea, and bromophenol blue). Isoelectrofocusing (first dimension) was performed at RT using ReadyStrip™ IPG (11 cm, pH 3-10; Bio-Rad, Hercules, CA, USA) at 300 V for 2 h linearly, 500 V for 2 h linearly, 1000 V for 2 h linearly, 8000 V for 8 h linearly, and 8000 V for 10 h rapidly. Then, the isolectrofocusing strips were equilibrated in 50 mM Tris-HCl (pH 6.8) containing 6 M urea, 1% (w/v) sodium dodecyl sulfate (SDS), 30% (v/v) glycerol, and 0.5% DTT for 10 min and again for another 10 min in 50 mM Tris-HCl (pH 6.8) containing 6 M urea, 1% (w/v) sodium dodecyl sulfate (SDS), 30% (v/v) glycerol, and 4.5% iodoacetamide (IA). The electrophoresis step (second dimension) was performed using 12% acrylamide precast Criterion gels. The gels were then placed in fixing solution (7% acetic acid, 10% methanol) and incubated for 45 min. Afterward, the gels were stained for 2 h with SYPRO Ruby Gel Stain (Bio-Rad, Hercules, CA, USA) and finally de-stained overnight in deionized water. 2D gels were then blotted onto nitrocellulose membranes (Bio-Rad, Hercules, CA, USA) and after 1 h incubation with 3% albumin in T-TBS to block unspecific binding sites, HNE-protein adducts were detected by 2 h staining at RT with a primary anti-HNE antibody (1:1000, NB-10063093 Novus Biologicals, Centennial, CO 80112, USA), followed by three washes with T-TBS and further incubation at RT for 1 h with a secondary alkaline phosphatase-conjugated anti-goat antibody (1:5000, Sigma-Aldrich, St. Louis, MO, USA). After extensively washes with T-TBS, the membranes were developed using a 5-bromo-4-chloro-3-indolyl phosphate/nitroblue tetrazolium solution (BCIP/NBT) prepared in alkaline phosphatase buffer solution. Redox proteomics gels and blots images were acquired by the Chemi-Doc MP System (Bio-Rad, Hercules, CA, USA) in TIFF format. Gels and blots (n = 12 each) were evaluated by PD-Quest 2D Analysis (7.2.0 version; Bio-Rad, Hercules, CA, USA). In brief, a master gel was chosen to identify spots on gels and blots, then spots from gels and blots were matched for normalization according to the total spot density. Gel-to-blot analysis was initially performed by manual matching of common spots and then, after obtaining a significant number of spots, by automated matching. Automated matching is created on user-defined parameters for spot detection based on the faintest spot, the largest spot, and the largest spot cluster that appear in the master gel. This process generates approximately 400-800 spots per gel or blot. Only proteins showing computer-determined significant (p < 0.05) differential levels between the groups (CTRL-M and HSV1-M) were considered for MS analysis. To define significant differences in proteins levels, analysis set manager software built in the PD-Quest software was used. The number of pixels that appear in a protein spot were calculated by the software, consistent with an increase or a decrease in protein levels. The image analysis was conducted in parallel on blots and on expression gels, and the results were then compared by the software to normalize HNE modification to expression value for each spot detected and matched. Data obtained by PD-QUEST software were analyzed using Student's t-test. Significance was accepted if the p value was <0.05. Significantly different protein spots identified by the comparison between CTRL-M and HSV1-M were removed from the gels and individually transferred into microcentrifuge tubes for trypsin digestion. By using DTT and IA, disulphide bonds were reduced and cap and then gel plugs were incubated overnight at 37 • C with shaking in a trypsin solution. Tryptic peptides solutions were resuspended in water and stored at −80 • C for MS/MS analysis. High-resolution HPLC-ESI-MS/MS experiments were performed using an Ultimate 3000 RSLC nano system coupled to an LTQ Orbitrap ELITE apparatus (Thermo Fisher Scientific, Waltham, MA, USA). Zorbax 300 SB-C18 (3. 50 µL/min for a total of 65 min. MS spectra were collected with 120,000 resolution and m/z range from 350 to 2000. In data-dependent acquisition mode, the five most intense multiply charged ions were selected and fragmented in ion trap by using CID 35% normalized collision energy. Tuning parameters were: capillary temperature 300 • C, source voltage 4.0 kV. MS/MS data were analyzed by Proteome Discoverer software (version 1.4.1.14, Thermo Fisher Scientific, Waltham, MA, USA), based on SEQUEST HT cluster as search engine against UniProtKB mouse database (released on 28 of February 2017, Mus musculus 16,839 entries). Search parameters were: 10 ppm tolerance for precursor ions and 0.5 Da for product ions; 2 missed cleavage. Carbamydomethylation of cysteine (+57.02 Da) and oxidation of methionine (+15.99Da) were taken into consideration as fixed and variable modification, respectively. Protein characterization was set with the identification of a minimum of two peptides per protein and two unique peptides by applying the high-confidence filter. To validate the results obtained by MS, the identity of two specific proteins, GRP78-and CRMP2-HNE adducts, was evaluated using immunochemical methods. Samples (500 µg of proteins) were incubated overnight with specific antibodies (anti-GRP78 SC-13539, Santa Cruz Biotechnology, Inc., Heidelberg, Germany or anti-CRMP2 #9393, Cell Signaling Technology, Danvers, MA, USA) diluted in IP buffer (10 mM Tris, pH 7.6, 140 mM NaCl, 0.5% NP40 including protease inhibitors). The neoformed immunocomplexes were then incubated 2 h with Protein G beads (Sigma-Aldrich, St Louis, MO, USA) and then rinsed three times with RIA buffer (10 mM Tris, pH 7.6; 140 mM NaCl; 1% NP40). Proteins were separated by SDS-PAGE, and then the gels were immunoblotted onto nitrocellulose membranes (Bio-Rad, Hercules CA, USA). The membranes were incubated with the anti-HNE primary antibody followed by 1 h incubation with a peroxidase-conjugated secondary antibody (ImmunStar 1:10000, Bio-Rad, Hercules, CA, USA) or peroxidase-conjugated-Protein G (1:5000 Sigma-Aldrich, St Louis, MO, USA #18-161). Chemiluminescence was detected with Clarity enhanced chemiluminescence (ECL) substrate (Bio-Rad, Hercules, CA, USA, #1705061), and membranes acquired with Chemi-Doc MP image system (Bio-Rad, Hercules, CA, USA) and analyzed using Image Lab software (Bio-Rad, Hercules, CA, USA). The IP results were normalized on the total amount of the proteins of interest. Mouse perfusion was performed as described [25] , and brains were post-fixed for 24 h in 4% paraformaldehyde (PFA), followed by cryoprotection for 24 h in 30% sucrose in PBS. Coronal brain sections (40-µm-thick) were obtained using a vibratome (VT1000 S, Leica Microsystems, Wetzlar, Germany). Free-floating sections were permeabilized 15 min at RT with 0.5% Triton X-100 in PBS, blocked 2h at RT with 10% horse serum (HS), 0.2% Triton X-100 in PBS, and then incubated overnight at 4 • C with specific primary antibodies (anti-GRP78 SC-1051 Santa Cruz Biotechnology, Heidelberg, Germany; anti-CRMP2 9393 Cell Signaling Technology, Danvers, MA, USA) diluted in antibody buffer (5% HS, 0.2% Triton X-100 in PBS). Sections were rinsed with PBS and then incubated 90 min at RT with secondary antibodies coupled to Alexa Fluor ® 488 or Alexa Fluor ® 546 dyes (Thermo Fisher Scientific, Waltham, MA, USA) diluted 1:1000 in antibody buffer. After extensively washing with PBS, the sections were incubated 20 min at RT with 4 ,6 -diamidine-2 -phenylindole dihydrochloride (DAPI, Thermo Fisher Scientific, Waltham, MA, USA) and finally mounted on coverslips with Fluoromount G (Thermo Fisher Scientific, Waltham, MA, USA). The stained sections were imaged with a confocal laser-scanning microscope (Leica SP5, Leica Microsystems, Wetzlar, Germany) under sequential mode, to avoid crosstalk between channels. Images were acquired at 20X or 63X magnification. Image analysis was performed by the Imaris Suite 7.4 software (Bitplane A.G., Zurich, Switzerland). Quantification of immunofluorescence was carried out by drawing regions of interest (ROIs) in specific brain areas of the slices (cortex) and quantifying the mean fluorescence intensity in the ROIs. For the production of figures, processing was done by using the Adobe Photoshop CS6 software (Adobe Systems Incorporated, San Josè, CA, USA). Statistical comparisons were performed with GraphPad 6.0 (Prism) by using Student's t-test. Data are presented as mean ± standard error of the mean (SEM). The level of significance was set at 0.05. All the experiments were performed by operators blind to the study conditions. We previously demonstrated that multiple HSV-1 reactivations reaching the brain induced hallmarks of neurodegeneration (i.e., accumulation of Aβ and phospho-tau and neuroinflammation) as well as cognitive deficits in mice [25] , thus resembling the occurrence of an AD-like phenotype. Here, we investigated whether recurrent HSV-1 infections in mouse brain may induce also a progressive increment of oxidative stress conditions, as evidenced during the development of AD [33] . Specifically, we decided to investigate in this study the presence of oxidative damage in brain tissue of mice showing the most significant accumulation of AD hallmarks [25] , which are a subset of mice that underwent seven cycles of TS-induced virus reactivation over their life and were sacrificed at 13 months of age (named hereafter HSV1-M, see experimental schedule in Figure 1 ). Hence, we first investigated the levels of oxidative stress markers in the brains of HSV1-M and matched mock-infected mice (named hereafter CTRL-M). In particular, we evaluated by western blot (WB) the level of oxidative modification of proteins indexed by the levels of HNE, 3-NT, and CP in cortical lysates from six HSV1-M and six CTRL-M. Results in Figure 2A and in Supplementary Figure S1 show that the levels of HNE (p = 0.01), 3-NT (p = 0.03), and CP (p = 0.03) were significantly increased in HSV1-M with respect to CTRL-M, indicating that recurrent HSV-1 infection triggers oxidative damage to proteins in mouse brains. We then employed redox proteomic analysis focusing on HNE modifications, since HNE levels and the levels of HNE-protein adducts are known to be increased in AD brains [6, 34, 35] . Specifically, we used a redox proteomic approach to identify protein targets of HNE peroxidation in cortices from the mice previously analyzed by WB for the levels of oxidatively modified proteins (see Figure 2 ). By comparing the densitometric intensities of individual spots of each HSV1-M sample with those measured in CTRL-M samples, we identified 12 HNE-modified spots whose levels were modulated (both upwards and downwards) in the cortex of HSV-1-M (see highlighted and numbered spots in Figure 2C ). These spots were further analyzed by HPLC-ESI-MS/MS to identify the corresponding proteins. Specifically, in HSV1-M cortices, we found increased levels with respect to controls, ranging from 1.5-to 10-fold, of the HNE-protein adducts for dihydropteridine reductase (Dhpr), phosphoglycerate mutase 1 (Pgam1), translationally controlled tumor protein (Tpt1), occludin (Ocln), SEC23-interacting protein (Sec23ip), U7 snRNA-associated Sm-like protein LSm11 (Lsm11), osocitrate dehydrogenase [NAD] subunit alpha (Idh3a), ATP synthase subunit alpha (Atp5a1), glucose-regulated protein 78 (Hspa5, GRP78), and apoptosis-resistant E3 ubiquitin protein ligase 1 (Arel) ( Table 1) . Decreased levels, between 10-to 45-fold, were found for dihydropyrimidinase-related protein 2 (Dpysl2, CRMP2), myosin-10 (Myh10), and guanine nucleotide-binding protein G(I)/G(S)/G(T) subunit beta-1 (Gnb1). We then employed redox proteomic analysis focusing on HNE modifications, since HNE levels and the levels of HNE-protein adducts are known to be increased in AD brains [6, 34, 35] . Specifically, we used a redox proteomic approach to identify protein targets of HNE peroxidation in cortices from the mice previously analyzed by WB for the levels of oxidatively modified proteins (see Figure 2 ). By comparing the densitometric intensities of individual spots of each HSV1-M sample with those measured in CTRL-M samples, we identified 12 HNE-modified spots whose levels were modulated (both upwards and downwards) in the cortex of HSV-1-M (see highlighted and numbered spots in Figure 2C ). The bioinformatics analysis of protein interaction networks and of biological processes for HNE-modified proteins in the comparison between HSV1-M and CTRL-M brains was built using the STRING software 11.0 [36, 37] . This analysis allowed us to detect the biological repercussions of the aberrant HNE modification of the MS/MS-identified proteins. The analysis of the networks produced by the aberrantly oxidized proteins in the comparison between HSV1-M and CTRL-M ( Figure 3A) showed a higher number of edges than expected (5 vs. 4) , implying that the HNE-modified proteins interacted more strongly than what would be expected for a random set of proteins of similar size, drawn from the genome. Interaction networks are displayed in a table ( Figure 3B ) indicating interaction nodes and score and in an image showing shared biological processes ( Figure 3C ) which are fully elucidated in table D of Figure 3 . Intriguingly, the majority of the proteins have catalytic, binding, or hydrolase activities, which are exploited in important cellular processes, such as energy metabolism (Dhpr, Pgam1, Idh3, and Atp5a1), protein folding (GRP78, Sec23ip), degradation processes (Arel1, Lsm11), and cell structure (Dpyl2, Gnb1, Myh10, Ocln, Tpt1), suggesting that their oxidative modification may deeply affect brain physiology, as previously demonstrated in AD [5, 6] . binding, or hydrolase activities, which are exploited in important cellular processes, such as energy metabolism (Dhpr, Pgam1, Idh3, and Atp5a1), protein folding (GRP78, Sec23ip), degradation processes (Arel1, Lsm11), and cell structure (Dpyl2, Gnb1, Myh10, Ocln, Tpt1), suggesting that their oxidative modification may deeply affect brain physiology, as previously demonstrated in AD [5, 6] . We next focused our attention on GRP78 (Hspa5), also known as BiP, that is a key protein in the UPR, a pathway involved in endoplasmic reticulum (ER) stress, and whose alterations are also linked to viral infection [30, 38] and AD pathogenesis [39] . First, we validated the redox proteomic results by checking the levels of GRP78-HNE adducts in mouse cortical homogenates through an immunoprecipitation assay. In particular, by the aid of specific antibodies, we immunoprecipitated GRP78 and revealed by western blot the HNE-modified protein. Representative results in Figure 4A show that the amount of the immunoprecipitated HNE-GRP78 was significantly higher in HSV1-M cortices as compared with those of CTRL-M (p = 0.0175). In addition, WB analysis of tissue homogenates revealed that GRP78 expression was similar in cortices from HSV1-M compared to those of matched CTRL-M ( Figure 4B , full length WB showed in Figure S2 ). Accordingly, IF analysis of brain slices did not reveal any difference in the expression levels of GRP78 between HSV1-M and CTRL-M ( Figure 4C ). We next focused our attention on GRP78 (Hspa5), also known as BiP, that is a key protein in the UPR, a pathway involved in endoplasmic reticulum (ER) stress, and whose alterations are also linked to viral infection [30, 38] and AD pathogenesis [39] . First, we validated the redox proteomic results by checking the levels of GRP78-HNE adducts in mouse cortical homogenates through an immunoprecipitation assay. In particular, by the aid of specific antibodies, we immunoprecipitated GRP78 and revealed by western blot the HNE-modified protein. Representative results in Figure 4A show that the amount of the immunoprecipitated HNE-GRP78 was significantly higher in HSV1-M cortices as compared with those of CTRL-M (p = 0.0175). In addition, WB analysis of tissue homogenates revealed that GRP78 expression was similar in cortices from HSV1-M compared to those of matched CTRL-M ( Figure 4B , full length WB showed in Figure S2 ). Accordingly, IF analysis of brain slices did not reveal any difference in the expression levels of GRP78 between HSV1-M and CTRL-M ( Figure 4C ). Next, since oxidative damage leads to GRP78 dysfunction and consequent aberrant activation of UPR, as previously observed in down syndrome (DS) [40, 41] , we investigated the activation of the three major arms of UPR pathways mediated by PERK, IRE, and ATF6. Specifically, we investigated changes in the levels of phospho-PERK (pT981, p-PERK) and phospho-IRE (pS724, p-IRE) and in the levels of ATF6 expression in cortical homogenates from HSV1-M and CTRL-M. Results in Figure 5 A,B show that PERK activation (i.e., levels of pT981) was significantly increased in HSV1-M with respect to CTRL-M (p = 0.0038), whereas IRE phosphorylation was unaffected (p = 0.504). In addition, the levels of total ATF6 (90 Kd) were found unchanged following recurrent HSV-1 infection, whereas its active cleaved form (50 Kd) was found significantly increased in HSV1-M with respect to CTRL-M (p = 0.02). Full length WBs are shown in Figure S3 . Overall, these data indicate that recurrent HSV-1 infection overactivates UPR in mouse brains, likely through GRP78 oxidation. Next, since oxidative damage leads to GRP78 dysfunction and consequent aberrant activation of UPR, as previously observed in down syndrome (DS) [40, 41] , we investigated the activation of the three major arms of UPR pathways mediated by PERK, IRE, and ATF6. Specifically, we investigated changes in the levels of phospho-PERK (pT981, p-PERK) and phospho-IRE (pS724, p-IRE) and in the levels of ATF6 expression in cortical homogenates from HSV1-M and CTRL-M. Results in Figure 5A ,B show that PERK activation (i.e., levels of pT981) was significantly increased in HSV1-M with respect to CTRL-M (p = 0.0038), whereas IRE phosphorylation was unaffected (p = 0.504). In addition, the levels of total ATF6 (90 Kd) were found unchanged following recurrent HSV-1 infection, whereas its active cleaved form (50 Kd) was found significantly increased in HSV1-M with respect to CTRL-M (p = 0.02). Full length WBs are shown in Figure S3 . Overall, these data indicate that recurrent HSV-1 infection overactivates UPR in mouse brains, likely through GRP78 oxidation. Next, we focused our attention on Dpysl2, also known as CRMP2, since its oxidative modifications in mouse brains seemed to be prevented by multiple HSV-1 reactivation (10-fold decrease vs. CTRL-M, see Figure 4D and Table 1 ). CRMP2 is a ubiquitous protein in the brain and it is involved in microtubule stabilization, with a consequent regulatory role in cytoskeleton dynamics, neurite outgrowth, vesicle trafficking, and synaptic transmission in the developing brain [42] . Its function in the adult brain is not fully understood, but CRMP2 has been associated with several neuropathologic conditions including AD, where it has been found collapsed in neurofibrillary tangles upon hyperphosphorylation on different sites, including Thr-514 [43] . Thus, we firstly validated the redox proteomic data by an immunoprecipitation assay performed with cortical homogenates from n = 3 HSV1-M and n = 3 matched control mice. Results in Figure 6A show a decreasing trend in the amount of the immunoprecipitated HNE-CRMP2 in HSV1-M cortices as compared with CTRL-M ones. Then, we evaluated CRMP2 expression levels in mouse cortices by western blot and immunofluorescence analyses, finding that multiple HSV-1 reactivations did not affect them ( Figure 6B,C) . Finally, we investigated CRMP2 phosphorylation level in T514 (p-CRMP2): we found that it was significantly decreased in HSV1-M cortices as compared to CTRL-M ones (p = 0.0108, Figure 6B ). Full length WBs are shown in Figure S2 . Altogether these data suggest that multiple virus reactivations preserve CRMP2 function in mouse brain. Next, we focused our attention on Dpysl2, also known as CRMP2, since its oxidative modifications in mouse brains seemed to be prevented by multiple HSV-1 reactivation (10-fold decrease vs. CTRL-M, see Figure 4D and Table 1 ). CRMP2 is a ubiquitous protein in the brain and it is involved in microtubule stabilization, with a consequent regulatory role in cytoskeleton dynamics, neurite outgrowth, vesicle trafficking, and synaptic transmission in the developing brain [42] . Its function in the adult brain is not fully understood, but CRMP2 has been associated with several neuropathologic conditions including AD, where it has been found collapsed in neurofibrillary tangles upon hyperphosphorylation on different sites, including Thr-514 [43] . Thus, we firstly validated the redox proteomic data by an immunoprecipitation assay performed with cortical homogenates from n = 3 HSV1-M and n = 3 matched control mice. Results in Figure 6A show a decreasing trend in the amount of the immunoprecipitated HNE-CRMP2 in HSV1-M cortices as compared with CTRL-M ones. Then, we evaluated CRMP2 expression levels in mouse cortices by western blot and immunofluorescence analyses, finding that multiple HSV-1 reactivations did not affect them ( Figure 6B,C) . Finally, we investigated CRMP2 phosphorylation level in T514 (p-CRMP2): we found that it was significantly decreased in HSV1-M cortices as compared to CTRL-M ones (p = 0.0108, Figure 6B ). Full length WBs are shown in Figure S2 . Altogether these data suggest that multiple virus reactivations preserve CRMP2 function in mouse brain. showing the levels of CRMP2 phosphorylation at threonine 514 (p-CRMP2) and total CRMP2 expression in cortical homogenates from CTRL-M and HSV1-M. Tubulin expression level was used as a sample loading control. Densitometric analysis of immunoreactive signals normalized to CRMP2 (for p-CRMP2) or tubulin (for CRMP2) are shown in the graphs. Values represent the normalized fold changes in protein levels from n = 6 HSV-M with respect to n = 6 CTRL-M (mean ± SEM; *p < 0.05). C: Confocal immunofluorescence analysis of coronal brain slices from HSV1-M and CTRL-M (n = 3). Panels show representative neurons of the cortex that were immunostained for CRMP2 (green) and NeuN (red). Cell nuclei were stained with DAPI (blue). Bar graphs show mean CRMP2 fluorescence intensity quantified in the frontal cortex and expressed as fold change with respect to CTRL-M (mean ± SEM). A growing body of evidence supports the role of HSV-1 in AD pathogenesis and progression [44] . We previously demonstrated that recurrent HSV-1 infection in mice triggers the accumulation of the main AD pathognomic hallmarks, including cognitive dysfunction [25] . Since oxidative stress has been implicated in many neurodegenerative disorders, including AD [45] [46] [47] [48] , and oxidative stress conditions are associated with HSV-1 infection in host cells [27, 29] , in the present study, we investigated whether recurrent HSV-1 infection triggers oxidative stress in mouse brain, thus affecting protein function and contributing to neurodegeneration. We revealed that multiple HSV-1 reactivations cause oxidative modification in proteins and lipids (HNE, 3-NT, and CP) in the mouse cortex, resembling what occurs in the progression of AD [49, 50] . In particular, studies from Butterfield's group showed that several proteins in the brain of mild cognitive impairment (MCI) and AD subjects are oxidatively modified compared with controls, and that these irreversible modifications interfere with their function [51, 52] . Indeed, these authors identified oxidatively modified proteins involved in specific cellular pathways, including energy metabolism, axonal A growing body of evidence supports the role of HSV-1 in AD pathogenesis and progression [44] . We previously demonstrated that recurrent HSV-1 infection in mice triggers the accumulation of the main AD pathognomic hallmarks, including cognitive dysfunction [25] . Since oxidative stress has been implicated in many neurodegenerative disorders, including AD [45] [46] [47] [48] , and oxidative stress conditions are associated with HSV-1 infection in host cells [27, 29] , in the present study, we investigated whether recurrent HSV-1 infection triggers oxidative stress in mouse brain, thus affecting protein function and contributing to neurodegeneration. We revealed that multiple HSV-1 reactivations cause oxidative modification in proteins and lipids (HNE, 3-NT, and CP) in the mouse cortex, resembling what occurs in the progression of AD [49, 50] . In particular, studies from Butterfield's group showed that several proteins in the brain of mild cognitive impairment (MCI) and AD subjects are oxidatively modified compared with controls, and that these irreversible modifications interfere with their function [51, 52] . Indeed, these authors identified oxidatively modified proteins involved in specific cellular pathways, including energy metabolism, axonal integrity, chaperone machinery, and antioxidant systems, that correlate with the pathology of MCI and AD [52, 53] . Along this line, our redox proteomic approach allowed us to identify 13 HNE adduct proteins profoundly modulated by recurrent HSV-1 infection (Table 1) : the majority of these proteins are involved in cellular processes such as energy metabolism (Dhpr, Pgam1, Idh3, and Atp5a1), protein folding (GRP78, Sec23ip), degradation process (Arel1, Lsm11), and cell structure (Dpyl2, Gnb1, Myh10, Ocln, Tpt1). These pathways have been linked to AD onset and progression [54] . Hence, our data support the hypothesis that oxidative stress conditions induced by HSV-1 reactivation can concur to the virus-induced neurodegeneration. Among the identified HNE-protein adducts that we found increased in HSV1-M cortices, we focused our attention on GRP78 (Hspa5: 5.6-fold oxidation, see Figure 2D and Table 1 ) that is a master protein in UPR response, a pathway involved in ER stress. In particular, GRP78 is a molecular chaperone member of the heat shock protein 70 (Hsp70) family of proteins, which transiently binds proteins into the ER to facilitate their folding, assembly, and transport [55] . Normally, GRP78 binds the intracellular domains of the three main UPR stress sensors ATF6, PERK, and IRE1, thus keeping them in an inactive state. Under ER stress, GRP78 acts as a chaperone to assist the folding of proteins accumulated in the ER and thus releases the UPR stress sensors. These, in turn, following phosphorylative events (e.g., for PERK and IRE1) as well as proteolytic ones (e.g., for ATF6), activate intracellular pathways to restore protein homeostasis or, when the stress persists, to induce apoptotic cell death. Specifically, activated PERK (p-PERK) phosphorylates the alpha subunit of eukaryotic translation initiation factor (eIF2α), which in turn inhibits protein translation to reduce ER protein burden; activated ATF6 (i.e., the cleaved form) translocates to the nucleus to upregulate the transcription of chaperone proteins, including GRP78; activated IRE1 acts as an endoribonuclease to upregulate genes responsible for protein folding and ER-associated degradation (ERAD) [56] . Intriguingly, among these three arms of the UPR, PERK is the most commonly linked to neurodegenerative disease, and increased p-PERK and p-eIF2α associated with pathology hallmarks have been found in AD brain tissues [39, 57, 58] . Indeed, ER stress and UPR are reportedly associated with AD pathogenesis, due to the accumulation of misfolded proteins in the brain, including Aβ and phospho-tau [39, 57, 59, 60] , and their chronic activation has been recently suggested as a key feature of this disease [61] . Interestingly, some papers showed that GRP78 is involved in Aβ clearance [62] [63] [64] [65] . For instance, GRP78 was demonstrated to interact with APP, reducing Aβ40 and Aβ42 secretion. This may influence APP metabolism by impairing its access to βand γ-secretases within the ER/Golgi or by facilitating its correct folding [55] . On the other hand, Aβ oligomers have been shown to induce cell death by eliciting ER stress, also though GRP78 overexpression [66] , suggesting the existence of a complex scenario where ER stress and UPR may play a key role in neurodegeneration. Similarly, controversial data associated GRP78 expression/activation with tau phosphorylation. Some in vitro data reported that overexpressed GRP78 is able to bind GSK-3β, enhancing its kinase activity and thus causing an increase in tau phosphorylation [67, 68] . On the contrary, other studies linked increased levels of tau phosphorylation and oligomer formation to aberrant UPR activation and GRP78 downregulation [48, 69] . Di Domenico and colleagues demonstrated that oxidative damage, including that caused by lipid peroxidation, may affect GRP78 capability to bind misfolded proteins, providing a plausible link between deficit of molecular chaperones clearance activity, accumulation of misfolded proteins, risk of cognitive decline, and neurodegeneration [70] . Accordingly, our results, showing an accumulation of HNE-modified GRP78 in those HSV-1-M that displayed the main accumulation of neurodegenerative hallmarks [25] , suggest that the oxidation-induced impairment of GRP78 activity may play a role in the virus-induced neurodegenerative effects. Along this line, we also provided evidence of UPR activation by measuring the levels of PERK phosphorylation as well as at those of ATF6 cleavage. This may actually be due both to virus-induced accumulation of misfolded proteins and to its replication cycles in the brain. Indeed, it is well known that the rapid synthesis of viral proteins can elicit the UPR in host cells as a first defense against the infection [71] , and several viruses have evolved mechanisms to modulate UPR signaling for their own benefit. For instance, HSV1 is known to elicit PERK activation in the early stage of the infection, followed by both virus-mediated eIF2α dephosphorylation and degradation of cellular mRNA to allow the switch on of the translation machinery for viral protein synthesis [72] . In the same directions are other in vitro studies reporting a complex cross talk between HSV-1 infection and UPR activation. For instance, GRP78 expression was found both up-and downregulated during HSV-1 infection, mainly depending on HSV-1 strain and stage of the infection [73] . Other authors reported that the virus disarms the UPR in the early phase of the infection, activating only the ATF6 pathway and leaving the GRP78 levels unchanged and IRE1 inactive, to prevent further activation of ERAD that may be detrimental to virus replication [74] . This is in line with our data showing activation of ATF6 but not of IRE1. However, some authors reported IRE1 activation in HSV-1-infected cells, even though they documented that the RNAse and the kinase domains of this protein exert opposite roles in HSV-1 replication, inhibiting or favoring virus replication, respectively [75] . As a matter of fact, the net results of UPR activation upon virus entry may depend on the ability of the virus to exploit UPR pathways to support the correct folding of viral proteins and prevent their degradation through ERAD [74] . Thus, also our data may result from a complex balance between the effects of repeated HSV-1 replication in the brain (including the accumulation of AD-like misfolded proteins) and the host innate response to the infection. Actually, we found both positively and negatively modulated HNE adducts, suggesting that the presence of the virus may also preserve the oxidative status of some proteins. This seems the case of CRMP2 (Dpysl2), whose oxidation was found preserved following multiple cycles of virus replication in the brain (−10-fold oxidation in HSV1-M vs. CTRL-M, see Figure 2D and Table 1 ). CRMP2 is known to bind tubulin heterodimers and to act as a regulator of the microtubule network. Interestingly, post-translational modifications including SUMOylation, O-GlcNAcylation, oxidation, and phosphorylation strictly regulate these functions [31] . Furthermore, according to computational analyses of CRMP structure [76] , phosphorylation events may act sinergically with a redox switch at specific thiols to regulate the protein activity. This is in line with previous data reporting that upon oxidation on Cys504, a disulphide-linked CRMP2 homodimer forms a transient complex with thirodoxin1, allowing the activation of glycogen synthase kinase-3β (GSK-3β) and in turn the phosphorylation at Thr514 [77] . Indeed, among the kinases responsible of CRMP-2 phosphorylation, there are GSK-3β and cyclin-dependent kinase-5 (CDK-5), both kinases being involved also in tau phosphorylation. Specifically, GSK-3β phosphorylates CRMP2 on Thr-509 and Thr-514 [43] , upon CDK-5-mediated phosphorylation on Ser-522 [43, 78] . Interestingly CRMP2 hyperphosphorylation at these sites has been documented in AD brains, underlining its strict analogies to tau [42, 79] . Our results show that CRMP2 phosphorylation levels at Thr-514 were decreased in HSV1-M compared to CTRL-M. However, compelling evidence, including ours, demonstrated that the virus promotes GSK3β activation and tau hyperphosphorylation in in vitro and in vivo models of infection [21, 22, 80] . In addition, we previously showed the presence of an aggregated form of phosphorylated tau in cortical homogenates from mice which underwent repeated cycles of virus reactivation [25] . Thus, it is possible to speculate that the virus preferentially affects tau protein functions. In addition, CRMP2, beside its link with tau, exerts also different functions, being involved in vesicle transport, endosomal-lysosomal trafficking, and autophagy [42] . In particular, CRMP2 functions as a cargo adaptor to the molecular motors kinesin-1 [81] (for anterograde transport) and dynein [82] (for retrograde transport). It is well known that HSV-1 exploits the host cytoskeleton (actin and microtubule networks) to travel along the cytoplasm, either from the membrane to the cell body (retrograde transport, exploited upon virus entry to reach the nucleus for productive infection or latency establishment) or from the cell body to the plasmatic membrane (anterograde transport, exploited for the release of virus progeny outside the infected cells) [83] . Thus, it is likely that the virus preserves CRMP2 function to assure its travelling along axons at the early stages of virus-driven AD onset. Once AD-like features are built, it is possible that also CRMP2 may undergo the well-known post-translational modifications that affect its function. Such a possibility is further supported by previous studies demonstrating the increased oxidation of CRMP2 in non-HSV1-infected AD brain [84] [85] [86] . Further experiments are required to clarify CRMP2 role in HSV-1 infection in the brain. Overall, our redox proteomic data suggest that recurrent HSV-1 infections perturb the host cellular environment, in particular UPR, cytoskeleton network, and the clearance machinery for oxidized proteins and DNA. These events may contribute to the progressive accumulation of neurodegenerative hallmarks, supporting the role of HSV-1 in AD pathogenesis. Supplementary Materials: The following are available online at http://www.mdpi.com/2076-2607/8/7/972/s1, Figure S1 : Densitometric analysis of total and protein-specific levels of protein 3-nitrotyrosine (3 NT) and protein carbonyls (PC) in the brain of mice representing an in vivo model of recurrent HSV-1 reactivation; Figure S2 : Full-length images of western blots shown in Figures 4 and 6; Figure S3 : Full-length images of western blots shown in Figure 5 . The amyloid hypothesis of Alzheimer's disease at 25 years Infectious agents and neurodegeneration Herpes simplex virus type 1 and other pathogens are key causative factors in sporadic Alzheimer's disease Apolipoprotein E and oxidative stress in brain with relevance to Alzheimer's disease Oxidative stress, dysfunctional glucose metabolism and Alzheimer disease Role of 4-hydroxy-2-nonenal (HNE) in the pathogenesis of alzheimer disease and other selected age-related neurodegenerative disorders. Free Radic Herpes simplex virus type 1 in brain and risk of Alzheimer's disease Herpes simplex virus type 1, apolipoprotein E and Alzheimer' disease Multiscale analysis of independent Alzheimer's cohorts finds disruption of molecular, genetic, and clinical networks by human herpesvirus Seropositivity to herpes simplex virus antibodies and risk of Alzheimer's disease: A population-based cohort study Reactivated herpes simplex infection increases the risk of Alzheimer's disease. Alzheimer's Dement Increase in the IgG avidity index due to herpes simplex virus type 1 reactivation and its relationship with cognitive function in amnestic mild cognitive impairment and Alzheimer's disease Titers of herpes simplex virus type 1 antibodies positively correlate with grey matter volumes in Alzheimer's disease. J. Alzheimer's Dis A genetic signature including apolipoprotein Eepsilon4 potentiates the risk of herpes simplex-associated Alzheimer's disease. Alzheimer's Dement Interaction between APOE4 and herpes simplex virus type 1 in Alzheimer's disease. Alzheimer's Dement Herpes simplex virus interferes with amyloid precursor protein processing Herpes simplex virus infection causes cellular beta-amyloid accumulation and secretase upregulation APP processing induced by herpes simplex virus type 1 (HSV-1) yields several APP fragments in human and rat neuronal cells HSV-1 promotes Ca2+ -mediated APP phosphorylation and Abeta accumulation in rat cortical neurons The viral protein corona directs viral pathogenesis and amyloid aggregation Herpes simplex virus type 1 induces nuclear accumulation of hyperphosphorylated tau in neuronal cells Herpes simplex virus type-1 infection induces synaptic dysfunction in cultured cortical neurons via GSK-3 activation and intraneuronal amyloid-beta protein accumulation Herpes simplex virus type 1 infection in neurons leads to production and nuclear localization of APP intracellular domain (AICD): Implications for Alzheimer's disease pathogenesis Herpes simplex virus-type1 (HSV-1) impairs DNA repair in cortical neurons Recurrent herpes simplex virus-1 infection induces hallmarks of neurodegeneration and cognitive deficits in mice Herpes simplex virus type-1 infection impairs adult hippocampal neurogenesis via amyloid-beta protein accumulation Imbalance in corneal redox state during herpes simplex virus 1-induced keratitis in rabbits. Effectiveness of exogenous glutathione supply Herpes simplex virus type 1 infection induces oxidative stress and the release of bioactive lipid peroxidation by-products in mouse P19N neural cell cultures Role of oxidative damage in the pathogenesis of viral infections of the nervous system GRP78: A cell's response to stress Dysregulation of CRMP2 post-translational modifications drive its pathological functions Rapid in vivo reactivation of herpes simplex virus in latently infected murine ganglionic neurons after transient hyperthermia Oxidative stress and neurotoxicity Oxidative modification of lipoic acid by HNE in Alzheimer disease brain Elevated protein-bound levels of the lipid peroxidation product, 4-hydroxy-2-nonenal, in brain from persons with mild cognitive impairment STRING 8-a global view on proteins and their functional interactions in 630 organisms The STRING database in 2017: Quality-controlled protein-protein association networks, made broadly accessible COVID-19 spike-host cell receptor GRP78 binding site prediction The unfolded protein response is activated in Alzheimer's disease Aberrant protein phosphorylation in Alzheimer disease brain disturbs pro-survival and cell death pathways Early and selective activation and subsequent alterations to the unfolded protein response in down syndrome mouse models Collapsin response mediator protein-2: An emerging pathologic feature and therapeutic target for neurodisease indications Molecular mechanism of axon guidance mediated by phosphorylation of CRMP2 Herpes simplex virus-1 in the brain: The dark side of a sneaky infection Oxidative stress in neurodegenerative diseases Brain regional correspondence between Alzheimer's disease histopathology and biomarkers of protein oxidation Roles of 3-nitrotyrosine-and 4-hydroxynonenal-modified brain proteins in the progression and pathogenesis of Alzheimer's disease The role of oxidative stress in neurodegenerative diseases Oxidative stress in the progression of Alzheimer disease in the frontal cortex Role of oxidative stress in the progression of Alzheimer's disease. J. Alzheimer's Dis Oxidatively modified proteins in Alzheimer's disease (AD), mild cognitive impairment and animal models of AD: Role of Abeta in pathogenesis Redox proteomics identification of 4-hydroxynonenal-modified brain proteins in Alzheimer's disease: Role of lipid peroxidation in Alzheimer's disease pathogenesis Redox proteomic analysis of carbonylated brain proteins in mild cognitive impairment and early Alzheimer's disease An investigation of the molecular mechanisms engaged before and after the development of Alzheimer disease neuropathology in Down syndrome: A proteomics approach. Free Radic The chaperone BiP/GRP78 binds to amyloid precursor protein and decreases Abeta40 and Abeta42 secretion The UPR(ER): Sensor and coordinator of organismal homeostasis The unfolded protein response is activated in pretangle neurons in Alzheimer's disease hippocampus PERK-opathies: An endoplasmic reticulum stress mechanism underlying neurodegeneration The unfolded protein response is activated in disease-affected brain regions in progressive supranuclear palsy and Alzheimer's disease Expression of heat shock proteins in Alzheimer's disease Emerging roles of ER stress in the etiology and pathogenesis of Alzheimer's disease Bip/GRP78-induced production of cytokines and uptake of amyloid-beta(1-42) peptide in microglia Microglial activation and amyloid-beta clearance induced by exogenous heat-shock proteins Possible involvement of ER chaperone Grp78 on reduced formation of amyloid-beta deposits Microorganisms 2020 ER-stress in Alzheimer's disease: Turning the scale? Intraneuronal amyloid beta oligomers cause cell death via endoplasmic reticulum stress, endosomal/lysosomal leakage, and mitochondrial dysfunction in vivo Bip enhanced the association of GSK-3beta with tau during ER stress both in vivo and in vitro SIL1 rescued bip elevation-related tau hyperphosphorylation in ER stress AAV2-mediated GRP78 transfer alleviates retinal neuronal injury by downregulating ER stress and tau oligomer formation Impairment of proteostasis network in down syndrome prior to the development of Alzheimer's disease neuropathology: Redox proteomics analysis of human brain Endoplasmic reticulum stress responses Herpes simplex virus 1 infection activates the endoplasmic reticulum resident kinase PERK and mediates eIF-2alpha dephosphorylation by the gamma(1)34.5 protein Changes in BiP (GRP78) levels upon HSV-1 infection are strain dependent Herpes simplex virus-1 disarms the unfolded protein response in the early stages of infection Opposite roles of RNase and kinase activities of inositol-requiring enzyme 1 (IRE1) on HSV-1 replication Molecular dynamics simulations and in vitro analysis of the CRMP2 thiol switch Thioredoxin mediates oxidation-dependent phosphorylation of CRMP2 and growth cone collapse Increased CRMP2 phosphorylation is observed in Alzheimer's disease; does this tell us anything about disease development? Collapsin response mediator protein-2 (CRMP2) is a plausible etiological factor and potential therapeutic target in Alzheimer's disease: Comparison and contrast with microtubule-associated protein Tau Inflammatory and neurodegeneration markers during asymptomatic HSV-1 reactivation Tubulin and CRMP-2 complex is transported via Kinesin-1 CRMP-2 directly binds to cytoplasmic dynein and interferes with its activity Herpes simplex viruses Proteomic identification of oxidatively modified proteins in Alzheimer's disease brain. Part II: Dihydropyrimidinase-related protein 2, alpha-enolase and heat shock cognate 71 Redox proteomics identification of oxidized proteins in Alzheimer's disease hippocampus and cerebellum: An approach to understand pathological and biochemical alterations in AD Redox proteomics analysis of HNE-modified proteins in Down syndrome brain: Clues for understanding the development of Alzheimer disease. Free Radic The authors declare no conflict of interest. The funder had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.