Zhihui Yin
Xinyue Geng
Zhengyi Zhang
Ying Wang
Xiaoyan Gao- School of Chinese Materia Medica, Beijing University of Chinese Medicine, Beijing, China
Neuronal mitochondrial oxidative stress induced by β-amyloid (Aβ) is an early event of Alzheimer’s disease (AD). Emerging evidence has shown that antioxidant therapy represents a promising therapeutic strategy for the treatment of AD. In this study, we investigated the antioxidant activity of rhein against Aβ1-42 oligomer-induced mitochondrial oxidative stress in primary neurons and proposed a potential antioxidant pathway involved. The results suggested that rhein significantly reduced reactive oxygen species (ROS) level, reversed the depletion of mitochondrial membrane potential, and protected neurons from oxidative stress-associated apoptosis. Moreover, further study indicated that rhein activated mitochondrial biogenesis accompanied by increased cytochrome C oxidase (CytOx) and superoxide dismutase (SOD) activities. CytOx on the respiratory chain inhibited the production of ROS from electron leakage and SOD helped to eliminate excess ROS. Finally, western blot analysis confirmed that rhein remarkedly increased the protein expression of peroxisome proliferator-activated receptor gamma coactivator 1-alpha (PGC-1α) together with its upstream deacetylase sirtuin 1 (SIRT1), and activated downstream transcription factor nuclear respiratory factor 1, promoting mitochondrial biogenesis. In conclusion, our results demonstrate that rhein activates mitochondrial biogenesis regulated by the SIRT1/PGC-1α pathway as an antioxidant defense system against Aβ1-42 oligomer-induced oxidative stress. These findings broaden our knowledge of improving mitochondrial biogenesis as an approach for relieving neuronal oxidative stress in AD.
Introduction
Alzheimer’s disease (AD), the most common neurodegenerative disease, is characterized by memory loss and cognitive dysfunction (Alzheimer’s Association, 2021). Mounting evidence indicates neuronal oxidative stress induced by β-amyloid (Aβ) is considered to be a critical pathological event in AD, although the etiology of AD has not been fully elucidated (Lin and Beal, 2006; Borger et al., 2011; Swerdlow, 2018). Aβ destroys the electron respiratory chain in the mitochondrial inner membrane by reducing the activities of enzymes in mitochondrial respiratory chain complexes, especially cytochrome C oxidase (CytOx) in complex IV, leading to electron leakage and excess production of reactive oxygen species (ROS) (Atamna and Frey, 2004; Battogtokh et al., 2018; Misrani et al., 2021). Excess ROS aggravates mitochondrial dysfunction, which opens the mitochondrial permeability transition pore (mPTP), accelerates the release of cytochrome c (cyto c) into the cytosol, and subsequently induces the caspase-related apoptosis cascade, triggering the apoptosis of neurons (Mehan et al., 2011; Xiao et al., 2019). Meanwhile, ROS upregulates the expression of amyloid precursor protein and β-site APP-cleaving enzyme 1 by activating the c-Jun N-terminal kinase pathway, facilitates the generation and accumulation of Aβ, and finally accelerates the progression of AD (Chami and Checler, 2012; Tamagno et al., 2012; Zhao and Zhao, 2013; Ganguly et al., 2017). Therefore, relieving oxidative stress is an effective therapeutic strategy for AD.
Mitochondrial biogenesis as a cellular endogenous antioxidant defense system is actively mobilized under pathological conditions to overcome intracellular oxidative stress (Du et al., 2017; Panes et al., 2020). The sirtuin 1 (SIRT1)/peroxisome proliferator-activated receptor-γ coactivator 1-α (PGC-1α) pathway is pivotal for mitochondrial biogenesis (Menzies et al., 2013). SIRT1 is a nicotinamide adenine dinucleotide-dependent deacetylase, which is responsible for deacetylation and activation of the mitochondrial regulatory factor PGC-1α. Activated PGC-1α targets downstream transcription factors, such as nuclear respiratory factor 1 (NRF1) and NRF-2, which upregulate mitochondrial transcription factor A (TFAM). TFAM initiates the transcription and replication of mitochondrial DNA (mtDNA), and ultimately, newly healthy mitochondria are generated for the self-repair of mitochondria and synthesis of antioxidant enzymes (Canto et al., 2009; Kume et al., 2010; Aquilano et al., 2013). The newly generated healthy mitochondria can repair the damaged mitochondrial respiratory chain via mitochondrial fusion through mitochondrial dynamics, inhibiting the production of mitochondria-derived ROS (Liao et al., 2020; Zeng et al., 2021). In addition, synthesis of the mitochondrial antioxidant enzymes, mainly superoxide dismutase (SOD), is increased by generating new mitochondria, thus reducing ROS (St-Pierre et al., 2006; Dhar et al., 2008; Piantadosi and Suliman, 2008; Yan et al., 2016). Unfortunately, under the oxidative stress induced by Aβ, the SIRT1/PGC-1α pathway is suppressed (Gao et al., 2011; Shaerzadeh et al., 2014). Therefore, the activation of SIRT1/PGC-1α-regulated mitochondrial biogenesis as an antioxidant defense system to maintain intracellular mitochondrial homeostasis is essential for relieving the oxidative stress of neurons (Katsouri et al., 2016).
1,8-dihydroxyanthraquinone derivatives are a class of natural compounds with antioxidant activity (Wang et al., 2012; Lu et al., 2015; Markovic et al., 2016; Wang et al., 2016). Representative compounds including rhein, emodin, aloe-emodin, chrysophanol, and physcion are considered candidate compounds for antioxidant therapy. Rhein can protect IEC-6 cells from oxidative damage by inhibiting hydrogen peroxide (H2O2)-induced ROS, caspase 3, and even apoptosis (Zhuang et al., 2019). Aloe-emodin improves cell viability by reducing the levels of nitric oxide (NO) and ROS in PC12 cells induced by H2O2 (Tao et al., 2014). Physcion exhibits neuroprotection effect in SH-SY5Y cells by reducing Aβ-induced oxidative stress and the influx of Ca2+ (Ho et al., 2015). In vivo studies showed that chrysophanol enhanced the activities of SOD and manganese SOD, and inhibited the production of ROS, ultimately relieving oxidative stress injury in the brain of a focal cerebral ischemia/reperfusion mouse model (Zhao et al., 2018). Emodin downregulated 4-hydroxynonenal, a cerebral oxidative damage indicator, increased the expression levels of SOD1 and catalase, and improved the antioxidant capacity of APP/PS1 mice (Li et al., 2021). Therefore, screening of natural compounds with antioxidant activity from anthraquinones is of great clinical significance for relieving oxidative stress and protecting neurons in AD.
As one of the in vitro evaluation methods for AD, the cell model is very important for the screening of pharmacodynamic (PD) activity. To date, the screening platform for the evaluation of antioxidants against oxidative stress in AD is still dominated by immortalized cell lines, such as PC12 and SH-SY5Y cells (Lopes et al., 2017; Wang and Xu, 2019; Sotolongo et al., 2020). Nevertheless, after multiple rounds of proliferation and passaging, some traits of immortalized cell lines are quite different from those of living cells in vivo. As a consequence, the positive PD results screened from a cell model cannot be replicated when applied to an in vivo model. In contrast, the primary neurons obtained from living organisms are closer to the neurons in vivo. In addition, recent studies have suggested that Aβ1-42 oligomers are the pathogens of early AD (Haass and Selkoe, 2007; Zahs and Ashe, 2013; Fontana et al., 2020). Thus, establishing an Aβ1-42 oligomer-induced primary neuron model similar to the pathological environment of AD will greatly improve the screening effectiveness of candidate compounds.
In this study, an AD-like neuronal oxidative stress model was established on primary neurons by optimizing the aggregation states of Aβ1-42. Thereafter, the antioxidant activities of five anthraquinones were evaluated by intracellular ROS and mitochondrial membrane potential. As a result, rhein exhibited excellent antioxidant activity and inhibited neuronal apoptosis. Further study was to explore the antioxidant mechanism of rhein. Finally, our results demonstrate that rhein activates mitochondrial biogenesis regulated by the SIRT1/PGC-1α pathway against Aβ1-42 oligomer-induced mitochondrial oxidative stress. Taken together, rhein may present a promising candidate molecule for relieving neuronal oxidative stress as well as oxidative stress-associated neurodegenerative disorders.
Materials and Methods
Materials and Reagents
Aβ1-42 (CAS: 107761-42-2 for lyophilized Aβ1-42) was obtained from Nanjing Peptide Biotech Co., Ltd. (Nanjing, China). Rhein (MW 284.22 of purity >98%, CAS: 478-43-3), 1,1,1,3,3,3-hexafluoro-2-propanol (HFIP, CAS: 920-66-1), thioflavin T (ThT, CAS: 2390-54-7), and phosphotungstic acid (CAS: 12501-23-4) were obtained from Shanghai Aladdin Biochemical Technology Co., Ltd. (Shanghai, China). Emodin (MW 270.24 of purity >98%, CAS: 518-82-1), aloe-emodin (MW 270.24 of purity >97%, CAS: 481-72-1), chrysophanol (MW 254.24 of purity >98%, CAS: 481-74-3), and physcion (MW 284.26 of purity >98%, CAS: 521-61-9) were obtained from Shanghai yuanye Bio-Technology Co., Ltd. (Shanghai, China). B27 (Cat: 17504044) and penicillin-streptomycin (PS, Cat: 10378016) were purchased from Gibco Corporation (Grand Island, NY, USA). Dulbecco’s modified Eagle medium (DMEM, Cat: C11995500), Neurobasal medium (Cat: 21103049), and l-glutamine (Cat: A2916801) were purchased from Invitrogen Corporation (Carlsbad, CA, USA). Poly-l-lysine (Cat: P2100), 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT, Cat: M8180), and CytOx activity detection kit (Cat: BC0945) were obtained from Beijing Solarbio Science & Technology Co., Ltd. (Beijing, China). Bicinchoninic acid (BCA) protein quantification kit (Cat: P0010), ROS assay kit (Cat: S0033), Mitochondrial membrane potential detection kit (Cat: C2006), SOD activity detection kit (Cat: S0101), Annexin V-FITC apoptosis detection kit (Cat: C1062), and cytosol/mitochondria fraction isolation kit (Cat: C3601) were obtained from Beyotime Biotechnology (Shanghai, China). Primary antibodies, MAP2 rabbit polyclonal antibody (Cat: 17490-1-AP), PGC-1α monoclonal antibody (Cat: 66369-1-Ig), NRF1 polyclonal antibody (Cat: 12482-1-AP), caspase 3 polyclonal antibody (Cat: 19677-1-AP), cyto c polyclonal antibody (Cat: 10993-1-AP), and β-Actin monoclonal antibody (Cat: 66009-1-Ig) were obtained from Proteintech Group, Inc. (Rosemont, IL, USA). SIRT1 monoclonal antibody (Cat: 9475S) was purchased from Cell Signaling Technology Inc. (Danvers, MA, USA). Secondary antibodies, goat anti-mouse IgG (Cat: SA00001-1) and goat anti-rabbit IgG (Cat: SA00001-2) were obtained from Proteintech Group, Inc. All other chemicals used were of the highest grade commercially available.
Aβ1-42 Preparation
Aβ1-42 monomers were prepared as described previously (Sinha et al., 2011). The lyophilized Aβ1-42 peptides were dissolved in HFIP at a concentration of 1 mg/ml and sonicated for 10 min in an ice-bath, followed by incubation with shaking at 4°C for 3 h. Subsequently, the resulting solution was evaporated under a gentle stream of N2 gas to remove HFIP, and peptide films were formed. Thereafter, the peptide films were dissolved in 10 mM sodium hydroxide alone with sonication for 1 min and centrifuged at 4°C (16000 × g, 10 min). Finally, the supernatant was collected and quantified using a NanoDrop One spectrometer (Thermo Fisher Scientific, Waltham, MA, United States) at 280 nm using an extinction coefficient of 1490 M−1 cm−1.
To prepare different aggregation stages of Aβ1-42, Aβ1-42 monomers were diluted in phosphate-buffered solution (PBS, 10 mM, pH 7.4) at a final peptide concentration of 50 μM and incubated with continuous orbital shaking at 150 rpm. At different incubation time points, Aβ1-42 solution was collected for corresponding experiments.
Circular Dichroism Spectra Assay
The prepared Aβ1-42 was diluted in PBS at a final concentration of 25 μM. 100 μL sample was added to a 0.5 mm path length quartz cell. The spectrum was collected with an average of three scans at a speed of 1 nm/s over the wavelength range from 195 to 260 nm using a Chirascan spectrometer (Applied Photophysics Ltd., Surrey, BA, United Kingdom).
ThT Fluorescence Assay
The kinetic process of Aβ1-42 aggregation was monitored by ThT, the fluorescence intensity of which is dependent on the formation of β-sheet. At different incubation times, Aβ1-42 solution was taken for fluorescence assay at the concentration of 1 μM for Aβ1-42 and 10 μM for ThT, respectively. The fluorescence was measured using a LS45 fluorescence spectrometer (PerkinElmer Inc., Waltham, MA, United States). The excitation wavelength and emission wavelength were set at 440 and 480 nm, respectively.
Transmission Electron Microscopy
The prepared Aβ1-42 (20 μL) was dropped on a carbon-coated copper grid. After 1 h, the redundant solution was removed with filter paper, and prepared Aβ1-42 was stained with 20 μL of 1% phosphotungstic acid for 90 s. Then, the transmission electron microscopy (TEM) images were taken using a JEM-1230 transmission electron microscope (JEOL Ltd., Akishima-shi, Japan) with an acceleration voltage of 100 kV.
Primary Neuron Cultures
Primary neurons were extracted and isolated from neonatal SD rats as described previously (Dong et al., 2018). Neonatal SD rats were obtained from SPF Biotechnology Co., Ltd. (Beijing, China). Primary neurons were seeding in poly-l-lysine precoated 96-well plates or 24-well plates with a cell density at 5 × 105 cells/mL and cultured in neurobasal medium supplemented with 2% B27 supplements, 1% l-glutamine, and 1% PS at 37°C in a humidified atmosphere containing 5% CO2. After culture for 7 days, the primary neurons can be used for experiments.
As illustrated in Supplementary Figure 1, the purity of primary neurons was identified to be above 95% by the immunofluorescence assay labelled with anti-MAP2 antibodies, and the purity of the primary neurons met the requirements for subsequent experiments.
Cell Treatments
To examine the cytotoxicity of Aβ1-42 in different aggregation, primary neurons were incubated with DMEM containing 5 μM Aβ1-42 aggregates. After incubation for 24 h, the cell viability was measured. To examine the antioxidant effects of the five anthraquinones against Aβ1-42 oligomer-induced oxidative stress, primary neurons which were incubated with DMEM containing 5 μM Aβ1-42 and different concentrations of anthraquinones (1, 5 and 10 μM as low, medium, and high doses, respectively) were the treated groups. After incubation for 24 h, treated groups were used for further experiments. Meanwhile, untreated primary neurons were the control group, and primary neurons treated with 5 μM Aβ1-42 oligomers alone were the Aβ1-42 group.
Cell Viability Assay
Cell viability was measured by the MTT assay. After the cell treatment, the medium in 96-well plates was discarded. Then, 100 μL DMEM containing MTT at a final concentration of 0.5 mg/ml was added into 96-well plates. After incubation for another 4 h, the solution was discarded followed by adding 100 μL DMSO. Finally, the absorbance at 490 nm was measured using an Epoch Microplate Spectrophotometer (BioTek Instruments Inc., Winooski, VT, United States).
Intracellular ROS Measurement
ROS in primary neurons was detected using a DCFH-DA probe. DCFH-DA at a concentration of 10 μM was added to 96-well plates. After incubation for 20 min at 37°C, primary neurons were washed three times to fully remove excessive DCFH-DA. Finally, fresh culture medium was added and fluorescence images were taken under a fluorescence microscope (Nikon Instruments Inc., Melville, NY, United States).
Mitochondrial Membrane Potential (ΔΨm) Measurement
ΔΨm in primary neurons was detected using a JC-1 probe. Mitochondria were stained with a JC-1 staining kit according to the manufacturer’s instructions. Briefly, JC-1 staining solution was added to 24-well plates. After incubation for 20 min at 37°C, primary neurons were washed twice to fully remove excessive JC-1. Finally, fresh culture medium was added and fluorescence images were taken under a fluorescence microscope.
Molecular Docking
To study the structure-activity relationship of the five anthraquinones, rhein, emodin, aloe-emodin, chrysophanol, and physcion were selected as ligand molecules, and docked with the receptor of SIRT1 (PDB code: 4I5I), respectively. The 3D structures of the five ligand molecules were acquired from the Pubchem Compound database (https://pubchem.ncbi.nlm.nih.gov/). The X-ray crystal structure of human Sirt1 was obtained from the Protein Data Bank (https://www.rcsb.org/). Molecular docking was performed by AutoDock Vina 1.1.2 and AutoDock Tools (The Scripps Research Institute, La Jolla, CA, United States). The size of the grid box in AutoDock Vina was kept as 40 × 40 × 40 for X, Y, and Z. The final model was determined based on the binding energy and molecular interaction. Docked complexes were analyzed and figures were exported using PyMOL 2.3.4 (Schrödinger, Inc., New York, NY, United States).
Measurement of CytOx Activity
CytOx activity was measured using a CytOx activity detection kit according to the manufacturer’s instructions. Briefly, primary neurons after treatment were washed once with cold PBS and lysed in an extracting solution. The lysate was centrifuged at 4°C (600 × g, 10 min), and the supernatant was transferred to another tube for centrifugation at 4°C (11000 × g, 15 min). Then, the sediment was collected and extracted by extracting solution. After sufficient ultrasonication, each sample was quantified by a BCA protein quantification kit. Finally, 20 μL sample, and 200 μL working solution were added to 96-well plates and the absorbance at 550 nm was measured using an Epoch Microplate Spectrophotometer.
Measurement of SOD Activity
SOD activity was measured using a SOD activity detection kit according to the manufacturer’s instructions. Briefly, primary neurons after treatment were washed once with cold PBS and lysed in a SOD preparation solution. The lysate was centrifuged at 4°C (16000 × g, 5 min), and the supernatant was collected. Each sample was quantified by a BCA protein quantification kit. Then, 20 μL protein sample, 160 μL working solution, and 20 μL reaction working solution were added to 96-well plates. After incubation for 30 min at 37°C, the absorbance at 450 nm was measured using an Epoch Microplate Spectrophotometer.
Annexin V-FITC/Propidium Iodide Staining Assay
Apoptosis was evaluated using an Annexin V-FITC apoptosis detection kit according to the manufacturer’s instructions. In brief, primary neurons after treatment were washed once with cold PBS. Then, 195 μL Annexin V-FITC binding buffer, 5 μL Annexin V- FITC, and 10 μL propidium iodide (PI) were sequentially added to 96-well plates. After incubation for 15 min in the dark, fresh culture medium was added and fluorescence images were taken under a fluorescence microscope.
Western Blot
Primary neurons after treatment were washed once with cold PBS and lysed in an ice-cold lysis buffer containing 20 mM Tris-HCl pH 7.5, 150 mM NaCl, 1% Triton X-100, 1 mM PMSF and other inhibitors including sodium pyrophosphate, β-glycerophosphate, EDTA, Na3VO4 and leupeptin. The lysate was centrifuged at 4°C (16000 × g, 5 min), and the supernatant was collected. Cytosolic fraction for cytosolic cyto c was prepared using a cytosol/mitochondria fraction isolation kit. After quantification of protein concentration, equivalent amounts of sample were separated by an SDS-PAGE and then transferred onto polyvinylidene fluoride membranes. After being blocked with 5% non-fat milk for 2 h at room temperature, the membranes were incubated with primary antibodies of interest overnight at 4°C and followed by the incubation with secondary antibodies for 2 h at room temperature. The immunoblots were developed with ECL reagents and visualized by a ChemiDoc™ MP Imaging System (Bio-Rad Laboratories, Inc., Hercules, CA, United States). The results were normalized to β-Actin and analyzed using ImageJ software (National Institutes of Health, Bethesda, MD, United States).
Statistical Analysis
All data were presented as means ± standard deviation (SD) for at least three independent experiments. Statistical significance was analyzed using one-way analysis of variance (ANOVA), followed by the least significant difference (LSD) test for multiple comparisons. Levels of significance were indicated as follows: *p < 0.05; **p < 0.01; ***p < 0.001. All statistical analyses were performed using SPSS Statistics 25.0 software (IBM Corporation, Armonk, NY, United States).
Results
Aβ1-42 Oligomer-Induced Oxidative Stress Model Established Based on Primary Neurons
To establish an in vitro model of AD-like oxidative stress, the preparation of Aβ1-42 was optimized. First of all, the changes in the secondary structure of prepared Aβ1-42 with different incubation times were detected by circular dichroism (CD) spectra. As shown in Figure 1A, Aβ was in a monomeric state at 0 h and was structured as a random coil with a negative absorption peak at 198 nm. With the extension of incubation time, the absorption at 198 nm gradually weakened, and a new absorption peak appeared at 215 nm at 6 h. Its characteristic β-sheet structure at 215 nm revealed that Aβ changed from random coil in the monomer state to β-sheet (i.e., gradually aggregated). With increasing incubation times of 24, 48, and 96 h, the absorption at 215 nm gradually enhanced as well. As shown in Figure 1B, the secondary structure illustrated that the proportions of β-sheets at 6, 24, 48, and 96 h were 32, 44, 50, and 51%, respectively, implying that Aβ gradually aggregated and tended to be stable. The aggregation process of Aβ1-42 was further investigated by ThT fluorescence assay. As shown in Figure 1C, the fluorescence intensity of ThT rapidly increased from 6 to 36 h, indicating that Aβ1-42 was in the exponential growth phase. After 48 h, the curve grew slowly and gradually reached the plateau phase. These results demonstrated that the prepared Aβ1-42 depicted the typical aggregation kinetic process of monomer-oligomer-protofibril-fibril.
FIGURE 1. The cytotoxicity of Aβ1-42 in different aggregation stages in primary neurons. Aβ1-42 monomers were incubated by continuous orbital shaking at 150 rpm and 37°C. (A) The changes in the secondary structure of prepared Aβ1-42 in different aggregation stages were detected by circular dichroism (CD) spectra. (B) The proportions of the secondary structure were further analyzed. (C) The aggregation process of Aβ1-42 was investigated by the thioflavin T (ThT) assay (n = 4) and the cytotoxicity of Aβ1-42 in different aggregation stages in primary neurons was determined by the MTT assay (n = 5). (D) Representative transmission electron microscope image of Aβ1-42 incubated for 24 h at 37°C. Scale bars: 200 nm. Data are presented as the mean ± standard deviation (SD). *p < 0.05, **p < 0.01 and ***p < 0.001 compared with the proportions of β-sheets at 0 h (one-way ANOVA).
To examine the cytotoxicity of Aβ1-42 in different aggregation stages in primary neurons, Aβ1-42 aggregates prepared at different incubation times in the final concentration of 5 μM were used to evaluate the cell inhibition rate by the MTT assay. The results are shown in Figure 1C. The Aβ1-42 monomer did not cause distinct toxicity, and the cell inhibition rate was lower than 10%. With the extension of incubation time, the toxicity of Aβ1-42 significantly increased. After incubation for 6, 12, and 24 h, the toxicity of Aβ1-42 continued to increase, with inhibition rates of 41, 42, and 51%, respectively. At the plateau phase, the toxicity began to decrease. At 48, 60, 72, and 96 h, the inhibition rates of neurons were 27, 20, 15, and 9%, respectively. According to the abovementioned results, Aβ1-42 at 24 h caused maximum damage to the neurons. The morphology of the Aβ1-42 aggregates at 24 h was further characterized by TEM. After phosphotungstic acid staining, Aβ1-42 had a width of ∼5 nm and a length of <200 nm in the TEM image, suggesting that it was in the oligomeric form (Figure 1D). Taken together, the Aβ1-42 oligomers obtained by incubation for 24 h were the most toxic, which could be used to establish an AD oxidative stress cell model.
After determining the preparation method for toxic Aβ, primary neurons were stimulated by Aβ1-42 oligomers to induce oxidative stress. The intracellular ROS level was detected by the DCFH-DA probe. After being oxidized by ROS, DCFH-DA transformed to DCF, a fluorescent substance, and the fluorescence intensity of DCF was used to evaluate intracellular oxidative stress. As illustrated in Figures 2A,B, compared with the control group, the fluorescence intensity of DCF in the Aβ1-42 group increased to 203%, indicating that Aβ1-42 oligomers induced high ROS levels in primary neurons. Because intracellular ROS was derived from mitochondrial dysfunction, the intracellular ΔΨm was detected with a JC-1 probe to evaluate mitochondrial function. As shown in Figures 2C,D, the Aβ1-42 group exhibited more JC-1 monomers (green fluorescence) and fewer JC-1 aggregates (red fluorescence) than the control group. The fluorescence intensity of the red/green ratio decreased to 32% after being normalized to the control group, suggesting that a lower ΔΨm was due to mitochondrial damage. These results demonstrated that Aβ1-42 oligomers damaged mitochondrial function and caused oxidative stress in primary neurons.
FIGURE 2. Aβ1-42 oligomer-induced oxidative stress model on primary neurons. Untreated primary neurons were the control group, and primary neurons treated with 5 μM Aβ1-42 oligomers alone were the Aβ1-42 group. (A) Representative fluorescence images of intracellular ROS. (B) Quantification of intracellular ROS by the fluorescence intensity of DCF and normalization to the control group (n = 7). (C) Representative fluorescence images of ΔΨm. (D) Quantification analysis of ΔΨm by the fluorescence intensity of red/green ratio with JC-1 and normalization to the control group (n = 9). Scale bars: 50 μm. Data are presented as the mean ± standard deviation (SD). ***p < 0.001 compared with the Aβ1-42 group (one-way ANOVA).
Evaluation of the Antioxidative Activities of the Five Anthraquinones in Alleviating Aβ1-42 Oligomer-Induced Oxidative Stress in Primary Neurons
After establishing an oxidative stress neuron model induced by Aβ1-42 oligomers, we used this model to investigate the antioxidant activities of five anthraquinones. The chemical structures of the five anthraquinones are displayed in Figure 3A. They shared 1,8-dihydroxy anthraquinone as the same structure, yet they were substituted with different functional groups at the 3- or 6-position. Thereafter, the effects of the five anthraquinones on intracellular ROS level were determined. As shown in Figure 3B and Supplementary Figure 2, compared with the Aβ1-42 group, rhein treatment at doses of 1, 5, and 10 μM inhibited intracellular ROS levels by 182, 145, and 122%, respectively, manifesting an obvious dose dependency. The 5 μM rhein was sufficient to produce a significant difference, while the 10 μM rhein reduced ROS to the level of the control group. Emodin was also able to lower intracellular ROS levels. Emodin treatment at doses of 1, 5, and 10 μM reduced ROS levels by 189, 172, and 149%, respectively. Nevertheless, a significant difference was only found in the 10 μM emodin group. Aloe-emodin treatment at doses of 1, 5, and 10 μM reduced intracellular ROS levels by 196, 180, and 159%, respectively. Similar to aloe-emodin, the ROS scavenging effects of chrysophanol at doses of 1, 5, and 10 μM were 215, 181, and 162%, respectively. However, none of the three doses of physcion were able to effectively reduce intracellular ROS levels. The abovementioned results manifested that rhein, emodin, aloe-emodin, and chrysophanol could alleviate intracellular oxidative stress to varying degrees, with the exception of physcion.
FIGURE 3. The antioxidant effects of the five anthraquinones on primary neurons. Primary neurons were incubated with 5 μM Aβ1-42 oligomers and different concentrations of anthraquinones (1, 5 and 10 μM) for 24 h at 37°C, respectively. Untreated primary neurons were the control group, and primary neurons treated with 5 μM Aβ1-42 oligomers alone were the Aβ1-42 group. (A) The chemical structures of rhein, emodin, aloe-emodin, chrysophanol, and physcion. (B) Quantification of intracellular ROS by the fluorescence intensity of DCF and normalization to the control group (n = 7). (C) Quantification analysis of ΔΨm by the fluorescence intensity of red/green ratio with JC-1 and normalization to the control group (n = 9). Data are presented as the mean ± standard deviation (SD). *p < 0.05, **p < 0.01, ***p < 0.001 compared with the Aβ1-42 group; ###p < 0.001, compared with the 10 μM rhein group (one-way ANOVA).
The difference in antioxidant activity among the five anthraquinones was analyzed under the same condition, by comparing emodin, aloe-emodin, chrysophanol, and physcion with rhein at the same dose of 10 μM. Statistical analyses and the tendency of ROS to alleviate oxidative stress showed that the comparison for antioxidant activities of the five anthraquinones was rhein > emodin > aloe-emodin > chrysophanol > physcion.
JC-1 was used as an indicator to detect the intracellular ΔΨm. Figure 3C and Supplementary Figure 3 illustrated the results of ΔΨm. With the exception of physcion, which did not alleviate oxidative stress, the other four anthraquinones restored ΔΨm to varying degrees. The 1 μM rhein slightly increased the ΔΨm and the difference was insignificant. The 10 μM rhein increased the ΔΨm to 89%, and different doses presented a distinct dose-dependent increase. The effects of emodin and aloe-emodin were similar, also manifesting a dose-dependent increase with the strongest effects occurring at the dose of 10 μM (71 and 60%, respectively). The effect of the 5 μM chrysophanol on ΔΨm was stronger than that of 1 μM but was not significantly different from that of 10 μM.
Likewise, the intensity difference of the recovery effect among the five anthraquinones on ΔΨm was analyzed by comparing emodin, aloe-emodin, chrysophanol, and physcion with rhein at the same dose of 10 μM. Statistical analyses showed that the comparison for the effects of the five anthraquinones was rhein > emodin > aloe-emodin > chrysophanol > physcion.
Molecular Docking Simulation and Structure-activity Relationship Analysis
SIRT1-dependent deacetylation of PGC-1α promotes mitochondrial biogenesis, which can repair the function of the neuronal mitochondrial respiratory chain and increase the level of antioxidant enzymes, thereby synergically alleviating Aβ-induced oxidative stress (Lagouge et al., 2006; Li et al., 2020; Lu et al., 2020). Under the oxidative stress in AD, the expression levels of both SIRT1 and PGC-1α are reduced, leading to impaired mitochondrial biogenesis. Hence, oxidative stress in neurons cannot be alleviated (Qin et al., 2009; Sheng et al., 2012). Here, to compare the impact of the structural differences of the five anthraquinones on their antioxidant activity, SIRT1 was selected as the potential receptor protein to study the structure-activity relationship by molecular docking.
First of all, SIRT1 was taken as the receptor protein (PDB code: 4I5I) to analyze its binding capacity to anthraquinones. As shown in Figures 4A–E, all the five anthraquinones were directly bound to SIRT1 and occupied its active site, indicating that anthraquinones could interact with SIRT1. As shown in Table 1, the binding energy between rhein, emodin, aloe-emodin, chrysophanol, and physcion and SIRT1 was −9.4, −9.1, −8.8, −9.0, and −8.8 kcal/mol, respectively. Among them, the binding capacity between rhein and SIRT1 was the strongest, followed by emodin and chrysophanol, while the binding capacity of aloe-emodin was relatively weaker and equivalent to that of physcion.
FIGURE 4. The molecular simulation between five anthraquinones and SIRT1 (PDB code: 4I5I). (A-E) Rhein, emodin, aloe-emodin, chrysophanol, and physcion in the active site of SIRT1. (F-J) The binding modes of the five ligands with receptor SIRT1. Green dotted line: H-bonds; red arc: hydrophobic interaction.
TABLE 1. Binding affinity of the five anthraquinones with SIRT1.
Subsequently, the effect of substituents of anthraquinones on the binding capacity to SIRT1 was further compared by analyzing the binding mode. As shown in Figures 4F–J and Table 1, hydrogen bonds and hydrophobic interactions existed between SIRT1 and the five anthraquinones, which jointly determined the binding stability of the receptor-ligand. In terms of hydrophobicity, since the five anthraquinones shared the same hydrophobic structure of 1,8-dihydroxyanthraquinone, they all formed corresponding hydrophobic interactions with different amino acid residues on SIRT1, and different substituents had little effect on the number of residues of hydrophobic interactions. In contrast, the substituents had a noticeable effect on the number of hydrogen bonds. Rhein formed six hydrogen bonds with residues Arg274, Ser442, Asn346, and Leu443 in the active site of SIRT1. Among them, the carboxyl group at position-3 contributed three hydrogen bonds. Emodin formed three hydrogen bonds with His363, Gln345, and Val412, and the phenolic hydroxyl group at position-6 contributed one hydrogen bond. Aloe-emodin formed four hydrogen bonds with Asn346, Gln345, His363, and Val412. Among them, the hydroxymethyl group at position-3 contributed one hydrogen bond. Due to the lack of hydrogen bond donor or acceptor on the methyl group at position-3, chrysophanol did not form a stable hydrogen bond with SIRT1. Thus, chrysophanol was bound to SIRT1 only through the hydrophobic interactions. The electron distribution in the oxygen atom environment of methoxyl group at position-6 in physcion was in the homogeneous distribution, which impeded the formation of hydrogen bonds with residues in the active site; so, it did not contribute any hydrogen bonds.
The abovementioned results suggested that the polar substituents of anthraquinones (i.e., carboxyl, phenolic hydroxyl, and hydroxymethyl) were dominant groups and enhanced the binding capacity or binding stability with SIRT1. This finding was consistent with the experimental results on antioxidant activity.
Rhein Inhibited Cell Apoptosis Induced by Mitochondrial Oxidative Stress
Mitochondrial oxidative stress can induce neuronal oxidative damage and apoptosis. Encouraged by the excellent antioxidant activity of rhein, we further confirmed the role of rhein in neuroprotection by alleviating mitochondrial oxidative stress. Annexin V-FITC/PI double staining assay was performed to evaluate the effect of rhein on neuronal apoptosis. As shown in Figure 5A, compared with the control group, the green and red fluorescence intensity in the Aβ1-42 group was greatly enhanced, indicating that the neuronal cells were undergoing early and late apoptosis, respectively. Interestingly, compared with the Aβ1-42 group, rhein reversed apoptosis. The 10 μM rhein significantly reduced the intracellular green and red fluorescence intensity. The abovementioned results proved that rhein could inhibit the apoptosis of primary neurons by relieving oxidative stress.
FIGURE 5. Rhein inhibited mitochondrial oxidative stress-associated apoptosis in primary neurons induced by Aβ1-42 oligomers. Primary neurons were incubated with 5 μM Aβ1-42 oligomers and rhein at different doses (1, 5 and 10 μM) for 24 h at 37°C, respectively. Untreated primary neurons were the control group, and primary neurons treated with 5 μM Aβ1-42 oligomers alone were the Aβ1-42 group. (A) Representative images of Annexin V-FITC/propidium iodide (PI) double staining. Annexin V-FITC stains early apoptotic cells; PI stains late apoptotic cells. Scale bars: 50 μm. (B) Representative western blot images of cytosolic cytochrome c (cyto c). (C) Relative expression level of cytosolic Cyto c and normalization to β-Actin (n = 3). (D) Representative western blot images of cleaved caspase 3. (E) Relative expression level of cleaved caspase 3 and normalization to β-Actin (n = 3). Data are presented as the mean ± standard deviation (SD). **p < 0.01 and ***p < 0.001 compared with the Aβ1-42 group (one-way ANOVA).
Oxidative stress can damage the ultrastructural integrity of mitochondria. Cyto c, an apoptotic factor, is released from the inner mitochondrial membrane to activate procaspase 3 into cleaved caspase 3, which induces the apoptosis cascade. This is a critical pathological process for the loss of neuronal synapses (Ye et al., 2013). To investigate the mechanism underlying anti-apoptosis by rhein, the expression levels of cyto c in the cytosol and cleaved caspase 3 in primary neurons were evaluated by western blotting. The results are shown in Figures 5B,C. Because Aβ1-42 induced oxidative damage in mitochondria, the expression of cytosolic cyto c in the Aβ1-42 group significantly increased. Impressively, after treatment with rhein, the level of overexpressed cytosolic cyto c decreased, indicating that rhein reduced the release of cyto c from mitochondria. The expression level of cleaved caspase 3 is shown in Figures 5D,E. Compared with the control group, the level of cleaved caspase 3 in the Aβ1-42 group was significantly increased, suggesting that the apoptosis cascade was activated. By contrast, rhein inhibited the level of overexpressed cleaved caspase 3. As shown by the above results, rhein inhibited the release of cyto c from mitochondria by relieving oxidative damage in mitochondria to prevent the activation of the apoptosis cascade and ultimately avoid neuronal apoptosis.
Rhein Increased the Activities of Mitochondrial CytOx and SOD
Considering that oxidative stress can be effectively alleviated through mitochondrial biogenesis, we focused on the effects of rhein on enzymes related to mitochondrial biogenesis, such as CytOx in the mitochondrial respiratory complex IV and antioxidant enzyme SOD. First, the activity of CytOx was evaluated. As shown in Figure 6A, compared with the control group, Aβ1-42 damaged the mitochondrial electron transport chain and diminished the activity of CytOx. After treatment with rhein, the activity of CytOx increased. The effect of the 10 μM rhein group was close to that of the control group. Then, the antioxidant capacity of mitochondria was evaluated by detecting SOD activity, a mitochondria-related antioxidant enzyme. As shown in Figure 6B, the activity of SOD decreased in the Aβ1-42 group; thus, the antioxidant capacity of primary neurons was reduced as well. Notably, rhein improved the activity of SOD. These results revealed that rhein not only repaired the mitochondrial electron transport chain to inhibit the production of ROS but also improved the activity of mitochondrial antioxidant enzymes by activating the mitochondrial antioxidant defense system, thus enhancing the elimination of intracellular ROS, which might be related to the improvement of mitochondrial biogenesis.
FIGURE 6. Rhein increased the activities of mitochondrial cytochrome C oxidase (CytOx) and superoxide dismutase (SOD). Primary neurons were incubated with 5 μM Aβ1-42 oligomers and rhein at different doses (1, 5 and 10 μM) for 24 h at 37°C, respectively. Untreated primary neurons were the control group, and primary neurons treated with 5 μM Aβ1-42 oligomers alone were the Aβ1-42 group. The effects of rhein on the activities of (A) CytOx and (B) SOD (n = 3). Data are presented as the mean ± standard deviation (SD). *p < 0.05, **p < 0.01 and ***p < 0.001 compared with the Aβ1-42 group (one-way ANOVA).
Rhein Improved Mitochondrial Biogenesis by Activating the SIRT1/PGC-1α Pathway
The SIRT1/PGC-1α pathway is critical for mitochondrial biogenesis. To further validate the regulatory mechanism of rhein for mitochondrial biogenesis, the expression levels of SIRT1 and PGC-1α were detected by western blotting. Compared with the control group, the expression level of SIRT1 was significantly decreased in the Aβ1-42 group. As expected, rhein effectively reversed the decreased expression of SIRT1 (Figures 7A,B). Also, Aβ1-42 reduced the expression of PGC-1α, indicating that SIRT1/PGC-1α was inhibited in the Aβ-burdened neuronal model (Figures 7C,D). After treatment with rhein, the expression level of PGC-1α was significantly increased. Hence, rhein was able to activate the SIRT1/PGC-1α pathway. In addition, the expression level of NRF1, a downstream transcription factor related to mitochondrial biogenesis, was also detected. The results are shown in Figures 7C,E. In the Aβ1-42 group, because of the dysfunction of mitochondrial biogenesis, the expression level of NRF1 decreased. Rhein induced the expression of NRF1 by activating the SIRT1/PGC-1α pathway. It could be concluded that rhein could restore the function of mitochondrial biogenesis to alleviate oxidative stress in primary neurons by increasing the expression of SIRT1, PGC-1α, and NRF1.
FIGURE 7. Rhein enhanced mitochondrial biogenesis by activating the SIRT1/PGC-1α pathway. Primary neurons were incubated with 5 μM Aβ1-42 oligomers and rhein at different doses (1, 5 and 10 μM) for 24 h at 37°C, respectively. Untreated primary neurons were the control group, and primary neurons treated with 5 μM Aβ1-42 oligomers alone were the Aβ1-42 group. (A) Representative western blot images of SIRT1. (B) Relative expression level of SIRT1 and normalization to β-Actin (n = 3). (C) Representative western blotting images of PGC-1α and NRF1. (D) Relative expression level of PGC-1α and normalization to β-Actin (n = 3). (E) Relative expression level of NRF1 and normalization to β-Actin (n = 3). Data are presented as the group mean ± standard deviation (SD). *p < 0.05, **p < 0.01 and ***p < 0.001 compared with the Aβ1-42 group (one-way ANOVA).
Discussion
Neuronal oxidative stress induced by Aβ is an important pathological feature of AD. High ROS level derived from mitochondrial induces neuronal apoptosis (Lustbader et al., 2004; Petersen et al., 2008). Meanwhile, ROS accelerates the progression of AD by increasing the production and accumulation of Aβ (Takuma et al., 2009; Borger et al., 2013). Thus, introducing antioxidants to relieve oxidative stress is an effective therapeutic strategy for AD. 1,8-dihydroxyanthraquinone derivatives including rhein, emodin, aloe-emodin, chrysophanol, and physcion have garnered widespread attention as antioxidative components. Thus, screening of antioxidants from anthraquinones is extremely promising for the antioxidant therapy of AD. Herein, a comparative analysis of the antioxidant activity of the five anthraquinones was carried out using an Aβ1-42 oligomer-induced oxidative stress model based on primary cultured neurons. Among all the five anthraquinones, rhein possessed excellent antioxidant activity and was selected to explore its mechanism of alleviating oxidative stress.
With the progression of AD, accumulated Aβ monomers gradually aggregate and form oligomers, protofibrils, fibrils, and finally senile plaques (Kung, 2012; Wang et al., 2015). The latest β-amyloid cascade hypothesis suggests that Aβ in the oligomeric form is the earliest culprit of AD (Lacor et al., 2004; Staderini et al., 2015; Lee et al., 2017). To establish a drug screening platform that is similar to the in vivo pathological environment of oxidative stress in AD, the preparation of toxic Aβ was examined firstly on primary neurons. By combining CD with the ThT assay, the kinetic process of Aβ aggregation was described with β-sheets as the characteristic structure. Then, the correlation between the aggregation form of Aβ and cytotoxicity was analyzed. The results indicated Aβ oligomers incubated for 24 h caused the most damage to the primary neurons, while Aβ monomers or fibrils prepared in vitro were less toxic. TEM analysis further confirmed the formation of Aβ oligomers. These findings are consistent with recent researches (Wang-Dietrich et al., 2013; Shea et al., 2019). After incubation with Aβ1-42 oligomers, the primary neurons exhibited high ROS level and depletion of ΔΨm, verifying that a mitochondrial oxidative stress cell model was successfully established.
For comparative analysis of antioxidant activities of five anthraquinones, we used aforementioned Aβ1-42 oligomer-induced primary neurons as an in vitro oxidative stress model. The results suggested that rhein, emodin, aloe-emodin, and chrysophanol reduced intracellular ROS levels to varying degrees in a dose-dependent manner, with the exception of physcion. As expected, rhein, emodin, aloe-emodin, and chrysophanol effectively restored ΔΨm, except physcion. These results indicate that the comparison for antioxidant activities of the five anthraquinones is: rhein > emodin > aloe-emodin > chrysophanol > physcion. To theoretically elucidate the effects of the structural differences of the five anthraquinones on the antioxidant activity, molecular docking was performed to compare the binding mode between anthraquinones and the receptor protein SIRT1. The results showed that rhein, emodin, aloe-emodin, chrysophanol, and physcion were directly bound to the active site of SIRT1. Among them, rhein showed the relatively lower binding energy of −9.4, indicating a better binding capacity with SIRT1. Furthermore, the binding mode showed that the carboxyl group on rhein formed more hydrogen bonds with amino acid residues and increased the stability of the rhein-SIRT1 complex. These findings might partly explain why rhein exhibited excellent antioxidant activity against oxidative stress induced by Aβ1-42 oligomers.
Ongoing oxidative stress can trigger apoptosis cascade and induce neuronal apoptosis. When neuronal mitochondria undergo oxidative damage, the apoptosis factor, cyto c, is released from the opened mPTP into the cytosol. Then, procaspase 3 is activated into cleaved caspase 3 by cyto c, which induces the downstream apoptosis cascade and neuronal apoptosis irreversibly (Bonda et al., 2014). Annexin V-FITC/PI staining showed that the number of neuronal cells with early apoptosis (green fluorescence) and late apoptosis (red fluorescence) was significantly increased in the Aβ1-42 group, indicating neuronal apoptosis occurred (Figure 5A). After treatment with rhein, the number of apoptotic neurons was effectively reduced. According to the western blot analysis, there was the overexpression of cytosolic cyto c and cleaved caspase 3 in the Aβ1-42 group. Interestingly, the levels of cytosolic cyto c and cleaved caspase 3 were reversed after rhein treatment. Thus, rhein reduced the release of cyto c from mitochondria, inhibiting the apoptosis cascade and ultimately protecting neurons from apoptosis.
Neuronal oxidative stress is closely related to mitochondrial dysfunction. Mitochondrial biogenesis is an intracellular antioxidative defense system, in which cells maintain healthy mitochondrial function by producing new mitochondria. Newly generated mitochondria repair the respiratory chain complexes by mitochondrial fusion and further prevent electron leakage, thereby inhibiting the continuous production of ROS. In addition, the increased synthesis of antioxidant enzymes in mitochondria can effectively remove the accumulated ROS to maintain mitochondrial redox homeostasis (Venditti et al., 2013). Under the oxidative stress induced by Aβ, mitochondrial biogenesis is blocked. The activities of enzymes in mitochondrial respiratory chain complexes and antioxidant enzymes decrease, and the antioxidant defense system is also destroyed. Indeed, in this study, it was found that the activity of CytOx, an enzyme in respiratory chain complex IV, was significantly reduced in the mitochondria of primary neurons damaged by Aβ1-42 oligomers (Atamna and Frey, 2004; Atamna et al., 2009). However, the activity of CytOx increased after treatment with rhein, indicating that the activity of enzymes in the respiratory chain was improved, which helped to reduce ROS production. Due to the impairment of mitochondrial function caused by Aβ1-42 oligomers, the activity of SOD, the main enzyme of mitochondrial antioxidants, was greatly reduced, making mitochondria unable to resist the existing oxidative stress. By contrast, SOD activity was significantly increased under the treatment of rhein, so that mitochondria could effectively remove the existing ROS. The results showed that rhein played a positive role in regulating both the enzymes in mitochondrial respiratory chain complexes and antioxidant enzymes, suggesting that it may be beneficial to improve mitochondrial biogenesis.
The SIRT1/PGC-1α pathway plays a crucial role in mitochondrial biogenesis. SIRT1, a deacetylase, is responsible for the deacetylation and activation of PGC-1α. PGC-1α, an important regulator of mitochondrial biogenesis, activates the downstream transcription factors NRF-1, NRF-2, and TFAM, and initiates the transcription and replication of mtDNA. Therefore, the expression of SIRT1 and PGC-1α is closely related to the self-repair of mitochondria and synthesis of antioxidant enzymes (St-Pierre et al., 2006; Kang et al., 2017). To explore the mechanism underlying the activation of mitochondrial biogenesis by rhein, the expression of related proteins was evaluated by western blot analysis. The results confirmed that mitochondrial biogenesis was impaired by Aβ1-42 oligomers accompanied by the reduced expression levels of SIRT1, PGC-1α, and downstream transcription factor NRF1. As a result, after treatment with rhein, the expression levels of SIRT1 and PGC-1α were significantly increased, indicating that the SIRT1/PGC-1α pathway was activated. Meanwhile, the increased expression of NRF1 confirmed that rhein was involved in the regulation of mitochondrial biogenesis. Hence, the activation of SIRT1/PGC-1α-mediated mitochondrial biogenesis by rhein might be a key mechanism for triggering the mitochondrial antioxidant defense system (Figure 8).
FIGURE 8. Schematic illustration of mitochondrial endogenous antioxidant defense system triggered by rhein. 1). Rhein activated the SIRT1/PGC-1α pathway by increasing the expression levels of SIRT1 and PGC-1α. 2). The PGC-1α in the nucleus activated the downstream transcription factor NRF1, and healthy mitochondria were generated by promoted mitochondrial biogenesis. 3). Newly generated healthy mitochondria restored the decreased activity of CytOx in mitochondrial respiratory chain complex IV induced by Aβ1-42 oligomers, inhibiting the production of ROS. 4). mitochondrial SOD effectively facilitated the elimination of excess ROS. 5). Rhein alleviated the mitochondrial oxidative damage, inhibited the caspase 3-related apoptosis cascade by reducing the release of cyto c from mitochondria, and ultimately, protected neurons from apoptosis.
In summary, we used an AD-like neuronal oxidative stress model induced by Aβ1-42 oligomers to screen the antioxidant activities of the five anthraquinones, including rhein, emodin, aloe-emodin, chrysophanol, and physcion. According to the results of intracellular ROS and ΔΨm, rhein exhibited outstanding antioxidant activity and inhibited oxidative stress-associated neuronal apoptosis. More importantly, rhein activated mitochondrial biogenesis as an endogenous antioxidant defense system against Aβ1-42 oligomer-induced oxidative stress. That is to say, CytOx in the respiratory chain complex IV inhibited the production of ROS from electron leakage and SOD helped to eliminate superfluous ROS by promoted mitochondrial biogenesis. Western blot analysis further confirmed that the SIRT1/PGC-1α pathway activated by rhein was a potential antioxidant pathway involved. Taken together, our results provide evidence that rhein activates mitochondrial biogenesis regulated by the SIRT1/PGC-1α pathway as an antioxidant defense system against Aβ1-42 oligomer-induced oxidative stress. This study contributes to the fundamental research for the antioxidant therapy of AD.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.
Author Contributions
XGa designed the research. ZY and XGe performed the experiments. XGe and ZZ analysed the data. ZY and YW wrote the paper. XGa and YW reviewed and edited the draft.
Funding
This work was supported by the National Natural Science Foundation of China (81973467/H2803) and the Fundamental Research Funds for the Central Universities (2020-JYB-ZDGG-033).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2021.746711/full#supplementary-material
Abbreviations
ΔΨm, mitochondrial membrane potential; Aβ, β-amyloid; AD, Alzheimer’s disease; BCA, bicinchoninic acid; CD, circular dichroism; CytOx, cytochrome C oxidase; cyto c, cytochrome c; H2O2, hydrogen peroxide; mPTP, mitochondrial permeability transition pore; mtDNA, mitochondrial DNA; MTT, 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide; NO, nitric oxide; NRF1, nuclear respiratory factor 1; PD, pharmacodynamic; PDB, Protein Data Bank; PGC-1α, peroxisome proliferator-activated receptor gamma coactivator 1-alpha; PI: propidium iodide; ROS, reactive oxygen species; SIRT1, sirtuin 1; SOD, superoxide dismutase; TEM, transmission electron microscopy; TFAM, transcription factor A; ThT, thioflavin T.
References
Alzheimer’s Association (2021). 2021 Alzheimer's Disease Facts and Figures. Alzheimers Dement 17 (3), 327–406. doi:10.1002/alz.12328
Aquilano, K., Baldelli, S., Pagliei, B., and Ciriolo, M. R. (2013). Extranuclear Localization of SIRT1 and PGC-1α: an Insight into Possible Roles in Diseases Associated with Mitochondrial Dysfunction. Curr. Mol. Med. 13 (1), 140–154. doi:10.2174/156652413804486241
Atamna, H., and Frey, W. H. (2004). A Role for Heme in Alzheimer's Disease: Heme Binds Amyloid Beta and Has Altered Metabolism. Proc. Natl. Acad. Sci. U. S. A. 101 (30), 11153–11158. doi:10.1073/pnas.0404349101
Atamna, H., Frey, W. H., and Ko, N. (2009). Human and Rodent Amyloid-Beta Peptides Differentially Bind Heme: Relevance to the Human Susceptibility to Alzheimer's Disease. Arch. Biochem. Biophys. 487 (1), 59–65. doi:10.1016/j.abb.2009.05.003
Battogtokh, G., Choi, Y. S., Kang, D. S., Park, S. J., Shim, M. S., Huh, K. M., et al. (2018). Mitochondria-targeting Drug Conjugates for Cytotoxic, Anti-oxidizing and Sensing Purposes: Current Strategies and Future Perspectives. Acta Pharm. Sin. B 8 (6), 862–880. doi:10.1016/j.apsb.2018.05.006
Bonda, D. J., Wang, X., Lee, H. G., Smith, M. A., Perry, G., and Zhu, X. (2014). Neuronal Failure in Alzheimer's Disease: A View through the Oxidative Stress Looking-Glass. Neurosci. Bull. 30 (2), 243–252. doi:10.1007/s12264-013-1424-x
Borger, E., Aitken, L., Muirhead, K. E., Allen, Z. E., Ainge, J. A., Conway, S. J., et al. (2011). Mitochondrial β-amyloid in Alzheimer's Disease. Biochem. Soc. Trans. 39, 868–873. doi:10.1042/bst0390868
Borger, E., Aitken, L., Du, H., Zhang, W., Gunn-Moore, F. J., and Yan, S. S. (2013). Is Amyloid Binding Alcohol Dehydrogenase a Drug Target for Treating Alzheimer's Disease? Curr. Alzheimer Res. 10 (1), 21–29. doi:10.2174/156720513804871381
Cantó, C., Gerhart-Hines, Z., Feige, J. N., Lagouge, M., Noriega, L., Milne, J. C., et al. (2009). AMPK Regulates Energy Expenditure by Modulating NAD+ Metabolism and SIRT1 Activity. Nature 458 (7241), 1056–1060. doi:10.1038/nature07813
Chami, L., and Checler, F. (2012). BACE1 Is at the Crossroad of a Toxic Vicious Cycle Involving Cellular Stress and β-amyloid Production in Alzheimer's Disease. Mol. Neurodegener. 7, 52. doi:10.1186/1750-1326-7-52
Dhar, S. S., Ongwijitwat, S., and Wong-Riley, M. T. (2008). Nuclear Respiratory Factor 1 Regulates All Ten Nuclear-Encoded Subunits of Cytochrome C Oxidase in Neurons. J. Biol. Chem. 283 (6), 3120–3129. doi:10.1074/jbc.M707587200
Dong, Y. T., Cao, K., Tan, L. C., Wang, X. L., Qi, X. L., Xiao, Y., et al. (2018). Stimulation of SIRT1 Attenuates the Level of Oxidative Stress in the Brains of APP/PS1 Double Transgenic Mice and in Primary Neurons Exposed to Oligomers of the Amyloid-β Peptide. J. Alzheimers Dis. 63 (1), 283–301. doi:10.3233/jad-171020
Du, F., Yu, Q., Yan, S., Hu, G., Lue, L. F., Walker, D. G., et al. (2017). PINK1 Signalling Rescues Amyloid Pathology and Mitochondrial Dysfunction in Alzheimer's Disease. Brain 140, 3233–3251. doi:10.1093/brain/awx258
Fontana, I. C., Zimmer, A. R., Rocha, A. S., Gosmann, G., Souza, D. O., Lourenco, M. V., et al. (2020). Amyloid-β Oligomers in Cellular Models of Alzheimer's Disease. J. Neurochem. 155 (4), 348–369. doi:10.1111/jnc.15030
Ganguly, G., Chakrabarti, S., Chatterjee, U., and Saso, L. (2017). Proteinopathy, Oxidative Stress and Mitochondrial Dysfunction: Cross Talk in Alzheimer's Disease and Parkinson's Disease. Drug Des. Devel Ther. 11, 797–810. doi:10.2147/dddt.S130514
Gao, Z., Zhang, J., Kheterpal, I., Kennedy, N., Davis, R. J., and Ye, J. (2011). Sirtuin 1 (SIRT1) Protein Degradation in Response to Persistent C-Jun N-Terminal Kinase 1 (JNK1) Activation Contributes to Hepatic Steatosis in Obesity. J. Biol. Chem. 286 (25), 22227–22234. doi:10.1074/jbc.M111.228874
Haass, C., and Selkoe, D. J. (2007). Soluble Protein Oligomers in Neurodegeneration: Lessons from the Alzheimer's Amyloid Beta-Peptide. Nat. Rev. Mol. Cel Biol. 8 (2), 101–112. doi:10.1038/nrm2101
Ho, S.-L., Poon, C.-Y., Lin, C., Yan, T., Kwong, D., Yung, K., et al. (2015). Inhibition of β-Amyloid Aggregation by Albiflorin, Aloeemodin and Neohesperidin and Their Neuroprotective Effect on Primary Hippocampal Cells against β-Amyloid Induced Toxicity. Car 12 (5), 424–433. doi:10.2174/1567205012666150504144919
Kang, H., Khang, R., Ham, S., Jeong, G. R., Kim, H., Jo, M., et al. (2017). Activation of the ATF2/CREB-PGC-1α Pathway by Metformin Leads to Dopaminergic Neuroprotection. Oncotarget 8 (30), 48603–48618. doi:10.18632/oncotarget.18122
Katsouri, L., Lim, Y. M., Blondrath, K., Eleftheriadou, I., Lombardero, L., Birch, A. M., et al. (2016). PPARγ-coactivator-1α Gene Transfer Reduces Neuronal Loss and Amyloid-β Generation by Reducing β-secretase in an Alzheimer's Disease Model. Proc. Natl. Acad. Sci. USA 113 (43), 12292–12297. doi:10.1073/pnas.1606171113
Kume, S., Uzu, T., Horiike, K., Chin-Kanasaki, M., Isshiki, K., Araki, S., et al. (2010). Calorie Restriction Enhances Cell Adaptation to Hypoxia through Sirt1-dependent Mitochondrial Autophagy in Mouse Aged Kidney. J. Clin. Invest. 120 (4), 1043–1055. doi:10.1172/jci41376
Kung, H. F. (2012). The β-Amyloid Hypothesis in Alzheimer's Disease: Seeing Is Believing. ACS Med. Chem. Lett. 3 (4), 265–267. doi:10.1021/ml300058m
Lacor, P. N., Buniel, M. C., Chang, L., Fernandez, S. J., Gong, Y., Viola, K. L., et al. (2004). Synaptic Targeting by Alzheimer's-Related Amyloid Beta Oligomers. J. Neurosci. 24 (45), 10191–10200. doi:10.1523/jneurosci.3432-04.2004
Lagouge, M., Argmann, C., Gerhart-Hines, Z., Meziane, H., Lerin, C., Daussin, F., et al. (2006). Resveratrol Improves Mitochondrial Function and Protects against Metabolic Disease by Activating SIRT1 and PGC-1alpha. Cell 127 (6), 1109–1122. doi:10.1016/j.cell.2006.11.013
Lee, S. J., Nam, E., Lee, H. J., Savelieff, M. G., and Lim, M. H. (2017). Towards an Understanding of Amyloid-β Oligomers: Characterization, Toxicity Mechanisms, and Inhibitors. Chem. Soc. Rev. 46 (2), 310–323. doi:10.1039/c6cs00731g
Li, B., Liu, J., Gu, G., Han, X., Zhang, Q., and Zhang, W. (2020). Impact of Neural Stem Cell-Derived Extracellular Vesicles on Mitochondrial Dysfunction, Sirtuin 1 Level, and Synaptic Deficits in Alzheimer's Disease. J. Neurochem. 154 (5), 502–518. doi:10.1111/jnc.15001
Li, Z., Bi, H., Jiang, H., Song, J., Meng, Q., Zhang, Y., et al. (2021). Neuroprotective Effect of Emodin against Alzheimer's Disease via Nrf2 Signaling in U251 Cells and APP/PS1 Mice. Mol. Med. Rep. 23 (2), 108. doi:10.3892/mmr.2020.11747
Liao, L. X., Wang, J. K., Wan, Y. J., Liu, Y., Dong, X., Tu, P. F., et al. (2020). Protosappanin A Maintains Neuronal Mitochondrial Homeostasis through Promoting Autophagic Degradation of Bax. ACS Chem. Neurosci. 11 (24), 4223–4230. doi:10.1021/acschemneuro.0c00488
Lin, M. T., and Beal, M. F. (2006). Mitochondrial Dysfunction and Oxidative Stress in Neurodegenerative Diseases. Nature 443 (7113), 787–795. doi:10.1038/nature05292
Lopes, F. M., da Motta, L. L., De Bastiani, M. A., Pfaffenseller, B., Aguiar, B. W., de Souza, L. F., et al. (2017). RA Differentiation Enhances Dopaminergic Features, Changes Redox Parameters, and Increases Dopamine Transporter Dependency in 6-Hydroxydopamine-Induced Neurotoxicity in SH-SY5Y Cells. Neurotox. Res. 31 (4), 545–559. doi:10.1007/s12640-016-9699-0
Lu, K., Zhang, C., Wu, W., Zhou, M., Tang, Y., and Peng, Y. (2015). Rhubarb Extract Has a Protective Role against Radiation-Induced Brain Injury and Neuronal Cell Apoptosis. Mol. Med. Rep. 12 (2), 2689–2694. doi:10.3892/mmr.2015.3693
Lu, J., Huang, Q., Zhang, D., Lan, T., Zhang, Y., Tang, X., et al. (2020). The Protective Effect of DiDang Tang against AlCl3-Induced Oxidative Stress and Apoptosis in PC12 Cells through the Activation of SIRT1-Mediated Akt/Nrf2/HO-1 Pathway. Front. Pharmacol. 11, 466. doi:10.3389/fphar.2020.00466
Lustbader, J. W., Cirilli, M., Lin, C., Xu, H. W., Takuma, K., Wang, N., et al. (2004). ABAD Directly Links Abeta to Mitochondrial Toxicity in Alzheimer's Disease. Science 304 (5669), 448–452. doi:10.1126/science.1091230
Marković, Z., Jeremić, S., Dimitrić Marković, J., Stanojević Pirković, M., and Amić, D. (2016). Influence of Structural Characteristics of Substituents on the Antioxidant Activity of Some Anthraquinone Derivatives. Comput. Theor. Chem. 1077, 25–31. doi:10.1016/j.comptc.2015.10.004
Mehan, S., Meena, H., Sharma, D., and Sankhla, R. (2011). JNK: a Stress-Activated Protein Kinase Therapeutic Strategies and Involvement in Alzheimer's and Various Neurodegenerative Abnormalities. J. Mol. Neurosci. 43 (3), 376–390. doi:10.1007/s12031-010-9454-6
Menzies, K. J., Singh, K., Saleem, A., and Hood, D. A. (2013). Sirtuin 1-mediated Effects of Exercise and Resveratrol on Mitochondrial Biogenesis. J. Biol. Chem. 288 (10), 6968–6979. doi:10.1074/jbc.M112.431155
Misrani, A., Tabassum, S., and Yang, L. (2021). Mitochondrial Dysfunction and Oxidative Stress in Alzheimer's Disease. Front. Aging Neurosci. 13, 617588. doi:10.3389/fnagi.2021.617588
Panes, J. D., Godoy, P. A., Silva-Grecchi, T., Celis, M. T., Ramirez-Molina, O., Gavilan, J., et al. (2020). Changes in PGC-1α/SIRT1 Signaling Impact on Mitochondrial Homeostasis in Amyloid-Beta Peptide Toxicity Model. Front. Pharmacol. 11, 709. doi:10.3389/fphar.2020.00709
Hansson Petersen, C. A., Alikhani, N., Behbahani, H., Wiehager, B., Pavlov, P. F., Alafuzoff, I., et al. (2008). The Amyloid Beta-Peptide Is Imported into Mitochondria via the TOM Import Machinery and Localized to Mitochondrial Cristae. Proc. Natl. Acad. Sci. U. S. A. 105 (35), 13145–13150. doi:10.1073/pnas.0806192105
Piantadosi, C. A., and Suliman, H. B. (2008). Transcriptional Regulation of SDHa Flavoprotein by Nuclear Respiratory Factor-1 Prevents Pseudo-hypoxia in Aerobic Cardiac Cells. J. Biol. Chem. 283 (16), 10967–10977. doi:10.1074/jbc.M709741200
Qin, W., Haroutunian, V., Katsel, P., Cardozo, C. P., Ho, L., Buxbaum, J. D., et al. (2009). PGC-1alpha Expression Decreases in the Alzheimer Disease Brain as a Function of Dementia. Arch. Neurol. 66 (3), 352–361. doi:10.1001/archneurol.2008.588
Shaerzadeh, F., Motamedi, F., and Khodagholi, F. (2014). Inhibition of Akt Phosphorylation Diminishes Mitochondrial Biogenesis Regulators, Tricarboxylic Acid Cycle Activity and Exacerbates Recognition Memory Deficit in Rat Model of Alzheimer's Disease. Cell. Mol. Neurobiol. 34 (8), 1223–1233. doi:10.1007/s10571-014-0099-9
Shea, D., Hsu, C. C., Bi, T. M., Paranjapye, N., Childers, M. C., Cochran, J., et al. (2019). α-Sheet Secondary Structure in Amyloid β-peptide Drives Aggregation and Toxicity in Alzheimer's Disease. Proc. Natl. Acad. Sci. U. S. A. 116 (18), 8895–8900. doi:10.1073/pnas.1820585116
Sheng, B., Wang, X., Su, B., Lee, H. G., Casadesus, G., Perry, G., et al. (2012). Impaired Mitochondrial Biogenesis Contributes to Mitochondrial Dysfunction in Alzheimer's Disease. J. Neurochem. 120 (3), 419–429. doi:10.1111/j.1471-4159.2011.07581.x
Sinha, S., Lopes, D. H., Du, Z., Pang, E. S., Shanmugam, A., Lomakin, A., et al. (2011). Lysine-specific Molecular Tweezers Are Broad-Spectrum Inhibitors of Assembly and Toxicity of Amyloid Proteins. J. Am. Chem. Soc. 133 (42), 16958–16969. doi:10.1021/ja206279b
Sotolongo, K., Ghiso, J., and Rostagno, A. (2020). Nrf2 Activation through the PI3K/GSK-3 axis Protects Neuronal Cells from Aβ-Mediated Oxidative and Metabolic Damage. Alzheimers Res. Ther. 12 (1), 13. doi:10.1186/s13195-019-0578-9
Staderini, M., Martín, M. A., Bolognesi, M. L., and Menéndez, J. C. (2015). Imaging of β-amyloid Plaques by Near Infrared Fluorescent Tracers: a New Frontier for Chemical Neuroscience. Chem. Soc. Rev. 44 (7), 1807–1819. doi:10.1039/c4cs00337c
St-Pierre, J., Drori, S., Uldry, M., Silvaggi, J. M., Rhee, J., Jäger, S., et al. (2006). Suppression of Reactive Oxygen Species and Neurodegeneration by the PGC-1 Transcriptional Coactivators. Cell 127 (2), 397–408. doi:10.1016/j.cell.2006.09.024
Swerdlow, R. H. (2018). Mitochondria and Mitochondrial Cascades in Alzheimer's Disease. J. Alzheimers Dis. 62 (3), 1403–1416. doi:10.3233/JAD-170585
Takuma, K., Fang, F., Zhang, W., Yan, S., Fukuzaki, E., Du, H., et al. (2009). RAGE-mediated Signaling Contributes to Intraneuronal Transport of Amyloid-Beta and Neuronal Dysfunction. Proc. Natl. Acad. Sci. U. S. A. 106 (47), 20021–20026. doi:10.1073/pnas.0905686106
Tamagno, E., Guglielmotto, M., Monteleone, D., and Tabaton, M. (2012). Amyloid-β Production: Major Link between Oxidative Stress and BACE1. Neurotox. Res. 22 (3), 208–219. doi:10.1007/s12640-011-9283-6
Tao, L., Xie, J., Wang, Y., Wang, S., Wu, S., Wang, Q., et al. (2014). Protective Effects of Aloe-Emodin on Scopolamine-Induced Memory Impairment in Mice and H₂O₂-induced Cytotoxicity in PC12 Cells. Bioorg. Med. Chem. Lett. 24 (23), 5385–5389. doi:10.1016/j.bmcl.2014.10.049
Venditti, P., Di Stefano, L., and Di Meo, S. (2013). Mitochondrial Metabolism of Reactive Oxygen Species. Mitochondrion 13 (2), 71–82. doi:10.1016/j.mito.2013.01.008
Wang, M., and Xu, T. (2019). Methyl B12 Protects PC12 Cells against Cytotoxicity Induced by Aβ25-35. J. Cel. Biochem. 120 (7), 11921–11930. doi:10.1002/jcb.28475
Wang, Y., Huang, X., Liang, Q. H., Fan, R., Qin, F., Guo, Y., et al. (2012). A Strategy for Detecting Absorbed Bioactive Compounds for Quality Control in the Water Extract of Rhubarb by Ultra Performance Liquid Chromatography with Photodiode Array Detector. Chin. J. Integr. Med. 18 (9), 690–698. doi:10.1007/s11655-012-1053-7
Wang, J., Zhao, C., Zhao, A., Li, M., Ren, J., and Qu, X. (2015). New Insights in Amyloid Beta Interactions with Human Telomerase. J. Am. Chem. Soc. 137 (3), 1213–1219. doi:10.1021/ja511030s
Wang, Y., Fan, X., Tang, T., Fan, R., Zhang, C., Huang, Z., et al. (2016). Rhein and Rhubarb Similarly Protect the Blood-Brain Barrier after Experimental Traumatic Brain Injury via Gp91phox Subunit of NADPH oxidase/ROS/ERK/MMP-9 Signaling Pathway. Sci. Rep. 6, 37098. doi:10.1038/srep37098
Wang-Dietrich, L., Funke, S. A., Kühbach, K., Wang, K., Besmehn, A., Willbold, S., et al. (2013). The Amyloid-β Oligomer Count in Cerebrospinal Fluid Is a Biomarker for Alzheimer's Disease. J. Alzheimers Dis. 34 (4), 985–994. doi:10.3233/jad-122047
Xiao, X., Chen, Q., Zhu, X., and Wang, Y. (2019). ABAD/17β-HSD10 Reduction Contributes to the Protective Mechanism of Huperzine a on the Cerebral Mitochondrial Function in APP/PS1 Mice. Neurobiol. Aging 81, 77–87. doi:10.1016/j.neurobiolaging.2019.05.016
Yan, S., Du, F., Wu, L., Zhang, Z., Zhong, C., Yu, Q., et al. (2016). F1F0 ATP Synthase-Cyclophilin D Interaction Contributes to Diabetes-Induced Synaptic Dysfunction and Cognitive Decline. Diabetes 65 (11), 3482–3494. doi:10.2337/db16-0556
Ye, J., Liu, Z., Wei, J., Lu, L., Huang, Y., Luo, L., et al. (2013). Protective Effect of SIRT1 on Toxicity of Microglial-Derived Factors Induced by LPS to PC12 Cells via the P53-caspase-3-dependent Apoptotic Pathway. Neurosci. Lett. 553, 72–77. doi:10.1016/j.neulet.2013.08.020
Zahs, K. R., and Ashe, K. H. (2013). β-Amyloid Oligomers in Aging and Alzheimer's Disease. Front. Aging Neurosci. 5, 28. doi:10.3389/fnagi.2013.00028
Zeng, K.-W., Wang, J.-K., Wang, L.-C., Guo, Q., Liu, T.-T., Wang, F.-J., et al. (2021). Small Molecule Induces Mitochondrial Fusion for Neuroprotection via Targeting CK2 without Affecting its Conventional Kinase Activity. Sig Transduct Target. Ther. 6 (1), 71. doi:10.1038/s41392-020-00447-6
Zhao, Y., and Zhao, B. (2013). Oxidative Stress and the Pathogenesis of Alzheimer's Disease. Oxid. Med. Cell Longev. 2013, 1–10. doi:10.1155/2013/316523
Zhao, Y., Huang, Y., Fang, Y., Zhao, H., Shi, W., Li, J., et al. (2018). Chrysophanol Attenuates Nitrosative/oxidative Stress Injury in a Mouse Model of Focal Cerebral Ischemia/reperfusion. J. Pharmacol. Sci. 138 (1), 16–22. doi:10.1016/j.jphs.2018.08.002
Keywords: β-amyloid, oxidative stress, rhein, mitochondrial biogenesis, sirtuin 1, peroxisome proliferator-activated receptor gamma coactivator 1-alpha
Citation: Yin Z, Geng X, Zhang Z, Wang Y and Gao X (2021) Rhein Relieves Oxidative Stress in an Aβ1-42 Oligomer-Burdened Neuron Model by Activating the SIRT1/PGC-1α-Regulated Mitochondrial Biogenesis. Front. Pharmacol. 12:746711. doi: 10.3389/fphar.2021.746711
Received: 24 July 2021; Accepted: 31 August 2021;
Published: 10 September 2021.
Edited by:
Fushun Wang, Nanjing University of Chinese Medicine, ChinaReviewed by:
Ke-Wu Zeng, Peking University Health Science Centre, ChinaMiao Qu, Capital Medical University, China
Copyright © 2021 Yin, Geng, Zhang, Wang and Gao. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ying Wang, wangy174@126.com; Xiaoyan Gao, gaoxiaoyan@bucm.edu.cn