Cytoprotective effect of selective small-molecule caspase inhibitors against staurosporine-induced apoptosis.

Caspases are currently known as the central executioners of the apoptotic pathways. Inhibition of apoptosis and promotion of normal cell survival by caspase inhibitors would be a tremendous benefit for reducing the side effects of cancer therapy and for control of neurodegenerative disorders such as Parkinson's, Alzheimer's, and Huntington's diseases. The objective of this study was to discover small-molecule caspase inhibitors with which to achieve cytoprotective effect. We completed the high-throughput screening of Bionet's 37,500-compound library (Key Organics Limited, Camelford, Cornwall, UK) against caspase-1, -3, and -9 and successfully identified 43 initial hit compounds. The 43 hit compounds were further tested for cytoprotective activity against staurosporine-induced cell death in NIH3T3 cells. Nineteen compounds were found to have significant cytoprotective effects in cell viability assays. One of the compounds, RBC1023, was demonstrated to protect NIH3T3 cells from staurosporine-induced caspase-3 cleavage and activation. RBC1023 was also shown to protect against staurosporine-induced impairment of mitochondrial membrane potential. DNA microarray analysis demonstrated that staurosporine treatment induced broad global gene expression alterations, and RBC1023 co-treatment significantly restored these changes, especially of the genes that are related to cell growth and survival signaling such as Egr1, Cdc25c, cdkn3, Rhob, Nek2, and Taok1. Collectively, RBC1023 protects NIH3T3 cells against staurosporine-induced apoptosis via inhibiting caspase activity, restoring mitochondrial membrane potential, and possibly upregulating some cell survival-related gene expressions and pathways.

Introduction

Apoptosis is a genetically programmed, morphologically distinct form of cell death that can be triggered by many chemicals in a variety of cellular events, such as calcium influx, oxidative stress, cytoskeletal interference, inhibitors of protein synthesis, membrane disruption, and DNA disruption induced by radiation.1 Caspases, currently known as the central executioners of the apoptotic pathways, are cysteine proteases that cleave substrates after aspartate residues,2–4 and at least 14 different subtypes have been identified in mammalian cells. Two pathways of caspase activation for apoptosis induction have been identified. The first one starts at death receptors such as Fas.5 In the second and more common pathway, diverse proapoptotic signals converge at the mitochondrial level, inducing the translocation of cytochrome c from mitochondria to cytoplasm, thereby triggering the cascade of caspase activation.6,7 Caspase-3 is one of the most important and active executors in the apoptosis pathway.

Caspases are important targets of drug discovery8 because inhibition of apoptosis is a means of promoting cell survival. The significance of caspases in cancer therapy is implicated in two directions: while the activation of caspases leading to the cancer cell-specific apoptosis is critical to numerous and diverse cancer treatments, caspase inhibition in nonmalignant cells is a defined target for cytoprotection and reducing organ-specific toxicity during cancer therapy.9 The toxicity of normal cells exerted by chemotherapy and radiation therapy is often a limiting factor for optimal therapeutic dosage. Cytoprotection of nonmalignant cells by caspase inhibition would allow a higher upper limit of dose escalation of existing therapeutics to achieve maximal efficacy. Caspase inhibition has also been suggested for ultraviolet protection therapy to alleviate side effects such as ultraviolet-induced skin damage, inflammation, and aging.10–13 Furthermore, caspase inhibition may provide prevention or treatment against acute episodes such as stroke, traumatic brain injury, spinal cord injury, or myocardial infarction. Neuron-degenerative diseases, treatment of which may benefit from caspase inhibition, include amyotrophic lateral sclerosis, Parkinson’s, Alzheimer’s, and Huntington’s.14

To discover new caspase inhibitors, we screened Bionet’s 37,500-compound library against human caspase-1, -3, and -9 with the fluorogenic substrates Ac-YVAD-AMC, Ac-DEVD-AMC, and Ac-LEHD-AMC, and identified 43 initial hit compounds. Using staurosporine-induced apoptosis assay in NIH3T3 cell line, we discovered 19 small molecules with cell protective activities, in which compound RBC1023, a selective caspase-3 inhibitor, had shown a broad protection spectrum across many different cell types. Further mechanistic studies indicated that RBC1023 protected NIH3T3 cells from staurosporine-induced caspase cleavage and activation and restored staurosporine-induced loss of mitochondrial membrane potential. Moreover, we investigated the changes in gene expression in NIH3T3 cells treated with staurosporine in the presence or absence of RBC1023 using a DNA microarray analysis to understand the molecular mechanisms of the cytoprotective effect of RBC1023.

Materials and methods

Chemicals and reagents

Bionet’s 37,500-compound library was obtained from Key Organics Limited (Camelford, Cornwall, UK). Purified cas-pase-1, -3, and -9 enzymes; fluorogenic caspase substrates Ac-YVAD-AMC, Ac-DEVD-AMC, and Ac-LEHD-AMC; and control caspase inhibitor Ac-DEVD-CHO, Ac-LEHD-CHO, and Ac-YVAD-CHO were all obtained from Biomol International (Enzo Life Sciences, Inc., Farmingdale, NY, USA). Staurosporine was obtained from Sigma-Aldrich (St Louis, MO, USA). Cell Titer-Glo® Luminescent Cell Viability Assay Kit and Apo-ONE® Homogeneous Caspase-3/7 Assay Kit were obtained from Promega Corporation (Fitchburg, WI, USA). MTT (3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide) Cell Proliferation Assay Kit was obtained from Cayman Chemical (Ann Arbor, MI, USA). JC-10 Mitochondrial Membrane Potential Assay Kit was obtained from Abcam (Cambridge, UK). Anti-caspase-3 polyclonal antibody (catalog#9665, 1:2000 dilutions) was obtained from Cell Signaling Technology (Danvers, MA, USA). Anti-α-tubulin monoclonal antibody (catalog#T6074, 1:5000 dilutions) was obtained from Sigma-Aldrich. The secondary antibodies Alexa Fluor® 633 goat anti-rabbit immunoglobulin (Ig)G (heavy and light chains [H+L]) and Alexa Fluor® 532 goat anti-mouse IgG (H+L) were obtained from Thermo Fisher Scientific (Life Technologies, Waltham, MA, USA). Other chemicals were of reagent grade and were obtained from Sigma-Aldrich or VWR International (Radnor, PA, USA).

High-throughput screening (HTS) and hit confirmation

The primary HTS against caspase-1, -3, and -9 enzymes was performed using the Bionet compound library. In brief, 5 μL of reaction buffer (50 mM 4-(2-hydroxyethyl)piperazine-1-ethanesulfonic acid (HEPES) pH 7.4, 100 mM NaCl, 0.1% 3-[(3-cholamidopropyl) dimethylammonio]-1-propanesulfonate (CHAPS), 1 mM ethylenediaminetetraacetic acid (EDTA), and 3 mM dithiothreitol) with 1X enzyme was added to each well of Genetix 7020 low-volume 384-well plates (Molecular Devices LLC, Sunnyvale, CA, USA). Five nanoliters of 10 mM of each compound from the Bionet 37,500-compound library was delivered to each well of the Genetix plates (final compound concentration 10 μM) by acoustic droplet ejection liquid handler (Echo 550, Labcyte Inc., Sunnyvale, CA, USA). The plates were shaken to mix for 1 minute. One microliter of positive control inhibitors (Ac-YVAD-CHO, Ac-DEVD-CHO, and Ac-LEHD-CHO) was added to the control wells. Ten nanoliters of 5 mM of caspase substrates was delivered to the plates by Echo 550 to start the reaction. The fluorescent signal was measured by EnVision® 2104 Multilabel Reader (PerkinElmer Inc., Waltham, MA, USA) at an excitation wavelength of 355 nm and an emission wavelength of 460 nm at the 2-hour time point. The fluorescent signal generated was proportional to the caspase-1, -3, and -9 activities. The initial hits from the HTS were cherry-picked and further tested in a ten-dose half maximal inhibitory concentration (IC-50)-mode against caspase-1, -3, and -9.

Cells, cell culture, and compound treatment

NIH3T3 mouse fibroblast, PC12 rat pheochromocytoma cells, and U87MG human glioblastoma cells were obtained from American Type Culture Collection (ATCC) (Manassas, VA, USA). SH-SY5Y cell line was kindly provided by Dr Honglin Zhou at the University of Pennsylvania (Philadelphia, PA, USA). NIH3T3 cells were grown in Dulbecco’s Modified Eagle’s Medium, PC12 cells were grown in RPMI 1640 (Roswell Park Memorial Institute medium 1640), U87MG cells were grown in Modified Eagle’s Medium, and SH-SY5Y cells were grown in a 1:1 mixture of Modified Eagle’s Medium and F12 media. All culture media were supplemented with 10% heat-inactivated fetal bovine serum, 100 μg/mL penicillin, and 100 μg/mL streptomycin. Cultures were maintained at 37°C in a humidified atmosphere of 5% CO2 and 95% air. For most experiments, unless specified, 5,000/well of logarithmically growing cells were seeded in each of the 384-well plates containing proper medium plus 1% fetal bovine serum with or without compound. Dimethyl sulfoxide ([DMSO] #0.3%) was added to all the control cultures.

Cell Titer-Glo® Luminescent Cell Viability Assay and adenosine triphosphate (ATP) measurement

NIH3T3, PC12, U87MG, or SH-SY5Y cells were seeded overnight in triplicates in 384-well plates. The cells were incubated with compounds for 24 hours. A volume of Cell Titer-Glo® reagent equal to the volume of cell culture medium was added to each well. Contents were mixed for 2 minutes on an orbital shaker to induce cell lysis. The plates were incubated at room temperature for 10 minutes to stabilize luminescent signal. Luminescence was recorded by the EnVision® 2104 Multilabel Reader. The numbers of viable cells in culture were determined based on quantitation of the ATP present in each well. For quantifying cellular ATP contents in NIH3T3 cells, an ATP standard curve was generated starting at 1 μM in serial ten-fold dilutions in culture medium. The luminescence signals were converted into ATP concentration based on the ATP standard curve.

Apo-ONE® Homogeneous Caspase-3/7 Assay

NIH3T3 cells were seeded overnight in triplicates in 384-well plates. The cells were incubated with compounds for 6 hours. A volume of Apo-ONE® reagent equal to the volume of cell culture medium was added to each well. The contents were gently mixed using a plate shaker at 300 rpm for 30 seconds. The plates were incubated at room temperature for 30 minutes. Upon the cleavage of caspase substrate Z-DEVD-rhodamine 110, bis-(N-CBZ-L-aspartyl-L-glutamyl-L-valyl-L-aspartic acid amide) by caspase-3/7 activity, the fluorescence of each well was measured at an excitation wavelength of 485 nm and an emission wavelength of 535 nm using the EnVision® 2104 Multilabel Reader. The amount of fluorescence generated was proportional to the amount of caspase-3/7 cleavage activity present in the sample.

MTT cell proliferation assay

NIH3T3 cells were seeded in triplicates in 96-well plates at a density of 2×104 cells/well in 100 μL of culture medium with or without compounds to be tested. The cells were cultured in a CO2 incubator at 37°C for 24 hours. Ten microliters of MTT reagent was added to each well and mixed gently for 1 minute. The cells were incubated for 4 hours at 37°C. After incubation, the formazan produced in the cells appeared as dark crystals in the bottom of the wells. The culture medium from each well was aspirated carefully and 100 μL of crystal-dissolving solution was added to each well to dissolve the formazan crystals. The absorbance of each sample was quantified at 590 nm using the EnVision® 2104 microplate reader. The absorbance intensity was proportional to the amount of viable cells.

JC-10 mitochondrial membrane potential assay

Mitochondrial membrane potential (Δψm) was evaluated by cationic dye JC-10. In normal cells, JC-10 aggregates in mitochondria, fluorescing red. In apoptotic cells, JC-10 accrues in the cytosol, as a green fluorescing monomer. Briefly, NIH3T3 cells were plated in triplicates in growth medium in 384-well plates overnight. The cells were treated with or without compounds for 6 hours. JC-10 dye-loading solution was added and the plates were incubated in a 37°C, 5% CO2 incubator for 30 minutes. Assay buffer was added into the dye-loading plates and the fluorescence intensities were monitored at excitation/emission =485/535 nm and 485/590 nm for ratio analysis in the EnVision® 2104 Multilabel microplate reader.

Western blot analysis

NIH3T3 cells were seeded in six-well plates (0.8×106 cells/well) and treated with or without compounds for 6 hours. After the treatment, cells were washed with cold phosphate-buffered saline and lysed in cell lysis buffer (150 mM NaCl, 1% Triton X-100, 50 mM Tris pH 7.5, proteinase inhibitor cocktail, 0.1 mM phenylmethylsulfonyl fluoride). Protein concentrations of the cell lysates were measured using Bio-Rad protein assay dye reagent concentrate (Bio-Rad Laboratories Inc., Hercules, CA, USA). About 30 μg of total protein for each sample was separated under reducing conditions by sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) and transferred to a nitrocellulose membrane (Life Technologies). The blots were blocked in 3% nonfat milk for 2 hours at room temperature with gentle shaking. This was followed by incubation with the respective antibody overnight at 4°C and then incubation with Alexa Fluor® 633 goat anti-rabbit IgG or Alexa Fluor® 532 goat anti-mouse IgG secondary antibody for 1 hour. Fluorescence signals of the immune complexes were scanned with a Typhoon 8610 variable mode imager (GE Healthcare Bio-Sciences Corp, Piscataway, NJ, USA).

DNA microarray analysis of gene expression profile

NIH3T3 cells were seeded in 10 cm dishes (2×106 cells/dish) and treated with DMSO or testing compounds. The cells were pelleted down by centrifuging and washed by cold phosphate-buffered saline. The cell pellets were frozen using dry ice, and submitted to the Penn Molecular Profiling Facility (University of Pennsylvania School of Medicine) for RNA extraction and transcript profiling. Microarray analysis was conducted in duplicate. RNA was extracted with TRIzol® by Life Technologies (Thermo Fisher Scientific) and purified using the Qiagen RNeasy Kit (Qiagen NV, Venlo, the Netherlands); total RNA yields ranged from 15.5 to 53.4 μg per sample with A260/280 of 2.07–2.18. Transcript profiling was conducted as described in the Affymetrix GeneChip Expression Analysis Technical Manual (Affymetrix Inc., Santa Clara, CA, USA).15 Briefly, 100 ng of total RNA was converted to first-strand complementary DNA (cDNA) using reverse transcriptase primed by poly (T) and random oligomers that incorporated the T7 promoter sequence. Second-strand cDNA synthesis was followed by in vitro transcription with T7 RNA polymerase for linear amplification of each transcript, and the resulting cRNA was converted to cDNA, fragmented, and biotinylated by terminal transferase end labeling. Labeled cDNA was added to Affymetrix hybridization cocktails, heated at 99°C for 5 minutes, and hybridized for 16 hours at 45°C to GeneChip Mouse Gene 1.0ST Arrays (Affymetrix Inc.). The microarrays were then washed at low (6× SSPE [sodium chloride, sodium phosphate monobasic, and ethylenediaminetetraacetic acid]) and high (100 mM MES [2-[N-morpholino] ethanesulfonic acid, tris base, ethylenediaminetetraacetic acid, and sodium dodecyl sulfate], 0.1 M NaCl) stringency, and stained with streptavidin/phycoerythrin. Fluorescence was amplified by adding biotinylated anti-streptavidin and an additional aliquot of streptavidin/phycoerythrin stain. A confocal scanner was used to collect the fluorescence signal after excitation at 570 nm. The intensity of the fluorescence signal was proportional to the levels of gene expression.

Statistical analysis

Statistical analysis in this study was performed using the unpaired Student’s t-test or one-way analysis of variance (ANOVA). The values are represented as mean ± standard deviation. A P-value <0.05 was considered to be statistically significant.

Results

HTS and hit confirmation

Using the Bionet 37,500-compound library, we successfully completed HTS against human caspase-1, -3, and -9. The screening was of high quality, with a robust Z-factor >0.6 (Table 1).16 By using simple signal-to-noise analysis, compounds with response signal >3 times the standard deviation of the noise were designated as “hits”. Totals of 453 compounds, 727 compounds, and 1,149 compounds were selected as initial hits for caspase-1, -3, and -9, respectively, which represent the hit rates of 1.20%, 1.94%, and 3.06%, respectively (Table 1). The hits were further confirmed in ten-dose IC50 mode in the presence of three different reducing agents (3 mM dithiothreitol, 5 mM β-mercaptoethanol, 5 mM cysteine) to eliminate any reactive electrophilic compounds.17 There were 143, 178, and 273 confirmed hits that inhibited caspase-1, -3, or -9, respectively. This provided overall hit rates of 0.38%, 0.47%, and 0.73%, respectively, which are typical for biochemical HTS against a diversity library. We further eliminated some compounds from these confirmed hits using three stringent selection processes: compounds with >20% inhibition on cytochrome P450 (Cyp) enzymes such as Cyp2C19, Cyp1A2, Cyp2A6, Cyp2C9, Cyp3A4, or Cyp2D6 at 10 μM (data not shown); compounds containing undesirable structures and metabolically fragile functionality (such as reactive chlorides, Michael acceptors, thiols, esters, quaternary salts, and detergents); and compounds showing significant toxicity against mouse embryonic fibroblast cells (data not shown). A final set of 43 compounds was selected for further testing in a cytoprotection assay (Table 2).

Table 1

Primary high-throughput screening (HTS) of caspase-1, -3, and -9

Target	Caspase-1	Caspase-3	Caspase-9
Compounds screened	37,500	37,500	37,500
HTS hits	453	727	1,149
Z-factor	>0.6	>0.6	>0.6
Primary hit rate	1.21%	1.94%	3.06%
Confirmation under three different reducing agents	143	178	273
Confirmation rate	0.38%	0.47%	0.73%

Table 2

Caspase inhibition activity IC50 (μM) of 43 hit compounds from the in vitro caspase-1, -3, and -9 screening assays

Compound #	Caspase-1	Caspase-3	Caspase-9
RBC1001	1.84	8.36	1.84
RBC1002	0.36	>20	0.36
RBC1003	>20	>20	0.14
RBC1004	14.25	>20	0.65
RBC1005	>20	>20	0.55
RBC1006	>20	>20	1.08
RBC1007	>20	>20	0.35
RBC1008	>20	>20	0.15
RBC1009	>20	>20	0.06
RBC1010	>20	>20	0.72
RBC1011	>20	>20	3.25
RBC1012	>20	>20	0.40
RBC1013	>20	>20	0.50
RBC1014	>20	>20	1.20
RBC1015	>20	1.13	>20
RBC1016	>20	>20	1.38
RBC1017	>20	>20	3.31
RBC1018	12.90	>20	0.60
RBC1019	12.08	>20	4.77
RBC1020	>20	>20	1.38
RBC1021	>20	1.41	>20
RBC1022	>20	>20	2.66
RBC1023	>20	0.38	>20
RBC1024	>20	4.84	3.05
RBC1025	12.00	>20	0.56
RBC1026	>20	>20	0.61
RBC1027	>20	>20	0.06
RBC1028	>20	>20	2.82
RBC1029	>20	>20	0.15
RBC1030	17.06	>20	1.80
RBC1031	>20	>20	3.94
RBC1032	>20	8.08	0.03
RBC1033	14.72	>20	2.34
RBC1034	>20	8.94	3.42
RBC1035	>20	3.05	>20
RBC1036	14.58	>20	0.15
RBC1037	>20	1.37	>20
RBC1038	14.02	>20	0.39
RBC1039	>20	>20	0.51
RBC1040	>20	>20	1.76
RBC1041	2.47	>20	5.93
RBC1042	>20	7.74	2.51
RBC1043	>20	0.40	>20

Abbreviation: IC50, half maximal inhibitory concentration.

Cytoprotective activities of the caspase inhibitors

Staurosporine is a broad-spectrum inhibitor of protein kinases, and has been widely used in the induction of apoptosis in diverse cellular models.18–22 Staurosporine preferentially activates the mitochondrial apoptotic pathway relying on caspase activation and causes cell death. To confirm the caspase inhibition of the hit compounds in cell-based assay, NIH3T3 cells were treated with 20 μM of each of the 43 hit compounds in the presence or absence of 1 μM staurosporine for 24 hours. Then, cell viability was assessed by Cell Titer-Glo Luminescent Cell Viability Assay. As shown in Figure 1, while staurosporine treatment alone caused over 90% cell death, 19 compounds showed significant protection against staurosporine-induced cell death in NIH3T3 cells (P<0.01). Several compounds that are toxic to cells alone were removed from the hit list (data not shown). Compound RBC1020 (a caspase-9 selective inhibitor) and RBC1023 (a caspase-3 selective inhibitor) with no cytotoxicity alone demonstrated 7.5-fold and 5.9-fold increases respectively of cell viability compared to cells treated with staurosporine alone, indicating a strong cytoprotective activity of these two compounds. Caspase inhibitors could have potential application in treating neurodegenerative diseases such as Parkinson’s and Alzheimer’s diseases, which are usually associated with neuron cell death. Therefore, we extended the study to further examine the cytoprotective activity of the leads RBC1020 and RBC1023 in several neural origin cell types. PC12, U87MG, and SH-SY5Y cells were treated with DMSO, staurosporine alone, RBC1020 or RBC1023 alone, or RBC1020 or RBC1023 plus staurosporine for 24 hours and then the cell viability was examined. As shown in Figure 2, quantitative measurement of cell viability showed that RBC1023 increased the percentage of staurosporine-treated viable cells from 4% to 66% in PC12 cells, from 8% to 71% in U87MG cells, and from 6% to 47% in SH-SY5Y cells. Interestingly, RBC1020 had no or minimum cytoprotection in these three cell lines. It is still unknown why RBC1020 and RBC1023 showed different effects in different cell lines. Additional studies are needed to investigate the mechanism that causes this difference. These results demonstrated that RBC1023, a selective caspase-3 inhibitor, exhibits strong cytoprotective effect against staurosporine-induced cell death in multiple cell lines.

Figure 1

Screening results for cytoprotection of 43 hit compounds against staurosporine (S)-induced cell death in NIH3T3 cells.

Notes: (A) NIH3T3 cells were treated with dimethyl sulfoxide, 20 μM of each of the 43 hit compounds, 1 μM of S, or 20 μM of each hit compound plus 1 μM of S for 24 hours. The cells were measured for viability by Cell Titer-Glo® luminescent assay (Promega Corporation, Fitchburg, WI, USA). The cell viability fold change between the combination group and S treatment group is presented. (B) Nineteen compounds were shown to have significant protection against S-induced cell death. The P-values between the two indicated groups were calculated by Student’s t-test. P<0.01 demonstrates a highly significant difference between the indicated groups.

Abbreviation: vs, versus.

Figure 2

RBC1023 protects STAU-induced cell death in PC-12, U87MG, and SH-SY5Y cells.

Notes: PC-12 (A), U87MG (B), and SH-SY5Y (C) cells were treated with DMSO, 20 μM of RBC1023, 1 μM of STAU, or 20 μM of RBC1023 plus 1 μM of STAU, 20 μM of RBC1020 and 20 μM of RBC1020 plus 1 μM of STAU for 24 hours. The cells were measured for viability by Cell Titer-Glo® luminescent assay (Promega Corporation, Fitchburg, WI, USA). *P<0.01 (a highly significant difference between the indicated groups) by Student’s t-test.

Abbreviations: DMSO, dimethyl sulfoxide; STAU, staurosporine.

Dose-dependent and time-dependent study of RBC1023 against staurosporine-induced cell death

Since the staurosporine/NIH3T3 apoptosis induction system is a very well-established model system for apoptosis-related mechanism studies, we next focused on the cytoprotection mechanism studies in the NIH3T3 cell line. We first investigated the protective effect of different concentrations of RBC1023 against different concentrations of staurosporine-induced cell death. NIH3T3 cells were treated with different dose combinations of RBC1023 (50 μM to 1 μM, three-fold dilution, ten doses) and staurosporine (0.768 μM or 2 μM) for 24 hours before the cell viability was measured. As shown in Figure 3A, RBC1023 started to show significant cytoprotection against the staurosporine-induced cell death at about 2 μM, when the concentration of staurosporine was at 0.768 μM. The cytoprotective effect by RBC1023 increased in a dose-dependent manner. To understand if the pretreatment or posttreatment of the caspase inhibitor would affect its cytoprotective effect, RBC1023 was added to cells at 1 hour before staurosporine treatment or 1 hour post-staurosporine treatment. The cell viability was measured after a total of 24 hours’ incubation. Interestingly, there were no significant differences of cytoprotective activity by RBC1023 between the 1-hour pretreatment and and 1-hour posttreatment group with staurosporine (Figure 3B).

Figure 3

Dose-response protection and treatment timing study of RBC1023 against STAU-induced cell death in NIH3T3 cells.

Notes: (A) NIH3T3 cells were treated with a combination of various concentrations of RBC1023 (50 μM to 1 μM, three-fold dilution, ten doses) and STAU (0.768 μM or 2 μM) for 24 hours. (B) NIH3T3 cells were treated with DMSO; 20 μM of RBC1023; 1 μM of STAU; or 20 μM of RBC1023 and 1 μM of staurosporine simultaneously; pretreated with 20 μM of RBC1023 for 1 hour and then treated with 1 μM of staurosporine; or pretreated with 1 μM of staurosporine for 1 hour and then treated with 20 μM of RBC1023. The cultures were then incubated for a total of 24 hours. The cell viability was measured by Cell Titer-Glo® assay (Promega Corporation, Fitchburg, WI, USA). *P<0.01 (a highly significant difference between the indicated groups) by Student’s t-test.

Abbreviations: DMSO, dimethyl sulfoxide; EC50, half maximal effective concentration; STAU, staurosporine.

RBC1023 blocks staurosporine-induced caspase activation in cell-based assays

To confirm RBC1023 as a caspase inhibitor in cell-based assays, we examined the caspase activity/activation-inducing effect in NIH3T3 cells that were treated with staurosporine in the presence or absence of RBC1023. As shown in Figure 4A, 6-hour staurosporine treatment induced significant caspase activation in NIH3T3 cells, as evaluated by Apo-ONE® Homogeneous Caspase-3/7 Assay using Z-DEVD-rhodamine 110 as a caspase substrate. Treatment of NIH3T3 cells with staurosporine for 6 hours also resulted in cleavage/activation of caspase-3, as was evident by the appearance of a 17 kD proteolytic product of caspase-3 as determined by Western blot analysis (Figure 4B). RBC1023 dramatically inhibited staurosporine-induced caspase-3 cleavage and activation (Figure 4B). These results confirmed that RBC1023 is a caspase inhibitor in cells.

Figure 4

RBC1023 inhibits STAU-induced caspase-3 cleavage and activation in NIH3T3 cells.

Notes: NIH3T3 cells were treated with DMSO, 20 μM of RBC1023, 1 μM of STAU, or 20 μM of RBC1023 plus 1 μM of STAU for 6 hours. Caspase activity was measured by Apo-ONE® Homogeneous Caspase-3/7 Assay (Promega Corporation, Fitchburg, WI, USA) (A), and processing of caspase-3 was detected by Western blotting with anti-caspase-3 antibody (B). Positions of the intact and cleaved forms of capase-3 are indicated by arrows. α-tubulin was used as a control for protein loading. One-way analysis of variance indicates: P>0.05 (DMSO versus [vs] RBC1023), P>0.05 (DMSO vs RBC1023 + STAU), and P<0.01 (STAU vs RBC1023 + STAU).

Abbreviations: DMSO, dimethyl sulfoxide; STAU, staurosporine.

Figure 6

Hierarchical clustering of 200 upregulated and downregulated genes from NIH3T3 samples treated with DMSO, RBC1023, STAU, or RBC1023 plus STAU.

Notes: NIH3T3 cells were treated with DMSO, 20 μM of RBC1023 (c23), 1 μM of STAU, or 20 μM of RBC1023 (c23) plus 1 μM of STAU for 20 hours. Total RNA was isolated and converted to complementaryDNA by reverse transcriptase reaction. The cDNA products were fragmented and hybridized to Affymetrix GeneChip Mouse Gene 1.0ST (Affymetrix Inc., Santa Clara, CA, USA).15 Clustering of 200 upregulated and downregulated genes that exhibit significant differences in all groups was determined by one-way analysis of variance. The clustering was based on fold changes within each paired data set. The color scale of the log2 ratios is shown at the bottom.

Abbreviations: DMSO, dimethyl sulfoxide; STAU, staurosporine.

RBC1023 restores NIH3T3 cells from the staurosporine-induced loss of cell viability and cellular ATP contents

To investigate if RBC1023 can rescue staurosporine-induced decrease of cell viability, we performed the MTT assay, in which the tetrazolium dye, MTT, is reduced to formazan by intracellular nicotinamide adenine dinucleotide phosphate (NADPH) oxidoreductases, which correlates with cellular metabolic activity and the number of viable cells present. As shown in Figure 5A, staurosporine treatment for 24 hours reduced the cell viability dramatically. RBC1023 co-treatment significantly rescued the loss of cell viability induced by staurosporine in NIH3T3 cells. We also performed ATP measurement by generating an ATP standard curve using Cell Titer-Glo® luminescent assay. As shown in Figure 5B, RBC1023 clearly restored the loss of ATP contents in cells treated with staurosporine. One-way ANOVA of the MTT and ATP assay data suggested that RBC1023 induces normalization or partial recovery of cell viability and ATP contents. Our results here further confirmed the cytoprotective effect of RBC1023 demonstrated by Cell Titer-Glo® assay.

Figure 5

RBC1023 restores NIH3T3 cells from STAU-induced loss of cell viability, ATP contents, and disruption of mitochondrial membrane potential.

Notes: NIH3T3 cells were treated with DMSO, 20 μM of RBC1023, 1 μM of STAU, or 20 μM of RBC1023 plus 1 μM of STAU for 24 hours (A and B), or 6 hours (C). The cell viability was measured by MTT assay (A) and ATP content assay (B), and the mitochondrial membrane potential was measured by JC-10 mitochondrial membrane potential assay (C). One-way analysis of variance indicates: (A) P>0.05 (DMSO versus [vs] RBC1023), P<0.01 (DMSO vs RBC1023 + STAU), and P<0.01 (STAU vs RBC1023 + STAU); (B) P<0.01 (DMSO vs RBC1023), P<0.001 (DMSO vs RBC1023 + STAU), and P<0.001 (STAU vs RBC1023 + STAU); (C) P>0.05 (DMSO vs RBC1023), P>0.05 (DMSO vs RBC1023 + STAU), and P<0.001 (STAU vs RBC1023 + STAU).

Abbreviations: ATP, adenosine triphosphate; DMSO, dimethyl sulfoxide; MTT, (3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide); STAU, staurosporine; Δψm, mitochondrial membrane potential.

RBC1023 restores NIH3T3 cells from staurosporine-induced disruption of mitochondrial membrane potential

Disruption of mitochondrial membrane potential (Δψm), which is observed at an early stage of apoptosis in many cells, results in opening of the permeability transition pores (PTPs), causing a local disruption of the outer mitochondrial membrane.23,24 As a consequence, the soluble inter-membrane proteins including cytochrome c are released, and these proteins contribute to the activation of caspases for chromatin condensation and DNA fragmentation.7 Based on this theory, we next assessed if RBC1023 restores NIH3T3 cells from staurosporine-induced loss of mitochondrial membrane potential. The mitochondrial membrane potential changes were measured by the cationic, lipophilic JC-10 dye. As shown in Figure 5C, NIH3T3 cells that were treated with 1 μM of staurosporine for 6 hours significantly lost mitochondrial membrane potential, whereas co-incubating the cells with RBC1023 and staurosporine apparently prevented staurosporine-induced mitochondrial membrane impairment. One-way ANOVA of the data also supports this statement. Our findings suggest that RBC1023 mediates mitochondrial protection in addition to its direct caspase inhibitory effect.

Investigation of gene profile for the cytoprotective effect of RBC1023 by DNA microarray analysis

To further explore the potential molecular mechanisms of the cytoprotection effect by RBC1023, we performed DNA microarray analysis using Affymetrix GeneChip Mouse Gene 1.0ST, which covers 26,166 coding transcripts. NIH3T3 cells were treated by one of four conditions as follows: DMSO, RBC1023, staurosporine, or RBC1023 plus staurosporine. For each condition, the NIH3T3 cells were exposed to the indicated treatment for 20 hours. RNA was harvested from the cells treated under each condition and then converted to cDNA, which was used to query the genes on the Affymetrix GeneChip Mouse Gene 1.0ST Array. The general microarray results showed that staurosporine treatment dramatically changes the global gene expression profiling, and the changes caused by treatment of staurosporine were significantly restored by the co-treatment of RBC1023 in NIH3T3 cells. The gene expressions between each group, analyzed using supervised clustering of the 200 upregulated and downregulated messenger (m)RNAs that exhibited the most significant differences in all groups, were determined by one-way ANOVA. As shown in the heat map (Figure 6), the control and RBC1023 alone treatment samples showed similar clustering of these specific gene expressions. However, there was a dramatic increase or decrease in the expression of many genes in the staurosporine treatment samples. Importantly, the expression levels of these downregulated or upregulated genes caused by treatment of staurosporine were significantly restored by the co-treatment of RBC1023 in NIH3T3 cells. Molecular functions of the top 200 up- and downregulated gene transcripts are listed in Table S1. The gene list includes some apoptosis and cellular signaling-related genes such as Htra2 (Omi), caspase-12, Etv1, Egr1, Nek7, Upk1a, Skap2, and Tuba1b (Table S1). Furthermore, we summarized the top 50 genes upregulated (ranked by fold-change) in the RBC1023-plus-staurosporine treatment group in comparison with the staurosporine-alone treatment group, in which displays ~2.5- to ~4.9-fold increases can be seen in Table 3. Interestingly, no known proapoptotic genes can be seen in this list (Table 3). We hypothesize that the upregulation of some cell survival genes may be partially responsible for the protective effect of RBC1023 against staurosporine-induced cell death.

Table 3

Top 50 genes upregulated (RBC1023 + STAU group versus STAU group)

Gene symbol	Fold-change (RBC1023 + STAU vs STAU)	P-value (RBC1023 + STAU vs STAU)	Function
Ly6a	4.94653	5.35E-06	Lymphocyte antigen 6 complex
Lrrn4cl	4.09678	2.17E-05	Function unknown
Egr1	4.0647	4.13E-05	Differentiation and mitogenesis
Plat	3.90703	0.00025859	Plasminogen activator
Bdkrb1	3.89931	2.63E-06	Inflammatory responses
Rhob	3.5974	0.000108495	Ras signaling cascade
Nek2	3.42607	6.47E-06	Serine/threonine-protein kinase that is involved in mitotic regulation
Serpinb1a	3.41146	3.16E-05	Regulates the activity of the neutrophil proteases elastase
Taok1	3.35134	0.00379721	Serine/threonine-protein kinase
Bhlhe40	3.29921	0.000687826	Cell differentiation
Angptl2	3.28147	6.37E-05	Growth factors
4930547N16Rik	3.24155	7.27E-06	Function unknown
Plxna2	3.20986	3.55E-06	Nervous system development
Aurka	3.17042	0.00014278	Serine/threonine-protein kinase 6
Tceal1	3.14666	1.88E-05	Transcriptional regulation
Il18	3.07789	8.13E-05	A proinflammatory cytokine
Tacr2	3.04967	5.85E-05	Receptors for tachykinins
Ccnb1	3.00447	1.19E-05	Control of the cell cycle at the G2/M (mitosis) transition
Hist1h1c	2.9718	2.15E-05	Chromatin remodeling, nucleosome spacing and DNA methylation
Mfsd6	2.96238	4.17E-06	Function unknown
Cdc25c	2.95929	8.77E-06	Cell mitotic control
Arhgef12	2.94683	0.00060089	Stimulate Rho-dependent signals
Tnfaip6	2.92751	2.48E-05	Involved in cell-cell and cell-matrix interactions
Depdc1a	2.92024	5.13E-05	Transcriptional regulation
Cep55	2.83691	4.68E-06	Mitotic exit and cytokinesis
Gas2l3	2.81642	8.16E-07	Promote and stabilize the formation of the actin and microtubule network
Ttll7	2.80163	3.77E-05	Neurite growth
Mastl	2.79123	3.11E-06	A regulator of mitosis entry and maintenance
Cdkn3	2.75997	4.81E-06	Cell cycle regulation
Kif2c	2.73935	4.82E-05	Regulates the turnover of microtubules during mitosis
Anxa8	2.72379	0.000119099	Involved in the blood coagulation cascade
Cyp39a1	2.71735	9.70E-06	Involved in drug metabolism and synthesis of cholesterol, steroids
Bend6	2.69474	0.000311133	Function unknown
Hyls1	2.69324	1.54E-06	Required for the formation of cilia
6720463M24Rik	2.66567	7.35E-06	Function unknown
Cenpe	2.65857	8.09E-07	Essential for the maintenance of chromosomal stability
Adamtsl5	2.65327	2.33E-06	Function unknown
Tlr4	2.64933	4.06E-05	Mediate the innate immune response to LPS
Dusp6	2.61857	0.000175186	Inactivates MAP kinases
Kif20b	2.60411	0.000105112	Required for completion of cytokinesis
Stxbp4	2.60053	0.0106744	Plays a role in the translocation of transport vesicles
Aldh3a1	2.59577	5.93E-05	The metabolism of corticosteroids, biogenic amines, neurotransmitters, and lipid
Prc1	2.58871	1.85E-05	Involved in cytokinesis
Ndc80	2.58587	0.000245121	Required for chromosome segregation and spindle checkpoint activity
Glmn	2.58467	0.00018307	Essential for normal development of the vasculature
Plk1	2.54604	2.21E-05	Regulators of cell cycle progression, mitosis, cytokinesis, and the DNA damage
Atf7ip2	2.54547	1.85E-05	Modulates transcription regulation and chromatin formation
Mir15b	2.53671	0.00251676	Involved in post-transcriptional regulation of gene expression
Suox	2.53666	9.13E-05	Catalyzes the oxidation of sulfite to sulfate
Serpinb9b	2.53635	0.000123141	Inhibits the activity of the effector molecule granzyme

Abbreviation: STAU, staurosporine.

Discussion

In this study, we identified 19 caspase inhibitors that showed cytoprotection against staurosporine-induced cell death by screening Bionet’s 37,500-compound library against caspase-1, -3, and -9, and then through multiselective processes. RBC1023, a selective caspase-3 inhibitor, showed dose-dependent cytoprotection against staurosporine-induced cell death in different types of cell lines. We also confirmed that RBC1023 protects NIH3T3 cells from the staurosporine-induced caspase-3 cleavage and activation. These results indicate that reduced apoptotic cell death and increased cell proliferation are attributed to the inhibition of caspase activation by RBC1023. Interestingly, RBC1023 protected against cell death even up to 1 hour after staurosporine treatment.

Mitochondria play a central role in apoptosis,25 and there are reports that demonstrate the critical role of mitochondria in cytoprotection.26,27 We evaluated the possible correlation between the protection by RBC1023 and the mitochondrial function. First, our MTT assay results demonstrated that RBC1023 co-treatment was able to rescue the staurosporine-triggered loss of cell viability, suggesting RBC1023 restored the loss of the enzyme activity in mitochondria that reduces MTT during staurosporine treatment. Secondly, we found that co-treatment with RBC1023 and staurosporine resulted in a significant increase of cellular ATP content in comparison with the staurosporine treatment group. Our results suggest that RBC1023 restored the loss of ATP production during the staurosporine treatment. Furthermore, our results indicated that RBC1023 restored staurosporine-induced disruption of mitochondrial membrane potential. It is well known that mitochondrial dysfunction is the primary cause of staurosporine-induced apoptosis. A critical factor mediating mitochondrial dysfunction is the opening of mitochondrial PTP (mPTP). The opening of the mPTP can lead to a bioenergetic, biosynthetic, and redox crisis in a cell that can directly threaten the survival of the cell.28 When the mPTP is open, the mitochondrial inner membrane becomes permeable to protons, which then lead to the uncoupling of the electron respiratory chain and the collapse of membrane potential, which in turn leads to a cessation of ATP generation in mitochondria.28,29 In the RBC1023-pretreated NIH3T3 cells, the staurosporine-induced loss of mitochondrial membrane potential and decline in ATP levels was alleviated, supporting the notion that the cytoprotection of RBC1023 is, in part, due to the prevention of mitochondrial dysfunction. Upon activation of mPTP, functional breakdown and morphological disintegration of mitochondria occur, thus initiating cell death.30 Another threatening consequence of the altered mitochondrial permeability is the release of apoptogenic proteins from the mitochondrial inter-membrane space into the cytosol.28,29 Cytochrome c is associated with the inner mitochondrial membrane and serves as an essential component of the electron transfer chain. With opening of the mPTP and translocation of cytochrome c into the cytosol, mitochondrial function is compromised. However, in this study, the release of cytochrome c from the mitochondrial matrix into the cytosol by staurosporine was not significantly inhibited by RBC1023 (data not shown). This was possibly because caspase-3 is activated downstream of cytochrome c release and apoptosome formation in the mitochondrial pathway. It has been reported that B-cell lymphoma 2 protects endothelial cells against gamma radiation-induced caspase activation via a Raf-MEK-ERK cell survival signaling pathway that is independent of cytochrome c release.31 It is possible that RBC1023 also exerts its cytoprotective effect via some cell survival signaling pathways that are independent of the mechanism of the blockage of cytochrome c release. Taken together, we conclude that RBC1023, in addition to its direct caspase inhibitory effect, has some mitochondrial protective effect, which is at least partially responsible for its cytoprotection. For our next study, the mitochondrial oxygen consumption will be measured to directly evaluate mitochondrial function and cellular respiratory activity. Whether RBC1023’s mitochondrial protective effect is related to its caspase inhibitory effect could be an interesting subject for future study.

The DNA microarray analysis in the present study demonstrates that RBC1023 can recover the loss of expression of many genes induced by staurosporine treatment, correlating with its cytoprotection. The top 50 genes that were upregulated in the RBC1023-plus-staurosporine treatment group compared to the staurosporine-alone treatment group include several functional categories, such as cell differentiation and growth genes (Egr1, Bhlhe40, and Angptl2), cell cycle regulation genes (Ccnb1, Cdc25c, Cep55, Mastl, Cdkn3, Kif2c, and Plk1), cellular signaling genes (Rhob, Nek2, Taok1, Aurka, and Arhgef12), and transcription regulation genes (Tceal1, Depdc1a, and Atf7ip2), all of which are mainly responsible for cell growth and differentiation. While further studies need to be done to verify the roles of these genes by real time-quantitative polymerase chain reaction and/or Western blot, our current findings clearly suggest that upregulation of these pro-growth and pro-survival genes by RBC1023 plays a critical role in its cytoprotective effect.

RBC1023 is a selective, reversible, small-molecule caspase-3 inhibitor, which has multiple advantages compared to the known caspase inhibitors that are being developed. For example, the caspase inhibitors from Pfizer’s IDN6556 series (Pfizer, Inc., New York, NY, USA)32 and Vertex’s 166 series (Vertex Pharmaceuticals, Boston, MA USA)33 are irreversible, peptide-mimic and Pan-caspase inhibitors. Merck’s caspase-3 inhibitor M867 (Merck & Co., Inc., White House Station, NJ, USA) is selective and reversible, but a peptide mimic.34,35 Modified peptides generally do not make good therapeutic reagents, due to problems such as proteolysis and poor cellular permeability;36 therefore, selective small-molecule caspase inhibitors offer significant potential for protection of numerous cell types. Indeed, our preliminary results have shown that RBC1023 protects primary normal human oral keratinocytes from radiation-induced cell death using radiation colony formation assay (data not shown). During head and neck cancer therapy, oral mucositis caused by normal (nonmalignant) cell apoptosis is recognized as a major side effect of the therapy. Therefore, the potential use of cytoprotectants, such as the caspase inhibitors, for prevention of severe mucositis during head and neck cancer therapy could present a unique opportunity for effective therapy. In addition, caspases have also been demonstrated to be involved in some non-apoptotic functions, such as the regulatory events that are important for neural functions including axon pruning and synapse elimination, which are necessary to refine mature neuronal circuits.37 Small-molecule caspase inhibitor drugs may therefore have potential therapeutic or preventive efficacies against neurodegenerative diseases.

Conclusion

Taken together, our experimental study results indicate that RBC1023 protects NIH3T3 cells against staurosporine-induced apoptosis via inhibiting caspase activity, restoring mitochondrial membrane potential, and possibly upregulating some cell survival-related gene expressions.

Table S1

Molecular functions of the top 200 upregulated and downregulated clustering genes in the heat map

Ref seq	Gene symbol	P-value (observed group)	Molecular function
1. NM_010655	Kpna2	5.05E-08	Nuclear protein import
2. NM_173752	1110067D22Rik	6.85E-08	Function unknown
3. NM_025415	Cks2	2.15E-07	Binds to cyclin dependent kinases
4. NM_010655	Kpna2	2.90E-07	Nuclear protein import
5. NM_001081127	Adamts14	5.96E-07	Metallopeptidase
6. NM_007960	Etv1	6.49E-07	Transcriptional activator
7. NM_019752	Htra2	8.23E-07	Caspase activation
8. NM_175149	2310022B05Rik	9.86E-07	Function unknown
9. NM_011158	Prkar2b	1.19E-06	cAMP signaling
10. NM_029074	Tmem188	1.44E-06	Protein phosphatase
11. –		1.52E-06	N/A
12. NM_028607	Bloc1s2	1.57E-06	Biogenesis of lysosome
13. NM_010358	Gstm1	2.17E-06	Glutathione S-transferases
14. NM_172757	Heatr3	2.51E-06	Function unknown
15. NM_008036	Fosb	3.03E-06	Form the transcription factor complex AP-1
16. NM_145692	Ppp1r42	3.03E-06	Regulates phosphatase activity
17. NM_009178	St3gal4	3.24E-06	Protein glycosylation
18. NM_198092	Usp2	3.36E-06	Hydrolase that deubiquitinates polyubiquitinated target proteins
19. ENSMUST00000000001	Gnai3	3.69E-06	Modulator in various transmembrane signaling systems
20. NM_013496	Crabp1	3.79E-06	Involved in retinoic acid-mediated differentiation and proliferation processes
21. –		4.39E-06	N/A
22. NM_022995	Pmepa1	4.56E-06	Involved in down-regulation of the androgen receptor
23. NM_001048192	Agbl5	4.73E-06	Mediates tubulin deglutamylation
24. NM_011654	Tuba1b	4.92E-06	Interaction with PI3KR1
25. NM_029556	Clybl	4.95E-06	Citrate lyase
26. NM_178942	Il17rc	4.99E-06	Cytokine signaling
27. NM_008495	Lgals1	5.13E-06	Regulates cell proliferation
28. NM_009672	Anp32a	5.18E-06	Cell proliferation and differentiation
29. NM_033595	Pcdhga12	5.32E-06	Involved in the establishment and maintenance of specific neuronal connections in the brain
30. NM_007951	Erh	5.36E-06	Involved in the establishment of specific neuronal connections in the brain
31. ENSMUST00000066352	Ptrh1	5.47E-06	tRNA hydrolase
32. NM_010918	Nktr	5.68E-06	Involved in the function of NK cells
33. NM_013726	Dbf4	5.81E-06	Plays a central role in DNA replication and cell proliferation
34. NM_031408	Gigyf1	5.85E-06	Regulates tyrosine kinase receptor signaling
35. NM_011056	Pde4d	5.86E-06	Degrade cAMP
36. ENSMUST00000029316	Exosc8	6.09E-06	Participates in cellular RNA processing and degradation events
37. NM_019511	Ramp3	6.55E-06	Receptor for adrenomedullin
38. NR_033813	Dlk1	6.72E-06	Plays a role in neuroendocrine differentiation
39. –		6.73E-06	N/A
40. NM_012043	Islr	6.93E-06	Function unknown
41. NM_027030	Dcps	7.14E-06	Decapping scavenger enzyme
42. –		7.55E-06	N/A
43. NM_172574	Pqlc3	7.74E-06	Protein coding
44. NM_025415	Cks2	8.20E-06	Involved in cyclin dependent kinase function
45. NM_010171	F3	8.68E-06	Coagulation factor III
46. NM_026268	Dusp6	8.95E-06	Inactivates MAP kinases
47. NM_016762	Matn2	9.13E-06	Formation of filamentous networks
48. NM_144829	Aarsd1	1.01E-05	tRNA synthetase
49. –		1.04E-05	N/A
50. NM_027859	Rnf215	1.05E-05	Ring finger protein
51. NM_009154	Sema5a	1.05E-05	Involved in axonal guidance during neural development
52. NM_001033165	Dnajc25	1.07E-05	Function unknown
53. NM_007476	Arf1	1.19E-05	Vesicular trafficking
54. NM_029701	Spcs3	1.21E-05	Removes signal peptides from nascent proteins
55. AK139984	Gm10616	1.21E-05	N/A
56. NM_026815	Upk1a	1.22E-05	Mediates signal transduction events
57. NM_172563	Hlf	1.22E-05	Function unknown
58. NM_009193	Slbp	1.28E-05	Involved in the histone pre-mRNA processing
59. NM_007960	Etv1	1.28E-05	Transcriptional activator
60. NM_001081636	Ccnd3	1.29E-05	Cell cycle regulation
61. NM_028770	Krt80	1.30E-05	Major structural fibers of cells
62. NM_016750	H2afz	1.31E-05	Nucleosome structure protein
63. –		1.37E-05	N/A
64. NM_011366	Sorbs3	1.39E-05	Activation of JNK
65. NR_028560	Scarna17	1.41E-05	Function unknown
66. NR_028560	Scarna17	1.41E-05	Function unknown
67. NM_008008	Fgf7	1.42E-05	Regulation of embryonic development, cell proliferation and cell differentiation
68. NM_027911	Raver1	1.44E-05	Modulate regulated alternative splicing events
69. NM_019840	Pde4b	1.44E-05	Hydrolyzes cAMP
70. NM_011873	Dazap2	1.46E-05	Cell signaling and transcription regulation
71. NM_025982	Tspan31	1.46E-05	Regulation of cell development, activation, growth and motility
72. NM_010655	Kpna2	1.50E-05	Nuclear protein import
73. NM_008891	Pnn	1.50E-05	Transcriptional activator
74. NM_008193	Guk1	1.51E-05	Essential for recycling GMP
75. NM_198658	Nr1h5	1.54E-05	Function unknown
76. NM_001242364	Fam160a2	1.56E-05	Promotes vesicle trafficking
77. NM_021605	Nek7	1.60E-05	Mitotic cell cycle progression
78. NM_009387	Tk1	1.62E-05	Thymidine kinase
79. NM_025281	Lyar	1.64E-05	Function unknown
80. NM_172458	Zfp871	1.64E-05	Function unknown
81. NM_177390	Myo1d	1.68E-05	Involved in intracellular movements
82. NM_153195	Fbxo7	1.71E-05	A negative regulator of NF-kappaB signaling
83. NM_144549	Trib1	1.77E-05	Interacts with MAPK kinases and regulates activation of MAP kinases
84. NM_023117	Cdc25b	1.80E-05	Dephosphorylates CDK1 and stimulates its kinase activity
85. NR_028480	Gm12238	1.84E-05	Function unknown
86. NM_011664	Ubb	1.86E-05	Ubiquitination and protein degradation
87. NM_008037	Fosl2	1.88E-05	Forming the transcription factor complex AP-1
88. –		1.88E-05	N/A
89. NM_022312	Tnr	1.88E-05	Mediates inhibition of cell adhesion and neurite outgrowth
90. NM_011693	Vcam1	1.89E-05	Mediates leukocyte-endothelial cell adhesion and signal transduction
91. NM_145706	Nup43	1.89E-05	Required for microtubule attachment, mitotic progression and chromosome segregation
92. NM_012003	Cops7a	1.90E-05	Involved in various cellular and developmental processes
93. NM_145823	Pitpnc1	1.91E-05	Involved in cell signaling and lipid metabolism
94. NM_008801	Pde6d	1.94E-05	Involved in cellular signal transduction
95. NR_028528	Snord57	1.95E-05	Small nucleolar RNA
96. NM_008086	Gas1	1.96E-05	Involved in growth suppression
97. NR_004445	Snord22	2.01E-05	Small nucleolar RNA
98. NM_001113486		2.01E-05	N/A
99. NM_026667	Fam114a1	2.06E-05	Plays a role in neuronal cell development
100. NM_001166065	Gcnt4	2.09E-05	Glycosyltransferase
101. NM_025415	Cks2	2.11E-05	Binds to cyclin dependent kinases
102. NM_001081229	Tsc22d2	2.12E-05	Function unknown
103. NM_011923	Angptl2	2.12E-05	Induces sprouting in endothelial cells
104. NM_011664	Ubb	2.15E-05	Ubiquitination and protein degradation
105. ENSMUST00000105363	Gamt	2.16E-05	Methyltransferase
106. NM_011654	Tuba1b	2.16E-05	Interaction with PI3KR1
107. NM_001159288	Morc2a	2.18E-05	Transcriptional repressor
108. NM_010560	Il6st	2.30E-05	Signal-transducing molecule
109. –		2.31E-05	N/A
110. NM_178403	Pus7	2.32E-05	Catalyzes pseudouridylation
111. NM_001030307	Dkc1	2.34E-05	Required for ribosome biogenesis and telomere maintenance
112. NM_008197	H1f0	2.35E-05	Required for the condensation of nucleosome chains into higher-order structures
113. NM_028156	1700001P01Rik	2.40E-05	Function unknown
114. NM_008871	Serpine1	2.45E-05	Serine protease inhibitor
115. L08651	Rpl29	2.46E-05	Ribosomal protein
116. –		2.48E-05	N/A
117. NM_001198985	Tjp2	2.50E-05	Plays a role in tight junctions and adherens junctions
118. NM_027349	Rbm25	2.55E-05	A regulator of alternative pre-mRNA splicing
119. NM_001110310	Snx12	2.66E-05	Involved in several stages of intracellular trafficking
120. NM_178615	Rgmb	2.79E-05	Promotes neuronal adhesion
121. NM_010234	Fos	2.79E-05	Forms the transcription factor complex AP-1
122. NM_011812	Fbln5	2.82E-05	Plays a role in vascular development and remodeling
123. –		2.83E-05	N/A
124. NM_001013762	Gm4952	2.86E-05	Function unknown
125. NM_001190466	Dact1	2.87E-05	Involved in regulation of intracellular signaling pathways
126. NM_029999	Lbh	2.94E-05	Involved in MAP kinase signaling pathway
127. NM_010026	Asap1	2.94E-05	Involved in the differentiation of fibroblasts into adipocytes
128. –		2.98E-05	N/A
129. ENSMUST00000012348	Gstm2	2.99E-05	Glutathione S-transferase
130. NM_001080129	Tmpo	2.99E-05	Helps maintain the structural organization of the nuclear envelope
131. NM_011898	Spry4	3.06E-05	Suppresses the insulin receptor and EGFR-transduced MAPK signaling pathway
132. NM_026115	Hat1	3.17E-05	Plays a role in DNA repair
133. NM_178591	Nrg1	3.18E-05	Plays critical roles in the growth and development of multiple organ systems
134. NR_046144	Rn5s20	3.19E-05	Required for transcription
135. NM_145354	Nsun2	3.20E-05	RNA methyltransferase
136. NM_001111059	Cd34	3.28E-05	Plays a role in the attachment of stem cells to stromal cells
137. NM_001204233	Spp1	3.29E-05	Involved in enhancing production of interferon-gamma
138. NM_026864	Rasl11a	3.30E-05	Ras signaling
139. NM_133714	Kansl2	3.32E-05	Involved in acetylation of nucleosomal histone H4
140. NM_021389	Sh3kbp1	3.36E-05	Involved in apoptosis, cytoskeletal rearrangement, cell adhesion
141. –		3.37E-05	N/A
142. NM_009609	Actg1	3.38E-05	Involved in various types of cell motility
143. NM_011859	Osr1	3.38E-05	Transcription factor
144. NM_010256	Gart	3.43E-05	Required for De novo purine biosynthesis
145. NM_011435	Sod3	3.44E-05	Antioxidant enzymes
146. –		3.46E-05	N/A
147. –		3.47E-05	N/A
148. NM_178604	Txnl4a	3.48E-05	Essential role in pre-mRNA splicing
149. BC089624	Rpl7a	3.53E-05	Catalyzes protein synthesis
150. NM_008047	Fstl1	3.54E-05	Modulates the action of some growth factors
151. NM_207161	BC048355	3.58E-05	Function unknown
152. NM_009626	Adh7	3.70E-05	Required for synthesis of retinoic acid
153. NM_172301	Ccnb1	3.71E-05	A regulatory protein involved in mitosis
154. NM_013543	H2-Ke6	3.73E-05	Plays a role in biosynthesis of fatty acids in mitochondria
155. NM_001081277	Ak5	3.74E-05	Involved in regulating the adenine nucleotide composition
156. NM_181390	Mustn1	3.79E-05	Involved in the development and regeneration of the musculoskeletal system
157. NM_009498	Vamp3	3.84E-05	Involved in vesicular transport from the late endosomes to the trans-Golgi network
158. NM_027722	Nudt4	3.86E-05	Plays a role in signal transduction
159. NM_022324	Sdf2l1	3.87E-05	Function unknown
160. NM_007621	Cbr2	3.88E-05	Plays a role in the uronate cycle of glucose metabolism
161. NM_026369	Arpc5	4.07E-05	Involved in regulation of actin polymerization
162. AK079043	Gm9908	4.09E-05	Function unknown
163. NM_011415	Snai2	4.10E-05	Transcriptional repressor
164. NM_008027	Flot1	4.13E-05	Involved in vesicular trafficking and signal transduction
165. NM_011452	Serpinb9b	4.14E-05	Granzyme B inhibitor
166. NM_001100454	Wfikkn1	4.18E-05	Protease inhibitor
167. NM_016750	H2afz	4.22E-05	Plays a central role in transcription regulation, DNA repair, DNA replication and chromosomal stability
168. NM_028626	Mcee	4.29E-05	Function unknown
169. NM_010729	Loxl1	4.32E-05	Plays roles in developmental regulation, senescence, tumor suppression
170. NR_034052	Snora2b	4.33E-05	Small nucleolar RNA
171. NM_146151	Tesk2	4.35E-05	Plays an important role in spermatogenesis
172. NM_019571	Tspan5	4.36E-05	Mediates signal transduction events
173. –		4.44E-05	N/A
174. NM_011297	Rps24	4.48E-05	Required for processing of pre-rRNA
175. NR_004432	Rnu12	4.49E-05	RNA gene
176. NM_007436	Aldh3a1	4.50E-05	Involved in the metabolism of corticosteroids, biogenic amines, neurotransmitters, and lipid peroxidation
177. NM_009504	Vdr	4.50E-05	Nuclear hormone receptor
178. NM_029645	Gatc	4.59E-05	Allows the formation of correctly charged Gln-tRNA in the mitochondria
179. NM_145360	Idi1	4.59E-05	Required for the synthesis of farnesyl diphosphate and cholesterol
180. NM_053273	Ttyh2	4.59E-05	Involved in cell proliferation and cell aggregation
181. NM_173756	Lin52	4.61E-05	Represses cell cycle-dependent genes
182. ENSMUST00000064795	Egr1	4.61E-05	Activates the transcription of target genes whose products are required for mitogenesis and differentiation
183. NM_010308	Gnao1	4.64E-05	Transducers in transmembrane signaling systems
184. NM_008304	Sdc2	4.66E-05	Mediates cell binding, cell signaling, and cytoskeletal organization
185. NM_029809	2310014L17Rik	4.66E-05	Function unknown
186. NM_009830	Ccne2	4.68E-05	Essential for the control of the cell cycle at the late G1 and early S phase
187. ENSMUST00000067692	Rad51c	4.71E-05	Essential for the homologous recombination pathway of DNA repair
188. NM_010219	Fkbp4	4.76E-05	Plays a role in the intracellular trafficking
189. NM_172590	Wdr41	4.78E-05	Function unknown
190. NM_177370	Rhbdd3	4.81E-05	Function unknown
191. ENSMUST00000027009	Casp12	4.83E-05	Reduces activation of NF-kappa-B in response to TNF
192. NM_021504	Ngly1	4.90E-05	Plays a role in the proteasome-mediated degradation of misfolded glycoproteins
193. NM_146206	Tpcn2	4.98E-05	Involved in smooth muscle contraction
194. NM_172784	Lrp11	5.01E-05	Function unknown
195. NM_018773	Skap2	5.02E-05	Plays a role in src signaling pathway
196. NM_028234	Rbm33	5.06E-05	Function unknown
197. ENSMUST00000097772	Il1rl1	5.18E-05	Involved in helper T cell function
198. NM_023871	Set	5.24E-05	Involved in apoptosis, transcription, nucleosome assembly and histone chaperoning
199. NM_007413	Adora2b	5.29E-05	Involved in axon elongation
200. NM_007713	Clk3	5.30E-05	Dual specificity kinase acting on both serine/threonine and tyrosine-containing substrates

Abbreviation: N/A, not available.