GLI1 is a transcriptional effector at the terminal end of the Hedgehog signaling (Hh) pathway and is tightly regulated during embryonic development and tissue patterning/differentiation. GLI1 has low-level expression in differentiated tissues, however, in certain cancers, aberrant activation of GLI1 has been linked to the promotion of numerous hallmarks of cancer, such as proliferation, survival, angiogenesis, metastasis, metabolic rewiring, and chemotherapeutic resistance. All of these are driven, in part, by GLI1’s role in regulating cell cycle, DNA replication and DNA damage repair processes. The consequences of GLI1 oncogenic activity, specifically the activity surrounding DNA damage repair proteins, such as NBS1, and cell cycle proteins, such as CDK1, can be linked to tumorigenesis and chemoresistance. Therefore, understanding the underlying mechanisms driving GLI1 dysregulation can provide prognostic and diagnostic biomarkers to identify a patient population that would derive therapeutic benefit from either direct inhibition of GLI1 or targeted therapy towards proteins downstream of GLI1 regulation.

Pancreatic cancer is a highly lethal disease with obesity as one of the risk factors. Oncogenic KRAS mutations are prevalent in pancreatic cancer and can rewire lipid metabolism by altering fatty acid (FA) uptake, FA oxidation (FAO), and lipogenesis. Identification of the underlying mechanisms could lead to improved therapeutic strategies for treating KRAS-mutant pancreatic cancer. Here, we observed that KRASG12D upregulated the expression of SLC25A1, a citrate transporter that is a key metabolic switch to mediate FAO, fatty acid synthesis, glycolysis, and gluconeogenesis. In genetically engineered mouse models and human pancreatic cancer cells, KRASG12D induced SLC25A1 upregulation via GLI1, which directly stimulated SLC25A1 transcription by binding its promoter. The enhanced expression of SLC25A1 increased levels of cytosolic citrate, FAs, and key enzymes in lipid metabolism. In addition, a high-fat diet (HFD) further stimulated the KRASG12D-GLI1-SLC25A1 axis and the associated increase in citrate and FAs. Pharmacologic inhibition of SLC25A1 and upstream GLI1 significantly suppressed pancreatic tumorigenesis in KrasG12D/+ mice on a HFD. These results reveal a KRASG12D-GLI1-SLC25A1 regulatory axis, with SLC25A1 as an important node that regulates lipid metabolism during pancreatic tumorigenesis, thus indicating an intervention strategy for oncogenic KRAS-driven pancreatic cancer. Significance: Upregulation of SLC25A1 induced by KRASG12D-GLI1 signaling rewires lipid metabolism and is exacerbated by HFD to drive the development of pancreatic cancer, representing a targetable metabolic axis to suppress pancreatic tumorigenesis.

Resistance to radiation and chemotherapy in colorectal cancer (CRC) patients contribute significantly to refractory disease and disease progression. Herein, we provide mechanistic rationale for acquired or inherent chemotherapeutic resistance to the anti-tumor effects of 5-fluorouracil (5-FU) that is linked to oncogenic GLI1 transcription activity and NBS1 overexpression. Patients with high levels of GLI1 also expressed high levels of NBS1. Non-canonical activation of GLI1 is driven through oncogenic pathways in CRC, like the BRAF^V600E mutation. GLI1 was identified as a novel regulator of NBS1 and discovered that by knocking down GLI1 levels in vitro, diminished NBS1 expression, increased DNA damage/apoptosis, and re-sensitization of 5-FU resistant cancer to treatment was observed. Furthermore, a novel GLI1 inhibitor, SRI-38832, which exhibited pharmacokinetic properties suitable for in vivo testing, was identified. GLI1 inhibition in a murine BRAF^V600E variant xenograft model of CRC resulted in the same down-regulation of NBS1 observed in vitro as well as significant reduction of tumor growth/burden. GLI1 inhibition could therefore be a therapeutic option for 5-FU resistant and BRAF^V600E variant CRC patients.

Epigenetic dysregulation characterized by aberrant DNA hypermethylation is a hallmark of cancer, and it can be targeted by hypomethylating agents (HMAs). Recently, we described the superior therapeutic efficacy of a novel HMA, namely, NTX-301, when used as a monotherapy and in combination with venetoclax in the treatment of acute myeloid leukemia. Following a previous study, we further explored the therapeutic properties of NTX-301 based on experimental investigations and integrative data analyses. Comprehensive sensitivity profiling revealed that NTX-301 primarily exerted anticancer effects against blood cancers and exhibited improved potency against a wide range of solid cancers. Subsequent assays showed that the superior efficacy of NTX-301 depended on its strong effects on cell cycle arrest, apoptosis, and differentiation. Due to its superior efficacy, low doses of NTX-301 achieved sufficiently substantial tumor regression in vivo. Multiomics analyses revealed the mechanisms of action (MoAs) of NTX-301 and linked these MoAs to markers of sensitivity to NTX-301 and to the demethylation activity of NTX-301 with high concordance. In conclusion, our findings provide a rationale for currently ongoing clinical trials of NTX-301 and will help guide the development of novel therapeutic options for cancer patients.

As a transcription factor, GLI1 plays an important role in cell cycle regulation, DNA replication, and DNA damage responses. The aberrant activation of GLI1 has been associated with cancers such as glioma, osteosarcoma, and rhabdomyosarcoma. The binding of GLI1 to a specific DNA sequence was achieved by five tandem zinc finger motifs (Zif motifs) on the N-terminal part of the molecule. Here, we reported a novel homodimeric crystal structure of Zif1-2. These two Zif motifs are linearized. Namely, Zif2 does not bend and interact with Zif1 of the same molecule. Instead, Zif1 from one molecule interacts with Zif2 from another molecule. The dimer interface of Zif1-2 is unique and different from the conformation of Zif1-2 from the GLI1-DNA co-crystal structure. The dimeric conformation of Zif motifs could represent the native conformation of apo form GLI1 Zif motifs in the cell. The molecular dynamics simulation result of the homodimer, the in silico mutagenesis, and the predicted protease stability of these mutants using a large language model are also presented.