- Research article
- Open Access
Signaling dynamics of palmitate-induced ER stress responses mediated by ATF4 in HepG2 cells
BMC Systems Biology volume 7, Article number: 9 (2013)
Palmitic acid, the most common saturated free fatty acid, has been implicated in ER (endoplasmic reticulum) stress-mediated apoptosis. This lipoapotosis is dependent, in part, on the upregulation of the activating transcription factor-4 (ATF4). To better understand the mechanisms by which palmitate upregulates the expression level of ATF4, we integrated literature information on palmitate-induced ER stress signaling into a discrete dynamic model. The model provides an in silico framework that enables simulations and predictions. The model predictions were confirmed through further experiments in human hepatocellular carcinoma (HepG2) cells and the results were used to update the model and our current understanding of the signaling induced by palmitate.
The three key things from the in silico simulation and experimental results are: 1) palmitate induces different signaling pathways (PKR (double-stranded RNA-activated protein kinase), PERK (PKR-like ER kinase), PKA (cyclic AMP (cAMP)-dependent protein kinase A) in a time dependent-manner, 2) both ATF4 and CREB1 (cAMP-responsive element-binding protein 1) interact with the Atf4 promoter to contribute to a prolonged accumulation of ATF4, and 3) CREB1 is involved in ER-stress induced apoptosis upon palmitate treatment, by regulating ATF4 expression and possibly Ca2+ dependent-CaM (calmodulin) signaling pathway.
The in silico model helped to delineate the essential signaling pathways in palmitate-mediated apoptosis.
Elevated serum free fatty acids (FFA) and hepatic lipid accumulation in non-adipose tissues can lead to cellular dysfunction or cell death, due in part to the diversion of unoxidized FFAs to nonoxidative pathways, resulting in lipoapoptosis . Excess amounts of nonesterified FFAs that fail to convert to triglyceride in liver cells, enhance the risk for hepatocellular lipoapoptosis, a pathogenic feature observed in non-alcoholic steatohepatitis (NASH) . The mechanisms involved in FFA-induced toxicity have remained unresolved, however recent studies suggest that hepatic lipoapoptosis arises predominantly from FFA-induced lipotoxic stress of intracellular organelles, in particular the endoplasmic reticulum (ER) and mitochondria .
The ER is one of the largest organelles in the cells, and perturbing ER homeostasis or inducing ER stress has profound effects on cell survival. Cellular perturbations, such as alterations in calcium storage in the ER lumen or an imbalance in the luminal-oxidizing environment will cause ER stress. This stress is sensed by the cells through three ER transmembrane proteins, inositol requiring enzyme (IRE) 1α, PKR-like ER kinase (PERK), and activating transcription factor (ATF) 6α . They activate signaling processes to restore ER homeostasis, and are collectively termed the unfolded protein response (UPR). UPR signaling pathways coordinate cellular response by down-regulating protein translation, enhancing expression of ER chaperone proteins that promote protein refolding, and activating proteases involved in the degradation of misfolded proteins. When these corrective actions are insufficient to attenuate ER-stress, the UPR switches to a pro-apoptotic mode. However, the mechanisms regulating ER-stress-induced apoptosis have not been clearly defined.
Of the ER transmembrane proteins, PERK is involved in regulating both translation and transcription to return the folding demand and capacity to homeostasis. Like PKR (double-stranded RNA-activated protein kinase, also known as eIF2αk2), one of the eukaryotic translation inhibition factor 2α (eIF2α) kinases, PERK (eIF2ak3) also phosphorylates eIF2α at Ser 51, and phosphorylated eIF2α impedes global translation initiation to decrease the protein load in the ER . However, activated eIF2α paradoxically favors an increase in the translation of the activating transcription factor-4 (ATF4), a member of the activating transcription factor/the cAMP-responsive element-binding protein (ATF/CREB) family of basic zipper-containing proteins . ATF4 binds to cAMP response element (CRE; TGACGT(C/A)(G/A)) as a homodimer to induce transcription of many pro-survival genes that are involved in amino-acid metabolism, redox reactions, stress response and protein secretion [7, 8]. Under sustained ER-stress, ATF4 activates CHOP (C/EBP-homologous protein; also known as GADD153) transcription by binding to the CCAAT/enhancer binding protein (C/EBP)/ATF response element (CARE) sequences, likely as heterodimers with members of the C/EBP family transcription factor. CHOP in turn inhibits the expression of anti-apoptotic Bcl-2 (B cell lymphoma 2) protein [9, 10]. Concomitantly, CHOP forms a heteromeric complex with the phosphorylated c-Jun to bind the PUMA (p53-upregulated modulator of apoptosis) promoter and contributes to the upregulation of pro-apoptotic proteins .
Saturated long chain-FFAs, especially palmitic acid, have been implicated in ER stress-induced apoptosis in liver cells [12, 13]. Since both transcriptional regulation and translational control of ATF4 plays a central role on ER stress-induced apoptosis, recent studies have focused on the regulation of the ATF4 gene mediated by palmitate. In H4IIE liver cells, a rat liver hepatoma cell line, palmitate has been shown to induce ER stress, as evident by increased mRNA levels of ATF4 [12–14]. Similarly human hepatocarcinoma cell lines, such as HepG2 cells, respond also to palmitate by elevating ATF4 mRNA and inducing ER stress . However, with the highly dynamic nature of UPR signaling, it would be judicious to integrate signaling pathways which are involved in the regulation of ATF4 and further elucidate the dynamics of the main signaling axis for ER stress-induced apoptosis, PERK/eIF2α/ATF4. The complexity of the signaling pathways involved makes it difficult to experimentally test all possible interactions to determine which pathways in the signaling network are functional for a given condition. Therefore, we applied a discrete dynamic network model to integrate potential interactions and components of ER signaling and its feedbacks in HepG2 cells to gain a better understanding of how palmitate induces ER stress in liver cells. The predictions are confirmed with further experiments.
The conventional signaling pathway (PERK/eIF2α/ATF4) cannot explain the prolonged ATF4-dependent ER stress induced by palmitate in HepG2 cells
To evaluate the dynamics of the conventional ER signaling pathway, we measured the protein expression levels of PERK, eIF2α, and ATF4 at 3, 6, and 24 hrs upon palmitate treatment. As shown in Figure 1A, the phosphorylated protein levels of PERK at Ser 713 and eIF2α at Ser 51 were signficantly increased by palmitate treatment at 6, 24 hr and 3, 6 hr, respectively, while their total protein levels remained constant. Palmitate transiently increased the phosphorylation levels of eIF2α but induced a sustained increase in the protein expression level of ATF4. We note several discrepancies from the conventional signaling pathway. First, although eIF2α is activated transiently by palmitate, ATF4 continued to be strongly induced by palmitate, i.e. remained stably upregulated through 24 hrs (Figure 1B). In support of our results, a recent study also observed similar trends in a time-course analysis of eIF2α phosphorylation and ATF4 protein expression in HepG2 cells treated with thapsigargine, an ER calcium-ATPase pump inhibitor that induces calcium release from the ER . The eIF2α phosphorylation was transientlty activated by thapsigarine and reached a maxium at 30 min, while ATF4 induction continued to rise for up to 8 hrs. However, it remains unclear why ATF4 remained stably upregulated even though the level of eIF2α phosphorylation was only transiently activated.
In addition, our group recently found that cAMP promoted apoptosis in palmitate-treated HepG2 cells, suggesting that cAMP-CREB1 signaling may be involved in palmitate-mediated ER stress . Furthermore, it has been shown that the phosphorylation of CREB1 can be up-regulated by thapsigarine (ER stress inducer) in human glioma cells . Since ATF4 (also known as CREB2) belongs to a family of ATF/CREB proteins, it raises a possibility that crosstalk exists between these two proteins during palmitate-induced ER stress. Therefore, we assessed whether the phosphorylation of CREB1 in HepG2 cells is up-regulated upon palmitate treatment. As shown in Figure 2, the phosphorylation of CREB1 was significantly induced by palmitate at all time points (3, 6, 24 hr).
To resolve these discrepancies from the conventional pathway and incorporate our findings with the current knowledge, we collected and integrated information from the literature on palmitate-induced ER stress signaling, with an emphasis on how ATF4 could be activated by palmitate. We constructed a network of palmitate-induced signaling processes by collating the literature information (Figure 3). We formularized the signaling network into a discrete dynamic model. Our discrete dynamic modeling approach was previously proposed and successfully applied to dissect the insulin signaling pathways in the same human liver cell system . We applied this modeling and simulation method to the palmitate-treated cells (see Methods), and integrated current knowledge to generate a signaling network, which was probed in silico to provide insight into the regulatory mechanisms involved.
PKR pathway is essential for eIF2α phosphorylation in palmitate
Our simulation results suggest plausible dynamic profiles of the network upon palmitate-stimulation (Figures 4A, 4B, 4C). The simulations are based on current knowledge of the regulatory interactions between the components in the network. We initially assume that the activation steps (mostly phosphorylation/de-phosphorylation) of the different components are at similar time scales. As shown in Figure 4A, CREB1 phosphorylation level was increased by palmitate over the simulation time, which matched the results obtained by western blotting analysis shown in Figure 2. However, the simulation results show that eIF2α and ATF4 were not activated by palmitate treatment (Figure 4A), which is inconsistent with the experimental results of Figure 1. The in silico results suggest an inconsistency with the current knowledge of the palmitate-induced signaling processes mediated by eIF2α and ATF4.
Since PKR, another eIF2α kinase, has been shown to be involved in ATF4-dependent ER stress in human embryonic kidney cells , we questioned whether PKR is involved in eIF2α phosphorylalion to regulate the pamitate-mediated ATF4 protein expression in HepG2 cells. The experimental results suggest that PKR shows a prompt response at 3 h (Figure 4D), while PERK takes longer to be activated (6 h) (Figure 1). In addition, since PKA (cyclic AMP (cAMP)-dependent protein kinase A) regulates CREB1 phosphorylation both directly and indirectly (through PP2A (nuclear protein phosphatase 2A)) [21, 22], we measured the protein expression level of PKA upon palmitate treatment. The catalytic subunits of PKA, which plays an important role for PKA activation, were significantly upregulated at 24 hr (Additional file 1: Figure S1). Thus we concluded that palmitate initiates signaling at different times through PKR (3 hr), PERK (6 hr), and PKA (24 hr). We subsequently incorporated this new information into the discrete dynamic model, by imposing different delays in the palmitate-induced PERK and PKA activation. The resultant eIF2α profile (Figure 4B) matches the experimental observation (Figure 1) of a transient elevation in its phosphorylation level, which initiates the early ER stress response. The network model thus suggests an important role of the different response time of PKR vs. PERK in the regulation of the downstream effectors of palmitate. Thus we hypothize that the prompt response of PKR could initiate the signal that activates the eIF2α pathway.
CREB1 may be involved in ATF4-dependent ER stress
Another significant discrepancy between the in silico simulation and the experimental observation lies in the prolonged activation of ATF4. The experimental measurements (Figure 1) show that ATF4 level is higher (than control cells without palmitate treatment) at both 6 h and 24 h. Such a prolonged activation cannot be explained by the model simulation where the ATF4 level is reduced to lower than control at 24 h, although the initial upregualtion of ATF4 in response to the upstream eIF2α is captured by the model. The discrepancy suggests the downsteam feedback regulation of ATF4 in the current model is incorrect in our liver cell system. The downstream feedback regulation in the model is mediated by phosphoprotein phosphatase 1 (PP1), which is known as a major regulator of ATF4 . We measured the level of phosphorylated PP1 at different times upon PA treatment and found that, in contrast to current knowledge of the ATF4 pathway, there is no significant change on the activity of PP1 (Additional file 1: Figure S2) in our system that could affect the ATF4 level. The lack of involvement of PP1 explains in part the discrepancy between the current knowledge and our experimental observation of ATF4 activation. Indeed, when we updated our model with this new information of a constant PP1 level (i.e., to remove its impact on other components), the ATF4 profile is no longer “inhibited” in the simulation (Figure 4C). Nevertheless, the simulations are still unable to capture the prolonged activation of ATF4, because there is no other regulators that connect to ATF4 in the model to support its sustained activation after 6 h. This suggests that the current knowledge of the signaling process is incomplete and there should be other (currently unknown) regulatory relationship(s) in the network that could lead to the accumulation of ATF4 and finally lipoapoptosis.
CREB1 is of the same family as ATF4 and the phosphorylation of CREB1 was significantly enhanced upon palmitate treatment (Figure 2). Multiple CRE binding sites (TGACG or CGTCA) are identified on the ATF4 and CREB1 genes. Both CREB1 and ATF4 proteins could bind the putative CRE binding sites to enhance their gene expressions. Thus, we tested whether silencing either gene affects the protein expression level of ATF4 and CREB1. As shown in Figure 5A, CREB1 silencing significantly reduced the protein expression of ATF4, while ATF4 silencing did not significantly affect the protein expression level of CREB1. From these results, we hypothesized that both ATF4 and CREB1 regulate the gene expression of ATF4. This potential crosstalk (Figure 5B) between ATF4 and CREB1 pathway was integrated into our system model (Figure 3), and the simulation results (Figure 5C and 5D) predict a prolonged activation of ATF4 level, with either ATF4 binding to its own promoter, or CREB1 binding to the Atf4 promoter to induce ATF4 gene expression.
ATF4 protein binds to the CRE site on the Atf4 promoter in vitro
Since ChIP-chip data for pCREB1 are available, we determined whether the CREB1 protein interacts with the Atf4 promoter. The hepatocyte-specific ChIP-chip data for pCREB1 binding on Creb1 and Atf4 promoters were obtained from the CREB1 target gene database . A cut-off value of binding ratio ≥ 2 and confidence level p-value ≤ 0.001 was used to determine pCREB positive promoters (see ref  for details). ChIP-chip data showed that both ATF4 (binding ratio = 3.1, p-value = 4.1e-06) and CREB1 (binding ratio = 2.1, p-value = 8.3e-04) are positive for pCREB1 in hepatocytes, suggesting that the CREB1 protein directly binds on both the Creb1 and Atf4 promoters. In addition, to experimentally examine whether the ATF4 protein interacts with the Atf4 promoter, we designed an ATF4 probe (-175 to -147) which contains one CRE site (CGTCA; -164 to -160) and performed an EMSA assay. As shown in Figure 6, the ATF4 probe formed a DNA-protein complex with the nuclear extract from the palmitate-treated HepG2 cells (lane 2; marked as an arrow). The IgG antibody does not cause any significant change of the DNA-protein complex (lane 3), while the ATF4 antibody eliminated the complex (lane 4). Since the disappearance of a band by the addition of a specific antibody, but without a supershift, provides evidence of DNA-protein complex [23, 25], our results suggest that the ATF4 protein directly binds to the CRE site on the Atf4 promoter. Thus, the ChiP-chip data and our EMSA data confirmed our potential models in Figure 5C and 5D, supporting that both ATF4 and CREB1 proteins could regulate the ATF4 gene expression.
In-silico perturbation study suggests the essential signaling pathways for ATF4 and CREB1 activation mediated by palmitate
To identify the essential signaling pathways of ATF4 and CREB1 activated by palmitate, we simulated a model that incorporated our experimental findings and then evaluated the dynamics of eIF2α phosphorylation, ATF4 expression, and CREB1 phosphorylation, by systematically deleting specific pathways. We addressed whether PKA, p38, Ca2+ activation regulated by palmitate are essential for CREB1 phosphorylation. The in silico knock-out results (Additional file 1: Figure S3) showed that PKA-CREB1 and p38-CREB1 interactions did not alter the dynamics of the model. However, the deletion of Ca2+ signaling dramatically reduced CREB1 phosphorylation (marked by an arrow in Additional file 1: Figure S3). The model further showed that the PKA-, p38-, Ras-deficient model was able to capture the dynamic response of CREB1 and ATF4, as long as the Ca2+ signaling was able to induce CREB1 phosphorylation in response to the palmitate treatment. It has been shown that Ca2+ dependent-calmodulin (CaM) plays an important role in the activation of the CREB1 pathway [26, 27], and our simulation results suggest that the calcium signaling pathway should be an essential component in the signaling network in response to palmitate stimulation of the liver cells. In other words, calcium signaling is sufficient to activate CREB1 phosphorylation, without the other pathways, and to recapitulate the dynamic behavior observed experimentally.
Since PKR and PERK are two main upstream signaling components of ATF4, we further evaluated whether knock-out of either PKR or PERK perturbs the signaling dynamics of ATF4. Simulation results (Figure 7A) suggest that the differences between PKR and PEKR perturbations (inhibitions) show up at an early time step (within 10 time steps, which correspond to ~6 hr in experimental time, according to the experiments in Figure 1 and the simulation results in Figures 4 and 5). The ATF4 level under PKR-silencing decreased much more significantly than the level under PERK-silencing in the earlier time point (e.g. at 5 time-steps, captured at ~3 hr in experimental time). This subtle difference comes from the underlying mechanism of the regulatory network where PKR and PERK plays different roles (i.e. palmitate initiates the PKR-eIF2α-ATF4 pathway through PACT-PKR activation at an earlier time then PERK), and is captured by the model. Our experimental results confirm this difference (Figure 7B and C), that PKR-knockdown cells showed significantly lower ATF4 gene expression level at the earlier time point (3 hr), while PERK-knockdown started to reduce the ATF4 gene expression level at a later time (6 hr). These results are consistent with the in-silico simulations (Figure 7A), suggesting that PKR is responsible for the early activation of ATF4. In addition, since ER stress in human embryonic kidney cells has been shown to involve the activation of PKR by PACT (PKR activating protein) , we experimentally tested whether PACT silencing blocks ATF4 mRNA expression levels. Treating the cells with palmitate for 3 hr resulted in a significant decrease in the ATF4 gene expression (Figure 7C), suggesting that the PACT-PKR pathway is essential for the early ATF4 gene expression, and possibly for an early ER stress response. Based on our computational analysis and experimental results, we identified the essential pathways involved in activating ATF4 and CREB1 in HepG2 liver cells upon palmitate treatment (Figure 8).
Increased long-chain saturated fatty acids, a characteristic feature of obesity and NASH , can induce ER stress, activate the UPR, and lead to cell death in hepatocytes [12, 13]. Increasing number of studies implicate an involvement of palmitate in ER stress-induced apoptosis in liver cells [11–14, 29, 30], but the mechanism by which palmitate regulates the UPR signaling pathways is unclear. The present study sought to determine which signaling pathways are essential for the activation of ATF4, which plays a crucial role as a mediator of ER stress-induced apoptosis.
PKR, an eIF2α kinase, has recently emerged to be involved in ER stress-induced apoptosis by pharmacological ER inducers, thapsigargine and tunicamycin [20, 31, 32]. Lee et al. found that activated PKR phosphorylated the eIF2α in a PERK-independent manner when HEK293A cells were treated with thapsigargine . They showed that PKR was responsible for approximately 40% and PERK for over 50% of the total phosphorylated eIF2a at 12 h upon thapsigargin treatment. However, our study suggests that the contributions of PKR and PERK on eIF2α phosphorylation change in a time dependent-manner upon palmitate treatment. At an earlier time (3 hr), PKR is entirely responsible for eIF2α phosphorylation (Figure 4D), with PERK adding to this contribution, to the eIF2α phosphorylation, at 6 hr (Figure 1). The relative contributions of PKR and PERK on eIF2α phosphorylation would also depend on the types of inducers and cells [20, 33, 34]. It would be worthy in future to investigate the kinetics of PKR and PERK activations to obtain more information on the roles of PKR and PERK-dependent ER stress.
Along with the positive regulation of PKR and PERK on eIF2α phosphorylation, a negative feedback loop could relieve the translational inhibition during ER stress, through GADD34-mediated PP1 activation [35–37]. ATF4 has been shown to bind specifically to a conserved ATF site on the GADD34 promoter in a stress-dependent manner. GADD34 in turn binds and activates the PP1 protein, which leads to eIF2α dephosphorylation and the resumption of general translation. In contrast to the conventional mechanism, we found that palmitate did not promote PP1 phosphorylation and the negative feedback mechanism was not activated in the palmitate-induced ER stress in HepG2 cells. Blockage of the negative feedback by palmitate may promote ATF4 accumulation in the cells while maintaining a normal level of eIF2α phosphorylation. In keeping with the concept of our study, a recent report suggested that salubrinal (well known as an eIF2α inhibitor) interrupts the feedback mechanism by inhibiting GADD34/PP1 complex activity, thereby potentiating palmitate-induced ER stress in pancreatic β-cells . Similar to the synergistic roles of salubrinal and palmitate in β-cells, palmitate alone may be sufficient to induce ER stress-induced apoptosis in HepG2 cells [17, 39].
ATF4 is known to form homo and heterodimers for DNA binding . One study showed that ER stress induces a complex formation between endogenous ATF1/CREB1 and ATF4, subsequently binding to the ATF/CRE site of the Grp78 promoter . Along with the ChIP-chip data , one group recently suggested that CREB1 is able to bind the CRE site (CGTCA ; −921 to −917 ) on the Atf4 promoter and a constitutively active PKA catalytic subunit dramatically activated Atf4 promoter through PKA-CREB1 signaling . Consistent with this study, our experimental and computational data showed that PKA activated by palmitate at a later time (24 hr) upregulates CREB phosphorylation, subsequently inducing the ATF4 expression. In addition, our EMSA data suggested that ATF4 protein binds its own Atf4 promoter region containing the CRE site (CGTCA ; −164 to −160) and the positive-feedback stabilizes the ATF4 expression. Upon palmitate treatment, PKA or Ca2+-dependent CaM signaling pathways could enhance the interaction of CREB1 alone or the ATF4/CREB1 heterodimer with the Atf4 promoter to prolong ATF4 accumulation. Although these remain open questions, our findings suggested that both ATF4 and CREB1 binding on the Atf4 promoter plays an important role in prolonging ATF4 accumulation upon palmitate exposure.
In this study, our in silico model provides a framework to integrate regulatory information into a complex network and test the network for consistency by comparing the dynamic simulations with experimental measurements. Our simulations show inconsistencies between current knowledge of the network and our observation of palmitate-induced ER stress in the liver cells, which led us to explore the temporal response of the different pathways (e.g. PACT-PKR, PERK, PKA) to palmitate stimulation. This led to the identification of an essential feedback structure in the downstream ATF4/CREB1 regulation that differs from the conventional understanding of the mechanism. The study demonstrates an iterative learning process in which we begin by integrating the current knowledge to build a model to generate hypothesis, which is tested with experiments to obtain novel information. The new information is then incorporated into the network model to correct and update our understanding of the regulation. The updated model can be used to guide new experiments, thereby forming an iterative process that can systematically be applied to study biological processes. The discrete dynamic modeling approach developed in  for signaling network was expanded to account for transcriptional regulation by introducing delays into the regulation. New data and new information could be easily incorporated by adding delays (temporal information), altering gene activity states (perturbations with activator/inhibitor) or the network wiring (binding information), as demonstrated in this study of the regulatory network underlying palmitate-induced ER stress. This systems biology approach can be applied to other biological process to elucidate the molecular mechanism underlying the regulatory network by integrating experimental measurements and computational simulation.
Integrating experiments and computational simulations helped to identify several notable findings, summarized in Figure 8. First, the dynamics of the signaling profiles show that palmitate initiates the PKR-eIF2α-ATF4 pathway through PACT-PKR activation at an earlier time and PERK later helps to maintain eIF2α phosphorylation. Second, ATF4 and CREB1 bind the Atf4 promoter to contribute to a prolonged ATF4 accumulation and their feedbacks appear necessary for ER-stress induced apoptosis. Third, palmitate responds later to PKA activation and possibly along with the Ca2+ dependent-CaM signaling pathway to increase CREB1 phosphorylation. Thus this systems biology approach helped provide insight into the molecular mechanisms by which palmitate induces ER stress in liver cells.
The discrete dynamic modeling of biological network
We applied the approach of discrete dynamic modeling proposed in our previous study . The model is constructed based on the topology and the regulatory relationships within a given network, with basic assumption that the network architecture defines the major dynamic characteristics of the system. We associate each component (protein) in the network with a discrete variable with three potential states (0: lower than control, 1: the control state, 2: higher than control), and simulate the system with transition rules (shift-up or shift-down) depending on the regulatory relationships between the components (activation or inhibition). More specifically, if an activator is in a state 2, (for example, a kinase is activated), the state of its target gene will be shifted upwards, from 0 to 1, or 1 to 2, depending on the target gene’s current state. In contrast, if the state of an inhibitor is higher than control, its target will be shifted downwards in the next updating event. The state of a component will decay (back to control state 1) if its regulators can no longer maintained their active state. Given this transition rules and a specific initial state of the system, we can compute a series of “state-changes” for each component along discrete time steps, where the current state of the system (i.e. the state of each component in the current time-step) depends on the both the previous states (i.e. states in the previous time-steps) as well as the transition process which imposes transition rules based on regulatory relationships within the network. Since the reaction rates may be different from cell to cell even for the same interaction, we apply asynchronous updating of the state, which is realized immediately, rather than renewing every variable simultaneously at each time-step. Thus, the relative rates of the different reactions can be specified by the ordering of the update, which implies that, although the response may be similar, the rate of response varies from cell to cell, or, between different runs. We perform 5,000 independent runs in each simulation to mimic a cell population that is measured by western blotting or RT-PCR experiments. We impose a distribution of initial states (centered at the control state 1) and a randomization of the ordering of the updates to represent cell-cell variation. The dynamic model was implemented by custom MATLAB code.
Two novel properties are introduced to the original modeling approach in  to extend the methodology to address different time-frames and to refine the model with posterior information (novel experimental observations).
Delayed processes. To deal with certain regulatory interactions that are expected to respond much slower than other signaling processes, such as transcription, translation, and the slower response of PERK upon palmitate treatment, we introduce a delay on these interactions such that the current state of the targeted gene depends on not the last state but an earlier previous state of its activator/inhibitor. Therefore, different from the original state transition: S i = f(S i-1 ), in which the current system state S i depends on the system state of the last step S i-1 , our new approach does not apply this Markov assumption, and instead S i = f(S i-1, S i-2, … S 1 ), which could present a history-dependent dynamics. In practice, we implement the idea by simply introducing multiple intermediate (virtual) components in between the regulator and its targets to delay the signal transduction.
Encode posterior information into the model. Novel information, such as time-separation of the interactions, or the measurement of the activity of particular component in the network, could be incorporated to refine the model by adding constraints to the simulation. The constraints could be a delay of a certain interaction, or impose a constant value of activity on certain components (e.g. PP1) during the simulation.
The discrete dynamic model is based on the simple logic (activation/inhibition) in a regulatory network, which corresponds to our current (mostly qualitative) understanding in this biological system. We measure the average gene expression level in an “in silico” cell population that has varying initial states and order of reactions, to enhance the robustness of the simulation results. The Additional file 1: Figure S4 shows the increase in robustness of the model when a larger sample size is applied. There are variances in the different replicates of the simulation when the sample size is small, e.g. 50 runs (Additional file 1: Figure S4). These differences come from the randomness of the initiation (i.e. initial states) and the shuffling of the reactions, which are not related to the regulatory mechanism that is being modeled. By increasing the sampling size, the model becomes more robust. With 5000 runs, the different initial states or order of reaction have no effect on the averaged response curve, thus any differences that is observed between the different perturbations should reflect changes in the regulation imposed by the perturbation.
Cell culture and reagents
HepG2/C3A human hepatocellular carcinoma cells were cultured in Dulbecco’s modified Eagle’s medium (DMEM) supplemented with 10% fetal bovine serum and 1% of penicillin–streptomycin (penicillin: 10,000 U/ml, streptomycin: 10,000 μg/ml; Invitrogen) in a humidified incubator at 37°C and 5% CO2. Sodium palmitate was purchased from Sigma. In all experiments, palmitate (0.7 mM) was complexed to 2% (w/v) fatty acid free BSA (US Biologicals) dissolved in regular medium. For palmitate treatment, HepG2 cells were seeded in the 6-well plate and cultured until 90% confluence.
RNA interference and reverse transfection
SiRNAs targeting human PACT (5′-GAGAGAAUAUACUACAAUUTT-3′ and 5′-AAUUGUAGUAUAUUCUCUTT-3′), human PKR (5′-GGUGAAGGUAGAUCAAAGATT-3′ and 5′-UCUUUGAUCUACCUUCACCTT-3′), human PERK (5′-GUGACGAAAUGGAACAAGATT-3′ and 5′-UCUUGUUCCAUUUCGUCACTA-3′), human CREB1 (5′-CCGUAACUCUAGUACAGCUTT-3′ and 5′-AGCUGUACUAGAGUUACGGTG-3′) were purchased from Ambion. SiRNAs targeting human ATF4 (siGENOME SMARTpool #M-005125-005) was purchased from Dharmacon. The scrambled non-targeting siRNA (Ambion), used as a negative control, or the siRNA targeting proteins was diluted in Opti-MEM (Invitrogen) and then mixed with the transfection reagent, Lipofectamine RNAiMAX (Invitrogen). The mixture of siRNA and Lipofectamine RNAiMAX in Opti-MEM was then added into six well plates and incubated at room temperature for 20 min. 6 × 105 HepG2 cells were suspended in antibiotic-free medium, seeded into 6-well plates, and cultured for 24 hr. After transfection, the cells were incubated in regular medium or palmitate-containing medium and then collected for further analysis.
Real-time quantitative RT-PCR analysis
Total RNA was extracted from cells with the RNeasy mini kit (Qiagen). One microgram of total mRNA was reverse-transcribed using an iScript cDNA synthesis kit (Bio-RAD). The first-strand cDNA was used as a template. The primers used for quantitative RT-PCR analyses of human ATF4 (5′-TGGACTTCGAGCAAGAGATG-3′ and 5′-AGGAAGGAAGGCTGGAAGAG-3′) and human actin (5′-ACATCGCCCTGTGGATGACT-3′, and 5′-TCACTTGTGGCCCAGATAGG-3′) were synthesized by Eurofins MWG Operon. Amplifications of the cDNA templates were detected by SYBR Green Supermix (Bio-Rad) using RT-PCR Detection System (Bio-Rad) and the cycle threshold values were determined by the MyIQ software (Bio-Rad). Each sample was performed in triplicate and normalized to the actin expression levels.
Western blot analysis
The HepG2 cells were washed twice with cold PBS and treated with CelLyticM cell lysis buffer (Sigma-Aldrich) or RIPA buffer (50 mM Tris (pH 8.0), 150 mM NaCl, 1% NP-40, 0.5% sodium deoxycholate, 0.1% SDS) supplemented with protease inhibitor cocktail (Sigma-Aldrich) for 10 min on ice. The cell lysate was clarified by centrifugation at 13000 rpm for 10 min, and the supernatant was collected. Total protein levels were quantified by Bradford assay (Bio-Rad). Thirty micrograms of total protein was loaded onto 9% SDS-PAGE gel, transferred to nitrocellulose membranes, and probed with antibodies for target proteins: PERK (Cell Signaling), p-PERK (BioLegend), eIF2a (Cell signaling), p-eIF2a (Cell signaling) ATF4 (Santa Cruz Biotechnology), CREB (Santa Cruz Biotechnology), p-CREB (Cell signaling), PKR (Novus Biologicals), and p-PKR (Novus Biologicals). The image was analyzed using the Molecular Imager ChemiDoc XRS System from Bio-Rad.
Nuclear extract preparation and EMSA (electrophoretic mobility shift assay)
Nuclear extracts were prepared as described in the literature . 0.7 mM palmitate-treated HepG2 cells (24 hr) were washed with ice-cold PBS, resuspended in Buffer (10 mM HEPES (pH 8.0), 1.5 mM MgCl2, 10 mM KCl, 1 mM DTT, 50 mM NaF, 1 mM orthovanadate, 1 complete protease inhibitor mini tablet/10 mL), and left to swell on ice for 15 min. Cells were then lysed by forcing them through a 25-gauge hypodermic needle 6-8 times. Lysate was centrifuged at 4000 g for 5 min to pellet the crude nuclear fraction. The nuclear pellet was then resuspended in Buffer (20 mM HEPES (pH 8.0), 420 mM NaCl, 1.5 mM MgCl2, 25% glycerol, 0.2 mM EDTA, 1 mM DTT, 50 mM NaF, 1 mM orthovanadate, 1 complete protease inhibitor mini tablet /10 mL) and incubated on ice for 30 min with gentle agitation. The nuclear extract was spun down at 16000 g for 15 min to pellet the nuclei. The supernatant was saved and used for EMSA.
A synthetic oligo representing −175 to −147 (ACTCCTTTTCTCGTCACAGCTACGCCCT) of the ATF4 promoter was used for EMSA. The probe was biotinylated with Biotin 3′-end DNA Labeling Kit from Thermo Scientific. EMSA method was modified from the previous study . 5 μg of nuclear extract was incubated in EMSA buffer (10 mM HEPES (pH 7.9), 50 mM KCl, 1 mM MgCl2, 1 mM EDTA, 0.1% NP-40, 50 mM NaF, 1 mM Orthovanadate) with 2 μg of ATF4 antibody for 1 hr on ice. The biotinylated DNA probe was added to each mixture (with the DNA probe only as a negative control and the DNA probe with nuclear extract as positive control) and incubated for 30 min at room temperature. The reactions were loaded on a 4-20% TGX gel and separated by gel electrophoresis. The gel was then transferred on to a nylon membrane and UV crosslinked with a hand held crosslinker for 10 min at 254 nm. The probe was then visualized using the Chemiluminescencent Nucleic Acid Detection Module from Thermo Scientific.
Activating transcription factor-4
B cell lymphoma 2
CCAAT/enhancer binding protein (C/EBP)/ATF response element
cAMP-responsive element-binding protein 1
Electrophoretic mobility shift assay
Eukaryotic translation inhibition factor 2α
Free fatty acid
Inositol requiring enzyme 1
PKR activating protein
PKR-like ER kinase
cAMP-dependent protein kinase A
Double-stranded RNA-activated protein kinase
p53-upregulated modulator of apoptosis
Unfolded protein response.
Unger RH, Orci L: Lipoapoptosis: its mechanism and its diseases. Biochimica Et Biophysica Acta-Molecular and Cell Biology of Lipids. 2002, 1585 (2–3): 202-212.
Feldstein AE, Canbay A, Angulo P, Taniai M, Burgart LJ, Lindor KD, Gores GJ: Hepatocyte apoptosis and Fas expression are prominent features of human nonalcoholic steatohepatitis. Gastroenterology. 2003, 125 (2): 437-443. 10.1016/S0016-5085(03)00907-7.
Borradaile NM, Han X, Harp JD, Gale SE, Ory DS, Schaffer JE: Disruption of endoplasmic reticulum structure and integrity in lipotoxic cell death. J Lipid Res. 2006, 47 (12): 2726-2737. 10.1194/jlr.M600299-JLR200.
Kaufman RJ: Stress signaling from the lumen of the endoplasmic reticulum: coordination of gene transcriptional and translational controls (vol 13, pg 1211, 1999). Genes Dev. 1999, 13 (14): 1898-1898.
Harding HP, Zhang YH, Ron D: Protein translation and folding are coupled by an endoplasmic-reticulum-resident kinase (vol 397, pg 271, 1999). Nature. 1999, 398 (6722): 90-90. 10.1038/18057.
Karpinski BA, Morle GD, Huggenvik J, Uhler MD, Leiden JM: Molecular cloning of human CREB-2: an ATF/CREB transcription factor that can negatively regulate transcription from the cAMP response element. Proc Natl Acad Sci USA. 1992, 89 (11): 4820-4824. 10.1073/pnas.89.11.4820.
Vallejo M, Ron D, Miller CP, Habener JF: C/ATF, a member of the activating transcription factor family of DNA-binding proteins, dimerizes with CAAT/enhancer-binding proteins and directs their binding to cAMP response elements. Proc Natl Acad Sci USA. 1993, 90 (10): 4679-4683. 10.1073/pnas.90.10.4679.
Harding HP, Zhang YH, Zeng HQ, Novoa I, Lu PD, Calfon M, Sadri N, Yun C, Popko B, Paules R, et al: An integrated stress response regulates amino acid metabolism and resistance to oxidative stress. Mol Cell. 2003, 11 (3): 619-633. 10.1016/S1097-2765(03)00105-9.
McCullough KD, Martindale JL, Klotz LO, Aw TY, Holbrook NJ: Gadd153 sensitizes cells to endoplasmic reticulum stress by down-regulating Bc12 and perturbing the cellular redox state. Mol Cell Biol. 2001, 21 (4): 1249-1259. 10.1128/MCB.21.4.1249-1259.2001.
Kilberg MS, Shan J, Su N: ATF4-dependent transcription mediates signaling of amino acid limitation. Trends Endocrinol Metab. 2009, 20 (9): 436-443. 10.1016/j.tem.2009.05.008.
Cazanave SC, Elmi NA, Akazawa Y, Bronk SF, Mott JL, Gores GJ: CHOP and AP-1 cooperatively mediate PUMA expression during lipoapoptosis. Am J Physiol Gastrointest Liver Physiol. 2010, 299 (1): G236-G243. 10.1152/ajpgi.00091.2010.
Wei YR, Wang D, Topczewski F, Pagliassotti MJ: Saturated fatty acids induce endoplasmic reticulum stress and apoptosis independently of ceramide in liver cells. Am J Physiol Endocrinol Metab. 2006, 291 (2): E275-E281. 10.1152/ajpendo.00644.2005.
Wei Y, Wang D, Pagliassotti MJ: Saturated fatty acid-mediated endoplasmic reticulum stress and apoptosis are augmented by trans-10, cis-12-conjugated linoleic acid in liver cells. Mol Cell Biochem. 2007, 303 (1–2): 105-113.
Wei YR, Wang D, Gentile CL, Pagliassotti MJ: Reduced endoplasmic reticulum luminal calcium links saturated fatty acid-mediated endoplasmic reticulum stress and cell death in liver cells. Mol Cell Biochem. 2009, 331 (1–2): 31-40.
Das SK, Chu WS, Mondal AK, Sharma NK, Kern PA, Rasouli N, Elbein SC: Effect of pioglitazone treatment on endoplasmic reticulum stress response in human adipose and in palmitate-induced stress in human liver and adipose cell lines. Am J Physiol Endocrinol Metab. 2008, 295 (2): E393-E400. 10.1152/ajpendo.90355.2008.
Magne L, Blanc E, Legrand B, Lucas D, Barouki R, Rouach H, Garlatti M: ATF4 and the integrated stress response are induced by ethanol and cytochrome P450 2E1 in human hepatocytes. J Hepatol. 2011, 54 (4): 729-737. 10.1016/j.jhep.2010.07.023.
Zhang L, Seitz LC, Abramczyk AM, Chan C: Synergistic effect of cAMP and palmitate in promoting altered mitochondrial function and cell death in HepG2 cells. Exp Cell Res. 2010, 316 (5): 716-727. 10.1016/j.yexcr.2009.12.008.
Kim YH, Joo HS, Kim DS: Nitric oxide induction of IRE1-alpha-dependent CREB phosphorylation in human glioma cells. Nitric Oxide-Biology and Chemistry. 2010, 23 (2): 112-120. 10.1016/j.niox.2010.04.009.
Wu M, Yang XR, Chan C: A Dynamic Analysis of IRS-PKR Signaling in Liver Cells: A Discrete Modeling Approach. PLoS One. 2009, 4 (12): e8040-10.1371/journal.pone.0008040.
Lee ES, Yoon CH, Kim YS, Bae YS: The double-strand RNA-dependent protein kinase PKR plays a significant role in a sustained ER stress-induced apoptosis. FEBS Lett. 2007, 581 (22): 4325-4332. 10.1016/j.febslet.2007.08.001.
Gonzalez GA, Montminy MR: Cyclic AMP stimulates somatostatin gene transcription by phosphorylation of CREB at serine 133. Cell. 1989, 59 (4): 675-680. 10.1016/0092-8674(89)90013-5.
Wadzinski BE, Wheat WH, Jaspers S, Peruski LF, Lickteig RL, Johnson GL, Klemm DJ: Nuclear protein phosphatase 2A dephosphorylates protein kinase A-phosphorylated CREB and regulates CREB transcriptional stimulation. Mol Cell Biol. 1993, 13 (5): 2822-2834.
Luo J, Zhou W, Zhou X, Li D, Weng J, Yi Z, Cho SG, Li C, Yi T, Wu X, et al: Regulation of bone formation and remodeling by G-protein-coupled receptor 48. Development. 2009, 136 (16): 2747-2756. 10.1242/dev.033571.
Zhang XM, Odom DT, Koo SH, Conkright MD, Canettieri G, Best J, Chen HM, Jenner R, Herbolsheimer E, Jacobsen E, et al: Genome-wide analysis of cAMP-response element binding protein occupancy, phosphorylation, and target gene activation in human tissues. Proc Natl Acad Sci USA. 2005, 102 (12): 4459-4464. 10.1073/pnas.0501076102.
Mahishi L, Usdin K: NF-Y, AP2, Nrf1 and Sp1 regulate the fragile X-related gene 2 (FXR2). Biochem J. 2006, 400 (2): 327-335. 10.1042/BJ20060734.
Sun P, Enslen H, Myung PS, Maurer RA: Differential activation of CREB by Ca2+/calmodulin-dependent protein kinases type II and type IV involves phosphorylation of a site that negatively regulates activity. Genes Dev. 1994, 8 (21): 2527-2539. 10.1101/gad.8.21.2527.
Wu GY, Deisseroth K, Tsien RW: Activity-dependent CREB phosphorylation: convergence of a fast, sensitive calmodulin kinase pathway and a slow, less sensitive mitogen-activated protein kinase pathway. Proc Natl Acad Sci USA. 2001, 98 (5): 2808-2813. 10.1073/pnas.051634198.
Nehra V, Angulo P, Buchman AL, Lindor KD: Nutritional and metabolic considerations in the etiology of nonalcoholic steatohepatitis. Dig Dis Sci. 2001, 46 (11): 2347-2352. 10.1023/A:1012338828418.
Ibrahim SH, Akazawa Y, Cazanave SC, Bronk SF, Elmi NA, Werneburg NW, Billadeau DD, Gores GJ: Glycogen synthase kinase-3 (GSK-3) inhibition attenuates hepatocyte lipoapoptosis. J Hepatol. 2011, 54 (4): 765-772. 10.1016/j.jhep.2010.09.039.
Pfaffenbach KT, Gentile CL, Nivala AM, Wang D, Wei YR, Pagliassotti MJ: Linking endoplasmic reticulum stress to cell death in hepatocytes: roles of C/EBP homologous protein and chemical chaperones in palmitate-mediated cell death. Am J Physiol Endocrinol Metab. 2010, 298 (5): E1027-E1035. 10.1152/ajpendo.00642.2009.
Shimazawa M, Ito Y, Inokuchi Y, Hara H: Involvement of double-stranded RNA-dependent protein kinase in ER stress-induced retinal neuron damage. Invest Ophthalmol Vis Sci. 2007, 48 (8): 3729-3736. 10.1167/iovs.06-1122.
Singh M, Fowlkes V, Handy I, Patel CV, Patel RC: Essential Role of PACT-Mediated PKR Activation in Tunicamycin-induced Apoptosis. J Mol Biol. 2009, 385 (2): 457-468. 10.1016/j.jmb.2008.10.068.
Anderson BR, Karikó K, Weissman D: Nucleofection induces transient eIF2α phosphorylation by GCN2 and PERK. Gene Ther. 2012, 10.1038/gt.2012.5.
Brunelli C, Amici C, Angelini M, Fracassi C, Belardo G, Santoro MG: The non-steroidal anti-inflammatory drug indomethacin activates the eIF2α kinase PKR, causing a translational block in human colorectal cancer cells. Biochem J. 2012, 443 (2): 379-386. 10.1042/BJ20111236.
Ma YJ, Hendershot LM: Delineation of a negative feedback regulatory loop that controls protein translation during endoplasmic reticulum stress. J Biol Chem. 2003, 278 (37): 34864-34873. 10.1074/jbc.M301107200.
Brush MH, Weiser DC, Shenolikar S: Growth arrest and DNA damage-inducible protein GADD34 targets protein phosphatase 1 alpha to the endoplasmic reticulum and promotes dephosphorylation of the alpha subunit of eukaryotic translation initiation factor 2. Mol Cell Biol. 2003, 23 (4): 1292-1303. 10.1128/MCB.23.4.1292-1303.2003.
Novoa I, Zeng HQ, Harding HP, Ron D: Feedback inhibition of the unfolded protein response by GADD34-mediated dephosphorylation of eIF2 alpha. J Cell Biol. 2001, 153 (5): 1011-1021. 10.1083/jcb.153.5.1011.
Cnop M, Ladriere L, Hekerman P, Ortis F, Cardozo AK, Dogusan Z, Flamez D, Boyce M, Yuan JY, Eizirik DL: Selective inhibition of eukaryotic translation initiation factor 2 alpha dephosphorylation potentiates fatty acid-induced endoplasmic reticulum stress and causes pancreatic beta-cell dysfunction and apoptosis. J Biol Chem. 2007, 282 (6): 3989-3997.
Srivastava S, Chan C: Hydrogen peroxide and hydroxyl radicals mediate palmitate-induced cytotoxicity to hepatoma cells: relation to mitochondrial permeability transition. Free Radic Res. 2007, 41 (1): 38-49. 10.1080/10715760600943900.
Ameri K, Harris AL: Activating transcription factor 4. Int J Biochem Cell Biol. 2008, 40 (1): 14-21. 10.1016/j.biocel.2007.01.020.
Luo S, Baumeister P, Yang S, Abcouwer SF, Lee AS: Induction of Grp78/BiP by translational block: activation of the Grp78 promoter by ATF4 through and upstream ATF/CRE site independent of the endoplasmic reticulum stress elements. J Biol Chem. 2003, 278 (39): 37375-37385. 10.1074/jbc.M303619200.
Celis JE: Cell biology: a laboratory handbook. 2006, Amsterdam; Boston: Elsevier Academic, 3
Ma Y, Brewer JW, Diehl JA, Hendershot LM: Two distinct stress signaling pathways converge upon the CHOP promoter during the mammalian unfolded protein response. J Mol Biol. 2002, 318 (5): 1351-1365. 10.1016/S0022-2836(02)00234-6.
We thank the members of the Cellular and Biomolecular Laboratory.
This study was supported in part by the National Science Foundation (CBET 0941055) and the National Institutes of Health (R01GM079688 and 1R01GM089866).
The authors declare that they have no competing interests.
HC and LZ conducted all experiments except for EMSA. MW performed the computational simulations. RT carried out the EMSA experiments. AN participated in the analysis of ChIP-chip data. HC and MW drafted the manuscript. CC guided and directed the overall design and coordination of the study and edited the manuscript. All authors read and approved the final manuscript.
Hyunju Cho, Ming Wu, Linxia Zhang contributed equally to this work.
Electronic supplementary material
Additional file 1:Contains experimental and computational simulation results: Figure S1. The protein expression levels of PKAc (PKA catalytic subunits) upon palmitate treatment, Figure S2. The protein expression levels of PP1 upon palmitate treatment, Figure S3. In silico knock-out of signaling pathways for CREB1 activation, Figure S4. Simulation of the discrete dynamic model with different amount of independent sampling size. (DOCX 763 kb) (DOCX 763 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.