G6PD upregulation is correlated with VSMC phenotypic switching

To understand the role of G6PD in VSMC phenotypic switching, primary VSMCs from Sprague‒Dawley (SD) rat thoracoabdominal aortic appendages were cultured and stimulated with PDGF-BB (10 ng/mL), a growth factor known to promote VSMC switching in an in vitro model. As presented in Fig. 1A, the cells were cultured for different durations (0–48 h) in serum-free DMEM to gradually induce a proliferative/synthetic phenotype [25]. One of the crucial markers of VSMC phenotypic switching is the loss of SM22α and α-actin, with switched VSMCs displaying a proliferative/synthetic phenotype [32,33,34]. To test this phenomenon, the protein levels of SM22α, α-actin, PCNA, and G6PD were analyzed at different time points during PDGF-BB stimulation via Western blot analysis. The levels of SM22α and α-actin decreased over time under PDGF-BB stimulation, indicating a gradual switch in VSMCs from a contractile/quiescent phenotype to a synthetic/proliferating phenotype. The level of PCNA, a marker of proliferation, gradually increased at 6 h and peaked at 12 h. Moreover, the expression of G6PD also increased in a time-dependent manner with PDGF-BB stimulation, reaching its peak at 12 h (Fig. 1B). In addition, G6PD is found in diverse oligomeric states, including monomers, dimers, tetramers, and hexamers, but only the dimeric and tetrameric forms are catalytically active [35]. Our observations revealed significant increases in both G6PD dimerization (Additional file 1: Figure S1A) and activity (Additional file 1: Figure S1B) following stimulation with PDGF-BB for 12 h. In summary, G6PD was upregulated and activated in synthetic/proliferated VSMCs, suggesting that G6PD plays an indispensable role in VSMC phenotypic switching.

Fig. 1

G6PD upregulation promotes VSMC phenotypic switching, as determined by bioinformatics analysis of G6PD-interacting proteins. A Schematic diagram of the phenotypic switching of VSMCs induced by platelet-derived growth factor-BB (PDGF-BB) stimulation (20 ng/mL). Primary VSMCs were cultured from the thoracoabdominal aortas of male Sprague–Dawley rats, cultured for 3–5 generations and then stimulated with or without PDGF-BB. The whole-cell lysates were subjected to Western blot analysis. B Representative Western blot analysis and analysis of grayscale images of G6PD, PCNA, and VSMC contractile markers (α-actin and SM22α) in the lysates of VSMCs from A, which were treated with PDGF-BB for 0, 12, 24, 36 and 48 h, respectively. β-actin was used as an internal control. The data are presented as the relative fold change at 0 h (n = 3). C Representative Western blot analysis and analysis of grayscale images of G6PD, PCNA, α-actin, and SM22α in VSMCs stimulated with or without PDGF-BB in response to siG6PD. n = 3. D Venn diagrams displaying the number of unique and shared proteins identified by LC‒MS between control, PDGF-BB, and IgG VSMCs. E–G The functional categories based on gene ontology (GO) term enrichment and the signaling pathways based on KEGG enrichment analysis. biological process (E), cellular component (F), and molecular function (G). H Mass spectrometry analysis of G6PD and VDAC1 peptides after purification of G6PD-associated proteins. I, J Three peptides of VDAC1 identified by MS (I) and their positions in the 3D structure of VDAC1 (J). Statistical significance was determined using one-way ANOVA in (B) and two-tailed Student’s t tests in (C). *P < 0.05; **P < 0.01; ns, no significant difference

To determine the essential role of G6PD in PDGF-BB-induced VSMC phenotypic switching, three siRNAs were designed to reduce the expression of G6PD, and siRNA-241 had the greatest interference efficiency (Additional file 1: Figure S1C). Besides, the G6PD activity was tested, and the result showed that knockdown of G6PD reduced the enzyme activity of G6PD in VSMCs (Additional file 1: Figure S1D). Interestingly, the increase in PCNA and decrease in SM22α and α-actin were reversed by siG6PD in proliferative/synthetic VSMCs (Fig. 1C), indicating that G6PD knockdown reversed PDGF-BB-induced VSMC phenotypic switching. To corroborate the findings of G6PD silencing, we studied cultured proliferative/synthetic VSMCs treated with the G6PD inhibitor 6-aminonicotinamide (6AN) [36]. Consistent with the above results, pharmacological inhibition of G6PD also reduced the enzyme activity of G6PD in VSMCs (Additional file 1: Figure S1E). The increase in PCNA and decrease in the contractile markers SM22α and α-actin induced by PDGF-BB were blocked by 6-AN (Additional file 1: Figure S1F). Taken together, these results suggest that increased expression and activity of G6PD contribute to phenotypic switching of VSMCs, which is restrained by pharmacological inhibition or knockdown of G6PD.

VDAC1 was identified as a novel interacting partner of G6PD

To explore the molecular mechanisms by which G6PD regulates VSMC phenotypic switching, we sought to identify protein complexes that interact with G6PD following stimulation with PDGF-BB. Immunoprecipitation (IP) coupled with mass spectrometry (MS) analysis was performed to systematically analyze G6PD-interacting proteins in VSMCs with (+) or without (−) PDGF-BB stimulation. Rabbit IgG was used as a mock control to exclude nonspecific interactions. The immunoprecipitated proteins were detected using antibodies specific for G6PD or IgG by SDS-polyacrylamide gel electrophoresis and Coomassie Brilliant Blue staining (Additional file 1: Figure S1G).

A total of 526, 905, and 555 proteins were identified in the control, PDGF-BB, and IgG groups, respectively. Venn diagram analysis of the differentially expressed proteins (fold change > 2, P < 0.05) between PDGF-BB-treated VSMCs and potential target proteins revealed 360 candidate target proteins, including voltage-dependent anion-selective channel protein 1 (VDAC1) (Fig. 1D). VDAC1 is a multifunctional mitochondrial protein that regulates cellular metabolic and energy functions as well as apoptosis by interacting with a variety of proteins, prompting us to explore its roles in VSMC biology.

A protein‒protein interaction (PPI) regulatory network was constructed with the STRING database and visualized with Cytoscape (Additional file 1: Figure S1H). The CytoNCA plugin was used to calculate the nodes with the highest degree scores [37]. To visualize the data, an online platform (http://www.bioinformatics.com.cn) was used to construct heatmaps and Venn diagrams. Biological process enrichment analysis based on gene ontology (GO) terms was also conducted to explore the possible role of G6PD, which revealed that the proteins associated with the negative regulation of the apoptotic process were significantly enriched (Fig. 1E). Cellular component analysis revealed that the upregulated proteins were enriched in the cytoplasm, extracellular exosome, nucleus, membrane, and mitochondrion (Fig. 1F). The analysis of molecular functions indicated that the upregulated proteins were significantly enriched in poly(A) RNA binding, molecule binding, and ion channel binding (Fig. 1G). Furthermore, to investigate the involvement of development-related signaling pathways, KEGG analysis was conducted, and metabolic pathways were found to be the most prominent component among the enriched pathways (Additional file 1: Figure S1I).

Based on the above findings, VDAC1, a gated protein associated with apoptosis, attracted our attention. The analysis revealed VDAC1 as a major interacting partner of G6PD (Fig. 1H). Three VDAC1 peptides were identified by mass spectrometry (Fig. 1I) and are displayed in the 3D structure of VDAC1 (Fig. 1J), suggesting that G6PD and VDAC1 physically interact and participate in functional crosstalk.

G6PD translocates to the mitochondria and interacts with VDAC1, influencing mitochondrial function in proliferative VSMCs

To ascertain the potential interaction between VDAC1, a protein located on the outer mitochondrial membrane [38], and G6PD, we initially evaluated the presence of G6PD in the mitochondria. Proteins extracted from enriched mitochondrial and cytosolic fractions were subjected to Western blot analysis to detect the expression levels of G6PD. The results indicated an increase in G6PD expression in mitochondria following PDGF-BB stimulation (Fig. 2A). To confirm this result, fluorescence colocalization analysis was performed using LAS AF software (Leica Microsystems, Germany). In detail, mitochondria and G6PD were labeled with MitoTracker (red) and a monoclonal anti-G6PD antibody, respectively, followed by a secondary antibody conjugated with FITC (green). As expected, the level of G6PD in mitochondria increased in response to PDGF-BB stimulation (yellow) (Fig. 2B). In addition, we examined the effect of G6PD on mitochondrial function. A Seahorse extracellular flux analyzer was used to measure the oxygen consumption rate (OCR) and extracellular acidification rate (ECAR), which are indirect measures of mitochondrial OXPHOS and glycolysis, respectively. Oligomycin, FCCP (carbonyl cyanide-4-(trifluoromethoxy) phenylhydrazone), and a mixture of rotenone A/antimycin were injected sequentially to assess mitochondrial respiration and nonmitochondrial respiration. The results revealed that the optimal concentration of FCCP for treating VSMCs was 0.5 μM (Additional file 1: Figure S2). Proliferative/synthetic VSMCs displayed a greater maximum OCR and ECAR, suggesting a greater energy demand within these cells. Treatment with 6AN or siG6PD significantly decreased the OCR and ECAR, particularly in G6PD-knockdown cells, in which there was almost a complete decrease in the OCR and ECAR. However, the basal OCR did not significantly differ among the three groups, except for the siG6PD group (Fig. 2C, D). Taken together, these findings suggest that G6PD translocates to mitochondria and influences the function of mitochondria in proliferative/synthetic VSMCs.

Fig. 2

G6PD translocates to mitochondria and affects mitochondrial function. A Mitochondrial and cytoplasmic fractions were prepared from VSMCs treated with 20 ng/mL PDGF-BB for 12 h. Western blot analysis of the cytosolic and mitochondrial fractions was performed to evaluate the translocation of G6PD from the cytosolic compartment to the mitochondria. GAPDH and TOM40 were used as the cytosolic and mitochondrial loading controls, respectively. B Immunofluorescence analysis was carried out after 12 h of PDGF-BB stimulation, and 6AN and siG6PD were also used in subsequent experiments. Mitochondria were identified with TOM40, nuclei were stained with DAPI, and a G6PD monoclonal antibody was used to indicate endogenous G6PD. Scale bar = 50 µm. n = 3. C, D Seahorse metabolic flux analyses showing the traces and quantification of the basal or maximum mitochondrial oxygen consumption rate (OCR, C) and extracellular acidification rate (ECAR, D) in VSMCs treated with 6AN and siG6PD under PDGF-BB stimulation. n = 3. Statistical significance was determined using two-tailed Student’s t tests in (A) and one-way ANOVA in (B–D). *P < 0.05; **P < 0.01; ***P < 0.001

VDAC1 was reported to be expressed in both the mitochondrial and cytoplasmic membranes [39]. To investigate the effect of PDGF-BB stimulation on the localization of VDAC1, we utilized laser confocal microscopy to observe VDAC1 and MitoTracker Red to determine its localization in VSMCs. The results indicated that VDAC1 was localized to mitochondria, regardless of PDGF-BB stimulation (Additional file 1: Figure S3A). In addition, the VDAC1 protein expression underwent a minor change upon PDGF-BB stimulation (Additional file 1: Figure S3B). To further decipher whether VDAC1 is part of the protein complex that interacts with G6PD, coimmunoprecipitation analysis was conducted to assess the interactions between VDAC1 and G6PD. In the VSMC extracts, the G6PD-specific antibody precipitated VDAC1, while the VDAC1-specific antibody precipitated G6PD. Furthermore, the interaction between these two proteins was enhanced after PDGF-BB stimulation (Fig. 3A, B). Notably, the G6PD-specific antibody precipitated VDAC2 and VDAC3, but their interaction was not affected by PDGF-BB stimulation (Additional file 1: Figure S3C, D). In addition, to avoid interference from additional proteins, we synthesized the eukaryotic expression vectors pcDNA3.1-HA-G6PD and pcDNA3.1-Flag-VDAC1 and coexpressed them in 293A cells. Coimmunoprecipitation was subsequently performed using anti-HA or anti-Flag antibodies. Consistent with the findings in VSMCs, the interaction of G6PD with VDAC1 was previously observed (Fig. 3C, D). To determine whether G6PD could directly interact with VDAC1, we performed a GST pull-down assay. Encouragingly, VDAC1 directly interacted with G6PD (Fig. 3E). This hypothesis was further supported by fluorescence colocalization. However, the interaction was reduced when the G6PD concentration was suppressed or pharmacologically inhibited (Fig. 3F–H). To determine the regions within G6PD that are associated with VDAC1, we coexpressed the HA-tagged full length (FL) and truncated N-terminal domain (HA-G6PD-NTD, aa). 1–210), and C-terminal domain (HA-G6PD-CTD, aa. 121–515) of G6PD with Flag-tagged VDAC1 in HEK293T cells (Fig. 3I; Additional file 1: Figure S3E). Specifically, VDAC1 interacted with the wild-type or N-terminal domain but not with the C-terminal domain (Fig. 3J). The G6PD-NTD region was also found to predominantly contribute to this interaction in VSMCs. Collectively, these results provide ample evidence supporting the direct interaction between G6PD and VDAC1, with G6PD-NTD predominantly contributing to this interaction.

Fig. 3

G6PD interacts with VDAC1 under PDGF-BB stimulation. A, B Endogenous VDAC1–G6PD interactions in VSMCs detected by coimmunoprecipitation (IP) experiments. VSMCs were left unstimulated or stimulated with PDGF-BB for 12 h. Cell lysates were immunoprecipitated (IP) and immunoblotted (IB) with the indicated antibodies. C, D HEK293A cells were cotransfected with the Flag-VDAC1 plasmid and the HA-G6PD plasmid. The cell lysates were precipitated with an anti-HA or anti-Flag antibody, and the precipitates were analyzed by immunoblotting with an anti-Flag or anti-HA antibody, respectively. E GST pulldown analysis using recombinant glutathione S-transferase (GST)-G6PD and Flag-tagged VDAC1 (Flag-VDAC1). F The colocalization of G6PD and VDAC1 in the mitochondria of VSMCs stimulated with or without PDGF-BB. Treatment with 6AN and siG6PD was further observed. Scale bar = 10 µm. G Histogram showing the arbitrary intensities of G6PD and VDAC1 across the white arrow shown in the zoomed and merged images in F. H The Pearson correlation coefficient for signal colocalization was determined using Fiji software. I Schematic illustration of the G6PD domains used to evaluate the interaction with VDAC1. J HEK293A cells were cotransfected with the Flag-VDAC1 plasmid and the HA-G6PD-FL (FL), HA-G6PD-NTD (NTD), or HA-G6PD-CTD (CTD) plasmids. The cell lysates were immunoprecipitated with an anti-HA antibody, and the precipitates were analyzed by immunoblotting with an anti-Flag antibody. K VSMCs were transfected with HA-G6PD-FL (FL), HA-G6PD-NTD (NTD), or HA-G6PD-CTD (CTD) plasmids. The cell lysates were immunoprecipitated with an anti-HA antibody, and the precipitates were analyzed by immunoblotting with an anti-VDAC1 antibody

G6PD inhibition induces VSMC apoptosis and decreases cell viability

VDAC1 is a protein channel located in the outer mitochondrial membrane [40, 41]. In addition to its role as an ion channel, VDAC1 is a well-known regulator of cell apoptosis [42,43,44]. Therefore, our interest was sparked in investigating whether G6PD is involved in mediating apoptosis. We analyzed the expression of various apoptosis-related proteins (Caspase7, Caspase9, PARP, Bcl-2, Cleaved Caspase7, Cleaved Caspase9, and Cleaved PARP) in VSMCs. Synchronously, the levels of the proapoptotic proteins cleaved caspase 7, cleaved caspase 9, and cleaved PARP were decreased, while the level of the antiapoptotic protein Bcl-2 was increased in VSMCs stimulated with PDGF-BB. However, 6AN reversed the changes in the levels of these proteins (Fig. 4A). Apoptosis was detected using the Annexin V FITC Apoptosis Detection Kit from BD Biosciences. Flow cytometry analysis revealed that the apoptosis rate of the untreated VSMCs was approximately 7%. Treatment with PDGF-BB reduced the apoptosis rate of VSMCs to approximately 3.5%, while treatment with 6-AN restored the apoptosis rate to a level similar to that of the control group, which was approximately 9% (Fig. 4B). A Cell Counting Kit-8 (CCK8) assay revealed that VSMC viability was enhanced upon stimulation with PDGF-BB. However, the viability of the VSMCs treated with 6-AN was even lower than that of the control group (Fig. 4C). A Cell Counting Kit-8 (CCK8) assay demonstrated that VSMC viability increased in response to stimulation with PDGF-BB, while the viability of VSMCs treated with 6-AN was lower than that of the control group (Fig. 4D). These findings demonstrated that pharmacologically inhibiting G6PD can enhance VSMC apoptosis and decrease VSMC viability.

Fig. 4

Both pharmacological inhibition and knockdown of G6PD abolish PDGF-BB-induced VSMC apoptosis inhibition and survival. A Western blot analysis of caspase-9, cleaved caspase-9, caspase-7, cleaved caspase-7, PARP, cleaved PARP, and Bcl-2 in VSMCs. n = 3. B Cell apoptosis was evaluated by PE-conjugated Annexin V/7-AAD staining and quantified by flow cytometry. n = 3. C The viability of VSMCs was assessed by a CCK8 assay. n = 22. D Cell proliferation was measured by a cell counting experiment. n = 10. E Western blot analysis of caspase-9, cleaved caspase-9, caspase-7, cleaved caspase-7, PARP, cleaved PARP, and Bcl-2 in VSMCs. n = 3. F Cell apoptosis was evaluated by PE-conjugated Annexin V/7-AAD staining and quantified by flow cytometry. n = 3. G The viability of VSMCs was assessed by a CCK8 assay. n = 7. H Cell proliferation was measured by a cell counting experiment, n = 3. Statistical significance was determined using one-way ANOVA. *P < 0.05; **P < 0.01; ***P < 0.001

To exclude the potential off-target effects of pharmacological inhibitors, siRNA was used to knock down G6PD in proliferative/synthetic VSMCs. Consistent with the pharmacological results, G6PD knockdown increased the levels of the proapoptotic proteins cleaved caspase 7, cleaved caspase 9, and cleaved PARP and reduced the level of the antiapoptotic protein Bcl-2 (Fig. 4E), suggesting that G6PD promoted VSMC apoptosis in the presence or absence of PDGF-BB stimulation, and flow cytometry analysis promoted this conversion (Fig. 4F). Accordingly, CCK-8 and cell counting assays revealed that G6PD silencing suppressed the cell proliferation (Fig. 4G) and viability (Fig. 4H) induced by PDGF-BB. Overall, these findings provide evidence that pharmacological inhibition or silencing of G6PD suppresses VSMC survival by promoting apoptosis and inhibiting proliferation.

VDAC1 is necessary for G6PD-mediated antiapoptotic effects, and the NTD of G6PD plays a major role in this process

Given that VDAC1 regulates cell apoptosis by interacting with antiapoptotic (Bcl-2 and Bcl-xL) or proapoptotic (Bax, Bak, and Bim) proteins [41], in light of the antiapoptotic function of G6PD in VSMC phenotypic switching and its increased interaction with VDAC1 during the antiapoptotic process, we hypothesized that G6PD could inhibit VSMC apoptosis in a VDAC1-dependent manner. To test our hypothesis, we cotransfected VSMCs with siRNAs targeting VDAC1 and an adenovirus encoding the G6PD open reading frame (Ad-G6PD) and subsequently induced them with PDGF-BB for 12 h. Three siRNAs targeting VDAC1 (si-168, si-261, si-423) were designed, and siRNA-168 had the highest interference efficiency (Additional file 1: Figure S4A). The optimal virus concentrations for cell transfection were 7.5 × 107 PFU/mL for Ad-G6PD and 3.0 × 107 PFU/mL for the control adenovirus (Ad-con) (Additional file 1: Figure S4B).

TUNEL assays showed that silencing VDAC1 led to reduced apoptosis induced by PDGF-BB stimulation, but G6PD overexpression did not reverse this effect (Fig. 5A, B). Knockdown of VDAC1 blocked PDGF-BB-induced cell viability and proliferation, while G6PD overexpression did not reverse these effects, as demonstrated by the CCK-8 and cell counting assays (Fig. 5C, D). Coincidentally, G6PD-NTD was involved mainly in the inhibition of VSMC apoptosis (Additional file 1: Figure S4C; Fig. 5E) and promoted the survival of VSMCs (Fig. 5F). Overall, these findings suggest that VDAC1 is necessary for G6PD-mediated antiapoptotic effects on VSMC phenotypic switching and that the N-terminal domain of G6PD plays a major role in this process.

Fig. 5

The antiapoptotic effects of G6PD on PDGF-BB-induced VSMC proliferation and viability are mediated by VDAC1. The expression of VDAC1 was reduced by siRNA in VSMCs, and G6PD was subsequently overexpressed by an adenovirus encoding an HA tag. Finally, the cells were treated with PDGF-BB for 12 h. A, B TUNEL analysis was used to analyze cell apoptosis, and the TUNEL-positive cell ratio was quantified from 5 microscopic views in each experiment. n = 5. Scale bar = 100 µm. C The viability of VSMCs was assessed by a CCK8 assay. n = 5. D Cell proliferation was measured by a cell counting experiment. n = 6. Ad-Null and si-Con were used as negative controls. E G6PD-NTD virus infection inhibited VSMC apoptosis, and western blotting was performed three times; representative images are presented. F G6PD-NTD virus infection promoted VSMC survival. Statistical significance was determined using two-way ANOVA in B–D, F. *P < 0.05; **P < 0.01; ***P < 0.001; ns, no significant difference

G6PD blocks the interaction between VDAC1 and Bax by reducing the oligomerization of VDAC1

Previous studies have demonstrated that VDAC1 can facilitate apoptosis by interacting with the proapoptotic molecule Bax. The molecular mechanism involves Bax translocating from the cytoplasm to the mitochondria and forming a hetero-oligomer with VDAC1 in the outer mitochondrial membrane. This interaction promotes the release of cytochrome c (CytoC) or apoptosis-inducing factor (AIF) from the intermembrane space to the cytoplasm, ultimately leading to apoptosis [45,46,47]. Furthermore, VDAC1 can also bind to the antiapoptotic protein Bcl-2, thereby inhibiting apoptosis [48]. We investigated whether G6PD is involved in the apoptotic process mediated by Bax or Bcl-2. We employed immunoprecipitation to evaluate the interaction of VDAC1 with Bax or Bcl-2 in the presence or absence of PDGF-BB stimulation. The results indicated that the interaction of VDAC1 with Bax was reduced (Fig. 6A) in PDGF-BB-induced synthetic VSMCs, while the interaction with Bcl-2 remained unchanged (Fig. 6B). Thus, we postulate that G6PD may compete with Bax for binding to VDAC1.

Fig. 6

G6PD competes with Bax for binding VDAC1 by reducing its oligomerization. A, B Coimmunoprecipitation (IP) assays were used to detect the interaction between VDAC1-Bax (A) and VDAC1-Bcl-2 (B) in VSMCs stimulated with or without PDGF-BB for 12 h, respectively. C, D The interaction between VDAC1 and Bax in VSMCs was examined by IP after the administration of 6-AN (C) or siG6PD (D). E, F After the addition of 6-AN (E) or siG6PD (F), G6PD was overexpressed, and the interaction between VDAC1 and Bax in VSMCs was examined by IP. G–I Oligomerization of VDAC1 was detected by WB analysis with or without PDGF-BB stimulation for 12 h (G). After treatment with siG6PD (H) or 6AN (I), G6PD-overexpressing cells were used for further study

To test our hypothesis, we treated cells with 6-AN or siRNA to inhibit G6PD and subsequently monitored the interaction between VDAC1 and Bax. As anticipated, VDAC1 restored the interaction with Bax when G6PD was knocked down or pharmacologically inhibited (Fig. 6C, D). However, these effects were abolished when G6PD expression was restored via adenovirus (Fig. 6E, F). These results suggest that G6PD suppresses the binding of VDAC1 to Bax, thereby inhibiting VSMC apoptosis.

The process of apoptosis, triggered by the binding of Bax to VDAC1, is often linked to an increase in the oligomerization of VDAC1. On the other hand, a decrease in apoptosis is typically associated with a reduction in VDAC1 oligomerization [49, 50]. To investigate the impact of competitive binding between G6PD and Bax on VDAC1 oligomerization, we first aimed to determine whether VDAC1 oligomerization was reduced in PDGF-BB-induced antiapoptotic VSMCs (Fig. 6G). Pharmacological inhibition or knockdown of G6PD led to the restoration of VDAC1 oligomerization, whereas G6PD overexpression reversed this effect (Fig. 6H, I). Overall, the results suggest that G6PD competes with Bax for binding to VDAC1, leading to a reduction in VDAC1 oligomerization.

The VDAC1 oligomerization inhibitor VBIT-12 mimics the antiapoptotic function of G6PD

To investigate whether PDGF-BB-induced VSMC apoptosis inhibition was mediated by reduced VDAC1 oligomerization, we treated VSMCs with the VDAC1 oligomerization inhibitor VBIT-12 under PDGF-BB stimulation and then overexpressed Ad-G6PD [38]. VBIT-12 inhibited VDAC1 oligomerization and subsequent apoptosis and associated processes, such as AIF and CytoC release from the mitochondrial intermembrane space (Fig. 7A). Consistent with our hypothesis, similar to G6PD overexpression, VBIT-12 treatment led to a decrease in PDGF-BB-induced VSMC apoptosis, as indicated by the TUNEL assay. However, the effect was only marginally enhanced with additional G6PD overexpression (Fig. 7B). These findings suggest that the inhibition of VDAC1 oligomerization may involve a mechanism similar to that of G6PD-mediated VSMC antiapoptotic effects. Cell viability and proliferation were significantly increased in G6PD-overexpressing VSMCs, but these effects were not evident after further PDGF-BB stimulation. Interestingly, VBIT-12 had no significant effect on cell viability (Fig. 7C) or proliferation (Fig. 7D) mediated by PDGF-BB and G6PD overexpression.

Fig. 7

VBIT-12 mimicked the antiapoptotic effect of G6PD. A The molecular mechanism by which VBIT-12 inhibits VDAC1 oligomerization and apoptosis resistance. B TUNEL analysis was used to analyze cell apoptosis, and the TUNEL-positive cell ratio was quantified in 5 microscopic fields of view in each experiment. Scale bar = 100 µm. C The viability of VSMCs was assessed by a CCK8 assay. n = 7. D Cell proliferation was measured by a cell counting assay. n = 6. Ad-Null and DMSO were used as negative controls. Statistical significance was determined using two-tailed Student’s t tests (B) and two-way ANOVA (C, D); *P < 0.05; ***P < 0.001; ns, no significant difference

Overall, in synthetic VSMCs, G6PD plays a dual role in promoting cell survival. First, VDAC1 competes with Bax for binding, inhibiting its oligomerization and promoting antiapoptotic cell death. Second, G6PD promotes cell viability and proliferation, thus further promoting VSMC survival.

G6PD accelerates vascular neointimal hyperplasia by promoting proliferation and inhibiting VSMC apoptosis

The in vivo role of G6PD was investigated in a mouse model of neointimal hyperplasia induced by common carotid artery ligation. Two weeks after ligation, the mouse model was established. The expression of G6PD was frequently elevated in the ligated arteries exhibiting neointimal hyperplasia (Fig. 8A).

Fig. 8

Overexpression of G6PD exacerbates neointimal hyperplasia by decreasing VSMC apoptosis. A Representative cross-sections and quantification of G6PD levels following IF staining after ligation surgery of the carotid artery for 14 days. Scale bar = 100 µm. B Immunofluorescence of mCherry after 14 days of in situ delivery of Ad-Null or Ad-G6PD into ligated mouse common carotid arteries. Scale bar = 100 µm. C Hematoxylin and eosin (HE)-stained cross-sections were obtained, and the intimal/media ratio was measured 14 days after carotid ligation surgery. Scale bar = 100 µm. D Immunofluorescence and quantification of cleaved caspase 3 levels after 14 days of in situ delivery of Ad-Null or Ad-G6PD into ligated mouse common carotid arteries. Scale bar = 100 µm. E Immunofluorescence and quantification of cleaved caspase 9 levels after 14 days of in situ delivery of Ad-Null or Ad-G6PD into ligated mouse common carotid arteries. Scale bar = 100 µm. F Immunofluorescence and quantification of PCNA levels after 14 days of in situ delivery of Ad-Null or Ad-G6PD into ligated mouse common carotid arteries. Scale bar = 100 µm. G Immunofluorescence and quantification of TUNEL-positive cells after 14 days of in situ delivery of Ad-Null or Ad-G6PD into ligated mouse common carotid arteries. Scale bar = 100 µm. The data are shown as the mean ± SEM; n ≥ 5. Statistical significance was determined using two-tailed Student’s t tests (A, D–F) and one-way ANOVA (C). *P < 0.05; **P < 0.01; ***P < 0.001; ns, not significant

An adenovirus encoding G6PD bearing a mCherry (red) fluorescence protein (Ad-G6PD) was generated in situ in the ligated arteries, and an adenovirus‐encoding mCherry protein (Ad‐Null) was used as a negative control. Red fluorescence was detected in both arteries and whole blood on day 14, indicating successful in situ delivery (Fig. 8B). After 14 days of in situ delivery, G6PD overexpression exacerbated neointimal formation, as confirmed by an increased intima-to-media ratio (Fig. 8C). Similarly, the levels of the apoptosis-related proteins cleaved caspase 3 (Fig. 8D) and cleaved caspase 9 (Fig. 8E) decreased, while the levels of PCNA (Fig. 8F) and the percentage of TUNEL-positive cells (Fig. 8G) increased. Therefore, G6PD overexpression exacerbates neointimal hyperplasia by promoting proliferation and decreasing VSMC apoptosis.

In addition to increasing apoptosis inhibition and proliferation, the synthetic phenotype also enhances fibrosis and extracellular matrix secretion. IF staining of collagen1 and MMP2 could provide additional evidence that in vivo G6PD overexpression results in vascular dysfunction linked to stenosis or atherosclerosis, due to intimal tunica growth. Our findings demonstrated a notable increase in the extracellular matrix component collagen 1 (Additional file 1: Figure S5A) and the fibrosis marker MMP2 (Additional file 1: Figure S5B) in vascular tissues overexpressing G6PD compared to those in the Ad-Null group, supporting our initial hypothesis.