Understand KRAS and the Quest for Anti-Cancer Drugs

Design and development of dual targeting CAR protein for the development of CAR T-cell therapy against KRAS mutated pancreatic ductal adenocarcinoma using computational approaches

Mutant KRAS promotes the proliferation, metastasis, and aggressiveness of various cancers including pancreatic ductal adenocarcinoma (PDAC), non-small cell lung cancer (NSCLC), and colorectal adenocarcinoma (CRC) respectively. Mutant KRAS therapeutics are limited, while Sotorasib and Adagrasib were the only FDA-approved drugs for the treatment of KRASG12C mutated NSCLC. Chimeric antigen receptor (CAR) T-cell therapy has been emerged as an effective strategy against hematological malignancies and being extended towards solid cancers including PDAC. mesothelin (MSLN) and Carcinoembryonic Antigen (CEA) were reported to be highly overexpressed in KRAS-mutated PDAC. Meanwhile, in clinical trials, several CAR T-cell therapy studies are mainly focused towards these two cancer antigens in PDAC, however, the dual targeting of these two neoantigens is not reported. In the present study, we have designed and developed a novel dual-targeting CAR protein by employing various bioinformatics approaches such as functional analysis (antigenicity, allergenicity, antigen binding sites & signalling cascades), qualitative analysis (physicochemical, prediction, refinement & validation of 2D and 3D structures), molecular docking, and in silico cloning. Our results revealed that the designed CAR protein specifically binds with both MSLN & CEA with significant binding affinities, and was predicted to be stable & non-allergenic. Additionally, the protein-protein interaction network reveals the T-cell mediated antitumor responses of each domain in the designed CAR. Conclusively, we have designed and developed a dual targeting (MSLN & CEA) CAR protein towards KRAS-mutated PDAC using computational approaches. Alongside, we further recommend to engineer this designed CAR in T-cells and evaluating their therapeutic efficiency in in vitro and in vivo studies in the near future.

Regulation of bone homeostasis: signaling pathways and therapeutic targets

As a highly dynamic tissue, bone is continuously rebuilt throughout life. Both bone formation by osteoblasts and bone resorption by osteoclasts constitute bone reconstruction homeostasis. The equilibrium of bone homeostasis is governed by many complicated signaling pathways that weave together to form an intricate network. These pathways coordinate the meticulous processes of bone formation and resorption, ensuring the structural integrity and dynamic vitality of the skeletal system. Dysregulation of the bone homeostatic regulatory signaling network contributes to the development and progression of many skeletal diseases. Significantly, imbalanced bone homeostasis further disrupts the signaling network and triggers a cascade reaction that exacerbates disease progression and engenders a deleterious cycle. Here, we summarize the influence of signaling pathways on bone homeostasis, elucidating the interplay and crosstalk among them. Additionally, we review the mechanisms underpinning bone homeostatic imbalances across diverse disease landscapes, highlighting current and prospective therapeutic targets and clinical drugs. We hope that this review will contribute to a holistic understanding of the signaling pathways and molecular mechanisms sustaining bone homeostasis, which are promising to contribute to further research on bone homeostasis and shed light on the development of targeted drugs.

KRAS Mutation Subtypes and Their Association with Other Driver Mutations in Oncogenic Pathways

The KRAS mutation stands out as one of the most influential oncogenic mutations, which directly regulates the hallmark features of cancer and interacts with other cancer-causing driver mutations. However, there remains a lack of precise information on their cooccurrence with mutated variants of KRAS and any correlations between KRAS and other driver mutations. To enquire about this issue, we delved into cBioPortal, TCGA, UALCAN, and Uniport studies. We aimed to unravel the complexity of KRAS and its relationships with other driver mutations. We noticed that G12D and G12V are the prevalent mutated variants of KRAS and coexist with the TP53 mutation in PAAD and CRAD, while G12C and G12V coexist with LUAD. We also noticed similar observations in the case of PIK3CA and APC mutations in CRAD. At the transcript level, a positive correlation exists between KRAS and PIK3CA and between APC and KRAS in CRAD. The existence of the co-mutation of KRAS and other driver mutations could influence the signaling pathway in the neoplastic transformation. Moreover, it has immense prognostic and predictive implications, which could help in better therapeutic management to treat cancer.

Computational design and validation of effective siRNAs to silence oncogenic KRAS

Oncogenic KRAS mutations drive cancer progression in lung, colon, breast, and pancreatic ductal adenocarcinomas. Apart from the current strategies, such as KRAS upstream inhibitors, downstream effector inhibitors, interaction inhibitors, cell cycle inhibitors, and direct KRAS inhibitors, against KRAS-mutated cancers, the therapeutic small interfering RNAs (siRNAs) represent a promising alternative strategy that directly binds with the target mRNA and inhibits protein translation via mRNA degradation. Here, in the present study, we utilized various in silico approaches to design potential siRNA candidates against KRAS mRNA. We have predicted nearly 17 siRNAs against the KRAS mRNA, and further through various criteria, such as U, R, and A rules, GC%, secondary structure formation, mRNA-siRNA duplex stability, Tm (Cp), Tm (Conc), and inhibition efficiency, they have been filtered into 4 potential siRNAs namely siRNA8, siRNA11, siRNA12, and siRNA17. Further, the molecular docking analysis revealed that the siRNA8, siRNA11, siRNA12, and siRNA17 showed higher negative binding energies, such as - 379.13 kcal/mol, - 360.19 kcal/mol, - 288.47 kcal/mol, and - 329.76 kcal/mol, toward the human Argonaute2 protein (hAgo2) respectively. In addition, the normal mode analysis of the hAgo2-siRNAs complexes indicates the structural changes and deformation of the hAgo2 protein upon the binding of siRNA molecules in the dynamic environment which suggests that these siRNAs could be effective. Finally, we conclude that these 4 siRNAs have therapeutic potential against KRAS mRNA and also have to be studied in vitro and in vivo to evaluate their specificity toward mutant KRAS (not degrading wild-type KRAS). Also, the current challenges in the use of siRNA therapeutics could be overcome by the emerging siRNA delivery methods, such as Antibody-siRNA conjugates (ARCs) and Gelatin-Antibody Delivery System (GADS), in the near future and these siRNAs could be employed as potential therapeutic agents against KRAS-mutated cancers.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03767-w.

Potential biomarkers uncovered by bioinformatics analysis in sotorasib resistant-pancreatic ductal adenocarcinoma

Background

Mutant KRAS-induced tumorigenesis is prevalent in lung, colon, and pancreatic ductal adenocarcinomas. For the past 3 decades, KRAS mutants seem undruggable due to their high-affinity GTP-binding pocket and smooth surface. Structure-based drug design helped in the design and development of first-in-class KRAS G12C inhibitor sotorasib (AMG 510) which was then approved by the FDA. Recent reports state that AMG 510 is becoming resistant in non-small-cell lung cancer (NSCLC), pancreatic ductal adenocarcinoma (PDAC), and lung adenocarcinoma patients, and the crucial drivers involved in this resistance mechanism are unknown.

Methods

In recent years, RNA-sequencing (RNA-seq) data analysis has become a functional tool for profiling gene expression. The present study was designed to find the crucial biomarkers involved in the sotorasib (AMG 510) resistance in KRAS G12C-mutant MIA-PaCa2 cell pancreatic ductal adenocarcinoma cells. Initially, the GSE dataset was retrieved from NCBI GEO, pre-processed, and then subjected to differentially expressed gene (DEG) analysis using the limma package. Then the identified DEGs were subjected to protein-protein interaction (PPI) using the STRING database, followed by cluster analysis and hub gene analysis, which resulted in the identification of probable markers.

Results

Furthermore, the enrichment and survival analysis revealed that the small unit ribosomal protein (RP) RPS3 is the crucial biomarker of the AMG 510 resistance in KRAS G12C-mutant MIA-PaCa2 cell pancreatic ductal adenocarcinoma cells.

Conclusion

Finally, we conclude that RPS3 is a crucial biomarker in sotorasib resistance which evades apoptosis by MDM2/4 interaction. We also suggest that the combinatorial treatment of sotorasib and RNA polymerase I machinery inhibitors could be a possible strategy to overcome resistance and should be studied in in vitro and in vivo settings in near future.

Structural basis of the oncogenic KRAS mutant and GJ101 complex

KRAS mutations occur in a quarter of all human cancers. When activated in its GTP-bound form, RAS stimulates diverse cellular systems, such as cell division, differentiation, growth, and apoptosis through the activations of various signaling pathways, which include mitogen-activated protein kinase (MAPK), phosphoinositide 3 kinases (PI3K), and RAL-GEFs pathways. We found that GJ101 (⁶⁵LYDVA⁶⁹) binds directly to the KRAS mutant (G12V) and showed tumor-suppressive activity. In addition, the GJ101 peptide inhibited KRAS mutant as determined by a [α-³²P] guanosine triphosphate (GTP) binding assay and suppressed pancreatic cell line in a cell proliferation assay. Herein, the complex structure of KRAS and GJ101 was clarified by X-ray crystallography. Isothermal titration calorimetry showed that GJ101 binds highly with KRAS mutant and the complex structure of KRAS G12V^.GJ101 complex presented that the residue of Q61 directly interacted with L65 of GJ101. Overall, the results suggest GJ101 be considered a developmental starting point for KRAS G12V inhibitor.

Single-cell RNA-sequencing dissects cellular heterogeneity and identifies two tumor-suppressing immune cell subclusters in HPV-related cervical adenosquamous carcinoma

The intratumor heterogeneity of human papillomavirus (HPV)-related cervical cancer remains poorly defined. We performed single-cell RNA sequencing on 18 046 individual cells derived from two HPV-related cervical adenosquamous carcinoma samples to analyze the transcriptional heterogeneity of both epithelial and immune constituents, identifying seven epithelial (Epi1-7) and 11 immune subclusters. Based on expression of known cervical cancer markers, Epi1-2 primarily displayed features of adenocarcinoma, whereas Epi3-6 were instead characterized by features of squamous carcinoma. Our analyses also revealed that hypoxia and Kirsten rat sarcoma viral oncogene signaling were highly represented within Epi1; metabolic pathways mediating glycolysis and oxidative phosphorylation were enriched in Epi2-4; while Epi5 was enriched in p53 pathway components and features of epithelial-mesenchymal transition. Moreover, CD8⁺ FGFBP2⁺ T cells and FGFBP2⁺ natural killer cells were found to display high levels of cytotoxic effectors (GZMA, GZMB, GNLY, and PRF1) and low levels of inhibitory markers (PDCD1, TIGIT, and CTLA4), such that tumor infiltration by these populations was positively associated with survival in a cohort of n = 165 patients with HPV-related cervical cancer from The Cancer Genome Atlas database (p = 0.017 and 0.014, respectively). These results shed new light on the intratumor heterogeneity of HPV-related cervical adenosquamous carcinoma, which will help to refine diagnostic and treatment approaches.

At the Research Frontiers of Small GTPases

Small GTPases act as molecular switches in regulating a myriad of cellular signaling, cytoskeletal dynamics, vesicular trafficking, and membrane/organelle transport processes. Here, I provide an editorial overview of papers collected in this Special Issue on the "Regulation and Function of Small GTPases 2.0".

Rerouting the drug response: Overcoming metabolic adaptation in KRAS-mutant cancers

Mutations in guanosine triphosphatase KRAS are common in lung, colorectal, and pancreatic cancers. The constitutive activity of mutant KRAS and its downstream signaling pathways induces metabolic rewiring in tumor cells that can promote resistance to existing therapeutics. In this review, we discuss the metabolic pathways that are altered in response to treatment and those that can, in turn, alter treatment efficacy, as well as the role of metabolism in the tumor microenvironment (TME) in dictating the therapeutic response in KRAS-driven cancers. We highlight metabolic targets that may provide clinical opportunities to overcome therapeutic resistance and improve survival in patients with these aggressive cancers.

Cryptotanshinone suppresses tumorigenesis by inhibiting lipogenesis and promoting reactive oxygen species production in KRAS‑activated pancreatic cancer cells

Pyruvate dehydrogenase kinase 4 (PDK4) is an important regulator of energy metabolism. Previously, knockdown of PDK4 by specific small interfering RNAs (siRNAs) have been shown to suppress the expression of Kirsten rat sarcoma viral oncogene homolog (KRAS) and the growth of lung and colorectal cancer cells, indicating that PDK4 is an attractive target of cancer therapy by altering energy metabolism. The authors previously reported that a novel small molecule, cryptotanshinone (CPT), which inhibits PDK4 activity, suppresses the in vitro three-dimensional (3D)-spheroid formation and in vivo tumorigenesis of KRAS-activated human pancreatic and colorectal cancer cells. The present study investigated the molecular mechanism of CPT-induced tumor suppression via alteration of glutamine and lipid metabolism in human pancreatic and colon cancer cell lines with mutant and wild-type KRAS. The antitumor effect of CPT was more pronounced in the cancer cells containing mutant KRAS compared with those containing wild-type KRAS. CPT treatment decreased glutamine and lipid metabolism, affected redox regulation and increased reactive oxygen species (ROS) production in the pancreatic cancer cell line MIAPaCa-2 containing mutant KRAS. Suppression of activated KRAS by specific siRNAs decreased 3D-spheroid formation, the expression of acetyl-CoA carboxylase 1 and fatty acid synthase (FASN) and lipid synthesis. The suppression also reduced glutathione-SH/glutathione disulfide and increased the production of ROS. Knockdown of FASN suppressed lipid synthesis in MIAPaCa-2 cells, partially promoted ROS production and mildly suppressed 3D-spheroid formation. These results indicated that CPT reduced tumorigenesis by inhibiting lipid metabolism and promoting ROS production in a mutant KRAS-dependent manner. This PDK4 inhibitor could serve as a novel therapeutic drug for KRAS-driven intractable cancers via alteration of cell metabolism.

Aberrant activation of KRAS in mouse theca-interstitial cells results in female infertility

KRAS plays critical roles in regulating a range of normal cellular events as well as pathological processes in many tissues mediated through a variety of signaling pathways, including ERK1/2 and AKT signaling, in a cell-, context- and development-dependent manner. The in vivo function of KRAS and its downstream targets in gonadal steroidogenic cells for the development and homeostasis of reproductive functions remain to be determined. To understand the functions of KRAS signaling in gonadal theca and interstitial cells, we generated a Kras mutant (tKrasMT) mouse line that selectively expressed a constitutively active KrasG12D in these cells. KrasG12D expression in ovarian theca cells did not block follicle development to the preovulatory stage. However, tKrasMT females failed to ovulate and thus were infertile. The phosphorylated ERK1/2 and forkhead box O1 (FOXO1) and total FOXO1 protein levels were markedly reduced in tKrasMT theca cells. KrasG12D expression in theca cells also curtailed the phosphorylation of ERK1/2 and altered the expression of several ovulation-related genes in gonadotropin-primed granulosa cells. To uncover downstream targets of KRAS/FOXO1 signaling in theca cells, we found that the expression of bone morphogenic protein 7 (Bmp7), a theca-specific factor involved in ovulation, was significantly elevated in tKrasMT theca cells. Chromosome immunoprecipitation assays demonstrated that FOXO1 interacted with the Bmp7 promoter containing forkhead response elements and that the binding activity was attenuated in tKrasMT theca cells. Moreover, Foxo1 knockdown caused an elevation, whereas Foxo1 overexpression resulted in an inhibition of Bmp7 expression, suggesting that KRAS signaling regulates FOXO1 protein levels to control Bmp7 expression in theca cells. Thus, the anovulation phenotype observed in tKrasMT mice may be attributed to aberrant KRAS/FOXO1/BMP7 signaling in theca cells. Our work provides the first in vivo evidence that maintaining normal KRAS activity in ovarian theca cells is crucial for ovulation and female fertility.

Genetic landscape of pancreatic adenocarcinoma patients: a pilot study from Pakistan

BACKGROUND

Pancreatic adenocarcinoma is one of the most aggressive malignancies with extremely low survival rate. Studies have shown that the exploration of key genes can provide a basis for targeted treatment of these patients. The genomic architecture of the Pakistani pancreatic adenocarcinoma patients remains unexplored. Keeping the scenario in place, the current study aims to analyse 88 cancer related genes in Pakistani pancreatic adenocarcinoma patients in order to elucidate candidate gene(s) for targeted molecular therapy.

METHODS AND RESULTS

A total 18 patients were included in the study initially and FFPE tumor samples were obtained. After confirmation of diagnosis and appropriate tumor content, DNA was extracted. Based on the quality and quantity of the extracted DNA, six pancreatic adenocarcinoma tumor samples were selected. Following to this, all the samples were subjected to targeted sequencing (Axen Cancer Panel 1). Variant detection was done and clinical significance of identified variants was assessed using ClinVar database. Targeted sequencing of tumor samples revealed a total of 29 alterations in the coding region of various genes. Among these five pathogenic variants were found in KRAS, BRCA1, TP53 and APC genes.

CONCLUSION

This is the first study that explores genes involved in pancreatic adenocarcinoma from the Pakistani population. Results obtained from the pilot study can guide us about the key genetic players in the Pakistani pancreatic adenocarcinoma population. This can lead to our better understanding of the molecular targeted therapies for these patients and designing future researches on larger sample size.