ReCorDE: a framework for identifying drug classes targeting shared vulnerabilities with applications to synergistic drug discovery.

Cancer is typically treated with combinatorial therapy, and such combinations may be synergistic. However, discovery of these combinations has proven difficult as brute force combinatorial screening approaches are both logistically complex and resource-intensive. Therefore, computational approaches to augment synergistic drug discovery are of interest, but current approaches are limited by their dependencies on combinatorial drug screening training data or molecular profiling data. These dataset dependencies can limit the number and diversity of drugs for which these approaches can make inferences. Herein, we describe a novel computational framework, ReCorDE (Recurrent Correlation of Drugs with Enrichment), that uses publicly-available cell line-derived monotherapy cytotoxicity datasets to identify drug classes targeting shared vulnerabilities across multiple cancer lineages; and we show how these inferences can be used to augment synergistic drug combination discovery. Additionally, we demonstrate in preclinical models that a drug class combination predicted by ReCorDE to target shared vulnerabilities (PARP inhibitors and Aurora kinase inhibitors) exhibits class-class synergy across lineages. ReCorDE functions independently of combinatorial drug screening and molecular profiling data, using only extensive monotherapy cytotoxicity datasets as its input. This allows ReCorDE to make robust inferences for a large, diverse array of drugs. In conclusion, we have described a novel framework for the identification of drug classes targeting shared vulnerabilities using monotherapy cytotoxicity datasets, and we showed how these inferences can be used to aid discovery of novel synergistic drug combinations.

Identification of VWA5A as a novel biomarker for inhibiting metastasis in breast cancer by machine-learning based protein prioritization.

Distant metastasis is the leading cause of death in breast cancer (BC). The timing of distant metastasis differs according to subtypes of BCs and there is a need for identification of biomarkers for the prediction of early and late metastasis. To identify biomarker candidates whose abundance level can discriminate metastasis types, we performed a high-throughput proteomics assay using tissue samples from BCs with no metastasis, late metastasis, and early metastasis, processed data with machine learning-based feature selection, and found that low VWA5A could be responsible for shorter duration of metastasis-free interval. Low expression of VWA5A gene in METABRIC cohort was associated with poor survival in BCs, especially in hormone receptor (HR)-positive BCs. In-vitro experiments confirmed tumor suppressive effect of VWA5A on BCs in HR+ and triple-negative BC cell lines. We found that expression of VWA5A can be assessed by immunohistochemistry (IHC) on archival tissue samples. Decreasing nuclear expression of VWA5A was significantly associated with advanced T stage and lymphatic invasion in consecutive BCs of all subtypes. We discovered lower expression of VWA5A as the potential biomarker for metastasis-prone BCs, and our results support the clinical utility of VWA5A IHC, as an adjunctive tools for prognostication of BCs.

Machine Learning-Based Proteomics Reveals Ferroptosis in COPD Patient-Derived Airway Epithelial Cells Upon Smoking Exposure.

BACKGROUND: Proteomics and genomics studies have contributed to understanding the pathogenesis of chronic obstructive pulmonary disease (COPD), but previous studies have limitations. Here, using a machine learning (ML) algorithm, we attempted to identify pathways in cultured bronchial epithelial cells of COPD patients that were significantly affected when the cells were exposed to a cigarette smoke extract (CSE). METHODS: Small airway epithelial cells were collected from patients with COPD and those without COPD who underwent bronchoscopy. After expansion through primary cell culture, the cells were treated with or without CSEs, and the proteomics of the cells were analyzed by mass spectrometry. ML-based feature selection was used to determine the most distinctive patterns in the proteomes of COPD and non-COPD cells after exposure to smoke extract. Publicly available single-cell RNA sequencing data from patients with COPD (GSE136831) were used to analyze and validate our findings. RESULTS: Five patients with COPD and five without COPD were enrolled, and 7,953 proteins were detected. Ferroptosis was enriched in both COPD and non-COPD epithelial cells after their exposure to smoke extract. However, the ML-based analysis identified ferroptosis as the most dramatically different response between COPD and non-COPD epithelial cells, adjusted P value = 4.172 × 10-6, showing that epithelial cells from COPD patients are particularly vulnerable to the effects of smoke. Single-cell RNA sequencing data showed that in cells from COPD patients, ferroptosis is enriched in basal, goblet, and club cells in COPD but not in other cell types. CONCLUSION: Our ML-based feature selection from proteomic data reveals ferroptosis to be the most distinctive feature of cultured COPD epithelial cells compared to non-COPD epithelial cells upon exposure to smoke extract.

Characterization of radiation-resistance mechanism in Spirosoma montaniterrae DY10T in terms of transcriptional regulatory system.

To respond to the external environmental changes for survival, bacteria regulates expression of a number of genes including transcription factors (TFs). To characterize complex biological phenomena, a biological system-level approach is necessary. Here we utilized six computational biology methods to infer regulatory network and to characterize underlying biologically mechanisms relevant to radiation-resistance. In particular, we inferred gene regulatory network (GRN) and operons of radiation-resistance bacterium Spirosoma montaniterrae DY10[Formula: see text] and identified the major regulators for radiation-resistance. Our results showed that DNA repair and reactive oxygen species (ROS) scavenging mechanisms are key processes and Crp/Fnr family transcriptional regulator works as a master regulatory TF in early response to radiation.

Transplantation of neuron-inducing grafts embedding positively charged gold nanoparticles for the treatment of spinal cord injury.

In this study, we aimed to investigate the recovery after traumatic spinal cord injury (SCI) by inducing cellular differentiation of transplanted neural stem cells (NSCs) into neurons. We dissociated NSCs from the spinal cords of Fisher 344 rat embryos. An injectable gel crosslinked with glycol chitosan and oxidized hyaluronate was used as a vehicle for NSC transplantation. The gel graft containing the NSC and positively charged gold nanoparticles (pGNP) was implanted into spinal cord lesions in Sprague-Dawley rats (NSC-pGNP gel group). Cellular differentiation of grafted NSCs into neurons (stained with ß-tubulin III [also called Tuj1]) was significantly increased in the NSC-pGNP gel group (***p < 0.001) compared to those of two control groups (NSC and NSC gel groups) in the SCI conditions. The NSC-pGNP gel group showed the lowest differentiation into astrocytes (stained with glial fibrillary acidic protein). Regeneration of damaged axons (stained with biotinylated dextran amines) within the lesion was two-fold higher in the NSC-pGNP gel group than that in the NSC gel group. The highest locomotor scores were also found in the NSC-pGNP gel group. These outcomes suggest that neuron-inducing pGNP gel graft embedding embryonic spinal cord-derived NSCs can be a useful type of stem cell therapy after SCI.

DRIM: A Web-Based System for Investigating Drug Response at the Molecular Level by Condition-Specific Multi-Omics Data Integration.

Pharmacogenomics is the study of how genes affect a person's response to drugs. Thus, understanding the effect of drug at the molecular level can be helpful in both drug discovery and personalized medicine. Over the years, transcriptome data upon drug treatment has been collected and several databases compiled before drug treatment cancer cell multi-omics data with drug sensitivity (IC 50, AUC) or time-series transcriptomic data after drug treatment. However, analyzing transcriptome data upon drug treatment is challenging since more than 20,000 genes interact in complex ways. In addition, due to the difficulty of both time-series analysis and multi-omics integration, current methods can hardly perform analysis of databases with different data characteristics. One effective way is to interpret transcriptome data in terms of well-characterized biological pathways. Another way is to leverage state-of-the-art methods for multi-omics data integration. In this paper, we developed Drug Response analysis Integrating Multi-omics and time-series data (DRIM), an integrative multi-omics and time-series data analysis framework that identifies perturbed sub-pathways and regulation mechanisms upon drug treatment. The system takes drug name and cell line identification numbers or user's drug control/treat time-series gene expression data as input. Then, analysis of multi-omics data upon drug treatment is performed in two perspectives. For the multi-omics perspective analysis, IC 50-related multi-omics potential mediator genes are determined by embedding multi-omics data to gene-centric vector space using a tensor decomposition method and an autoencoder deep learning model. Then, perturbed pathway analysis of potential mediator genes is performed. For the time-series perspective analysis, time-varying perturbed sub-pathways upon drug treatment are constructed. Additionally, a network involving transcription factors (TFs), multi-omics potential mediator genes, and perturbed sub-pathways is constructed, and paths to perturbed pathways from TFs are determined by an influence maximization method. To demonstrate the utility of our system, we provide analysis results of sub-pathway regulatory mechanisms in breast cancer cell lines of different drug sensitivity. DRIM is available at: http://biohealth.snu.ac.kr/software/DRIM/.

Impact of mutations in DNA methylation modification genes on genome-wide methylation landscapes and downstream gene activations in pan-cancer.

BACKGROUND: In cancer, mutations of DNA methylation modification genes have crucial roles for epigenetic modifications genome-wide, which lead to the activation or suppression of important genes including tumor suppressor genes. Mutations on the epigenetic modifiers could affect the enzyme activity, which would result in the difference in genome-wide methylation profiles and, activation of downstream genes. Therefore, we investigated the effect of mutations on DNA methylation modification genes such as DNMT1, DNMT3A, MBD1, MBD4, TET1, TET2 and TET3 through a pan-cancer analysis. METHODS: First, we investigated the effect of mutations in DNA methylation modification genes on genome-wide methylation profiles. We collected 3,644 samples that have both of mRNA and methylation data from 12 major cancer types in The Cancer Genome Atlas (TCGA). The samples were divided into two groups according to the mutational signature. Differentially methylated regions (DMR) that overlapped with the promoter region were selected using minfi and differentially expressed genes (DEG) were identified using EBSeq. By integrating the DMR and DEG results, we constructed a comprehensive DNA methylome profiles on a pan-cancer scale. Second, we investigated the effect of DNA methylations in the promoter regions on downstream genes by comparing the two groups of samples in 11 cancer types. To investigate the effects of promoter methylation on downstream gene activations, we performed clustering analysis of DEGs. Among the DEGs, we selected highly correlated gene set that had differentially methylated promoter regions using graph based sub-network clustering methods. RESULTS: We chose an up-regulated DEGs cluster where had hypomethylated promoter in acute myeloid leukemia (LAML) and another down-regulated DEGs cluster where had hypermethylated promoter in colon adenocarcinoma (COAD). To rule out effects of gene regulation by transcription factor (TF), if differentially expressed TFs bound to the promoter of DEGs, that DEGs did not included to the gene set that effected by DNA methylation modifiers. Consequently, we identified 54 hypomethylated promoter DMR up-regulated DEGs in LAML and 45 hypermethylated promoter DMR down-regulated DEGs in COAD. CONCLUSIONS: Our study on DNA methylation modification genes in mutated vs. non-mutated groups could provide useful insight into the epigenetic regulation of DEGs in cancer.