Systems biology perspectives are crucial for understanding the pathophysiology of complex diseases, and therefore hold great promise for the discovery of novel treatment strategies. Drug combinations have been shown to improve durability and reduce resistance to available first-line therapies in a variety of cancers; however, traditional drug discovery approaches are prohibitively cost and labor-intensive to evaluate large-scale matrices of potential drug combinations. Computational methods are needed to efficiently model complex interactions of drug target pathways and identify mechanisms underlying drug combination synergy. In this study, we employ a computational approach, SynGeNet (Synergy from Gene expression and Network mining), which integrates transcriptomics-based connectivity mapping and network centrality analysis to analyze disease networks and predict drug combinations. As an exemplar of a disease in which combination therapies demonstrate efficacy in genomic-specific contexts, we investigate malignant melanoma. We employed SynGeNet to generate drug combination predictions for each of the four major genomic subtypes of melanoma (BRAF, NRAS, NF1, and triple wild type) using publicly available gene expression and mutation data. We validated synergistic drug combinations predicted by our method across all genomic subtypes using results from a high-throughput drug screening study across. Finally, we prospectively validated the drug combination for BRAF -mutant melanoma that was top ranked by our approach, vemurafenib (BRAF inhibitor) + tretinoin (retinoic acid receptor agonist), using both in vitro and in vivo models of BRAF -mutant melanoma and RNA-sequencing analysis of drug-treated melanoma cells to validate the predicted mechanisms. Our approach is applicable to a wide range of disease domains, and, importantly, can model disease-relevant protein subnetworks in precision medicine contexts.

Objective: Hepatocellular carcinoma (HCC) is frequently diagnosed in patients with late-stage disease who are ineligible for curative surgical therapies. The majority of patients become resistant to sorafenib, the only approved first-line therapy for advanced cancer, underscoring the need for newer, more effective drugs. The purpose of this study is to expedite identification of novel drugs against sorafenib resistant (SR)-HCC. Methods: We employed a transcriptomics-based drug repurposing method termed connectivity mapping using gene signatures from in vitro-derived SR Huh7 HCC cells. For proof of concept validation, we focused on drugs that were FDA-approved or under clinical investigation and prioritized two anti-neoplastic agents (dasatinib and fostamatinib) with targets associated with HCC. We also prospectively validated predicted gene expression changes in drug-treated SR Huh7 cells as well as identified and validated the targets of Fostamatinib in HCC. Results: Dasatinib specifically reduced the viability of SR-HCC cells that correlated with up-regulated activity of SRC family kinases, its targets, in our SR-HCC model. However, fostamatinib was able to inhibit both parental and SR HCC cells in vitro and in xenograft models. Ingenuity pathway analysis of fostamatinib gene expression signature from LINCS predicted JAK/STAT, PI3K/AKT, ERK/MAPK pathways as potential targets of fostamatinib that were validated by Western blot analysis. Fostamatinib treatment reversed the expression of genes that were deregulated in SR HCC. Conclusion: We provide proof of concept evidence for the validity of this drug repurposing approach for SR-HCC with implications for personalized medicine.

The emergence of drug resistance to traditional chemotherapy and newer targeted therapies in cancer patients is a major clinical challenge. Reactivation of the same or compensatory signaling pathways is a common class of drug resistance mechanisms. Employing drug combinations that inhibit multiple modules of reactivated signaling pathways is a promising strategy to overcome and prevent the onset of drug resistance. However, with thousands of available FDA-approved and investigational compounds, it is infeasible to experimentally screen millions of possible drug combinations with limited resources. Therefore, computational approaches are needed to constrain the search space and prioritize synergistic drug combinations for preclinical studies. In this study, we propose a novel approach for predicting drug combinations through investigating potential effects of drug targets on disease signaling network. We first construct a disease signaling network by integrating gene expression data with disease-associated driver genes. Individual drugs that can partially perturb the disease signaling network are then selected based on a drug-disease network "impact matrix", which is calculated using network diffusion distance from drug targets to signaling network elements. The selected drugs are subsequently clustered into communities (subgroups), which are proposed to share similar mechanisms of action. Finally, drug combinations are ranked according to maximal impact on signaling sub-networks from distinct mechanism-based communities. Our method is advantageous compared to other approaches in that it does not require large amounts drug dose response data, drug-induced "omics" profiles or clinical efficacy data, which are not often readily available. We validate our approach using a BRAF-mutant melanoma signaling network and combinatorial in vitro drug screening data, and report drug combinations with diverse mechanisms of action and opportunities for drug repositioning.

NRAS is a commonly mutated oncogene in melanoma and therapeutic targeting of NRAS has proven difficult. We characterized the expression and phenotypic functions of 5 recently discovered splice isoforms of NRAS in melanoma. Canonical NRAS (isoform-1) was expressed to the highest degree and its expression was significantly increased in melanoma metastases compared to primary lesions. Isoform-5 expression in metastases showed a significant, positive correlation with survival and tumors over-expressing isoform-5 had significantly decreased growth in a xenograft model. In contrast, over-expression of any isoform resulted in enhanced proliferation, and invasiveness was increased with over-expression of isoform-1 or isoform-2. Downstream signaling analysis indicated that the isoforms signaled differentially through the MAPK and PI3K pathways and A375 cells over-expressing isoform-2 or isoform-5 showed resistance to vemurafenib treatment in vitro. The NRAS isoforms appear to play varying roles in melanoma phenotype and could potentially serve as biomarkers for therapeutic response and disease prognosis.

Background: Malignant melanoma has a propensity for development of hepatic and pulmonary metastases. MicroRNAs (miRs) are small, non-coding RNA molecules containing about 22 nucleotides that mediate protein expression and can contribute to cancer progression. We aim to identify clinically useful differences in miR expression in metastatic melanoma tissue.Methods: RNA was extracted from formalin fixed, paraffin embedded samples of hepatic and pulmonary metastatic melanoma, benign nevi, and primary cutaneous melanoma. Assessment of miR expression was performed on purified RNA using the NanoString nCounter miRNA assay. miRs with >2 fold change in expression when compared to other tumor sites (p-value ≤0.05, modified t-test) were identified as dysregulated. Common gene targets were then identified among dysregulated miRs unique to each metastatic site.