1. Introduction
Consistently escalating morbidity and complications associated with diabetes have become a matter of growing concern. As of 2021, globally, 537 million (1 in 10) adults were living with diabetes, and these staggering figures are predicted to even rise up to 783 million by the 2045 [
1]. Although in the past century, the discovery of insulin has improved the therapeutic options for patients with diabetes, certain limiting factors such as insufficient efficacy, cost, weight gain, and life-threatening hypoglycemic episodes remain the prohibitive challenges that are commonly associated with use of insulin in humans [
2,
3].
Type 2 diabetes mellitus (T2DM) is a chronic metabolic disorder, accounting for 90–95% of all diabetic cases. Its pathophysiology involves a combined effect of long-term insulin resistance and β-cell impairment and consequently the reduced insulin secretion and hyperglycemia, leading ultimately to pancreatic β-cell dysfunction [
4]. Type 1 diabetes (T1D) is an autoimmune and inflammatory disease which is characterized by a loss of pancreatic β-cells. Its treatment is basically dependent on the exogenous insulin to maintain normoglycemia which requires continuous infusion or multiple daily injections, with no oral medications available. Moreover, treatment with exogenous insulin does not reverse the loss of β-cell mass and function [
5]. Although the therapeutic approaches aimed at improving insulin delivery and glucose uptake are necessary for the treatment of diabetes, they are not designed to target the disease at its root cause i.e., β-cell dysfunction or loss.
Verapamil is a calcium channel blocker (CCB) which was approved by the Food and Drug Administration (FDA) in the late 1990s for the treatment of hypertension, angina, and arrhythmia [
6]. Notably, it was found to lower the risk of developing new-onset diabetes [
7,
8,
9], as compared to other CCBs; and was also found to lower fasting glucose levels in patients with diabetes [
3,
10,
11,
12,
13,
14]. Verapamil is known to act by inhibiting the L-type calcium channels and reducing the cellular calcium influx, thereby inhibiting the transcription of thioredoxin-interacting protein (TXNIP) [
14,
15,
16] which would otherwise accumulate and increase oxidative stress related toxicity in the pancreas [
11,
17,
18]. While, the intracellular calcium is required to trigger insulin exocytosis, chronic increase in the calcium levels may result in β-cell impairment [
19]. This preventative mechanism of verapamil is highly attributed to TXNIP downregulation which leads to a decrease in β-cell apoptosis [
14,
15,
19].
Verapamil ameliorated β-cell survival and activity, enhanced insulin levels, improved glucose sensitivity and homeostasis, and reduced hypoglycemic episodes in both mouse model and T1D patients [
3,
11,
13,
14,
15]. Oral verapamil administered to patients with T1D lowered dependency on exogenous insulin and improved endogenous β-cell function, measured using the C-peptide area under the curve (AUC), as compared to placebo. Additionally, it decreased the frequency of hypoglycemic events [
3,
11,
13,
14,
15]. Similarly, in mice, oral verapamil lowered TXNIP expression, increased endogenous insulin levels, and reduced β-cell apoptosis, rescuing the mice from streptozotocin (STZ)-induced diabetes [
15,
20]. In human clinical trials, verapamil showed a reduced risk of new-onset diabetes, as compared to atenolol use [
8,
9].
These interesting studies led us to further investigate the mechanisms of the beneficial effects of verapamil on MIN6 mouse pancreatic β-cells and transgenic zebrafish models, specifically focusing on β-cells viability and function, following challenge with immune-metabolic stressors mimicking the T1D/T2D insults.
2. Materials and Methods
Cell Culture and Treatments
MIN6 cells, a mouse β-cell line, were generously gifted by Dr. Jun-ichi Miyazaki, Kumamoto University Medical School, Japan [
21]. The cells were cultured and maintained in Dulbecco’s modified Eagle’s medium (DMEM) which contained 25.2 mM glucose, supplemented with 15% heat-inactivated fetal bovine serum (FBS), 2 mM l-glutamine, 1 mM sodium pyruvate, 10 mM HEPES, 50 μg/mL penicillin, 100 μg/mL streptomycin, and 70 mM β-mercaptoethanol at 37°C, with 5% CO
2 (all reagents were purchased from Invitrogen, Waltham, MA, USA). The cells were re-fed every 2 to 3 days [
21]. For experimental assays, cells were maintained in the supplemented DMEM containing 5.6 mM glucose. All analyses using MIN6 cells were performed at a passage of <25.
A stock solution of verapamil hydrochloride (152-11-4, Calbiochem, Burlington, MA, USA) was freshly prepared in phosphate-buffered saline (PBS). The cells were treated with verapamil at different concentrations, depending on each experimental setting.
Cell Survival
Cells were seeded at a concentration of 5×104 cells/well in a 96-well plate (Costar, High Wycombe, UK). Based on the experimental set-up, cells were treated with verapamil and/or different stressors (Details included in Supplementary Methods), and cell survival was determined using an MTT (3-[4,5-dimethylthiazol-2-yl]-2,5 diphenyl tetrazolium bromide) assay (Trevigen, Minneapolis, MN, USA), following the manufacturer’s instructions. Plates were read using Synergy H4 Hybrid Microplate Reader (BioTek, Winooski, VT, USA), and the data were analyzed using Gen5 software.
Cell Proliferation
MIN6 β-cell proliferation was determined by monitoring growth, consecutively for 24h. To perform proliferation assays, cells were seeded at concentration of 5×104 cells/well in 24-well plates (Costar, High Wycombe, UK) and incubated overnight at 37°C in 5% CO2. Later, verapamil was added (50 μM), every 24h for 7 days, followed by trypsinization, harvesting, and cell counting using trypan blue dye exclusion assay.
Immunohistochemical Detection of Ki67
To test the effect of verapamil on MIN6 β-cell proliferation, immunohistochemical (IHC) analysis of Ki67 was performed, as described previously [
22]. Briefly, cells were seeded on coverslips in media supplemented with 2% bovine serum albumin (BSA) free fatty acid (FFA) (Cat# 820024, Sigma, USA) for 24h. Then, cells were either treated with verapamil (50 μM) for another 24h or left untreated to serve as control. Cells were stained using rabbit anti-mouse Ki67 primary antibody and goat anti-rabbit secondary antibody; actin filaments and nuclei were counterstained using Phalloidin-iFluor 594 reagent and mounting medium with DAPI (Details included in Supplementary Methods) and confocal images were obtained using inverted Zeiss LSM710 spectral confocal microscope (Carl Zeiss, Gottingen, Germany). All samples were analyzed using the same parameters and the resulting color markup for analysis of each sample was confirmed. Correlated total cell fluorescence (CTCF) was calculated from 10 different fields of each number (n) via ImageJ, using the following equation:
CTCF = Integrated Density - (Area of selected cell × Mean fluorescence of background readings) [
23].
Western Blotting
MIN6 cells were seeded in 6-well plates (Corning, Somerville, MA, USA), and incubated in media supplemented with 2% BSA free-FFA for 24h. Then, cells were either treated with verapamil (50 μM) and incubated for another 24h or left untreated (control). Later, total protein was extracted, and Western blotting was performed as previously described [
24,
25] (Details included in Supplementary Methods).
Cell Cycle Analysis
Cell cycle analysis was carried out using flow cytometry, based on propidium iodide (PI) staining of the DNA, as described before [
26]. Briefly, MIN6 cells were seeded in 12-well plates (Costar, High Wycombe, UK) in media supplemented with FFA-free 2% BSA for 24h. Then, cells were treated with verapamil (50 μM) for 1, 2, 3, or 4h, or left untreated (control). Following trypsinization, re-suspended cells were fixed with ice-cold 70% ethanol, harvested by centrifugation (800 ×g for 5 min), treated at room temperature (RT) with RNase (100 µg/mL) for 30 min, followed by incubation with PI (40 µg/mL) for another 30 min. DNA content and cell cycle were analyzed using BD FACSCant Flow Cytometer (BD Bioscience, San Jose, Ca, USA). Cells were gated and analyzed for doublet exclusion. A single cell population was acquired for analysis at 490 nm Ex and 630 nm Em wavelengths. Data were analyzed using BD FACSDivaTM Software 8 (BD Biosciences, USA).
Protection Assays
For protection assays, verapamil was used at concentrations of 1, 5, 10, and 50 μM. The protection assays were conducted for verapamil treatment in pretreatment, cotreatment, and pretreatment followed by cotreatment conditions, and challenged the cells with STZ and T1D-/T2D-cytomixes, as indicated elsewhere [
27,
28]. (Details included in Supplementary Methods).
Glucose-Stimulated Insulin Secretion (GSIS) Assay
A GSIS assay was performed using static incubation method, as described before [
29] (Details included in Supplementary Methods). Briefly, 1×10
5 cells/well were seeded in 24-well plates (Corning, Somerville, MA, USA) and incubated for 24h. Later, MIN6 cells were treated with verapamil (50 μM) for another 24h. Then, the cells were incubated with Krebs –Ringer HEPES buffer supplemented with different concentrations of glucose at 37°C. GSIS assay was terminated, the conditioned media were collected, and levels of secreted insulin were detected using ultrasensitive mouse insulin ELISA kit (Mercordia, Sylveniusgatan, Uppsala, Sweden). Total insulin content in cells were measured by adding cold acidified ethanol to each well and keeping it overnight at -80°C. After three freeze/thaw cycles, cells were scraped off and supernatants were collected by centrifugation, neutralized with 1M Tris pH 7.5 and insulin levels were measured using ELISA. The levels of total insulin content and secreted insulin were normalized by the level of total protein measured by Bradford assay against a standard curve.
Whole-Transcriptome Analysis
Whole-transcriptome analysis was conducted as described [
25] (Details included in Supplementary Methods). Briefly, total RNA was isolated from verapamil-treated (50 μM; 24h) and untreated control MIN6 cells using RNeasy kit (Qiagen, Hilden, Germany), following standard protocol. For whole-transcriptome sequencing, RNA (40 ng) was used to prepare the transcriptome libraries, using Truseq stranded mRNA kit (Illumina Inc., San Diego, CA, USA) and following the manufacturer’s protocol. The libraries were validated and quantified using bioanalyzer (Agilent Technologies, Santa Clara, CA, USA) and qubit fluorometer (Thermofisher Scientific, Santa Clara, CA, USA), respectively. Paired-end sequencing was carried out using Novaseq 6000 system (Illumina Inc. San Diego, CA, USA) and the BCL files were converted to Fastq using bcl2fastq v.2.20 tool. Quality control of Fastq files was performed using FastQC (v0.11.9). Trimmomatic (v0.39) was used for adaptor and quality trimming, and for removing extremely short reads and HISAT2 (v2.1.0) for data alignment. To enumerate the number of reads associated with the genes, htseq-count tool was used in HTSeq (v 0.9.1). Differential gene expression analysis was performed using Bioconductor package edgeR, using the default setting.
Proteome Sample Preparation
Proteome sample preparation was conducted based on the standard protocol [
31] (Details included in Supplementary Methods). Briefly, frozen pellets of MIN6 cells treated with verapamil (50 μM; 24h) and untreated cells (control) were lysed, reduced using reduction buffer at 37°C, and then digested by incubating with trypsin at 37°C overnight. Next, the samples were acidified, cleaned by using C-18 MACROSpin plates, and the digested peptides were eluted in 340 µl/well of C18 elution solution (0.1% v/v) TFA in 50% acetonitrile. Samples were dried, reconstituted (1 µg/µl) and analyzed using Q-Exactive liquid chromatography-electrospray ionization-tandem mass spectrometry (LC-ESI-MS/MS) coupled with an EASY-nLC™ 1200 nano-LC System, through an EASY-Spray Ion Source (Thermo Fisher Scientific, Waltham, MA, USA). Peptides were eluted, separated using a 170-minute run, and the eluent was ionized using Easy Spray nano ESI source operating in positive ion mode. Label-free quantitation was performed using Thermo Scientific Proteome Discoverer 2.4 software and SEQUEST
® HT search engine. The data were searched against a Mus musculus fastafile with a 1% false discovery rate (FDR) using Percolator. Further processing was performed using new Rt-Aligner and created Feature Mapper nodes for the untargeted label-free quantification workflow. FDRs for proteins, peptides, and peptide spectral matches were kept at 1%. All results were filtered by a q-value of <0.01 (equals an FDR of 1% on peptide level and a filter of minimum 2 unique peptides) Each experiment was performed in technical triplicates.
In vivo Study Using Transgenic Zebrafish Embryos Model
Zebrafish (
Danio rerio) were reared at 28°C under a 12-h light/dark cycle at the Dasman Diabetes Institute’s Animal Facility. The fish were bred as per the guidelines approved by the Animal Care Ethical Committee (RA-2019-005), in accordance with the National Institutes of Health (NIH) guide for the care and use of laboratory animals (NIH Publications No. 8023, revised 1978) and in compliance with the ARRIVE guidelines as well as the standard laboratory procedures for zebrafish [
33]. Ins:NfsB-mCherry transgenic zebrafish was a generous gift from Dr. Michael J. Parsons, University of California, Irvine, USA. In this ins:NfsB-mCherry transgenic model, the bacterial gene encoding nitro-reductase (NTR) enzyme (which converts a prodrug such as metronidazole (MTZ) to cytotoxins) was fused to mCherry fluorescent reporter, under the control of insulin proximal promoter [
34]. MTZ (M3761, Sigma, St. Louis, Mo, USA) was dissolved in E3 media as described [
35] and verapamil was dissolved in deionized water. Embryos were reared in E3 media at 28°C in dark and to study the β-cell protective effect of verapamil, larvae were monitored 3-6 days post fertilization (dpf) β-cell regeneration by detecting mCherry fluorescence intensity. The transgenic larvae were grouped as follows: group 1, control (untreated larvae); group 2, larvae treated with 10 µM verapamil; group 3, larvae treated with 10 µM verapamil for 24h, and then exposed to fresh E3 media containing 10 mM MTZ for 48h; and group 4, embryos treated with 10 mM MTZ at 4 dpf for 48h. Full body images (75× magnification; 200 ms exposure time) were acquired using Discovery V12 Stereo microscope (Zeiss, Jena, Germany) and Alexa Flour 594 red fluorescent emission filter.
Statistical Analysis
The data obtained were expressed as Mean ± SEM values from at least three individual experiments. The significance of group differences was calculated using unpaired Student’s t-test or one-/two-way Analysis of Variance (ANOVA) test, followed by Tukey’s multiple comparisons test using GraphPad Prism 9.0 (GraphPad Software, CA, USA). P-values <0.05 were considered statistically significant, and expressed as *P˂0.05, **P˂0.01, ***P˂0.001, or ****P˂0.0001.
4. Discussion
Insulin-producing pancreatic β-cells are the key players in glucose homeostasis and their loss is associated with the development of T1D and the progression of T2D [
4,
5]. Not surprisingly, residual β-cell function has been associated with reduced incidence of microvascular complications in both T1D and T2D [
41,
42]. Currently, promising antidiabetic therapeutic strategies aim at rescuing or preserving the β-cell function and insulin secretory capacity. Herein, we investigated the β-cell cytoprotective and insulinotropic effects of verapamil using T1D- and/or T2D-mimicking stressor challenges in MIN6 β-cells and transgenic zebrafish models. In our MIN6 T1D/T2D stressor models, verapamil treatment improved the β-cell viability/proliferation in a dose- and time-dependent manner, as supported by increased expression of Ki67/H3 proliferation markers and evidence of actively replicating β-cells in the G2/M phase. Moreover, in verapamil-treated β-cells, an increase in the insulin secretion was observed under high glucose culture conditions. Indeed, the normalized secreted insulin in verapamil-treated cells was significantly promoted only under the condition of high glucose (16.8 mM) stimulation, compared to the lower glucose concentrations tested in our model. It implies that β cells get hypersensitized by abnormally high glucose levels, a notion which has previously been associated with enhanced blood glucose control and reduction in exogenous insulin requirements [
11,
13]. Together, these observations lead to the conjecture that verapamil may rather have a more promising insulinotropic benefit under advanced hyperglycemic conditions.
Notably, while most other studies evaluated verapamil use in diabetic cohorts, we assessed its direct β-cell protective and insulinotropic effects using modified, more clinically relevant models in setting of verapamil pre-treatments or co-treatments, followed by diverse cytotoxic challenges including T1D-/T2D-cytomixes. The β-cell protective effect was further confirmed at the organismal level using MTZ-mediated β-toxicity in mCherry-reporter zebrafish embryos. We found a significantly enhanced β-cell viability in presence of verapamil in both
in vitro and
in vivo, establishing its preventative and β-protective benefits, in line with previous evidence that verapamil lowers the risk of new-onset T2D [
3,
7,
8,
11,
13,
14].
Our analyses of changes in the proteomic and transcriptomic landscapes in verapamil-treated MIN6 β-cells are of particular interest. It is noteworthy that the treatment with verapamil brought to light a wide array of proteins and genes that were differentially expressed in MIN6 β-cells. Of those, some are remarkable, owing to their involvement in fundamental pathways ultimately contributing toward maintaining glucose homeostasis and/or preserving pancreatic β-cells. In our transcriptomic analysis, we found
Ccnd1 upregulated at the mRNA level, and the role of this gene in inducing β-cell proliferation is well established [
43,
44]. In addition, through transcriptome analysis we also detected the lower expression of
Txnip gene in verapamil-treated, compared to untreated cells. TXNIP, a redox regulator, was identified as an attractive therapeutic target for restoring β-cells dysfunction, owing to its increased expression in the β-cells of diabetic patients [
17,
20]. TXNIP overexpression is known to play a key role in β-cell apoptosis, whereas its deficiency promotes the endogenous β-cell survival [
17,
20,
45]. Thioredoxin, a protein involved in the antioxidant defense system of β-cells, is inhibited by TXNIP, thereby promoting oxidative stress. Moreover, TXNIP is also associated with inflammasome activity, and its downregulation may have anti-inflammatory effects [
46].
In response to verapamil, the most dominantly upregulated factor, at both mRNA and protein levels, was the glucose-like peptide-1 (GLP-1)-induced incretin — called cholecystokinin (CCK). To the best of our knowledge, we are the first to show that verapamil induces expression of
CCK transcripts and proteins in MIN6 β-cells. CCK plays a role in preserving β-cell mass with increasing age and protecting them against apoptosis induced by STZ and β-toxic proinflammatory cytokines [
47]. In our
in vitro studies, verapamil induced β-cell protection against insults including STZ and T1D-/T2D-cytomixes, which could be due, in part, to upregulated CCK expression.
We also found upregulated levels of
Wnt4, which is a ligand heterogeneously expressed in β-cells and it plays a central role in activating calcium signaling and lowering blood glucose levels, along with regulating cellular redox homeostasis [
48]. Wnt4 is the most abundantly expressed protein in β-cells in high fat, high glucose, or insulin-resistant conditions; and it is suggested to play a pivotal role in β-cell proliferation or regulating and enhancing the insulin secretory response induced by the above-mentioned conditions [
49]. Kurita et al. suggested that
Wnt4 might possibly be involved in glucose-induced vesicle transport of insulin in MIN6 β-cells [
49]. Moreover,
Wnt4 overexpression was shown to upregulate genes controlling β-cell maturation and function and related with higher mitochondrial abundance [
48].
An increase in both mRNA and protein levels of CaMK4 upon verapamil treatment in our study is suggestive of its impact on MIN6 β-cell survival and proliferation. CaMK4, a serine/threonine kinase, has been previously shown to mediate
Irs2 expression in MIN6 cells and mouse islets stimulated by glucose,
via the activation of cAMP response element binding protein (CREB) [
50]. The present study, based on previous knowledge that
Irs2 knockout reduces β-cell mass and increases β-cell apoptosis, suggests that CaMK4 may regulate β-cell survival and proliferation, possibly
via the IRS2-dependent stimulation of cell proliferation and apoptosis inhibition.
In this study we found higher level of basal and maximal mitochondrial respiration in verapamil-treated MIN6 β-cells which could be attributed to higher
Wnt4 expression, and that, along with reduced TXNIP gene expression could possibly explain the ability of verapamil to render protection against stressors [
15]. Indeed, several studies showed that verapamil treatment reduced the pro-apoptotic TXNIP expression, thereby promoting the survival of β-cells and preventing diabetes [
15,
20]. Altogether, Wnt4 upregulation by verapamil treatment may have potential as a therapeutic target regarding T2D and, therefore, requires further in-depth studies.
Mitochondrial dysfunction is related to pathophysiology of diabetes, resulting in impaired ATP production and defective insulin secretion. VDAC, a key protein found in the outer mitochondrial membrane, regulates the movement of metabolites across the membrane [
51], and is a key element in apoptotic signaling [
52]. VDAC1 overexpression in various cells has been shown to induce apoptotic cell death [
52], and it increased the susceptibility to apoptosis in response to high glucose conditions [
53], suggesting its regulatory role as a gatekeeper in mitochondria-mediated apoptosis.
VDAC1 overexpression was found to be driven by impaired glycemic levels and its inhibition restored GSIS and prevented hyperglycemia in
db/db mice [
54]. Of note, we found that verapamil treatment downregulated protein expression of both VDAC isoforms (VDAC1/2); and considering that VDAC1 was upregulated by a state of glucotoxicity over a prolonged period in pre-diabetic patients [
53,
54,
55], verapamil could be an ideal therapy to prevent both the development and progression of diabetes. Taken together, based on growing numbers of studies, we suggest that the anti-diabetic benefit of verapamil could be attributed to: (i) decreased TXNIP mRNA expression [
15], (ii) CCK upregulation [
47], (iii) Wnt4 upregulation [
49], (iv) increased mRNA and protein levels of CaMK4 [
56], and/or (v) downregulation of VDAC 1/2 [
52].
Although our study provides valuable insights into the molecular and cellular mechanisms that are altered by verapamil, it involves certain limitations. Firstly, we have used a single cell-line (MIN6 cells) and.it would be interesting to also study effects of verapamil on other β cell lines and primary isolated islets. Next, the role of genes such as Wnt4 must be inferred with caution, owing to their heterogeneity and biphasic involvement in both canonical and non-canonical pathways, and mainly due to a lack of deeper understanding of their precise role in diabetes and related disorders. Further studies are required to elucidate the exact role of such factors.
Author Contributions
Conceptualization, Hossein Arefanian, Ashraf Al Madhoun, Fatema Al-Rashed, Fawaz Alzaid, Fatemah Bahman and Fahd Al-Mulla; Funding acquisition, Hossein Arefanian, Ashraf Al Madhoun, Rasheed Ahmad and Fahd Al-Mulla; Investigation, Hossein Arefanian, Ashraf Al Madhoun, Fatema Al-Rashed, Fawaz Alzaid, Fatemah Bahman, Rasheeba Nizam, Sardar Sindhu, Mohamed Abu-Farha, Rasheed Ahmad and Fahd Al-Mulla; Methodology, Hossein Arefanian, Ashraf Al Madhoun, Fatema Al-Rashed, Fawaz Alzaid, Fatemah Bahman, Rasheeba Nizam, Mohammed Alhusayan, Sumi John, Sindhu Jacob, Michayla Williams, Nermeen Abu Khalaf, Steve Shenouda, Shibu Joseph, Halemah AlSaeed, Shihab Kochumon, Mohamed Abu-Farha and Fahd Al-Mulla; Supervision, Hossein Arefanian, Ashraf Al Madhoun, Fatema Al-Rashed, Fawaz Alzaid, Mohamed Abu-Farha, Jehad Abubaker, Thangavel Alphonse Thanaraj, Rasheed Ahmad and Fahd Al-Mulla; Visualization, Hossein Arefanian, Ashraf Al Madhoun, Fatema Al-Rashed, Rasheeba Nizam, Mohammed Alhusayan, Sumi John, Michayla Williams, Nermeen Abu Khalaf, Steve Shenouda and Sardar Sindhu; Writing – original draft, Hossein Arefanian, Ashraf Al Madhoun, Fatema Al-Rashed, Fatemah Bahman, Rasheeba Nizam, Lubaina Koti, Sardar Sindhu, Mohamed Abu-Farha and Fahd Al-Mulla; Writing – review & editing, Hossein Arefanian, Ashraf Al Madhoun, Fatema Al-Rashed, Fawaz Alzaid, Anwar Mohammad, Lubaina Koti, Sardar Sindhu, Mohamed Abu-Farha, Jehad Abubaker, Thangavel Alphonse Thanaraj, Rasheed Ahmad and Fahd Al-Mulla.
Figure 1.
Proliferative effect of verapamil in MIN6 cells. (A) MIN6 cells were maintained in 5.6 mM glucose and were treated with different concentrations (1, 5, 10, and 50 μM) of verapamil for 24 hours. (B) Growth curve of MIN6 cells cultured in 5.6 mM glucose media and treated with verapamil (50 μM, Red line), in comparison with untreated cells (Blue line). (C,D) Ki67 analysis showing the influence of verapamil on cell proliferation. MIN6 cells were cultured in serum free conditions in the presence or absence of verapamil. Cells cytoskeletons were labeled with phalloidin (red) and the proliferation level was monitored with the expression of KI67 protein (green). DAPI was used to counterstain the nucleus. Scale bar = 50 μM, n= 2. The corrected total cell fluorescence (CTCF) was calculated from 10 different fields of each (n) number. (E,F) Western blot analysis of histone H3 expression corrected to β-actin with representative immunoblots, n=2. Flow cytometry analysis was conducted to investigate the influence of verapamil induction in a time-dependent manner. (G) Histograms of the cell cycle phases. n=4. Data are presented as mean ± SEM values and were analyzed using one-way ANOVA with Tukey’s multiple comparisons test. ns = non-significant, *p≤0.05, **p≤0.01, ***p≤0.001, ****p≤0.0001 versus control.
Figure 1.
Proliferative effect of verapamil in MIN6 cells. (A) MIN6 cells were maintained in 5.6 mM glucose and were treated with different concentrations (1, 5, 10, and 50 μM) of verapamil for 24 hours. (B) Growth curve of MIN6 cells cultured in 5.6 mM glucose media and treated with verapamil (50 μM, Red line), in comparison with untreated cells (Blue line). (C,D) Ki67 analysis showing the influence of verapamil on cell proliferation. MIN6 cells were cultured in serum free conditions in the presence or absence of verapamil. Cells cytoskeletons were labeled with phalloidin (red) and the proliferation level was monitored with the expression of KI67 protein (green). DAPI was used to counterstain the nucleus. Scale bar = 50 μM, n= 2. The corrected total cell fluorescence (CTCF) was calculated from 10 different fields of each (n) number. (E,F) Western blot analysis of histone H3 expression corrected to β-actin with representative immunoblots, n=2. Flow cytometry analysis was conducted to investigate the influence of verapamil induction in a time-dependent manner. (G) Histograms of the cell cycle phases. n=4. Data are presented as mean ± SEM values and were analyzed using one-way ANOVA with Tukey’s multiple comparisons test. ns = non-significant, *p≤0.05, **p≤0.01, ***p≤0.001, ****p≤0.0001 versus control.
Figure 2.
Protective effect of verapamil in MIN6 cells exposed to STZ, T1D-cytomix, or T2D-cytomix stressors. Bar graphs presenting MTT assay results as ratio of each treatment over viability of MIN6 cells maintained in 5.6 mM glucose. (A–C) Cells were stressed with streptozotocin (STZ). (A) Cells were pretreated with different concentrations (1, 5, 10, and 50 μM) verapamil for 24 hours, the media were changed, and cells were stressed by 3 mM of STZ for another 24 hours. (B) Cells were treated with a combination of STZ (3 mM) and verapamil (1, 5, 10, and 50 μM) for 24 hours. (C) Cells were pretreated with verapamil (1, 5, 10, and 50 μM) for 24 hours, then fresh media containing both STZ (3 mM) and verapamil (1, 5, 10, and 50 μM) were supplied to the cells for another 24 hours. (D–F) Cells were stressed with T1D-cytomix (IL-1β: 50 ng/ml, TNF-α: 50 ng/ml, and INF-γ: 100 ng/ml) in each setting independently. (D) Cells were pretreated with different concentration (1, 5, 10, and 50 μM) of verapamil for 24 hours, then media were changed, and cells were stressed with T1D-cytomix for another 24 hours. (E) Cells were treated with a combination of T1D-cytomix and verapamil (1, 5, 10, and 50 μM) for 24 hours. (F) Cells were pretreated with verapamil (1, 5, 10, and 50 μM) for 24 hours, then fresh media containing both T1D-cytomix and verapamil (1, 5, 10, and 50 μM) were supplied to the cells for another 24 hours. (G–I) Bar graphs presenting MTT assay results as ratio of each treatment over viability of cells stressed by T2D-cytomix (IL-1β: 50 ng/ml, TNF-α: 50 ng/ml, and Palmitic Acid: 500 μM) in each setting independently. (G) Cells were pretreated with different concentration (1, 5, 10, and 50 μM) of verapamil for 24 hours, then media were changed, and cells were stressed with T2D-cytomix for another 24 hours. (H) Cells were treated with a combination of T2D-cytomix and different concentration (1, 5, 10, and 50 μM) of verapamil for 24 hours. (I) Cells were pretreated with different concentration (1, 5, 10, and 50 μM) of verapamil for 24 hours, then fresh media containing both T2D-cytomix and different concentration (1, 5, 10, and 50 μM) of verapamil were supplied to the cells for another 24 hours. Each experiment was performed at least by 5 independent repeats. Difference between the groups was statistically analyzed by two-way analysis of variance (ANOVA) and p values < 0.05 considered as a significant difference. Data are presented as Mean ± SEM, * p-value< 0.05, ** p-value <0.01, *** p-value <0.001, and **** p-value <0.0001.
Figure 2.
Protective effect of verapamil in MIN6 cells exposed to STZ, T1D-cytomix, or T2D-cytomix stressors. Bar graphs presenting MTT assay results as ratio of each treatment over viability of MIN6 cells maintained in 5.6 mM glucose. (A–C) Cells were stressed with streptozotocin (STZ). (A) Cells were pretreated with different concentrations (1, 5, 10, and 50 μM) verapamil for 24 hours, the media were changed, and cells were stressed by 3 mM of STZ for another 24 hours. (B) Cells were treated with a combination of STZ (3 mM) and verapamil (1, 5, 10, and 50 μM) for 24 hours. (C) Cells were pretreated with verapamil (1, 5, 10, and 50 μM) for 24 hours, then fresh media containing both STZ (3 mM) and verapamil (1, 5, 10, and 50 μM) were supplied to the cells for another 24 hours. (D–F) Cells were stressed with T1D-cytomix (IL-1β: 50 ng/ml, TNF-α: 50 ng/ml, and INF-γ: 100 ng/ml) in each setting independently. (D) Cells were pretreated with different concentration (1, 5, 10, and 50 μM) of verapamil for 24 hours, then media were changed, and cells were stressed with T1D-cytomix for another 24 hours. (E) Cells were treated with a combination of T1D-cytomix and verapamil (1, 5, 10, and 50 μM) for 24 hours. (F) Cells were pretreated with verapamil (1, 5, 10, and 50 μM) for 24 hours, then fresh media containing both T1D-cytomix and verapamil (1, 5, 10, and 50 μM) were supplied to the cells for another 24 hours. (G–I) Bar graphs presenting MTT assay results as ratio of each treatment over viability of cells stressed by T2D-cytomix (IL-1β: 50 ng/ml, TNF-α: 50 ng/ml, and Palmitic Acid: 500 μM) in each setting independently. (G) Cells were pretreated with different concentration (1, 5, 10, and 50 μM) of verapamil for 24 hours, then media were changed, and cells were stressed with T2D-cytomix for another 24 hours. (H) Cells were treated with a combination of T2D-cytomix and different concentration (1, 5, 10, and 50 μM) of verapamil for 24 hours. (I) Cells were pretreated with different concentration (1, 5, 10, and 50 μM) of verapamil for 24 hours, then fresh media containing both T2D-cytomix and different concentration (1, 5, 10, and 50 μM) of verapamil were supplied to the cells for another 24 hours. Each experiment was performed at least by 5 independent repeats. Difference between the groups was statistically analyzed by two-way analysis of variance (ANOVA) and p values < 0.05 considered as a significant difference. Data are presented as Mean ± SEM, * p-value< 0.05, ** p-value <0.01, *** p-value <0.001, and **** p-value <0.0001.
Figure 3.
Functional effect of verapamil on the insulin content and glucose stimulated-insulin secretion (GSIS) rate of MIN6 cells. (A) Measured insulin content in MIN6 cells treated with 50 µM verapamil (Red) for 24 hours or left untreated (Blue). (B) Total secreted insulin by MIN6 cells and (C) Total insulin content in MIN6 cells, when stimulated with different concentrations of glucose and either treated with 50 µM verapamil (Red) for 24 hours or left untreated (Blue). (D) Percentage of secreted insulin over total insulin content of each condition. Each experiment was performed in 4 independent repeats. Difference between the groups was statistically analyzed using two-way analysis of variance (ANOVA) and p-values < 0.05 were considered statistically significant. Data are presented as Mean ± SEM; * p-value < 0.05; ** p-value < 0.01; *** p-value <0.001; and **** p-value <0.0001.
Figure 3.
Functional effect of verapamil on the insulin content and glucose stimulated-insulin secretion (GSIS) rate of MIN6 cells. (A) Measured insulin content in MIN6 cells treated with 50 µM verapamil (Red) for 24 hours or left untreated (Blue). (B) Total secreted insulin by MIN6 cells and (C) Total insulin content in MIN6 cells, when stimulated with different concentrations of glucose and either treated with 50 µM verapamil (Red) for 24 hours or left untreated (Blue). (D) Percentage of secreted insulin over total insulin content of each condition. Each experiment was performed in 4 independent repeats. Difference between the groups was statistically analyzed using two-way analysis of variance (ANOVA) and p-values < 0.05 were considered statistically significant. Data are presented as Mean ± SEM; * p-value < 0.05; ** p-value < 0.01; *** p-value <0.001; and **** p-value <0.0001.
Figure 4.
Verapamil pretreatment increases mitochondrial respiration in MIN6 cells treated with streptozotocin. MIN6 cells (4 × 104) were cultured in 5.6 mM glucose media and treated overnight with verapamil (50 µM) or were left untreated. Cells were subsequently treated with streptozotocin (STZ, 3 mM) or left untreated and subjected to metabolic flux analysis. (A) Metabolic flux analysis measured oxygen consumption as cells were treated with Oligomycin (1 µM), Carbonyl cyanide 4-(trifluoromethoxy)phenylhydrazone (FCCP, 2 µM), and Rotenone/Antimycin A (Rot/AA, 0.5 µM). (B) Basal respiration in metabolic flux analysis, point before Oli administration. (C) Maximal respiration in metabolic flux analysis, point before Rot/AA administration. (D) Metabolic flux analysis of cells pretreated with verapamil prior to STZ treatment. (E) Basal respiration and (F) maximal respiration in metabolic flux analysis of cells pretreated with verapamil. Difference between the groups was statistically analyzed using two-way analysis of variance (ANOVA) and p-values < 0.05 were considered statistically significant, n=4 per group, p-value < 0.05; ** p-value < 0.01; *** p-value <0.001; and **** p-value <0.0001.
Figure 4.
Verapamil pretreatment increases mitochondrial respiration in MIN6 cells treated with streptozotocin. MIN6 cells (4 × 104) were cultured in 5.6 mM glucose media and treated overnight with verapamil (50 µM) or were left untreated. Cells were subsequently treated with streptozotocin (STZ, 3 mM) or left untreated and subjected to metabolic flux analysis. (A) Metabolic flux analysis measured oxygen consumption as cells were treated with Oligomycin (1 µM), Carbonyl cyanide 4-(trifluoromethoxy)phenylhydrazone (FCCP, 2 µM), and Rotenone/Antimycin A (Rot/AA, 0.5 µM). (B) Basal respiration in metabolic flux analysis, point before Oli administration. (C) Maximal respiration in metabolic flux analysis, point before Rot/AA administration. (D) Metabolic flux analysis of cells pretreated with verapamil prior to STZ treatment. (E) Basal respiration and (F) maximal respiration in metabolic flux analysis of cells pretreated with verapamil. Difference between the groups was statistically analyzed using two-way analysis of variance (ANOVA) and p-values < 0.05 were considered statistically significant, n=4 per group, p-value < 0.05; ** p-value < 0.01; *** p-value <0.001; and **** p-value <0.0001.
Figure 5.
Transcriptomic and proteomic landscape of MIN6 cells treated with verapamil. (A) Out of a set size of 2131 genes (yellow) and proteins (pink), 2099 genes and 112 proteins were found to be differentially expressed. Of these, shared differential expressions of 32 genes and proteins were observed. (B) Volcano plot showing differentially expressed genes in MIN6 mouse beta cell line treated with verapamil (50 uM) for 24 hours (|FC| > 2, p < 0.05). (C,D) Perturbation vs. Over-representation Pathway Plot: Dots representing the top 10 impacted pathways are positioned by their p-values from two different analyses: an Impact analysis measuring total perturbation accumulation (pAcc) vs. a classical Over-representation analysis (pORA). Pathways with significant combined P-values are shown in red. The selected pancreatic secretion pathway is shown in yellow. The size of each dot denotes the total number of genes in the corresponding pathway. (E) Protein and mRNA levels of key genes found to be up- or down-regulated when exposed to verapamil treatment, along with their respective molecular functions shown in the y axis. Each cell depicts one independent experiment. The expression profiles of the respective genes are depicted by color gradient (p < 0.05, |log2FC| >1). (F) Western blot analysis of CCK expression normalized to β-actin expression; a representative immunoblot is shown.
Figure 5.
Transcriptomic and proteomic landscape of MIN6 cells treated with verapamil. (A) Out of a set size of 2131 genes (yellow) and proteins (pink), 2099 genes and 112 proteins were found to be differentially expressed. Of these, shared differential expressions of 32 genes and proteins were observed. (B) Volcano plot showing differentially expressed genes in MIN6 mouse beta cell line treated with verapamil (50 uM) for 24 hours (|FC| > 2, p < 0.05). (C,D) Perturbation vs. Over-representation Pathway Plot: Dots representing the top 10 impacted pathways are positioned by their p-values from two different analyses: an Impact analysis measuring total perturbation accumulation (pAcc) vs. a classical Over-representation analysis (pORA). Pathways with significant combined P-values are shown in red. The selected pancreatic secretion pathway is shown in yellow. The size of each dot denotes the total number of genes in the corresponding pathway. (E) Protein and mRNA levels of key genes found to be up- or down-regulated when exposed to verapamil treatment, along with their respective molecular functions shown in the y axis. Each cell depicts one independent experiment. The expression profiles of the respective genes are depicted by color gradient (p < 0.05, |log2FC| >1). (F) Western blot analysis of CCK expression normalized to β-actin expression; a representative immunoblot is shown.
Figure 6.
Verapamil pretreatment protects the pancreatic cells against MTZ-induced cytotoxic damage in zebrafish embryos. (A) Schematic timeline showing the experimental design and drug treatments. At 3 days post fertilization (dpf), four groups of embryos were treated as following: group 1, untreated embryos; group 2, embryos treated with 10 µM verapamil for 72 h; group 3, embryos treated with 10 µM verapamil for 24 h, followed by drug removal and administration of 10 mM MTZ for 48 h; group 4, embryos treated with 10 mM MTZ at 4 dpf for 48 h. Insulin-producing pancreatic cells co-expressing mCherry fluorescence reporter protein (ChFP) were imaged at 4 dpf (T0) and 6 dpf (T1). For each group, the differences in ChFP intensity between T0 and T1 were determined. (B) Representative images of pancreatic cells in each group at T1. Inserts depict the magnified ChFP area. Images were taken using Stereo discovery 1.2 ZIESS microcopy (C) Quantification of ChFP intensity in four groups at T1. No significant difference in fluorescence intensity was found between groups 1 and 2. However, the fluorescence intensity in group 3 which was pretreated with verapamil was significantly higher, as compared to group 3 which was not pretreated with verapamil. Experiments were done in triplicates (n = 20–30 embryos/group). Data represents the Mean ± SEM values. ns: non-significant.
Figure 6.
Verapamil pretreatment protects the pancreatic cells against MTZ-induced cytotoxic damage in zebrafish embryos. (A) Schematic timeline showing the experimental design and drug treatments. At 3 days post fertilization (dpf), four groups of embryos were treated as following: group 1, untreated embryos; group 2, embryos treated with 10 µM verapamil for 72 h; group 3, embryos treated with 10 µM verapamil for 24 h, followed by drug removal and administration of 10 mM MTZ for 48 h; group 4, embryos treated with 10 mM MTZ at 4 dpf for 48 h. Insulin-producing pancreatic cells co-expressing mCherry fluorescence reporter protein (ChFP) were imaged at 4 dpf (T0) and 6 dpf (T1). For each group, the differences in ChFP intensity between T0 and T1 were determined. (B) Representative images of pancreatic cells in each group at T1. Inserts depict the magnified ChFP area. Images were taken using Stereo discovery 1.2 ZIESS microcopy (C) Quantification of ChFP intensity in four groups at T1. No significant difference in fluorescence intensity was found between groups 1 and 2. However, the fluorescence intensity in group 3 which was pretreated with verapamil was significantly higher, as compared to group 3 which was not pretreated with verapamil. Experiments were done in triplicates (n = 20–30 embryos/group). Data represents the Mean ± SEM values. ns: non-significant.
Table 1.
Up-regulated and down-regulated 32 proteins and their respective genes identified by proteomics/transcriptomics analyses.
Table 1.
Up-regulated and down-regulated 32 proteins and their respective genes identified by proteomics/transcriptomics analyses.
Gene name |
logfc protein (mass spec) |
adjpv protein (mass spec) |
logfc mRNA (RNA-seq) |
adjpv mRNA (RNA-seq) |
Acss2 |
0.546956178 |
0.045430822 |
1.468439926 |
0.000001 |
Ap3m2 |
-1.258425153 |
0.004563029 |
-0.631083607 |
0.00535079 |
Atad1 |
-0.878321443 |
0.003214256 |
-0.69090951 |
0.001886914 |
Atl3 |
-0.775959726 |
0.01924926 |
-0.705494101 |
0.000782498 |
Atp2b1 |
-0.768567592 |
0.001599313 |
-0.903858193 |
0.00000179 |
Cadm1 |
6.64385619 |
0.000001 |
-1.111895426 |
0.000001 |
Camk4 |
0.618238656 |
0.025051947 |
1.667580225 |
0.000001 |
Cck |
1.343692069 |
0.000001 |
2.086382128 |
0.000001 |
Fdps |
1.236339539 |
0.000001 |
2.06424327 |
0.000001 |
Fn1 |
0.750177706 |
0.00019395 |
0.707236676 |
0.000001 |
Gabarap |
0.779049553 |
0.000817148 |
0.608883531 |
0.000001 |
Gc |
0.575312331 |
0.010226369 |
-0.668515132 |
0.000001 |
Gdap1 |
-0.744197163 |
0.03032308 |
-0.715727824 |
0.009402028 |
Hax1 |
6.64385619 |
0.000001 |
0.708603816 |
0.000001 |
Hectd1 |
1.205392513 |
0.000001 |
-0.766523599 |
0.004572659 |
Hmgcs1 |
1.608809243 |
0.000001 |
1.574403987 |
0.000001 |
Idi1 |
1.136191386 |
0.000001 |
1.063478538 |
0.00444946 |
Ins2 |
-0.706041021 |
0.003576738 |
-1.047050788 |
0.000001 |
Kif23 |
-0.899695094 |
0.029438363 |
-1.278078727 |
0.000001 |
Ldlr |
1.607862903 |
0.000001 |
1.459525662 |
0.000001 |
Mafa |
0.715454127 |
0.019937366 |
0.86170978 |
0.000001 |
Mid1ip1 |
0.82130204 |
0.000734422 |
0.902299589 |
0.000001 |
Nsg1 |
0.568518598 |
0.029438363 |
0.821403563 |
0.000001 |
Pam |
0.600269754 |
0.008748346 |
0.772163694 |
0.000001 |
Ppfia2 |
-2.53951953 |
0.000001 |
-1.455788602 |
0.000001 |
Prlr |
0.619178216 |
0.031249276 |
-0.938951004 |
0.003238752 |
Sdf2l1 |
0.606915942 |
0.007338519 |
1.333480725 |
0.000001 |
Slc25a1 |
-0.798366139 |
0.005415551 |
0.796013338 |
0.000001 |
Slc9a3r1 |
0.508935662 |
0.025006868 |
0.846755275 |
0.000001 |
Sqstm1 |
0.789103218 |
8.91352E-05 |
0.605062208 |
0.000001 |
Stc1 |
0.804053559 |
0.000223206 |
1.114392736 |
0.000001 |
Timm8b |
0.670840336 |
0.003900441 |
0.676846573 |
0.000001 |