Melatonin-Mediated Intracellular Insulin during 2-Deoxy-d-glucose Treatment Is Reduced through Autophagy and EDC3 Protein in Insulinoma INS-1E Cells
2-DG triggers glucose deprivation without altering other nutrients or metabolic pathways and then activates autophagy via activation of AMPK and endoplasmic reticulum (ER) stress. We investigated whether 2-DG reduced intracellular insulin increased by melatonin via autophagy/EDC3 in insulinoma INS-1E cells. p-AMPK and GRP78/BiP level were significantly increased by 2-DG in the presence/absence of melatonin, but IRE1α level was reduced in 2-DG treatment. Levels of p85α, p110, p-Akt (Ser473, Thr308), and p-mTOR (Ser2481) were also significantly reduced by 2-DG in the presence/absence of melatonin. Mn-SOD increased with 2-DG plus melatonin compared to groups treated with/without melatonin alone. Bcl-2 was decreased and Bax increased with 2-DG plus melatonin. LC3II level increased with 2-DG treatment in the presence/absence of melatonin. Intracellular insulin production increased in melatonin plus 2-DG but reduced in treatment with 2-DG with/without melatonin. EDC3 was increased by 2-DG in the presence/absence of melatonin. Rapamycin, an mTOR inhibitor, increased GRP78/BiP and EDC3 levels in a dose-dependent manner and subsequently resulted in a decrease in intracellular production of insulin. These results suggest that melatonin-mediated insulin synthesis during 2-DG treatment involves autophagy and EDC3 protein in rat insulinoma INS-1E cells and subsequently results in a decrease in intracellular production of insulin.
Pancreatic β-cells are the only source of the peptide hormone insulin, which is essential for stimulating glucose uptake in peripheral tissue and inhibiting glucose production in liver. Pancreatic β-cell dysfunction and death result from a combination of genetic predisposition and exposure to environmental risk factors and contribute to type 2 diabetes (T2D) . The endoplasmic reticulum (ER) contains the protein-folding machinery for secretory proteins and is therefore crucial for insulin biosynthesis. Disruption of β-cell ER function by elevated free fatty acid levels and insulin resistance can result in an imbalance in protein homeostasis and ER stress, which has been recognized as an important mechanism for T2D [2, 3]. Recent reports have shown that autophagy is activated in response to ER stress in pancreatic β-cells [3–5].
AMP-activated protein kinase (AMPK) is a central metabolic regulator which balances energy production through the AMP : ATP ratio in eukaryotic cells and functions as a therapeutic target for anticancer and metabolic diseases in diabetes mellitus and lipid metabolism [6–8]. AMPK appears to be involved in insulin signaling pathways by positively regulating Akt/Protein Kinase B (PKB) in endothelial cells and cardiomyocytes [9, 10]. However, AMPK activation has also been reported to negatively affect Akt/PKB and insulin production . Activation or phosphorylation of AMPK is associated with activation or phosphorylation of Akt/PKB and mammalian target of rapamycin (mTOR) for insulin synthesis or secretion in pancreatic β-cells [12, 13].
Melatonin can control insulin secretion both in vivo and in vitro via MT1 or MT2 receptors [14, 15], Gi protein , cyclic AMP-response element-binding protein (CREB) , and Gq proteins/phospholipase C/IP3 . Circadian rhythms are influenced by melatonin and also induce a phase shift in insulin secretion . Melatonin directly influences diabetes and associated metabolic diseases . The genes and proteins involved in direct regulation of insulin biosynthesis and secretion in pancreatic β-cells have been reported: IRE1α , Sirt1 , and Sirt4 . Recently, melatonin was shown to directly influence insulin biosynthesis and secretion under ER stress via the RNA-binding protein human antigen D (HuD) in rat pancreatic INS-1E cells, suggesting that nuclear HuD may be an important negative regulator in insulin production and secretion . However, an interaction between the enhancer of mRNA decapping complex 3 (EDC3)  and insulin production has not been reported.
The glucose analogue 2-deoxy-D-glucose (2-DG) is metabolized by hexokinase and acts as an inhibitor of glycolysis . 2-DG triggers glucose deprivation without altering other nutrients or metabolic pathways  and then activates autophagy via activation of AMPK  and endoplasmic reticulum (ER) stress . Therefore, 2-DG appears to be an ideal tool to understand the interactions between autophagy and ER stress. In this study, we provide the first demonstration that 2-DG reduces intracellular insulin which was increased by melatonin via autophagy and EDC3 in insulinoma INS-1E cells.
2. Materials and Methods
2.1. Cell Culture
INS-1E cells, a clonal pancreatic β-cell line, were obtained from Professor Claes B. Wollheim and were cultured in RPMI 1640 medium (Invitrogen, Carlsbad, CA, USA) containing 11 mM glucose supplemented with 10 mM HEPES (pH 7.3), 10% heat-inactivated fetal bovine serum (FBS; Invitrogen), 50 μM β-mercaptoethanol, 1 mM sodium pyruvate, 50 μg/mL penicillin, and 100 μg/mL streptomycin at 37°C with 5% CO2 in a humidified incubator.
2.2. Treatment with 2-DG, Melatonin, and Rapamycin
INS-1E cells were cultured in RPMI 1640 medium plus 2% heat-inactivated FBS in a 37°C and 5% CO2 incubator, with or without melatonin (Sigma-Aldrich, St. Louis, MO, USA) and/or 2-DG (5 mM) (Sigma-Aldrich) for 24 hr, or with rapamycin (20 to 80 nM) (Calbiochem, San Diego, MO, USA) for 24 hr.
2.3. Western Blot Analysis
Cells were harvested, washed twice with ice-cold phosphate buffered saline (PBS), and then resuspended in 20 mM Tris-HCl buffer (pH 7.4) containing protease inhibitors (0.1 mM phenylmethylsulfonyl fluoride, 5 μg/mL aprotinin, 5 μg/mL pepstatin A, and 1 μg/mL chymostatin) and phosphatase inhibitors (5 mM Na3VO4 and 5 mM NaF). Whole cell lysate was prepared using 20 strokes of a Dounce homogenizer, followed by centrifugation at 13,000 ×g for 20 min at 4°C. The protein concentration was determined using the BCA assay (Sigma). Proteins (40 μg) were separated by 12% sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) and the detection of insulin proteins was performed by 16.5% tricine SDS-PAGE. These gels transferred onto a polyvinylidene difluoride (PVDF) membrane. The membrane was incubated with antibodies (diluted as indicated in brackets) directed against the following proteins: p-AMPK and AMPK (1 : 500) (Santa Cruz Biotechnology); IRE1α (1 : 500) (Santa Cruz Biotechnology); p-PERK (1 : 1000) (Cell Signaling Technology, Beverly, MA, USA); GRP78/BiP (1 : 1000) (Cell Signaling Technology); p85α and p85β (1 : 500) (Santa Cruz Biotechnology); p110 (1 : 500) (Santa Cruz Biotechnology); p-Akt (T308, S473) and Akt (1 : 1000) (Cell Signaling Technology); p-mTOR (Ser2448, 2481) and mTOR (1 : 1000) (Cell Signaling Technology); Cu/Zn-SOD (1 : 500) (Santa Cruz Biotechnology); Mn-SOD (1 : 500) (Santa Cruz Biotechnology); catalase (1 : 500) (Santa Cruz Biotechnology); Bcl-2 (1 : 500) (Santa Cruz Biotechnology); Bax (1 : 500) (Santa Cruz Biotechnology); insulin (1 : 500) (Santa Cruz Biotechnology); EDC3 (1 : 500) (Santa Cruz Biotechnology); and Actin (1 : 1000) (Assay Designs, Ann Arbor, MI, USA). Immunoreactive proteins were visualized by exposure to X-ray film. Protein bands were analyzed by image-scanning, and optical density was measured using ImageJ analysis software (version 1.37, Wayne Rasband, NIH, Bethesda, MD, USA). The data were corrected for background subtraction and normalized by including Actin as an internal control.
2.4. Statistical Analysis
Significant differences were detected by ANOVA, followed by Tukey’s test for multiple comparisons. Analysis was performed using the Prism Graph Pad v4.0 (Graph Pad Software Inc., San Diego, CA, USA). Values are expressed as means ± SD of at least three separated experiments, in which case a representative experiment is depicted in the figures. values < 0.05 were considered statistically significant.
3.1. Expressions of p-AMPK, GRP78/BiP, IRE1α, and p-PERK Proteins
We investigated whether 2-DG reduced insulin synthesis which had been increased by melatonin via autophagy in rat insulinoma INS-1E cells. 2-DG activates autophagy via reactive oxygen species-mediated activation of AMPK  and ER stress . Therefore, we first investigated whether 2-DG and/or melatonin induced phosphorylation of AMPK and ER stress in rat insulinoma INS-1E cells. Phosphorylation of AMPK was significantly increased by treatment with 2-DG in the presence or absence of melatonin, compared to treatment with melatonin or FBS only (Figures 1(a) and 1(b)). ER stress marker GRP78/BiP protein was also significantly increased by treatment with 2-DG in the presence or absence of melatonin (Figures 1(c) and 1(d)), but expression of IRE1α protein was reduced (Figures 1(c) and 1(d)) and p-PERK protein expression was not significantly affected by 2-DG treatment (Figures 1(c) and 1(e)).
3.2. Expressions of p-PI3K, Akt/PKB, and mTOR Proteins
We determined how 2-DG and/or melatonin affected known regulators of insulin signaling, phosphoinositide 3-kinase (PI3K), Akt/PKB, and mTOR, in rat insulinoma INS-1E cells. Class IA PI3K is composed of a heterodimer of p110 catalytic subunit and p85 regulatory subunit . 2-DG treatment with or without melatonin caused a significant reduction in expression of PI3K subunit p85α and p110 and a significant increase in p85β protein expression (Figures 2(a)–2(d)). 2-DG also significantly decreased expression of p-Akt (Ser473, Thr308) (Figures 3(a)–3(c)) and p-mTOR (Ser2481) (Figures 3(d)–3(f)) proteins in the presence or absence of melatonin. These results indicated that ER stress in the presence of 2-DG or 2-DG-plus melatonin negatively regulates the PI3K/Akt/mTOR pathway.
3.3. Expressions of SOD, Catalase, Bcl-2, and Bax Proteins
Next, to characterize intracellular and mitochondrial conditions according to antioxidant and prooxidant proteins, we examined how 2-DG and/or melatonin affected superoxide dismutase (SOD), catalase, Bcl-2, and Bax proteins in rat insulinoma INS-1E cells. Cu/Zn-SOD was not significantly affected (Figures 4(a) and 4(b)), but Mn-SOD protein expression was increased by treatment with 2-DG plus melatonin compared to groups treated with/without melatonin alone or the nontreated control group (Figures 4(a) and 4(c)). Catalase protein expression was decreased by 2-DG treatment in the presence of melatonin (Figures 4(a) and 4(d)). Bcl-2 protein expression was decreased by treatment with 2-DG plus melatonin compared to treatment with/without melatonin (Figures 4(e) and 4(f)), and Bax protein expression was increased in cells treated with 2-DG and/or melatonin, compared to nontreated control group (Figures 4(e) and 4(g)). These results suggested that 2-DG or 2-DG-plus melatonin influenced mitochondria, thus resulting in a change of antioxidant and prooxidant protein levels.
3.4. Expressions of LC3II, Insulin, and EDC3 Proteins
Subsequently, we investigated whether 2-DG reduced insulin synthesis which had been increased by melatonin via EDC3 in rat insulinoma INS-1E cells. Expression of autophagy markers LC3 II was increased by 2-DG treatment in the presence or absence of melatonin (Figures 5(a) and 5(b)). Insulin synthesis was increased by melatonin treatment compared to nontreated melatonin group, but it was reduced by treatment with 2-DG (Figures 5(a) and 5(c)). Expression of EDC3, a negative regulator in insulin synthesis, was increased by 2-DG in the presence or absence of melatonin (Figures 5(a) and 5(d)), suggesting that the decrease in intracellular insulin synthesis was associated with an increase of autophagy and EDC3 protein expression.
3.5. The Relationship of Insulin Production with Autophagy and EDC3 Protein
Lastly, to confirm the relationship of insulin production with autophagy and EDC3 protein, we investigated whether the mTOR inhibitor, rapamycin, induced autophagy and expression of GRP78/BiP and EDC3 proteins in rat insulinoma INS-1E cells. Rapamycin treatment resulted in an increase in GRP78/BiP (Figures 6(a) and 6(c)) and EDC3 proteins in a dose-dependent manner (Figures 6(a) and 6(d)) and a subsequent decrease in insulin production (Figures 6(a) and 6(b)). These results suggested that melatonin-mediated insulin synthesis during 2-DG treatment involved autophagy-induced ER stress and EDC3 protein in rat insulinoma INS-1E cells, subsequently resulting in a decrease in insulin protein biosynthesis.
ER dysfunction has been implicated in insulin resistance and is an important cause of T2D [2–4]. Autophagy in ER stress-induced pancreatic β cells plays a role as an important regulator of insulin production and secretion [2–5, 13]. In addition, autophagy in pancreatic β cells is essential to maintain normal morphology, cell mass, and function of β-cells and is regarded as a crucial stress response factor to protect β-cells under insulin-resistant states [4, 5]. The glucose analog 2-DG acts as a glycolytic inhibitor, leading to ER stress and an unfolded protein response [26–28]. Although 2-DG-induced ER stress has been shown to activate autophagy, the precise mechanism is not fully understood. Here we demonstrated the relevance of ER stress-mediated autophagy caused by 2-DG for insulin synthesis in pancreatic β cells. Our previous report showed that ER stress caused by melatonin, especially in the presence of thapsigargin, decreased intracellular insulin biosynthesis and that extracellular secretion of insulin may be regulated by melatonin in rat insulinoma INS-1E cells . Here we found that melatonin reduced insulin production in the presence of 2-DG via autophagy-induced ER stress in rat insulinoma INS-1E cells (Figures 1(c)–1(f) and 6(a)–6(d)). mTOR inhibitor, rapamycin, induced autophagy and then subsequently decrease in insulin production (Figure 6). Therefore, autophagy can regulate insulin production at the cellular level through physiological and biochemical changes of ER homeostasis involving the unfolded protein response [30, 31].
GRP78/BiP is localized in the ER and its expression is increased by environmental stresses in many types of cells. Other studies have showed that conditions including glucose deprivation, 2-DG treatment, and hypoxia or increase of intracellular Ca2+ concentration induce GRP78/BiP expression [32–35]. Our previous research demonstrated that expression of the ER-stress protein GRP78/BiP is significantly decreased by tunicamycin/melatonin treatment and remains almost unperturbed with thapsigargin/melatonin treatment when compared to control/melatonin in rat pancreatic INS-1E cells . This suggests that GRP78/BiP protein can be regulated by several cellular stresses which perturb ER function and homeostasis, including tunicamycin, thapsigargin, and calcium ionophore A23187 . 2-DG is an especially potent inducer of GRP78/BiP protein expression . In this study we found that melatonin significantly regulated insulin production through the expression of GRP78/BiP protein which induced ER stress under conditions of 2-DG or 2-DG plus melatonin treatment in rat insulinoma INS-1E cells (Figures 1(c) and 1(f)). Therefore, GRP78/BiP accompanied with ER stress is a crucial element in insulin biosynthesis and secretion in rat pancreatic β-cells .
The interaction between 2-DG and the AMPK/mTOR signal pathway has not been reported in pancreatic β cells. Activation of AMPK and inactivation of mTOR signal pathway have only been studied under conditions of nutrient deprivation, including glucose and/or leucine starvation, in pancreatic beta-cells, rat skeletal muscle cells, and hepatocytes [12, 13, 39, 40]. Gleason et al.  reported that because glucose and amino acids stimulate insulin release from pancreatic β cells, AMPK represents an important factor for control of β-cell function. In addition, the activation/phosphorylation of AMPK in β-cells by low glucose correlates with inactivation of the mTOR pathway as a cellular sensor for nutritional conditions in β cells [12, 41]. The other condition, that is, lack of amino acid leucine, causes an increase in the activation/phosphorylation of AMPK and a decrease in the activation of mTOR pathway [12, 39, 40]. It is therefore unlikely that the ability of AMPK and the mTOR signal pathway to sense both glucose and amino acids plays a role in regulation of insulin synthesis or secretion [12, 41]. 2-DG induces glucose deprivation and results in the activation/phosphorylation of AMPK [26, 27]. In our study, phosphorylation of AMPK was significantly increased (Figures 1(a) and 1(b)) and phosphorylation of mTOR (Ser2481) was significantly decreased (Figures 3(d) and 3(f)) by 2-DG in the presence or absence of melatonin, showing that β-cells do indeed possess a mechanism by which glucose derivation can suppress the mTOR signal pathway via activation of AMPK. This result suggests that a connection between AMPK and mTOR in the β-cell is of particular relevance to insulin synthesis/secretion during 2-DG treatment. The finding that AMPK and the mTOR signal pathway can sense both glucose and amino acids may apply to 2-DG treatment conditions in the regulation of insulin synthesis or its secretion from β-cells.
The genes/proteins involved in insulin biosynthesis and secretion in pancreatic β-cells have been studied: IRE1α , Sirt1 , Sirt4 , and HuD . The positive regulators of insulin biosynthesis and secretion are IRE1α, Sirt1, and WSF1, whereas Eny2, Sirt4, and HuD are negative regulators. Our recent study suggests that nuclear HuD under ER stress conditions with thapsigargin and melatonin functions as an important negative regulator of insulin production from rat pancreatic INS-1E cells . However, our current study demonstrated that HuD under autophagy and ER stress with 2-DG and melatonin could not negatively regulate insulin production (data not shown). Instead, EDC3 was involved in insulin synthesis during 2-DG/melatonin or rapamycin treatment and subsequently induced a decrease in insulin protein (Figures 5 and 6). An interaction between EDC3 and insulin production has not yet reported . Our study is the first to suggest that EDC3 functions as a negative regulator of insulin production in rat pancreatic INS-1E cells.
Glucose deprivation by mechanisms including 2-DG induces the production of reactive oxygen species (ROS) to activate AMPK in pancreatic β cells . Hypoxia and hyperoxia trigger AMPK activation via mitochondrial ROS production . These studies support a role for AMPK as an important redox sensor through a mitochondrial ROS mechanism. We showed that 2-DG and/or melatonin affected SOD, catalase, Bcl-2, and Bax proteins in rat insulinoma INS-1E cells (Figure 4), suggesting that 2-DG or 2-DG plus melatonin influences antioxidant and prooxidant proteins through a mechanism involving mitochondrial ROS.
Class IA PI3K is composed of a p110 catalytic subunit and a p85 regulatory subunit and plays a pivotal role in insulin signaling [29, 44]. We found that 2-DG treatment significantly reduced expression of both PI3K subunits in the presence or absence of melatonin (Figure 2). This result indicated that the p85 proteins may bind and stabilize the p110 subunit .
We report that melatonin may directly regulate the intracellular biosynthesis of insulin in the presence of 2-DG in rat insulinoma INS-1E cells. Intracellular insulin production may be partially regulated by the expression of ER stress GRP78/BiP-induced autophagy and EDC3 protein under conditions of 2-DG treatment.
These results suggest that melatonin-mediated insulin synthesis during 2-DG treatment involves autophagy and EDC3 protein in rat insulinoma INS-1E cells and subsequently results in a decrease in intracellular production of insulin.
The authors declare that there is no conflict of interests regarding the publication of this paper.
This research was supported by the Leading Foreign Research Institute Recruitment Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Science, ICT & Future Planning (2010-00757).
S. L. Longnus, C. Ségalen, J. Giudicelli, M. P. Sajan, R. V. Farese, and E. Van Obberghen, “Insulin signalling downstream of protein kinase B is potentiated by 5′AMP-activated protein kinase in rat hearts in vivo,” Diabetologia, vol. 48, no. 12, pp. 2591–2601, 2005.View at: Publisher Site | Google Scholar
L. Bertrand, A. Ginion, C. Beauloye et al., “AMPK activation restores the stimulation of glucose uptake in an in vitro model of insulin-resistant cardiomyocytes via the activation of protein kinase B,” American Journal of Physiology—Heart and Circulatory Physiology, vol. 291, no. 1, pp. H239–H250, 2006.View at: Publisher Site | Google Scholar
E. Mühlbauer, E. Albrecht, I. Bazwinsky-Wutschke, and E. Peschke, “Melatonin influences insulin secretion primarily via MT(1) receptors in rat insulinoma cells (INS-1) and mouse pancreatic islets,” Journal of Pineal Research, vol. 52, no. 4, pp. 446–459, 2012.View at: Publisher Site | Google Scholar
E. Mühlbauer, E. Albrecht, K. Hofmann, I. Bazwinsky-Wutschke, and E. Peschke, “Melatonin inhibits insulin secretion in rat insulinoma β-cells (INS-1) heterologously expressing the human melatonin receptor isoform MT2,” Journal of Pineal Research, vol. 51, no. 3, pp. 361–372, 2011.View at: Publisher Site | Google Scholar
E. Peschke, E. Mühlbauer, U. Mußhoff, V. J. Csernus, E. Chankiewitz, and D. Peschke, “Receptor (MT1) mediated influence of melatonin on cAMP concentration and insulin secretion of rat insulinoma cells INS-1,” Journal of Pineal Research, vol. 33, no. 2, pp. 63–71, 2002.View at: Publisher Site | Google Scholar
I. Bazwinsky-Wutschke, S. Wolgast, E. Mühlbauer, E. Albrecht, and E. Peschke, “Phosphorylation of cyclic AMP-response element–binding protein (CREB) is influenced by melatonin treatment in pancreatic rat insulinoma β-cells (INS-1),” Journal of Pineal Research, vol. 53, no. 4, pp. 344–357, 2012.View at: Publisher Site | Google Scholar
C. B. Wollheim and T. J. Biden, “Second messenger function of inositol 1,4,5-trisphosphate. Early changes in inositol phosphates, cytosolic Ca2+, and in insulin release in carbamylcholine-stimulated RINm5F cells,” The Journal of Biological Chemistry, vol. 261, no. 18, pp. 8314–8319, 1986.View at: Google Scholar
G. Boden, J. Ruiz, J.-L. Urbain, and X. Chen, “Evidence for a circadian rhythm of insulin secretion,” American Journal of Physiology—Endocrinology and Metabolism, vol. 271, no. 2, part 1, pp. E246–E252, 1996.View at: Google Scholar
A. Sole and R. K. Crane, “The inhibition of brain hexokinase by adenosinediphosphate and sulfhydryl reagents,” Journal of Biological Chemistry, vol. 206, no. 2, pp. 925–936, 1954.View at: Google Scholar
W. Jiang, Z. Zhu, and H. J. Thompson, “Modulation of the activities of AMP-activated protein kinase, protein kinase B, and mammalian target of rapamycin by limiting energy availability with 2-deoxyglucose,” Molecular Carcinogenesis, vol. 47, no. 8, pp. 616–628, 2008.View at: Publisher Site | Google Scholar
G. Kuznetsov, K. T. Bush, P. L. Zhang, and S. K. Nigam, “Perturbations in maturation of secretory proteins and their association with endoplasmic reticulum chaperones in a cell culture model for epithelial ischemia,” Proceedings of the National Academy of Sciences of the United States of America, vol. 93, no. 16, pp. 8584–8589, 1996.View at: Publisher Site | Google Scholar
W. W. Li, Y. Hsiung, Y. Zhou, B. Roy, and A. S. Lee, “Induction of the mammalian GRP78/BiP gene by Ca2+ depletion and formation of aberrant proteins: activation of the conserved stress-inducible grp core promoter element by the human nuclear factor YY1,” Molecular and Cellular Biology, vol. 17, no. 1, pp. 54–60, 1997.View at: Publisher Site | Google Scholar
Z. Yu, H. Luo, W. Fu, and M. P. Mattson, “The endoplasmic reticulum stress-responsive protein GRP78 protects neurons against excitotoxicity and apoptosis: suppression of oxidative stress and stabilization of calcium homeostasis,” Experimental Neurology, vol. 155, no. 2, pp. 302–314, 1999.View at: Publisher Site | Google Scholar
L. Zhang, E. Lai, T. Teodoro, and A. Volchuk, “GRP78, but not protein-disulfide isomerase, partially reverses hyperglycemia-induced inhibition of insulin synthesis and secretion in pancreatic β-cells,” The Journal of Biological Chemistry, vol. 284, no. 8, pp. 5289–5298, 2009.View at: Publisher Site | Google Scholar
Y. Cai, G. A. Martens, S. A. Hinke, H. Heimberg, D. Pipeleers, and M. Van de Casteele, “Increased oxygen radical formation and mitochondrial dysfunction mediate beta cell apoptosis under conditions of AMP-activated protein kinase stimulation,” Free Radical Biology and Medicine, vol. 42, no. 1, pp. 64–78, 2007.View at: Publisher Site | Google Scholar