How can we kill cancer cells: Insights from the computational models of apoptosis

Disulfideptosis-associated lncRNAs reveal features of prognostic, immune escape, tumor mutation, and tumor malignant progression in renal clear cell carcinoma

PURPOSE

Investigating the role of lncRNAs associated with the latest cell death mode (Disulfideptosis) in renal clear cell carcinoma, as well as their correlation with tumor prognosis, immune escape, immune checkpoints, tumor mutational burden, and malignant tumor progression. Searching for potential biomarkers and targets for renal clear cell carcinoma.

METHODS

Downloaded the expression profile data and clinical data of 533 cases of renal clear cell carcinoma from the TCGA database, and randomly divided them into a test set (267 cases) and a validation set (266 cases). Based on previous research, 13 genes associated with Disulfideptosis were obtained. Using R software, lncRNAs with a differential expression that is related to the prognosis of renal clear cell carcinoma and associated with Disulfideptosis were screened out. After univariate Cox regression analysis, Lasso regression analysis, and multivariate Cox regression analysis, lncRNAs with independent predictive ability were obtained. A predictive risk model was established based on the risk scores. Verification was carried out between the obtained high-risk and low-risk groups and their subgroups (including Age, Gender, tumor mutational burden (TMB), tumor grading, and staging). Subsequently, a nomogram was established, and a calibration curve was generated for verification. Performed GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) functional enrichment analyses. Downloaded the values of Tumor Immune Dysfunction and Exclusion (TIDE) for all samples and calculated the difference between the high and low-risk groups. Selected human renal tumor cell lines (786-O, OS-RC-2, A-498, ACHN) and human renal cortex proximal tubule epithelial cell line (HK-2). The RNA expression levels of the above lncRNAs in each cell line were analyzed using RT-qPCR (Real-time Quantitative PCR Detecting System). Used siRNA (small interfering RNA) to knock down FAM225B in 786-O and OS-RC-2 cell lines, and then performed in vitro cell experiments to validate the functional characteristics of FAM225B.

RESULTS

Our constructed predictive model includes 5 lncRNAs with an independent predictive ability (FAM225B, ZNF503-AS1, SPINT1-AS1, WWC2-AS2, LINC01338), which can effectively distinguish between patients in high and low-risk groups and their subgroups. The 1, 3, and 5-year AUC (Area Under the ROC Curve) values of the established nomogram are 0.756, 0.752, and 0.781, respectively. The 5-year AUC value is higher compared to other clinical characteristics (Age: 0.598, Gender: 0.488, Grade: 0.680, Stage: 0.717). After the knockdown of FAM225B, the proliferation, migration, and invasion abilities of renal cancer cell lines OS-RC-2 and 786-O all decreased.

CONCLUSION

We have constructed and validated a prognostic model based on Disulfideptosis-associated lncRNAs. This model can effectively predict the high or low risk of patient prognosis and can distinguish the tumor cell mutational burden and immune escape capabilities among high-risk and low-risk patients. This predictive model can serve as an independent prognostic factor for renal clear cell carcinoma, providing a new direction for personalized treatment of patients with renal clear cell carcinoma.

Selective Cytotoxicity of 1,3,4-Thiadiazolium Mesoionic Derivatives on Hepatocarcinoma Cells (HepG2)

In this work, we evaluated the cytotoxicity of mesoionic 4-phenyl-5-(2-Y, 4-X or 4-X-cinnamoyl)-1,3,4-thiadiazolium-2-phenylamine chloride derivatives (MI-J: X=OH, Y=H; MI-D: X=NO2, Y=H; MI-4F: X=F, Y=H; MI-2,4diF: X=Y=F) on human hepatocellular carcinoma (HepG2), and non-tumor cells (rat hepatocytes) for comparison. MI-J, M-4F and MI-2,4diF reduced HepG2 viability by ~ 50% at 25 μM after 24-h treatment, whereas MI-D required a 50 μM concentration, as shown by 3-(4,5-dimethythiazol-2-yl)-2,5-diphenyltetrazolium bromide assays. The cytotoxicity was confirmed with lactate dehydrogenase assay, of which activity was increased by 55, 24 and 16% for MI-J, MI-4F and MI-2,4diF respectively (at 25 μM after 24 h). To identify the death pathway related to cytotoxicity, the HepG2 cells treated by mesoionic compounds were labeled with both annexin V and PI, and analyzed by flow cytometry. All compounds increased the number of doubly-stained cells at 25 μM after 24 h: by 76% for MI-J, 25% for MI-4F and MI-2,4diF, and 11% for MI-D. It was also verified that increased DNA fragmentation occurred upon MI-J, MI-4F and MI-2,4diF treatments (by 12%, 9% and 8%, respectively, at 25 μM after 24 h). These compounds were only weakly, or not at all, transported by the main multidrug transporters, P-glycoprotein, ABCG2 and MRP1, and were able to slightly inhibit their drug-transport activity. It may be concluded that 1,3,4-thiadiazolium compounds, especially the hydroxy derivative MI-J, constitute promising candidates for future investigations on in-vivo treatment of hepatocellular carcinoma.

Identification of proapoptopic, anti-inflammatory, anti- proliferative, anti-invasive and anti-angiogenic targets of essential oils in cardamom by dual reverse virtual screening and binding pose analysis

BACKGROUND

Cardamom (Elettaria cardamom), also known as "Queen of Spices", has been traditionally used as a culinary ingredient due to its pleasant aroma and taste. In addition to this role, studies on cardamom have demonstrated cancer chemopreventive potential in in vitro and in vivo systems. Nevertheless, the precise poly-pharmacological nature of naturally occurring chemo-preventive compounds in cardamom has still not been fully demystified.

METHODS

In this study, an effort has been made to identify the proapoptopic, anti-inflammatory, anti-proliferative, anti-invasive and anti-angiogenic targets of Cardamom's bioactive principles (eucalyptol, alpha-pinene, beta-pinene, d-limonene and geraniol) by employing a dual reverse virtual screening protocol. Experimentally proven target information of the bioactive principles was annotated from bioassay databases and compared with the virtually screened set of targets to evaluate the reliability of the computational identification. To study the molecular interaction pattern of the anti-tumor action, molecular docking simulation was performed with Auto Dock Pyrx. Interaction studies of binding pose of eucalyptol with Caspase 3 were conducted to obtain an insight into the interacting amino acids and their inter-molecular bondings.

RESULTS

A prioritized list of target proteins associated with multiple forms of cancer and ranked by their Fit Score (Pharm Mapper) and descending 3D score (Reverse Screen 3D) were obtained from the two independent inverse screening platforms. Molecular docking studies exploring the bioactive principle targeted action revealed that H- bonds and electrostatic interactions forms the chief contributing factor in inter-molecular interactions associated with anti-tumor activity. Eucalyptol binds to the Caspase 3 with a specific framework that is well-suited for nucleophilic attacks by polar residues inside the Caspase 3 catalytic site.

CONCLUSION

This study revealed vital information about the poly-pharmacological anti-tumor mode-of-action of essential oils in cardamom. In addition, a probabilistic set of anti-tumor targets for cardamom was generated, which can be further confirmed by in vivo and in vitro experiments.

Modeling heterogeneous responsiveness of intrinsic apoptosis pathway

BACKGROUND

Apoptosis is a cell suicide mechanism that enables multicellular organisms to maintain homeostasis and to eliminate individual cells that threaten the organism's survival. Dependent on the type of stimulus, apoptosis can be propagated by extrinsic pathway or intrinsic pathway. The comprehensive understanding of the molecular mechanism of apoptotic signaling allows for development of mathematical models, aiming to elucidate dynamical and systems properties of apoptotic signaling networks. There have been extensive efforts in modeling deterministic apoptosis network accounting for average behavior of a population of cells. Cellular networks, however, are inherently stochastic and significant cell-to-cell variability in apoptosis response has been observed at single cell level.

RESULTS

To address the inevitable randomness in the intrinsic apoptosis mechanism, we develop a theoretical and computational modeling framework of intrinsic apoptosis pathway at single-cell level, accounting for both deterministic and stochastic behavior. Our deterministic model, adapted from the well-accepted Fussenegger model, shows that an additional positive feedback between the executioner caspase and the initiator caspase plays a fundamental role in yielding the desired property of bistability. We then examine the impact of intrinsic fluctuations of biochemical reactions, viewed as intrinsic noise, and natural variation of protein concentrations, viewed as extrinsic noise, on behavior of the intrinsic apoptosis network. Histograms of the steady-state output at varying input levels show that the intrinsic noise could elicit a wider region of bistability over that of the deterministic model. However, the system stochasticity due to intrinsic fluctuations, such as the noise of steady-state response and the randomness of response delay, shows that the intrinsic noise in general is insufficient to produce significant cell-to-cell variations at physiologically relevant level of molecular numbers. Furthermore, the extrinsic noise represented by random variations of two key apoptotic proteins, namely Cytochrome C and inhibitor of apoptosis proteins (IAP), is modeled separately or in combination with intrinsic noise. The resultant stochasticity in the timing of intrinsic apoptosis response shows that the fluctuating protein variations can induce cell-to-cell stochastic variability at a quantitative level agreeing with experiments. Finally, simulations illustrate that the mean abundance of fluctuating IAP protein is positively correlated with the degree of cellular stochasticity of the intrinsic apoptosis pathway.

CONCLUSIONS

Our theoretical and computational study shows that the pronounced non-genetic heterogeneity in intrinsic apoptosis responses among individual cells plausibly arises from extrinsic rather than intrinsic origin of fluctuations. In addition, it predicts that the IAP protein could serve as a potential therapeutic target for suppression of the cell-to-cell variation in the intrinsic apoptosis responsiveness.

Monte carlo study elucidates the type 1/type 2 choice in apoptotic death signaling in healthy and cancer cells

Apoptotic cell death is coordinated through two distinct (type 1 and type 2) intracellular signaling pathways. How the type 1/type 2 choice is made remains a central problem in the biology of apoptosis and has implications for apoptosis related diseases and therapy. We study the problem of type 1/type 2 choice in silico utilizing a kinetic Monte Carlo model of cell death signaling. Our results show that the type 1/type 2 choice is linked to deterministic versus stochastic cell death activation, elucidating a unique regulatory control of the apoptotic pathways. Consistent with previous findings, our results indicate that caspase 8 activation level is a key regulator of the choice between deterministic type 1 and stochastic type 2 pathways, irrespective of cell types. Expression levels of signaling molecules downstream also regulate the type 1/type 2 choice. A simplified model of DISC clustering elucidates the mechanism of increased active caspase 8 generation and type 1 activation in cancer cells having increased sensitivity to death receptor activation. We demonstrate that rapid deterministic activation of the type 1 pathway can selectively target such cancer cells, especially if XIAP is also inhibited; while inherent cell-to-cell variability would allow normal cells stay protected.