|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
From the a Laboratory for Experimental Medicine and Endocrinology, University Hospital Gasthuisberg, Catholic University of Leuven, Herestraat 49, box 902, B-3000 Leuven, Belgium, d Centre for Biological Sequence Analysis, BioCentrum-DTU, Technical University of Denmark, Building 208, DK-2800 Lyngby, Denmark, e Laboratory for Experimental Medicine, Université Libre de Bruxelles, route de Lennik 808, box CP618, B-1070 Brussels, Belgium, g Laboratory of Neuroplasticity and Neuroproteomics and h Functional Genomics and Proteomics Research Unit, Catholic University of Leuven, Naamsestraat 59, B-3000 Leuven, Belgium, i ProMeta and k Laboratory of Biochemistry, University Hospital Gasthuisberg, Catholic University of Leuven, Herestraat 49, box 901, B-3000 Leuven, Belgium, and j ESAT-SDC, Department of Electrical Engineering, Catholic University of Leuven, Kasteelpark Arenberg 10, box 2446, B-3001 Heverlee, Belgium
| ABSTRACT |
|---|
|
|
|---|
0.01) by a combination of interleukin-1β and interferon-
, whereas only 42 and 23 proteins were altered by either of the cytokines alone, giving rise to 199 distinct differentially expressed spots. Identification of 141 of these by MALDI-TOF/TOF revealed proteins playing a role in insulin secretion, cytoskeleton organization, and protein and RNA metabolism as well as proteins associated with endoplasmic reticulum and oxidative stress/defense. We investigated the interactions of these proteins and discovered a significant interaction network (p < 1.27e–05) containing 42 of the identified proteins. This network analysis suggests that proteins of different pathways act coordinately in a β-cell dysfunction/apoptotic β-cell death interactome. In addition the data suggest a central role for chaperones and proteins playing a role in RNA metabolism. As many of these identified proteins are regulated at the protein level or undergo post-translational modifications, a proteomics approach, as performed in this study, is required to provide adequate insight into the mechanisms leading to β-cell dysfunction and apoptosis. The present findings may open new avenues for the understanding and prevention of β-cell loss in type 1 diabetes.
(IFN-
), and tumor necrosis factor
, secreted by these immune cells are one of the main causes of β-cell death (1, 2). β-Cells are highly specialized cells primarily oriented toward tightly regulated insulin synthesis and secretion. Because of the high demand of proinsulin/insulin processing, they possess a highly developed endoplasmic reticulum (ER) compartment. The correct folding of proinsulin in this ER compartment is mediated by a specific subset of proteins, the molecular chaperones. Altered expression of several of these chaperones in response to different stress situations has been associated with the pathogenesis of type 1 and even type 2 diabetes (3).
Not surprisingly, one of the major changes in β-cells upon exposure to inflammatory cytokines is an increased proinsulin/insulin ratio and loss of insulin secretion in response to glucose (4–6). Upon prolonged exposure of the β-cells to cytokines this functional impairment leads to apoptotic β-cell death as shown convincingly by in vitro exposure of human or rodent β-cells to cytokines (2, 7, 8). Many studies have been undertaken in an attempt to unravel the pathways involved in cytokine-induced β-cell destruction, revealing that stress response genes, either protective or deleterious for β-cell survival, are induced, whereas several genes related to differentiated β-cell functions are down-regulated (7, 9). Central factors in the regulation of these cytokine-induced genes are nuclear factor
B and signal transducer and activator of transcription 1 (STAT1) (7, 10, 11). The strength of the microarray technique used for the majority of these studies is that it allows analysis of a very large number of genes; however, the technique yields little information on the functional state of the cells as changes in transcript levels are not always reflected at the protein level (12). In addition, many forms of regulation take place at a post-transcriptional level, such as post-translational modifications (PTMs), proper protein folding, and selective proteolytic degradation (13, 14). It is thus crucial to complement cDNA array analysis with a global evaluation of protein expression. Previous proteomics studies have been performed on mechanisms of β-cell death mainly on whole rat islets. Most of these studies made use of classical two-dimensional gel electrophoresis (2DE) and focused on the effects of IL-1β alone on rat islets (15–17). Although these contributed greatly to our current understanding on the mechanisms involved in β-cell death, we believe that the use of other inflammatory stimuli, such as a combination of IL-1β and IFN-
, will add valuable data improving our knowledge on the mechanisms involved in cytokine-induced β-cell death. Indeed although exposure of purified human or rodent β-cells to IL-1β alone will induce β-cell dysfunction, it is not sufficient to induce apoptosis (2). Combining IL-1β with IFN-
will induce an extensive level of apoptosis, as reported in human and rodent islets as well as in INS-1E cells, even at relatively low cytokine concentrations. This suggests that IFN-
signal transduction synergizes with IL-1β signaling to trigger β-cell apoptosis (3).
With the recent advances in differential dye labeling and the introduction of the 2D DIGE technique (18) as well as the development of more sophisticated MS equipment for protein identification, an attractive technology for the quantitative analysis of differences in protein profiles has become available (19, 20). In the present project we investigated by 2D DIGE analysis the complex processes involved in cytokine-induced β-cell death, making use of the rat insulin-producing INS-1E cell line. These are well differentiated cells and have been found to be a stable and reliable β-cell model (21) with a global pattern of cytokine-induced gene expression similar to that observed in primary β-cells (9, 22). The effects of single cytokine treatments as well as a combined treatment with IL-1β and IFN-
were analyzed, each covering a broad pH range by using an overlap of two different IPG strips (pH 4–7 and 6–9). This allowed us to conclude that, at the concentrations presently used, IL-1β and IFN-
act synergistically in the induction of apoptosis, concurrent with major alterations in protein levels, as compared with single cytokine treatment alone. Interestingly the minor changes observed in single cytokine treatments point to a protection or adaptation of the INS-1E cells toward the treatment. For the combined treatment, where we were able to identify 119 differentially expressed proteins, the changes are consistent with an induction of the apoptotic cell type. Importantly many of the regulations observed occur at a post-translational level, for instance for the ER chaperones immunoglobulin heavy chain-binding protein/glucose-regulated protein, 78 kDa (BiP/GRP78) and oxygen-regulated protein, 150 kDa (ORP150) and for the heterogeneous nuclear ribonucleoproteins (hnRNPs). A protein interaction analysis of the differentially expressed proteins revealed a significant network (p < 1.27e–05) containing 42 of these proteins. These results suggest that many of the pathways involved act in a coordinate way with a central role for molecular chaperones and proteins involved in RNA metabolism.
| EXPERIMENTAL PROCEDURES |
|---|
|
|
|---|
Cell Death Analysis and Nitrite Determination—
The percentage of living, apoptotic, and necrotic cells was assessed by cell counting. INS-1E cells were cultured in 96-well plates (7500 cells/well). Typically after 24 h of control or cytokine-treatment (10 units/ml IL-1β, 500 units/ml IFN-
, or the combination of the two cytokines with the same concentrations) cells were incubated for 15 min with propidium iodide (10 µg/ml) and Hoechst HO 342 (20 µg/ml) (23). A minimum of 500 cells was counted in each experimental condition by two researchers counting independently of each other on a Zeiss Axiovert 40 CFL fluorescence microscope (Zeiss, Jena, Germany).
For nitrite determination, 50 µl of culture supernatant was collected after 24 h from the cell cultures for the proteomics experiments (2.5 x 106 cells/100-mm dish). Nitrite determination was performed using Griess reagent as described previously (24), and the concentration was determined spectrophotometrically at 546 nm.
Sample Preparation for 2D DIGE Analysis—
INS-1E cells were incubated for 1, 4, and 24 h with 10 units/ml recombinant human IL-1β (a kind gift from Dr. C. W. Reinolds, National Cancer Institute, National Institutes of Health, Bethesda, MD) and 500 units/ml recombinant rat IFN-
(R&D Systems, Minneapolis, MN). For the 24-h condition, cells were also incubated with individual cytokines using the same concentrations. For each condition, 12.5 x 106 INS-1E cells (2.5 x 106 cells/100-mm dish) were collected. For each time point or treatment condition quadruplicate experiments were performed, originating from four independent experiments. The cells were washed twice in phosphate-buffered saline and resuspended in 100 µl of lysis buffer (7 M urea, 2 M thiourea, 4% (w/v) CHAPS, 40 mM Tris base, 1% (w/v) DTT, and a mixture of protease inhibitors (Complete protease inhibitor, Roche Diagnostics)). Samples were sonicated (UP50H, Hielscher, Teltow, Germany) five times for 10 s after which they were centrifuged for 13 min at 20,000 x g (4 °C) to precipitate all non-soluble proteins and membrane fractions. The supernatant was desalted by dialysis (PlusOne Mini Dialysis kit, GE Healthcare), and the pH was adjusted to 8.5. Protein concentration was determined using the Bradford assay (25). 50 µg of protein sample was labeled with 200 pmol of Cy3 or Cy5 (two "forward" and two "reverse" gels were made at each time point), whereas 50 µg of pooled internal standard was labeled with 200 pmol of Cy2 (18). The internal standard consisted of a pool of all control and treated samples of all four replicate experiments. The labeling reaction was carried out for 30 min on ice and quenched with 10 mM lysine (15 min on ice). Labeled protein extracts were pooled, and sample loading buffer was added (7 M urea, 2 M thiourea, 4% (w/v) CHAPS, 0.5% (v/v) IPG buffer, and depending on the pH range of the strips, 1% (w/v) DTT for IPG strips pH 4–7 or 1.2% (v/v) Destreak (GE Healthcare) for pH 6–9).
2D DIGE and Gel Imaging—
IPG strips (24 cm) were rehydrated overnight in 450 µl of rehydration buffer (7 M urea, 2 M thiourea, 4% (w/v) CHAPS, 0.5% (v/v) IPG buffer, 0.05% (w/v) OrangeG, and either 1% (w/v) DTT (for IPG strips pH 4–7) or 1.2% (w/v) Destreak (for IPG strips pH 6–9)). The pooled samples containing sample loading buffer were loaded onto the rehydrated strips using anodic cup loading and separated according to their isoelectric point on an Ettan IPGphor II manifold. The complete process was tracked using the Ettan IPGphor control software (version 1.01.03) (GE Healthcare). The first dimension was ended when the current reached a stable phase (at
60 kV-h).
Prior to the second dimension, the strips were equilibrated during two intervals of 15 min each in an equilibration buffer (6 M urea, 30% (v/v) glycerol, 2% (w/v) SDS, and 50 mM Tris-HCl, pH 8.8) containing 1% (w/v) DTT in the first step and 4% (w/v) iodoacetamide and 0.02% bromphenol blue in the second step. Equilibrated strips were placed on top of 12.5% SDS-polyacrylamide gel and separated on an Ettan DaltSix system (GE Healthcare).
Scanning of the gels was performed using a Typhoon 9400 (GE Healthcare) at 100-µm pixel size. Prior to analysis with the DeCyderTM Version 6.5 software, gel images were cropped using ImageQuant TL. Spot detection and matching was performed automatically using the "Batch Processor" module of the DeCyder Version 6.5 software followed by careful manual rematching of wrongly matched spots or unmatched spots. The 12 spot maps corresponding to the four gels from each time point and each pH range were used to calculate average abundance.
Spot Digestion and Protein Identification by MALDI-TOF/TOF Analysis—
For spot picking, two preparative gels for each pH range were run (350 µg of protein lysate each). First and second dimension runs were performed as described above except that CyDye labeling was omitted. Glass plates were pretreated with BindSilane, and two reference markers were applied to enable automatic picking. The gels were poststained using Deep Purple (GE Healthcare). Matching with the analytical gels was performed using the biological variation analysis module of the DeCyder Version 6.5 software. A pick list was generated and exported into the Spot Picker Version 1.20 software that controls the Ettan Spot Picker (GE Healthcare).
Spots were picked in MilliQ water, transferred to 100 µl of fixation solution (50% (v/v) methanol, 5% (v/v) acetic acid, and 45% (v/v) MilliQ) and rinsed three times with MilliQ water and three times with ACN (LC-MS quality, CHROMASOLV®, Sigma-Aldrich). The gels were hydrated in a 100 mM NH4HCO3 solution for 10 min followed by a dehydration step in 100% ACN for 10 min with vigorous vortexing. This step was repeated twice prior to dehydrating the gel pieces in a SpeedVac.
Gel pieces were rehydrated in digestion buffer (50 mM NH4HCO3 and 5 mM CaCl2) containing 5 ng/µl modified trypsin (Promega, Madison, WI) and incubated overnight at 37 °C. The resulting peptides were extracted out of the gel plugs in four steps: once with 50 mM NH4HCO3, twice with 50% (v/v) ACN and 5% (v/v) formic acid, and once with 95% ACN and 5% formic acid (30 min each). Supernatants were dried in a SpeedVac.
Upon concentrating and desalting the tryptic fragments using Millipore C18 ZipTips, the samples were mixed in a 1:1 (v/v) ratio with
-cyano-4-hydroxycinnamic acid matrix (saturated solution in 50% ACN and 2.5% TFA in HPLC water), spotted onto the MALDI target plate, and allowed to air dry. MS/MS analyses were performed on a 4800 MALDI-TOF/TOF instrument (Applied Biosystems, Foster City, CA). The instrument was calibrated with Applied Biosystems Calibration Mixture 1. Measurements were taken in the positive ion mode between 900 and 3000 m/z. Sequences were automatically acquired by scanning first in MS mode and selecting the 15 most intense ions for MS/MS using an exclusion list of peaks arising from tryptic autodigestion. Air was used as the collision gas, whereas the collision energy was adapted automatically. Data interpretation was carried out using the GPS Explorer software (Version 3.5), and database searching was carried out using the Mascot program (Version 2.0.00). Because all experiments were performed on rat INS-1E cells, MS/MS searches were conducted with the following settings: NCBI (taxonomy set on Rattus, 33,734 protein entries) and MSDB (Mass Spectrometry Protein Sequence Database) (taxonomy set on Rodentia, 88,439 protein entries) as database, MS/MS tolerance for precursor and fragment ions between 0.2 and 1 Da depending on the sample, methionine oxidation as variable modification, and carbamidomethylation of cysteine as fixed modification. As enzyme, trypsin was selected, and a maximum of one missed cleavage was allowed. Using these parameters the probability-based MOWSE (Molecular Weight Search) scores greater than the given cutoff value for MS/MS fragmentation data were taken as significant (p < 0.05).
For protein identifications where no hit was found in the rat databases, protein identity was based on comparison with the orthologous mouse sequence because the mouse genome is better annotated compared with the rat genome. In these cases the number of identical amino acids is indicated for each peptide (Supplemental Table 1). For all identified proteins, each sequenced peptide was individually aligned using BLAST (26). For protein identifications where all the individual peptides completely matched to more than one Swiss-Prot database accession number/protein name, both protein sequences were aligned using BLAST (26). If this alignment resulted in 100% sequence identity, thus pointing to a single protein present in the database under different names/accession numbers, the Swiss-Prot accession number of the entry with the best description is consistently given in the first place to eliminate redundancy. For completeness also the second accession number is reported (Supplemental Table 1). For peptides matching to different isoforms or to multiple members of a protein family, the following criteria were used for selecting which one to report. 1) If at least one of the identified peptides matched exclusively to a specific isoform or protein member, this protein isoform/member could be identified unambiguously. 2) The experimental molecular weight and pI obtained from the 2DE gel was compared with the theoretical molecular weight and pI of the different isoforms/protein members. In case no single protein could be ruled out based on these criteria or possibly more than one member of the protein group could be present in a single spot, the different names/accession numbers are reported, including the percent identity between them (Supplemental Table 1).
Constructing a Human Protein Interaction Network—
We constructed a human protein interaction network as reported recently (27, 28). Data were downloaded from MINT 48, BIND 49, IntAct 50, KEGG (Kyoto Encyclopedia of Genes and Genomes), and Reactome. To increase the coverage of interactions, interolog data 53 (the transfer of protein interactions between orthologous protein pairs in different organisms) were included. Interactions were transferred from 17 eukaryotic organisms and added to the network. Orthology was assigned using the Inparanoid database 55 using stringent thresholds. The statistical significance of the networks was estimated using a randomization scheme in which networks were generated from random input sets of the same size as in the present experiment. We performed 1000 randomizations from which we derived a probability distribution that was used to calculate the significance of the networks in this study (28).
Statistical Analysis—
For all experiments the paired Student's t test was used. To analyze differences in protein levels, the DeCyder Version 6.5 software was used, and a p value of less than or equal to 0.01 was considered statistically significant.
| RESULTS |
|---|
|
|
|---|
(500 units/ml), or the combination of both cytokines. In parallel to the cell culture experiments performed for proteomics analysis, we investigated the effects of the cytokines on cell viability and nitric oxide (NO) production. As reported previously (9), only a minor increase in the proportion of apoptotic cells was observed with IL-1β (from 1.88 ± 0.83 to 5.40 ± 2.09%, p < 0.005, n = 4) or IFN-
alone (from 1.88 ± 0.83 to 7.10 ± 2.09%, p < 0.05, n = 4), whereas a synergistic effect on apoptosis was seen in the combination-treated cells (from 1.88 ± 0.83 to 21.42 ± 4.93%, p < 0.0005, n = 4). No significant alterations in the proportion of necrotic cells were observed. In parallel with this increase in apoptosis, NO production increased from 29.31 ± 20.7 pmol of nitrite/104 cells in control cells to 624.38 pmol of nitrite/104 cells in IL-1β-exposed cells or to 809.85 ± 170.7 pmol of nitrite/104 cells in IL-1β + IFN-
-exposed cells (p < 0.001, n = 4). For IFN-
alone, no significant rise in NO production was detected (35.27 ± 12.9 pmol of nitrite/104 cells). These findings are in accordance with previously performed studies investigating alterations in mRNA levels in cytokine-induced INS-1E cells under similar conditions (9).
Differential Proteomic Profile of Cytokine-exposed INS-1E Cells—
Proteomic profiles of cytokine-exposed INS-1E cells were determined at two early time points (1 and 4 h) and at one late time point after cytokine induction (24 h) in an attempt to distinguish between early and late effects induced by IL-1β and IFN-
. In addition, after 24 h of incubation, also the effects of IL-1β and IFN-
alone were investigated. For each of these conditions quadruplicate experiments were performed at two different pH ranges, comparing cytokine-exposed with non-exposed control INS-1E cells. In the pH 4–7 range, 2328 ± 151 spots could be matched for IL-1β and IFN-
alone, and 2214 ± 164 spots could be matched for the combination of IL-1β and IFN-
. Similarly for the pH 6–9 range a total of 1821 ± 184 spots were matched for IL-1β alone and IFN-
alone, and 1641 ± 73 spots were matched for the combination. In the pH 4–7 range 31, 23, and 116 spots were observed as differentially expressed after 24-h incubation with IL-1β alone, IFN-
alone, or IL-1β + IFN-
, respectively (p
0.01) (Fig. 1A and Supplemental Fig. 1). In addition, only 1 and 19 spots were differentially expressed after IL-1β + IFN-
treatment at the early time points of 1 and 4 h of incubation, respectively. Because some of these spots were differential in different treatment conditions or at different time points, a total of 159 spots was observed to be differential. For the pH range 6–9, 11, 0, and 30 spots were observed as differentially expressed after treatment with IL-1β alone, IFN-
alone, and IL-1β + IFN-
, respectively (p
0.01, 24 h). For the early time points no spots were differentially expressed after 1 h of incubation, and only two spots were differentially expressed after 4 h of incubation with the combination treatment (p
0.01). Taken together, a total of 40 different differential spots were observed when analyzing this pH range (Fig. 1B and Supplemental Fig. 2).
|
. Of a total of 199 differentially expressed spots, we were able to identify 141 spots (70.9%) by MALDI-TOF/TOF, corresponding to 92 different proteins. Of note, 12 of these spots were found to contain two proteins. Therefore, it is at present not possible to conclude which protein in the spot actually changed expression upon cytokine treatment. Further analysis using narrow pH range zoom strips could tackle this problem. The identified proteins were grouped in different functional classes: 1) metabolism; 2) protein synthesis, modification, and secretion; 3) ionic channels, ion transporters, and related proteins; 4) processing of peptides and neurotransmitters; 5) cell adhesion, cytoskeleton, and related genes; 6) RNA synthesis/turnover; 7) cell cycle and related genes; 8) defense/repair; 9) apoptosis/endoplasmic reticulum stress-related; and 10) miscellaneous (Tables I–IV). The same functional classes were used in our previous microarray analysis (9). The complete list of identified proteins is shown in Tables I–IV and includes information on the -fold regulation, Swiss-Prot accession number, number of identified peptide sequences, number of unique peptides, and percent sequence coverage. Additional information is shown in Supplemental Table 1. Although the majority of protein identifications are based on multiple peptide sequences, three identifications are single peptide-based protein identifications. For these latter proteins the peptide sequence, precursor m/z and charge ratio, and the score/Expect value are shown in Table V, and the annotated MS/MS spectra are shown in Supplemental Fig. 3. The identified proteins play a role in many different pathways, such as insulin biosynthesis (summarized in Fig. 2), energy metabolism (summarized in Fig. 3), cytoskeleton reorganization (summarized in Fig. 4) (29–33), RNA metabolism (such as hnRNP H, hnRNP K, DEAD box protein 1 (DDX1), DEAD box protein 3 X-chromosomal (DD3X3)), oxidative stress (such as manganese-superoxide dismutase (SOD2), peroxiredoxin-1 (PRDX1), and peroxiredoxin-3 (PRDX3)), and molecular chaperones (BiP/GRP78, ORP150, heat shock 70-kDa protein 8 (HSPA8), and mortalin (HSPA9)).
|
|
|
|
|
|
|
|
|
0.01), whereas an isoform of PRDX1 (spot 1975, pH 6–9) was down-regulated (2.7-fold, p
0.01). Another group of proteins prone to PTM are the hnRNPs. Two different members have been identified (hnRNP H1 and hnRNP K), represented in a total of eight different spots; all were down-regulated (p
0.01). Especially hnRNP K seems to be prone to PTMs because six different spots contained this protein. These can be divided into two distinct groups of phosphorylated isoforms, with molecular mass of
70 and 51 kDa, respectively. Also of note, for PCSK1 and PCSK2, two and six isoforms were identified, respectively, all with the same molecular weight and all significantly down-regulated (p
0.01).
Upon exposure to individual cytokines, we were able to identify 38 (30 and 17 identifications for IL-1β and IFN-
, respectively) proteins of a total of 57 (40 and 23 spots for IL-1β and IFN-
, respectively) differentially expressed proteins (75% of the spots identified for IL-1β and 74% of the spots identified for IFN-
). These are listed in Tables III and IV, classified using the same functional groups as above. An important subset of proteins identified in IL-1β-exposed cells were proteins involved in metabolism, whereas the main group of proteins altered by IFN-
plays a role in defense and repair pathways. Of interest, more than half of the proteins identified in the single cytokine condition were not altered significantly in combination-exposed cells. In addition, more than 80 additional proteins were identified in combination-exposed cells compared with single cytokine-exposed cells, confirming the synergistic nature of the combined cytokine exposure.
Network Analysis—
Because the majority of changes in protein profile were observed by treatment with the combination of IL-1β and IFN-
, and this after 24 h of stimulation, concurrent with a major increase in apoptosis, this condition was chosen for further network analysis. Considering that most proteins do not act as single entities but work together in networks, we analyzed the interactions between the proteins identified in this work. This allowed us to examine whether the differentially expressed proteins form one or more functional pathways and to place them in a framework of cellular systems. For this purpose the identified proteins were placed in a protein interaction network recently developed by us (27, 28). This analysis revealed that 42 of the differentially expressed proteins identified constitute a significant interaction network (p < 1.27e–05). This suggests extensive cross-talk between the different proteins and the pathways in which they are involved (Fig. 6 and Supplemental Material 3) with some proteins such as the chaperones BiP/GRP78, HSPA8, and HSPA9 and the RNA synthesis/turnover proteins placed at the center of different networks.
|
| DISCUSSION |
|---|
|
|
|---|
, or IFN-
+ IL-1β treatment of INS-1E cells, a rat insulin-producing cell line used as a model for β-cell dysfunction and death in type 1 diabetes. By applying the 2D DIGE technique to investigate the differential protein expression profiles we identified relevant proteins with a high confidence rate not possible by classical 2DE. The importance of the differentially expressed proteins was further supported by the close linkage they appear to have in the predicted interaction network. Indeed this suggests that the proteins identified, although belonging to a broad range of different classes and functional pathways, interact with each other in close relationships, creating a functional network that may regulate β-cell dysfunction and eventually β-cell death.
The concentrations of cytokines used in the present study clearly act synergistically in inducing apoptotic cell death in INS-1E cells. Indeed only minor alterations in the proportion of apoptotic cells were observed upon exposure to either IL-1β or IFN-
alone, whereas the combined cytokines resulted in a dramatic increase in apoptosis. This is reflected in the observed changes in protein profiles where again a synergistic effect between IL-1β and IFN-
was observed with major changes in protein levels only occurring in the combined treatment and this after a long incubation of 24 h.
Because the primary goal of β-cells is to process and secrete high amounts of insulin, it is not surprising that different proteins playing a crucial role in insulin biosynthesis were altered in expression following cytokine treatment (Fig. 2). In accordance with the reported increase in the proinsulin/insulin ratio and a decreased insulin secretion upon cytokine treatment (4, 6, 34), a clear down-regulation of the different isoforms of PCSK1 and PCSK2, the enzymes responsible for conversion of proinsulin to insulin, was observed. Another interesting protein important in the processing of insulin is thioredoxin-like protein 1 (TXNL1). TXNL1 reduces the disulfide bridges of insulin, thereby initiating the degradation of insulin (35). The observed down-regulation of this protein is yet another detrimental effect of IL-1β and IFN-
on the insulin processing capacities in INS-1E cells.
In recent years, it has become clear that ER stress is implicated in β-cell death (3). One of the well known effects of cytokines on INS-1E cells is the formation of NO (2). Although inducible nitric-oxide synthase itself could not be detected in the 2D gels because of the high molecular weight of this protein, we observed two abundantly expressed ER chaperones, BiP/GRP78 and ORP150, that were clearly shifted toward more acidic forms suggestive for phosphorylation or ribosylation of these proteins. Although there are many studies implicating an important role of ER chaperones in β-cell death, very little information is available on the functional impact of their phosphorylation or ribosylation status. Overexpression of BiP/GRP78 has a protective effect on β-cell death (36), and the present observations on decreased BiP expression is in agreement with our previous findings on an IFN-
-induced decrease in BiP mRNA expression (37, 38). Furthermore phosphorylation of BiP results in a lower binding capacity for proteins (39, 40) that, together with the above described decrease in expression of non-phosphorylated BiP, may aggravate ER stress by interfering with the clearance of unfolded or misfolded proteins (41).
ORP150, another ER chaperone, also shifted toward expression of more acidic isoforms although less markedly than for BiP/GRP78. Increased ORP150 expression improves insulin resistance in animal models of type 2 diabetes and may have a protective role against β-cell death (42). In this context, it will be of interest to study the functional impact of the observed ORP150 PTM.
Another chaperone modified by cytokine treatment is HSPA8, which mainly resides in the cytosol. HSPA8 is an ATP-dependent chaperone (43) and has antiapoptotic effects (44). The observed down-regulation of this chaperone in all three different isoforms identified may hamper β-cell defenses. Opposite results have been described in neonatal rat islets where an increase in HSPA8 protein expression was observed upon treatment with IL-1β alone (45). Because in our experiments INS-1E cells were exposed to IL-1β + IFN-
, this points toward an aggravating effect of IFN-
on β-cell death via inhibition of molecular chaperones as shown previously in both INS-1E and primary β-cells (38). The fourth cytokine-modified member of the hsp70 family is HSPA9, a mitochondrial chaperone involved in protein folding and intracellular trafficking (46). HSPA9 has been implicated in the regulation of apoptosis through transcriptional inactivation of the tumor suppressor protein p53 (47).
Apart from NO formation, another detrimental effect of cytokine exposure to β-cells is the formation of reactive oxygen species (ROS). The observed induction of SOD2 points toward a protective response of the INS-1E cells against cytokine-induced ROS formation in agreement with published findings from microarray (9) or Western blot (10) analyses. However, a longer exposure of β-cells to cytokines will finally result in down-regulation of SOD2, further aggravating β-cell death (10). SOD2-mediated conversion of superoxide molecules will result in a high local production of H2O2. This hydrogen peroxide will be eliminated by reduction to H2O by another group of proteins, the peroxiredoxins (48), and we have presently observed an up-regulation of the cytoplasmic PRDX1. It is puzzling why the other isoform of PRDX1 was down-regulated, which was also the case for the mitochondrial PRDX3, which exerts a function similar to that of PRDX1 (49). PRDX1 and PRDX3 use thioredoxins as electron donors to reduce hydrogen peroxide, and two proteins containing thioredoxin domains were presently identified: the ER protein-disulfide isomerase A6/P5 (PDIA6) and the cytoplasmic TXNL1. Cytokines down-regulated both these proteins, which may impair the β-cell ROS scavenging machinery. As a whole, these observations suggest a complex regulation of the SOD2-peroxiredoxin-thioredoxin system by cytokines that may initially contribute to free radical scavenging but, in the case of protracted exposure to cytokines, may contribute to cell death (10).
Different members of two large groups of proteins playing crucial roles in RNA metabolism were identified in the present study, namely the hnRNPs and the DEAD box family proteins. The main functions of hnRNPs, a class of abundantly expressed proteins, are in transcription, processing of pre-mRNA, alternative splicing, and nucleocytoplasmic shuttling of proteins. hnRNPs themselves are prone to alternative splicing and PTM such as cleavage by caspases (50), altering their activation state. In turn, by having important functions in mRNA splicing (e.g. splicing of Bcl family members) and nucleocytoplasmic shuttling they play an important role in regulation of apoptosis (50). Especially hnRNP K (six spots, pH 4–7) seems to be prone to PTMs. The six different spots that have been identified can be divided into two distinct groups of modified isoforms. In addition, the location of the spots suggests a post-translational cleavage. Although this needs further investigation, cleavage by caspases, accompanied by subcellular translocation, may be involved in this process (50, 51). Two different members of the DEAD box protein family were also identified as altered by cytokine treatment, namely DDX1 and DDX3X. DEAD box proteins are implicated in diverse cellular processes such as RNA splicing, ribosome assembly, and translation initiation. Although the function of DDX1 and DDX3X is poorly understood, the increased expression of DDX3X correlates with its reported proapoptotic role in NIH-3T3 fibroblasts (52).
Compared with previous microarray data (7, 53, 54), few proteins, e.g. SOD2, PCSK1, PCSK2, and BiP/GRP78, were found to correlate with the present proteomics data. Microarray data of cytokine-treated rat or human β-cells (7, 9, 53, 55, 56) revealed many genes that were classified as playing a role in mononuclear cell infiltration/activation (e.g. chemokines), signal transduction (e.g. STAT1, IRF1, and c-jun), and dysfunction (e.g. insulin, pancreas/duodenum homeobox 1 (pdx-1) and insulin gene enhancer protein (isl-1)). Several of these genes were confirmed at the protein level by ELISA, Western blotting, or other methods (5, 57, 58). The majority of these proteins, however, were not detected in the present study. A weak point of 2DE in general is the limitation of protein separation for extreme molecular weight and pI, low abundance, or membrane proteins. This is probably the reason why we could not detect for instance inducible nitric-oxide synthase, chemokines, or transcription factors as modified by exposure to cytokines. It may be possible to overcome these problems, to a certain extent, by making use of narrower pH ranges or prefractionation techniques, but this demanding approach was beyond the scope of the present study. On the other hand we identified many new proteins, not previously identified by microarray analysis, that play a role in for instance metabolism, insulin biosynthesis, cytoskeleton organization, ER stress and defense/repair, and RNA metabolism. Many of these proteins are abundantly expressed and are undergoing PTMs, which could of course not be detected by microarray analysis. It is therefore beyond any doubt that the present data are to a great extent complementary to the previously performed microarray results and add new interesting insights into the mechanisms involved in cytokine-induced β-cell dysfunction and apoptotic cell death.
In addition, treatment with a combination of IL-1β and IFN-
clearly has much more pronounced effects on the global protein profile, as discussed above, as compared with the single cytokine treatments. Of note, although only 17 proteins were identified as changed by IFN-
alone, some of these changes are very pronounced and not observed in the combination-treated cells, such as interferon-inducible GTPase (9.32-fold induced), interferon-
-induced GTPase (9.82-fold induced), or guanylate nucleotide-binding protein 4 (3.70-fold induced). Interestingly these proteins are involved in cellular defense mechanisms. This suggests that treatment with IFN-
alone at the concentration used may activate different protective mechanisms resulting in an adaptation of the cells against the cytotoxic effect of IFN-
consistent with the finding that no induction of apoptosis when using the single cytokine is observed. Next to this, also an interesting picture is observed for IL-1β alone. Here most proteins that are affected are involved in β-cell function and not β-cell death. This is in concordance with previous observations that IL-1β plays a crucial role in inhibiting β-cell function (4, 5), supported by the intriguing recent observations that blocking IL-1β in type 2 diabetic patients partially restores β-cell function (59). The role of IL-1β in in vivo β-cell survival and death, however, is still controversial. Conflicting reports are present in the literature on whether IL-1β actually induces apoptosis in levels that are of relevant impact to affect β-cell mass in vivo or whether, under some conditions, it actually promotes β-cell survival (60, 61). Issues in the different reports are timing of exposure, dose of IL-1β used, and the presence of other cytokines that may amplify or inhibit IL-1β signaling (22, 60, 62).
Transferred to a more clinical situation where β-cells are under inflammatory attack by a combination of different cytokines the picture will be even more complex. Considering that we are dealing with a dynamic network of protein-to-protein interactions, it is most probable that perturbations of this network, induced by the local inflammatory environment, will result in β-cell dysfunction and eventually apoptotic cell death, leading to a diabetic phenotype. Molecular chaperones are suggested to play a crucial role in maintaining this dynamic network, and the present data suggest that the balance between the post-translationally modified and unmodified forms of the chaperones may be of relevance for the cellular outcome. Other central players are the proteins involved in RNA metabolism, such as the hnRNPs and DEAD box proteins. Use of the novel tools for network analysis (28) provides evidence for protein-to-protein interactions between these central players on the one hand and many of the proteins involved in energy metabolism and cytoskeleton organization on the other hand (Fig. 6) that will have a major impact on β-cell function. These findings lead us to hypothesize that, although upon an inflammatory attack the β-cell will eventually die from apoptosis, many important alterations are taking place that will affect the homeostasis of the cell state, resulting in a global β-cell dysfunction. This is also evidenced by the major detrimental effects on insulin processing and secretion, which are after all the primary tasks of this highly specialized cell. In conclusion, the present findings provide new insights into the protein pathways triggered in β-cells undergoing an inflammatory attack, pointing to a crucial role for chaperones and proteins involved in RNA processing. Together with previously published studies (9), a more comprehensive knowledge of the complexity of the networks involved will make it possible to develop new strategies for the treatment of type 1 diabetes.
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
Published, MCP Papers in Press, October 5, 2007, DOI 10.1074/mcp.M700085-MCP200
1 The abbreviations used are: IL-1β, interleukin-1β; BiP/GRP78, immunoglobulin heavy chain-binding protein/glucose-regulated protein, 78 kDa; DDX1, DEAD box protein 1; DDX3X, DEAD box protein 3 X-chromosomal; ER, endoplasmic reticulum; hnRNP, heterogeneous nuclear ribonucleoprotein; HSPA8, heat shock 70-kDa protein 8; HSPA9, mortalin; IFN-
, interferon-
; NO, nitric oxide; ORP150, oxygen-regulated protein, 150 kDa; PTM, post-translational modification; PRDX, peroxiredoxin; ROS, reactive oxygen species; SOD2, manganese-superoxide dismutase; STAT1, signal transducer and activator of transcription 1; 2D, two-dimensional; 2DE, two-dimensional gel electrophoresis; TXNL1, thioredoxin-like protein 1; BLAST, Basic Local Alignment Search Tool; NCBI, National Center for Biotechnology Information; iHOP, Information Hyperlinked over Proteins. ![]()
* This work was supported in part by the Catholic University of Leuven (Grants GOA 2004/10 and EF/05/007), the Flemish Research Foundation (Grants FWO G.0084.02 and G.0233.04, a postdoctoral fellowship to C. G. and G. V. d. B. and a clinical research fellowship to C. M.), the Belgium Program on Interuniversity Poles of Attraction initiated by the Belgian State (Grants IUAP P5/17 and P6/40), the European Union Project SAVEBETA (Contract 036903 in the Framework Program 6 of the European Community), and the Centre of Excellence SymBioSys (Research Council KULeuven Grant EF/05/007). The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact. ![]()
S The on-line version of this article (available at http://www.mcponline.org) contains supplemental material. ![]()
b Both authors contributed equally to this work. ![]()
f A research associate of the Fonds National de la Recherche Scientifique (Belgium). ![]()
c To whom correspondence should be addressed: Laboratory for Experimental Medicine and Endocrinology, UZ-Gasthuisberg, Onderwijs en Navorsing, Herestraat 49, bus 902, B-3000 Leuven, Belgium. Tel.: 32-16-34-61-63; Fax: 32-16-34-60-35; E-mail lut.overbergh{at}med.kuleuven.be
| REFERENCES |
|---|
|
|
|---|
B-dependent genes in primary rat pancreatic β-cells.
J. Biol. Chem. 276, 48879
–48886
B in the islet β-cell apoptosis induced by interleukin (IL)-1β.
J. Cell Physiol. 204, 124
–130[CrossRef][Medline]
-interferon signaling pathway at the level of signal transducer and activator of transcription-1 prevents immune destruction of β-cells.
Diabetes 54, 2396
–2403[CrossRef][Medline]
B activation.
Mol. Endocrinol. 20, 1867
–1879
potentiates endoplasmic reticulum stress-induced death by reducing pancreatic beta cell defence mechanisms.
Diabetologia 49, 1229
–1236[CrossRef][Medline]