Exploring subclass-specific therapeutic agents for hepatocellular carcinoma by informatics-guided drug screen

In silico approaches for drug repurposing in oncology: a scoping review

Introduction: Cancer refers to a group of diseases characterized by the uncontrolled growth and spread of abnormal cells in the body. Due to its complexity, it has been hard to find an ideal medicine to treat all cancer types, although there is an urgent need for it. However, the cost of developing a new drug is high and time-consuming. In this sense, drug repurposing (DR) can hasten drug discovery by giving existing drugs new disease indications. Many computational methods have been applied to achieve DR, but just a few have succeeded. Therefore, this review aims to show in silico DR approaches and the gap between these strategies and their ultimate application in oncology. Methods: The scoping review was conducted according to the Arksey and O'Malley framework and the Joanna Briggs Institute recommendations. Relevant studies were identified through electronic searching of PubMed/MEDLINE, Embase, Scopus, and Web of Science databases, as well as the grey literature. We included peer-reviewed research articles involving in silico strategies applied to drug repurposing in oncology, published between 1 January 2003, and 31 December 2021. Results: We identified 238 studies for inclusion in the review. Most studies revealed that the United States, India, China, South Korea, and Italy are top publishers. Regarding cancer types, breast cancer, lymphomas and leukemias, lung, colorectal, and prostate cancer are the top investigated. Additionally, most studies solely used computational methods, and just a few assessed more complex scientific models. Lastly, molecular modeling, which includes molecular docking and molecular dynamics simulations, was the most frequently used method, followed by signature-, Machine Learning-, and network-based strategies. Discussion: DR is a trending opportunity but still demands extensive testing to ensure its safety and efficacy for the new indications. Finally, implementing DR can be challenging due to various factors, including lack of quality data, patient populations, cost, intellectual property issues, market considerations, and regulatory requirements. Despite all the hurdles, DR remains an exciting strategy for identifying new treatments for numerous diseases, including cancer types, and giving patients faster access to new medications.

Multidimensional Transcriptomics Unveils RNF34 as a Prognostic Biomarker and Potential Indicator of Chemotherapy Sensitivity in Wilms' Tumour

Nephroblastoma, colloquially known as Wilms' tumour (WT), is the predominant malignant renal neoplasm arising in the paediatric population. Modern therapeutic approaches for WT incorporate a synergistic combination of surgical intervention, radiotherapy, and chemotherapy, which substantially ameliorate the overall patient survival rate. Despite this, the optimal sequence of chemotherapy and surgical intervention remains a matter of contention, with each strategy presenting its own strengths and weaknesses that could influence clinical decision-making. To make some headway on this clinical dilemma, we deployed a multidimensional transcriptomics integration approach by analysing bulk RNA sequencing data with 136 samples, as well as single-nucleus RNA sequencing (snRNA-seq) and paired spatial transcriptome sequencing (stRNA) data from 32 WT specimens. Our findings identified a distinct elevation of RNF34 expression within WT samples, which correlated with unfavourable prognostic outcomes. Leveraging the Genomics of Drug Sensitivity in Cancer (GDSC), we simultaneously revealed that patients with high expression of RNF34 have higher sensitivity to commonly used chemotherapy drugs for WT. Furthermore, our analysis of snRNA and stRNA data unveiled a reduced proportion of RNF34 expression in neoplastic cells after chemotherapy. Moreover, stRNA data delineated a significant association between a higher proportion of RNF34 expression in cancer cells and adverse features such as anaplastic histology and tumour recurrence. Intriguingly, we also observed a close association between elevated RNF34 expression and a characteristic exhausted tumour immune microenvironment. Collectively, our findings underscore the pivotal role of RNF34 in the prognostic prediction potential and treatment sensitivity of WT. This comprehensive analysis can potentially inform and refine clinical decision-making for WT patients and guide future studies towards the development of optimized, rational therapeutic strategies.

An integrated framework for prognosis prediction and drug response modeling in colorectal liver metastasis drug discovery

BACKGROUND

Colorectal cancer (CRC) is the third most prevalent cancer globally, and liver metastasis (CRLM) is the primary cause of death. Hence, it is essential to discover novel prognostic biomarkers and therapeutic drugs for CRLM.

METHODS

This study developed two liver metastasis-associated prognostic signatures based on differentially expressed genes (DEGs) in CRLM. Additionally, we employed an interpretable deep learning model utilizing drug sensitivity databases to identify potential therapeutic drugs for high-risk CRLM patients. Subsequently, in vitro and in vivo experiments were performed to verify the efficacy of these compounds.

RESULTS

These two prognostic models exhibited superior performance compared to previously reported ones. Obatoclax, a BCL-2 inhibitor, showed significant differential responses between high and low risk groups classified by prognostic models, and demonstrated remarkable effectiveness in both Transwell assay and CT26 colorectal liver metastasis mouse model.

CONCLUSIONS

This study highlights the significance of developing specialized prognostication approaches and investigating effective therapeutic drugs for patients with CRLM. The application of a deep learning drug response model provides a new drug discovery strategy for translational medicine in precision oncology.

Immune perturbation network identifies an EMT subtype with chromosomal instability and tumor immune-desert microenvironment

Most gastric cancer (GC) subtypes are identified through transcriptional profiling overlooking dynamic changes and interactions in gene expression. Based on the background network of global immune genes, we constructed sample-specific edge-perturbation matrices and identified four molecular network subtypes of GC (MNG). MNG-1 displayed the best prognosis and vigorous cell cycle activity. MNG-2 was enriched by immune-hot phenotype with the potential for immunotherapy response. MNG-3 and MNG-4 were identified with epithelial-mesenchymal transition (EMT) peculiarity and worse prognosis, termed EMT subtypes. MNG-3 was characterized by low mutational burden and stromal cells and considered a replica of previous subtypes associated with poor prognosis. Notably, MNG-4 was considered a previously undefined subtype with a dismal prognosis, characterized by chromosomal instability and immune-desert microenvironment. This subtype tended to metastasize and was resistant to respond to immunotherapy. Pharmacogenomics analysis showed three therapeutic agents (NVP-BEZ235, LY2606368, and rutin) were potential interventions for MNG-4.

Unraveling molecular networks in thymic epithelial tumors: deciphering the unique signatures

Thymic epithelial tumors (TETs) are a rare and diverse group of neoplasms characterized by distinct molecular signatures. This review delves into the complex molecular networks of TETs, highlighting key aspects such as chromosomal abnormalities, molecular subtypes, aberrant gene mutations and expressions, structural gene rearrangements, and epigenetic changes. Additionally, the influence of the dynamic tumor microenvironment on TET behavior and therapeutic responses is examined. A thorough understanding of these facets elucidates TET pathogenesis, offering avenues for enhancing diagnostic accuracy, refining prognostic assessments, and tailoring targeted therapeutic strategies. Our review underscores the importance of deciphering TETs' unique molecular signatures to advance personalized treatment paradigms and improve patient outcomes. We also discuss future research directions and anticipated challenges in this intriguing field.

Senescence risk score: a multifaceted prognostic tool predicting outcomes, stemness, and immune responses in colorectal cancer

Colorectal cancer (CRC) remains a primary cause of cancer mortality globally, necessitating precise prognostic indicators for effective clinical management. Our study introduces the Senescence Risk Score (SRRS), based on several senescence-related genes (SRGs), a potent prognostic tool designed to measure cellular senescence in CRC. The higher SRRS predicts a poorer prognosis, providing a novel and efficient approach to patient stratification. Notably, we found that SRRS correlates with methylation and mutation variations, and increased immune infiltration in the tumor microenvironment, thus revealing potential therapeutic targets. We also discovered an inverse relationship between SRRS and cell stemness, which could have significant implications for cancer treatment strategies. Utilizing bioinformatics resources and machine learning, we identified LIMK1 and WRN as key genes associated with SRRS, further enhancing its prognostic value. Importantly, the modulation of these genes significantly impacts cellular senescence, proliferation, and stemness in CRC cells. In summary, our development of SRRS offers a powerful tool for CRC prognosis and paves the way for novel therapeutic strategies, underscoring its potential in transforming CRC patient management.

Integrated multiomics analysis and machine learning refine molecular subtypes and prognosis for muscle-invasive urothelial cancer

Muscle-invasive urothelial cancer (MUC), characterized by high aggressiveness and significant heterogeneity, is currently lacking highly precise individualized treatment options. We used a computational pipeline to synthesize multiomics data from MUC patients using 10 clustering algorithms, which were then combined with 10 machine learning algorithms to identify molecular subgroups of high resolution and develop a robust consensus machine learning-driven signature (CMLS). Through multiomics clustering, we identified three cancer subtypes (CSs) that are related to prognosis, with CS2 exhibiting the most favorable prognostic outcome. Subsequent screening enabled identification of 12 hub genes that constitute a CMLS with robust predictive power for prognosis. The low-CMLS group exhibited a more favorable prognosis and greater responsiveness to immunotherapy and was more likely to exhibit the "hot tumor" phenotype. The high-CMLS group had a poor prognosis and lower likelihood of benefitting from immunotherapy, but dasatinib and romidepsin may serve as promising treatments for them. Comprehensive analysis of multiomics data can offer important insights and further refine the molecular classification of MUC. Identification of CMLS represents a valuable tool for early prediction of patient prognosis and for screening potential candidates likely to benefit from immunotherapy, with broad implications for clinical practice.

Prognosis and personalized treatment prediction in lung adenocarcinoma: An in silico and in vitro strategy adopting cuproptosis related lncRNA towards precision oncology

Background: There is a rapid increase in lung adenocarcinomas (LUAD), and studies suggest associations between cuproptosis and the occurrence of various types of tumors. However, it remains unclear whether cuproptosis plays a role in LUAD prognosis. Methods: Dataset of the TCGA-LUAD was treated as training cohort, while validation cohort consisted of the merged datasets of the GSE29013, GSE30219, GSE31210, GSE37745, and GSE50081. Ten studied cuproptosis-related genes (CRG) were used to generated CRG clusters and CRG cluster-related differential expressed gene (CRG-DEG) clusters. The differently expressed lncRNA that with prognosis ability between the CRG-DEG clusters were put into a LASSO regression for cuproptosis-related lncRNA signature (CRLncSig). Kaplan-Meier estimator, Cox model, receiver operating characteristic (ROC), time-dependent AUC (tAUC), principal component analysis (PCA), and nomogram predictor were further deployed to confirm the model's accuracy. We examined the model's connections with other forms of regulated cell death, including apoptosis, necroptosis, pyroptosis, and ferroptosis. The immunotherapy ability of the signature was demonstrated by applying eight mainstream immunoinformatic algorithms, TMB, TIDE, and immune checkpoints. We evaluated the potential drugs for high risk CRLncSig LUADs. Real-time PCR in human LUAD tissues were performed to verify the CRLncSig expression pattern, and the signature's pan-cancer's ability was also assessed. Results: A nine-lncRNA signature, CRLncSig, was built and demonstrated owning prognostic power by applied to the validation cohort. Each of the signature genes was confirmed differentially expressed in the real world by real-time PCR. The CRLncSig correlated with 2,469/3,681 (67.07%) apoptosis-related genes, 13/20 (65.00%) necroptosis-related genes, 35/50 (70.00%) pyroptosis-related genes, and 238/380 (62.63%) ferroptosis-related genes. Immunotherapy analysis suggested that CRLncSig correlated with immune status, and checkpoints, KIR2DL3, IL10, IL2, CD40LG, SELP, BTLA, and CD28, were linked closely to our signature and were potentially suitable for LUAD immunotherapy targets. For those high-risk patients, we found three agents, gemcitabine, daunorubicin, and nobiletin. Finally, we found some of the CRLncSig lncRNAs potentially play a vital role in some types of cancer and need more attention in further studies. Conclusion: The results of this study suggest our cuproptosis-related CRLncSig can help to determine the outcome of LUAD and the effectiveness of immunotherapy, as well as help to better select targets and therapeutic agents.

Gene interaction perturbation network deciphers a high-resolution taxonomy in colorectal cancer

Molecular subtypes of colorectal cancer (CRC) are currently identified via the snapshot transcriptional profiles, largely ignoring the dynamic changes of gene expressions. Conversely, biological networks remain relatively stable irrespective of time and condition. Here, we introduce an individual-specific gene interaction perturbation network-based (GIN) approach and identify six GIN subtypes (GINS1-6) with distinguishing features: (i) GINS1 (proliferative, 24%~34%), elevated proliferative activity, high tumor purity, immune-desert, PIK3CA mutations, and immunotherapeutic resistance; (ii) GINS2 (stromal-rich, 14%~22%), abundant fibroblasts, immune-suppressed, stem-cell-like, SMAD4 mutations, unfavorable prognosis, high potential of recurrence and metastasis, immunotherapeutic resistance, and sensitive to fluorouracil-based chemotherapy; (iii) GINS3 (KRAS-inactivated, 13%~20%), high tumor purity, immune-desert, activation of EGFR and ephrin receptors, chromosomal instability (CIN), fewer KRAS mutations, SMOC1 methylation, immunotherapeutic resistance, and sensitive to cetuximab and bevacizumab; (iv) GINS4 (mixed, 10%~19%), moderate level of stromal and immune activities, transit-amplifying-like, and TMEM106A methylation; (v) GINS5 (immune-activated, 12%~24%), stronger immune activation, plentiful tumor mutation and neoantigen burden, microsatellite instability and high CpG island methylator phenotype, BRAF mutations, favorable prognosis, and sensitive to immunotherapy and PARP inhibitors; (vi) GINS6, (metabolic, 5%~8%), accumulated fatty acids, enterocyte-like, and BMP activity. Overall, the novel high-resolution taxonomy derived from an interactome perspective could facilitate more effective management of CRC patients.

Integrative genomic analysis facilitates precision strategies for glioblastoma treatment

Glioblastoma (GBM) is the most common form of malignant primary brain tumor with a dismal prognosis. Currently, the standard treatments for GBM rarely achieve satisfactory results, which means that current treatments are not individualized and precise enough. In this study, a multiomics-based GBM classification was established and three subclasses (GPA, GPB, and GPC) were identified, which have different molecular features both in bulk samples and at single-cell resolution. A robust GBM poor prognostic signature (GPS) score model was then developed using machine learning method, manifesting an excellent ability to predict the survival of GBM. NVP−BEZ235, GDC−0980, dasatinib and XL765 were ultimately identified to have subclass-specific efficacy targeting patients with a high risk of poor prognosis. Furthermore, the GBM classification and GPS score model could be considered as potential biomarkers for immunotherapy response. In summary, an integrative genomic analysis was conducted to advance individual-based therapies in GBM.

•

A multiomics-based classification of GBM was established

•

Single-cell transcriptomic profiling of GBM subclasses was revealed using Scissor

•

A robust prognostic risk model was developed for GBM by machine learning method

•

Prediction of potential agents based on molecular and prognostic risk stratification

Subclass-specific Prognosis and Treatment Efficacy Inference in Head and Neck Squamous Carcinoma

Exploring the prognostic classification and biomarkers in Head and Neck Squamous Carcinoma (HNSC) is of great clinical significance. We hybridized three prominent strategies to comprehensively characterize the molecular features of HNSC. We constructed a 15-gene signature to predict patients death risk with an average AUC of 0.744 for 1-, 3-, and 5-year on TCGA-HNSC training set, and average AUCs of 0.636, 0.584, 0.755 in GSE65858, GSE-112026, CPTAC-HNSCC datasets, respectively. By combined with NMF clustering and consensus clustering of fraction of tumor immune cell infiltration (ICI) in the tumor microenvironment (TME), we captured a more refined biological characteristics of HNSC, and observed a prognosis heterogeneity in high tumor immunity patients. By matching tumor subset-specific expression signatures to drug-induced cell line expression profiles from large-scale pharmacogenomic databases in the OCTAD workspace, we identified a group of HNSC patients featured with poor prognosis and demonstrated that the individuals in this group are likely to receive increased drug sensitivity to reverse differentially expressed disease signature genes. This trend is especially highlighted among those with higher death risk and tumour immunity.

A survey of optimal strategy for signature-based drug repositioning and an application to liver cancer

Pharmacologic perturbation projects, such as Connectivity Map (CMap) and Library of Integrated Network-based Cellular Signatures (LINCS), have produced many perturbed expression data, providing enormous opportunities for computational therapeutic discovery. However, there is no consensus on which methodologies and parameters are the most optimal to conduct such analysis. Aiming to fill this gap, new benchmarking standards were developed to quantitatively evaluate drug retrieval performance. Investigations of potential factors influencing drug retrieval were conducted based on these standards. As a result, we determined an optimal approach for LINCS data-based therapeutic discovery. With this approach, homoharringtonine (HHT) was identified to be a candidate agent with potential therapeutic and preventive effects on liver cancer. The antitumor and antifibrotic activity of HHT was validated experimentally using subcutaneous xenograft tumor model and carbon tetrachloride (CCL₄)-induced liver fibrosis model, demonstrating the reliability of the prediction results. In summary, our findings will not only impact the future applications of LINCS data but also offer new opportunities for therapeutic intervention of liver cancer.

Molecular classification of hepatocellular carcinoma: prognostic importance and clinical applications

Hepatocellular carcinoma (HCC) is a lethal human malignancy with a very low overall and long-term survival rate. Poor prognostic outcomes are predominantly associated with HCC due to a huge landscape of heterogeneity found in the deadliest disease. However, molecular subtyping of HCC has significantly improved the knowledge of the underlying mechanisms that contribute towards the heterogeneity and progression of the disease. In this review, we have extensively summarized the current information available about molecular classification of HCC. This review can be of great significance for providing the insight information needed for development of novel, efficient and personalized therapeutic options for the treatment of HCC patients globally.

Exploring liver cancer biology through functional genetic screens

As the fourth leading cause of cancer-related death in the world, liver cancer poses a major threat to human health. Although a growing number of therapies have been approved for the treatment of hepatocellular carcinoma in the past few years, most of them only provide a limited survival benefit. Therefore, an urgent need exists to identify novel targetable vulnerabilities and powerful drug combinations for the treatment of liver cancer. The advent of functional genetic screening has contributed to the advancement of liver cancer biology, uncovering many novel genes involved in tumorigenesis and cancer progression in a high-throughput manner. In addition, this unbiased screening platform also provides an efficient tool for the exploration of the mechanisms involved in therapy resistance as well as identifying potential targets for therapy. In this Review, we describe how functional screens can help to deepen our understanding of liver cancer and guide the development of new therapeutic strategies.

The origin of human mutation in light of genomic data

Despite years of active research into the role of DNA repair and replication in mutagenesis, surprisingly little is known about the origin of spontaneous human mutation in the germ line. With the advent of high-throughput sequencing, genome-scale data have revealed statistical properties of mutagenesis in humans. These properties include variation of the mutation rate and spectrum along the genome at different scales in relation to epigenomic features and dependency on parental age. Moreover, mutations originated in mothers are less frequent than mutations originated in fathers and have a distinct genomic distribution. Statistical analyses that interpret these patterns in the context of known biochemistry can provide mechanistic models of mutagenesis in humans.

Mapping the landscape of synthetic lethal interactions in liver cancer

Almost all the current therapies against liver cancer are based on the "one size fits all" principle and offer only limited survival benefit. Fortunately, synthetic lethality (SL) may provide an alternate route towards individualized therapy in liver cancer. The concept that simultaneous losses of two genes are lethal to a cell while a single loss is non-lethal can be utilized to selectively eliminate tumors with genetic aberrations. Methods: To infer liver cancer-specific SL interactions, we propose a computational pipeline termed SiLi (statistical inference-based synthetic lethality identification) that incorporates five inference procedures. Based on large-scale sequencing datasets, SiLi analysis was performed to identify SL interactions in liver cancer. Results: By SiLi analysis, a total of 272 SL pairs were discerned, which included 209 unique target candidates. Among these, polo-like kinase 1 (PLK1) was considered to have considerable therapeutic potential. Further computational and experimental validation of the SL pair TP53-PLK1 demonstrated that inhibition of PLK1 could be a novel therapeutic strategy specifically targeting those patients with TP53-mutant liver tumors. Conclusions: In this study, we report a comprehensive analysis of synthetic lethal interactions of liver cancer. Our findings may open new possibilities for patient-tailored therapeutic interventions in liver cancer.