Lung cancer is one of the most important health risks worldwide for human. Non-small cell lung cancer (NSCLC) is the most common cause of premature death from malignant disease. This study provides in-depth insights from systems biology analyses to identify molecular to inform systemic drug targeting in NSCLC. Gene expression profiles from non small cell lung cancer were analyzed with genome-scale biomolecular networks (I,e., protein-protein interaction, transcriptional and post transcriptional regulatory networks). The aim of the study was to determine the pathways and expression profile of the genes to discover molecular signature at RNA and protein levels which could serve as potential drug targets for therapeutics innovation and the identification of novel targets. Eight proteins, six TFs and seven miRNAs came into prominence as potential drug targets. The differential expression profiles of these reporter biomolecules were cross-validated by independent RNA-Seq and miRNA-Seq. Risk discrimination performance of the reporter biomolecules NPR3, JUN, PPARG, TP53, CKMT1A, SP3 and TFAP2A were also evaluated. Total 213 drugs and 7 proteins were found for non small cell lung cancer through dgidb. Among these identified drugs seven drugs such as- Gemcitabine, Carboplatin, paclitaxel, Docetaxel, Crizotinib, Bevacizumab and Gemcitabine is used for NSCLC which is approved by National Cancer Institute. The molecular signatures and repurposed drugs presented here permit further attention for experimental studies which are offer significant potential as biomarkers and candidate therapeutics for precision medicine approaches to clinical management of NSCLC.

Background: Identification of cancer biomarkers that are differentially expressed (DE) under two biological conditions is an important task in many microarray studies. There exist several methods in the literature in this regards and most of these methods designed especially for unpaired samples, which does not satisfy the requirements of paired samples where the gene expressions are taken from the same patients before and after treatment. Furthermore, the traditional biomarker identification methods based on either p-values or fold change (FC) values. However, sometimes, p-value based results do not comply with FC based results due to the smaller variance of gene expressions. There are some methods that combine both p-values and FC values to solve this problem. But, these methods also show weak performance for small-sample case in presence of outlying expressions. To overcome this problem, in this paper an attempt is made to develop a hybrid robust SAM-FC approach by combining rank of FC values and rank of p-values based on SAM statistic using minimum β-divergence method, which is designed for paired samples. This method introduces a weight function known as β-weight function. This weight function produces larger weights corresponding to usual/normal expressions and smaller weights for unusual/outlying expressions. The β-weight function plays the significant role on the performance of the proposed method. Results: The proposed method uses β-weight function as a measure of outlier detection by setting β=0.2. We unify both classical and robust estimates using β-weight function such that maximum likelihood estimators (MLEs) are used in absence of outliers and minimum β-divergence estimators are used in presence of outliers to obtain reasonable p-values and FC values in the proposed method. We examined the performance of proposed method in a comparison of some popular methods (t-test, SAM, LIMMA, Wilcoxon, WAD, RP and FCROS) using both simulated and real gene expression profiles for both small-and large-sample cases. From the simulation and a real spike in data analysis results we observed that the proposed method outperforms other methods for small-sample case in presence of outliers and it keeps almost equal performance with other robust methods (Wilcoxon, RP and FCROS) otherwise. From a head-and-neck cancer (HNC) dataset the proposed method identified 2 genes (CYP3A4, NOVA1) that are significantly enriched in linoleic acid metabolism, drug metabolism, steroid hormone biosynthesis and metabolic pathways. The survival analysis through Kaplan-Meier curve revealed that combined effect of these 2 genes has prognostic capability and they might be promising biomarker of HNC. Moreover, we retrieved the 12 candidate drugs based on gene interaction from glad4u and drug bank databases. Conclusion The identified drugs showed statistical significance and critical role of the proteins indicate that these proteins might be therapeutic target in cancer. Thus, elucidating the associations between the drugs identified in the present study require further investigations.

Assessment of drugs toxicity and associated biomarker genes is one of the most important tasks in the pre-clinical phase of drug development pipeline as well as in the toxicogenomic studies. There are few statistical methods for the assessment of doses of drugs (DDs) toxicity and their associated biomarker genes. However, these methods consume more time for computation of the model parameters using the EM (Expectation-Maximization) based iterative approaches. To overcome this problem, in this paper, an attempt is made to propose an alternative approach based on hierarchical clustering (HC) for the same purpose. There are several types of HC approaches whose performance depends on different similarity/distance measures. Therefore, we explored suitable combinations of distance measures and HC methods based on Japanese Toxicogenomics Project (TGP) datasets for better clustering/co-clustering between DDs and genes as well as to detect toxic DDs and their associated biomarker genes. We observed that Word’s HC method with each of Euclidean, Manhattan and Minkowski distance measures produces better clustering/co-clustering results. For an example, in case of glutathione metabolism pathway (GMP) dataset LOC100359539/Rrm2, Gpx6, RGD1562107, Gstm4, Gstm3, G6pd, Gsta5, Gclc, Mgst2, Gsr, Gpx2, Gclm, Gstp1, LOC100912604/Srm, Gstm4, Odc1, Gsr, Gss are the biomarker genes and Acetaminophen_Middle, Acetaminophen_High, Methapyrilene_High, Nitrofurazone_High, Nitrofurazone_Middle, Isoniazid_Middle, Isoniazid_High are their regulatory (associated) DDs explored by our proposed co-clustering algorithm based on the distance and HC method combination Euclidean: Word. Similarly, for the PPAR signaling pathway (PPAR-SP) dataset Cpt1a, Cyp8b1, Cyp4a3, Ehhadh, Plin5, Plin2, Fabp3, Me1, Fabp5, LOC100910385, Cpt2, Acaa1a, Cyp4a1, LOC100365047, Cpt1a, LOC100365047, Angptl4, Aqp7, Cpt1c, Cpt1b, Me1 are the biomarker genes and Aspirin_Low, Aspirin_Middle, Aspirin_High, Benzbromarone_Middle, Benzbromarone_High, Clofibrate_Middle, Clofibrate_High, WY14643_Low, WY14643_High, WY14643_Middle, Gemfibrozil_Middle, Gemfibrozil_High are their regulatory DDs. These results are validated by the available literature and functional annotation.

Background and objectives: Colorectal cancer (CRC) is the second most common cause of cancer-related death in the world, but early diagnosis ameliorates the survival of CRC. This report directed to identify molecular biomarker signatures in CRC. Materials and Methods: We analyzed two microarray datasets (GSE35279 and GSE21815) from Gene Expression Omnibus (GEO) to identify mutual differentially expressed genes (DEGs). We integrated DEGs with protein-protein interaction and transcriptional/post-transcriptional regulatory networks to identify reporter signaling and regulatory molecules; utilized functional overrepresentation and pathway enrichment analyses to elucidate their roles in biological processes and molecular pathways; performed survival analyses to evaluate their prognostic performance; and applied drug repositioning analyses through Connectivity map (CMap) and geneXpharma tools to hypothesize possible drug candidates targeting reporter molecules. Results: A total of 727 up-regulated and 99 down-regulated DEGs were detected. The PI3K-Akt signaling, Wnt signaling, ECM-interaction, and cell cycle were identified as significantly enriched pathways. Ten hub proteins (ADNP, CCND1, CD44, CDK4, CEBPB, CENPA, CENPH, CENPN, MYC, and RFC2), 10 transcription factors (ETS1, ESR1, GATA1, GATA2, GATA3, AR, YBX1, FOXP3, E2F4, and PRDM14) and 2 miRNAs (miR-193b-3p and miR-615-3p) were detected as reporter molecules. The survival analyses through Kaplan Meier curves indicated remarkable performance of reporter molecules in estimation of survival probability in CRC patients. In addition, several drug candidates including anti-neoplastic and immunomodulating agents were repositioned. Conclusions: This study presents biomarker signatures at protein and RNA levels with prognostic capability in CRC. We think that the molecular signatures and candidate drugs presented in this study might be useful in future studies indenting development of accurate diagnostic and/or prognostic biomarker screens and efficient therapeutic strategies in CRC.