Targeted drug combination therapy design based on driver genes

Revealing cancer driver genes through integrative transcriptomic and epigenomic analyses with Moonlight

Cancer involves dynamic changes caused by (epi)genetic alterations such as mutations or abnormal DNA methylation patterns which occur in cancer driver genes. These driver genes are divided into oncogenes and tumor suppressors depending on their function and mechanism of action. Discovering driver genes in different cancer (sub)types is important not only for increasing current understanding of carcinogenesis but also from prognostic and therapeutic perspectives. We have previously developed a framework called Moonlight which uses a systems biology multi-omics approach for prediction of driver genes. Here, we present an important development in Moonlight2 by incorporating a DNA methylation layer which provides epigenetic evidence for deregulated expression profiles of driver genes. To this end, we present a novel functionality called Gene Methylation Analysis (GMA) which investigates abnormal DNA methylation patterns to predict driver genes. This is achieved by integrating the tool EpiMix which is designed to detect such aberrant DNA methylation patterns in a cohort of patients and further couples these patterns with gene expression changes. To showcase GMA, we applied it to three cancer (sub)types (basal-like breast cancer, lung adenocarcinoma, and thyroid carcinoma) where we discovered 33, 190, and 263 epigenetically driven genes, respectively. A subset of these driver genes had prognostic effects with expression levels significantly affecting survival of the patients. Moreover, a subset of the driver genes demonstrated therapeutic potential as drug targets. This study provides a framework for exploring the driving forces behind cancer and provides novel insights into the landscape of three cancer sub(types) by integrating gene expression and methylation data.

A co-formulation of interferons alpha2b and gamma distinctively targets cell cycle in the glioblastoma-derived cell line U-87MG

BACKGROUND

HeberFERON is a co-formulation of α2b and γ interferons, based on their synergism, which has shown its clinical superiority over individual interferons in basal cell carcinomas. In glioblastoma (GBM), HeberFERON has displayed promising preclinical and clinical results. This led us to design a microarray experiment aimed at identifying the molecular mechanisms involved in the distinctive effect of HeberFERON compared to the individual interferons in U-87MG model.

METHODS

Transcriptional expression profiling including a control (untreated) and three groups receiving α2b-interferon, γ-interferon and HeberFERON was performed using an Illumina HT-12 microarray platform. Unsupervised methods for gene and sample grouping, identification of differentially expressed genes, functional enrichment and network analysis computational biology methods were applied to identify distinctive transcription patterns of HeberFERON. Validation of most representative genes was performed by qPCR. For the cell cycle analysis of cells treated with HeberFERON for 24 h, 48 and 72 h we used flow cytometry.

RESULTS

The three treatments show different behavior based on the gene expression profiles. The enrichment analysis identified several mitotic cell cycle related events, in particular from prometaphase to anaphase, which are exclusively targeted by HeberFERON. The FOXM1 transcription factor network that is involved in several cell cycle phases and is highly expressed in GBMs, is significantly down regulated. Flow cytometry experiments corroborated the action of HeberFERON on the cell cycle in a dose and time dependent manner with a clear cellular arrest as of 24 h post-treatment. Despite the fact that p53 was not down-regulated, several genes involved in its regulatory activity were functionally enriched. Network analysis also revealed a strong relationship of p53 with genes targeted by HeberFERON. We propose a mechanistic model to explain this distinctive action, based on the simultaneous activation of PKR and ATF3, p53 phosphorylation changes, as well as its reduced MDM2 mediated ubiquitination and export from the nucleus to the cytoplasm. PLK1, AURKB, BIRC5 and CCNB1 genes, all regulated by FOXM1, also play central roles in this model. These and other interactions could explain a G2/M arrest and the effect of HeberFERON on the proliferation of U-87MG.

CONCLUSIONS

We proposed molecular mechanisms underlying the distinctive behavior of HeberFERON compared to the treatments with the individual interferons in U-87MG model, where cell cycle related events were highly relevant.

Cancer fitness genes: emerging therapeutic targets for metastasis

Development of cancer therapeutics has traditionally focused on targeting driver oncogenes. Such an approach is limited by toxicity to normal tissues and treatment resistance. A class of 'cancer fitness genes' with crucial roles in metastasis have been identified. Elevated or altered activities of these genes do not directly cause cancer; instead, they relieve the stresses that tumor cells encounter and help them adapt to a changing microenvironment, thus facilitating tumor progression and metastasis. Importantly, as normal cells do not experience high levels of stress under physiological conditions, targeting cancer fitness genes is less likely to cause toxicity to noncancerous tissues. Here, we summarize the key features and function of cancer fitness genes and discuss their therapeutic potential.

Exploiting the Molecular Basis of Oesophageal Cancer for Targeted Therapies and Biomarkers for Drug Response: Guiding Clinical Decision-Making

Worldwide, oesophageal cancer is the sixth leading cause of deaths related to cancer and represents a major health concern. Sub-Saharan Africa is one of the regions of the world with the highest incidence and mortality rates for oesophageal cancer and most of the cases of oesophageal cancer in this region are oesophageal squamous cell carcinoma (OSCC). The development and progression of OSCC is characterized by genomic changes which can be utilized as diagnostic or prognostic markers. These include changes in the expression of various genes involved in signaling pathways that regulate pathways that regulate processes that are related to the hallmarks of cancer, changes in the tumor mutational burden, changes in alternate splicing and changes in the expression of non-coding RNAs such as miRNA. These genomic changes give rise to characteristic profiles of altered proteins, transcriptomes, spliceosomes and genomes which can be used in clinical applications to monitor specific disease related parameters. Some of these profiles are characteristic of more aggressive forms of cancer or are indicative of treatment resistance or tumors that will be difficult to treat or require more specialized specific treatments. In Sub-Saharan region of Africa there is a high incidence of viral infections such as HPV and HIV, which are both risk factors for OSCC. The genomic changes that occur due to these infections can serve as diagnostic markers for OSCC related to viral infection. Clinically this is an important distinction as it influences treatment as well as disease progression and treatment monitoring practices. This underlines the importance of the characterization of the molecular landscape of OSCC in order to provide the best treatment, care, diagnosis and screening options for the management of OSCC.

Current progress and novel strategies that target CDK12 for drug discovery

CDK12 is a cyclin-dependent kinase that plays critical roles in DNA replication, transcription, mRNA splicing, and DNA damage repair. CDK12 genomic changes, including mutation, amplification, deletion, and fusion, lead to various cancers, such as colorectal cancer, gastric cancer, and ovarian cancer. An increasing number of CDK12 inhibitors have been reported since CDK12 was identified as a biomarker and cancer therapeutic target. A major challenge lies in that CDK12 and CDK13 share highly similar sequences, which leads to great difficulties in the development of highly selective CDK12 inhibitors. In recent years, great efforts were made in developing selective CDK12 blockers. Techniques including PROTAC and molecular glue degraders were also applied to facilitate their development. Also, the drug combination strategy of CDK12 small molecule inhibitors were studied. This review discusses the latest studies on CDK12 inhibitors and analyzes their structure-activity relationships, shedding light on their further development.

Anticancer drug resistance: An update and perspective

Driver mutations promote initiation and progression of cancer. Pharmacological treatment can inhibit the action of the mutant protein; however, drug resistance almost invariably emerges. Multiple studies revealed that cancer drug resistance is based upon a plethora of distinct mechanisms. Drug resistance mutations can occur in the same protein or in different proteins; as well as in the same pathway or in parallel pathways, bypassing the intercepted signaling. The dilemma that the clinical oncologist is facing is that not all the genomic alterations as well as alterations in the tumor microenvironment that facilitate cancer cell proliferation are known, and neither are the alterations that are likely to promote metastasis. For example, the common KRasG12C driver mutation emerges in different cancers. Most occur in NSCLC, but some occur, albeit to a lower extent, in colorectal cancer and pancreatic ductal carcinoma. The responses to KRasG12C inhibitors are variable and fall into three categories, (i) new point mutations in KRas, or multiple copies of KRAS G12C which lead to higher expression level of the mutant protein; (ii) mutations in genes other than KRAS; (iii) original cancer transitioning to other cancer(s). Resistance to adagrasib, an experimental antitumor agent exerting its cytotoxic effect as a covalent inhibitor of the G12C KRas, indicated that half of the cases present multiple KRas mutations as well as allele amplification. Redundant or parallel pathways included MET amplification; emerging driver mutations in NRAS, BRAF, MAP2K1, and RET; gene fusion events in ALK, RET, BRAF, RAF1, and FGFR3; and loss-of-function mutations in NF1 and PTEN tumor suppressors. In the current review we discuss the molecular mechanisms underlying drug resistance while focusing on those emerging to common targeted cancer drivers. We also address questions of why cancers with a common driver mutation are unlikely to evolve a common drug resistance mechanism, and whether one can predict the likely mechanisms that the tumor cell may develop. These vastly important and tantalizing questions in drug discovery, and broadly in precision medicine, are the focus of our present review. We end with our perspective, which calls for target combinations to be selected and prioritized with the help of the emerging massive compute power which enables artificial intelligence, and the increased gathering of data to overcome its insatiable needs.

Identification of a novel prognosis-associated ceRNA network in lung adenocarcinoma via bioinformatics analysis

Background

Lung adenocarcinoma (LUAD) is the most common subtype of nonsmall-cell lung cancer (NSCLC) and has a high incidence rate and mortality. The survival of LUAD patients has increased with the development of targeted therapeutics, but the prognosis of these patients is still poor. Long noncoding RNAs (lncRNAs) play an important role in the occurrence and development of LUAD. The purpose of this study was to identify novel abnormally regulated lncRNA–microRNA (miRNA)–messenger RNA (mRNA) competing endogenous RNA (ceRNA) networks that may suggest new therapeutic targets for LUAD or relate to LUAD prognosis.

Methods

We used the SBC human ceRNA array V1.0 to screen for differentially expressed (DE) lncRNAs and mRNAs in four paired LUAD samples. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed to annotate the DE lncRNAs and mRNAs. R bioinformatics packages, The Cancer Genome Atlas (TCGA) LUAD database, and Kaplan–Meier (KM) survival analysis tools were used to validate the microarray data and construct the lncRNA–miRNA–mRNA ceRNA regulatory network. Then, quantitative real-time PCR (qRT-PCR) was used to validate the DE lncRNAs in 7 LUAD cell lines.

Results

A total of 2819 DE lncRNAs and 2396 DE mRNAs (P < 0.05 and fold change ≥ 2 or ≤ 0.5) were identified in four paired LUAD tissue samples. In total, 255 of the DE lncRNAs were also identified in TCGA. The GO and KEGG analysis results suggested that the DE genes were most enriched in angiogenesis and cell proliferation, and were closely related to human cancers. Moreover, the differential expression of ENST00000609697, ENST00000602992, and NR_024321 was consistent with the microarray data, as determined by qRT-PCR validation in 7 LUAD cell lines; however, only ENST00000609697 was associated with the overall survival of LUAD patients (log-rank P = 0.029). Finally, through analysis of ENST00000609697 target genes, we identified the ENST00000609697–hsa-miR-6791-5p–RASL12 ceRNA network, which may play a tumor-suppressive role in LUAD.

Conclusion

ENST00000609697 was abnormally expressed in LUAD. Furthermore, downregulation of ENST00000609697 and its target gene RASL12 was associated with poor prognosis in LUAD. The ENST00000609697–hsa-miR-6791-5p–RASL12 axis may play a tumor-suppressive role. These results suggest new potential prognostic and therapeutic biomarkers for LUAD.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12938-021-00952-x.

Performance Assessment of the Network Reconstruction Approaches on Various Interactomes

Beyond the list of molecules, there is a necessity to collectively consider multiple sets of omic data and to reconstruct the connections between the molecules. Especially, pathway reconstruction is crucial to understanding disease biology because abnormal cellular signaling may be pathological. The main challenge is how to integrate the data together in an accurate way. In this study, we aim to comparatively analyze the performance of a set of network reconstruction algorithms on multiple reference interactomes. We first explored several human protein interactomes, including PathwayCommons, OmniPath, HIPPIE, iRefWeb, STRING, and ConsensusPathDB. The comparison is based on the coverage of each interactome in terms of cancer driver proteins, structural information of protein interactions, and the bias toward well-studied proteins. We next used these interactomes to evaluate the performance of network reconstruction algorithms including all-pair shortest path, heat diffusion with flux, personalized PageRank with flux, and prize-collecting Steiner forest (PCSF) approaches. Each approach has its own merits and weaknesses. Among them, PCSF had the most balanced performance in terms of precision and recall scores when 28 pathways from NetPath were reconstructed using the listed algorithms. Additionally, the reference interactome affects the performance of the network reconstruction approaches. The coverage and disease- or tissue-specificity of each interactome may vary, which may result in differences in the reconstructed networks.

Challenges Faced by Clinicians in the Personalized Treatment Planning: A Literature Review and the First Results of the Russian National Cancer Program

Advances in cancer molecular profiling have enabled the development of more effective approaches to the diagnosis and personalized treatment of tumors. However, treatment planning has become more labor intensive, requiring hours or even days of clinician effort to optimize an individual patient case in a trial-and-error manner. Lessons learned from the world cancer programs provide insights into ways to develop approaches for the treatment strategy definition which can be introduced into clinical practice. This article highlights the variety of breakthroughs in patients' cancer treatment and some challenges that this field faces now in Russia. In this report, we consider the key characteristics for planning an optimal clinical treatment regimen and which should be included in the algorithm of clinical decision support systems. We discuss the perspectives of implementing artificial intelligence-based systems in cancer treatment planning in Russia.

Identification of anticancer drug target genes using an outside competitive dynamics model on cancer signaling networks

Each cancer type has its own molecular signaling network. Analyzing the dynamics of molecular signaling networks can provide useful information for identifying drug target genes. In the present study, we consider an on-network dynamics model—the outside competitive dynamics model—wherein an inside leader and an opponent competitor outside the system have fixed and different states, and each normal agent adjusts its state according to a distributed consensus protocol. If any normal agent links to the external competitor, the state of each normal agent will converge to a stable value, indicating support to the leader against the impact of the competitor. We determined the total support of normal agents to each leader in various networks and observed that the total support correlates with hierarchical closeness, which identifies biomarker genes in a cancer signaling network. Of note, by experimenting on 17 cancer signaling networks from the KEGG database, we observed that 82% of the genes among the top 3 agents with the highest total support are anticancer drug target genes. This result outperforms those of four previous prediction methods of common cancer drug targets. Our study indicates that driver agents with high support from the other agents against the impact of the external opponent agent are most likely to be anticancer drug target genes.

A Methodological Framework to Discover Pharmacogenomic Interactions Based on Random Forests

The identification of genomic alterations in tumor tissues, including somatic mutations, deletions, and gene amplifications, produces large amounts of data, which can be correlated with a diversity of therapeutic responses. We aimed to provide a methodological framework to discover pharmacogenomic interactions based on Random Forests. We matched two databases from the Cancer Cell Line Encyclopaedia (CCLE) project, and the Genomics of Drug Sensitivity in Cancer (GDSC) project. For a total of 648 shared cell lines, we considered 48,270 gene alterations from CCLE as input features and the area under the dose-response curve (AUC) for 265 drugs from GDSC as the outcomes. A three-step reduction to 501 alterations was performed, selecting known driver genes and excluding very frequent/infrequent alterations and redundant ones. For each model, we used the concordance correlation coefficient (CCC) for assessing the predictive performance, and permutation importance for assessing the contribution of each alteration. In a reasonable computational time (56 min), we identified 12 compounds whose response was at least fairly sensitive (CCC > 20) to the alteration profiles. Some diversities were found in the sets of influential alterations, providing clues to discover significant drug-gene interactions. The proposed methodological framework can be helpful for mining pharmacogenomic interactions.

Bioinformatic analysis linking genomic defects to chemosensitivity and mechanism of action

A joint analysis of the NCI60 small molecule screening data, their genetically defective genes, and mechanisms of action (MOA) of FDA approved cancer drugs screened in the NCI60 is proposed for identifying links between chemosensitivity, genomic defects and MOA. Self-Organizing-Maps (SOMs) are used to organize the chemosensitivity data. Student's t-tests are used to identify SOM clusters with enhanced chemosensitivity for tumor cell lines with versus without genetically defective genes. Fisher's exact and chi-square tests are used to reveal instances where defective gene to chemosensitivity associations have enriched MOAs. The results of this analysis find a relatively small set of defective genes, inclusive of ABL1, AXL, BRAF, CDC25A, CDKN2A, IGF1R, KRAS, MECOM, MMP1, MYC, NOTCH1, NRAS, PIK3CG, PTK2, RPTOR, SPTBN1, STAT2, TNKS and ZHX2, as possible candidates for roles in chemosensitivity for compound MOAs that target primarily, but not exclusively, kinases, nucleic acid synthesis, protein synthesis, apoptosis and tubulin. These results find exploitable instances of enhanced chemosensitivity of compound MOA's for selected defective genes. Collectively these findings will advance the interpretation of pre-clinical screening data as well as contribute towards the goals of cancer drug discovery, development decision making, and explanation of drug mechanisms.

Dual inhibitors of histone deacetylases and other cancer-related targets: A pharmacological perspective

Epigenetic enzymes histone deacetylases (HDACs) are clinically validated anticancer drug targets which have been studied intensively in the past few decades. Although several drugs have been approved in this field, they are still limited to a subset of hematological malignancies (in particular T-cell lymphomas), with therapeutic potential not fully realized and the drug-resistance occurred after a certain period of use. To maximize the therapeutic potential of these classes of anticancer drugs, and to extend their application to solid tumors, numerous combination therapies containing an HDACi and an anticancer agent from other mechanisms are currently ongoing in clinical trials. Recently, dual targeting strategy comprising the HDACs component has emerged as an alternative approach for combination therapies. In this perspective, we intend to gather all HDACs-containing dual inhibitors related to cancer therapy published in literature since 2015, classify them into five categories based on targets' biological functions, and discuss the rationale why dual acting agents should work better than combinatorial therapies using two separate drugs. The article discusses the pharmacological aspects of these dual inhibitors, including in vitro biological activities, pharmacokinetic studies, in vivo efficacy studies, as well as available clinical trials. The review of the current status and advances should provide better analysis for future opportunities and challenges of this field.