Neighbours of cancer-related proteins have key influence on pathogenesis and could increase the drug target space for anticancer therapies

Cancer drug resistance as learning of signaling networks

Drug resistance is a major cause of tumor mortality. Signaling networks became useful tools for driving pharmacological interventions against cancer drug resistance. Signaling datasets now cover the entire human cell. Recently, network adaptation became understood as a learning process. We review rapidly increasing evidence showing that the development of cancer drug resistance can be described as learning of signaling networks. During drug adaptation, the network forgets drug-affected pathways by desensitization and relearns by strengthening alternative pathways. Thus, resistant cancer cells develop a drug resistance memory. We show that all key players of cellular learning (i.e., IDPs, protein translocation, microRNAs/lncRNAs, scaffolding proteins and epigenetic/chromatin memory) have important roles in the development of cancer drug resistance. Moreover, all of them are central components of the epithelial-mesenchymal transition leading to metastases and resistance. Phenotypic plasticity was recently listed as a hallmark of cancer. We review how network plasticity induces rare, pre-existent drug-resistant cells in the absence of drug treatment. Key network methods assessing the development of drug resistance and network pharmacological interventions against drug resistance are summarized. Finally, we highlight the class of cellular memory drugs affecting cellular learning and forgetting, and we summarize current challenges to prevent or break drug resistance using network models.

A bioinformatics approach to elucidate conserved genes and pathways in C. elegans as an animal model for cardiovascular research

Cardiovascular disease (CVD) is a collective term for disorders of the heart and blood vessels. The molecular events and biochemical pathways associated with CVD are difficult to study in clinical settings on patients and in vitro conditions. Animal models play a pivotal and indispensable role in CVD research. Caenorhabditis elegans, a nematode species, has emerged as a prominent experimental organism widely utilized in various biomedical research fields. However, the specific number of CVD-related genes and pathways within the C. elegans genome remains undisclosed to date, limiting its in-depth utilization for investigations. In the present study, we conducted a comprehensive analysis of genes and pathways related to CVD within the genomes of humans and C. elegans through a systematic bioinformatic approach. A total of 1113 genes in C. elegans orthologous to the most significant CVD-related genes in humans were identified, and the GO terms and pathways were compared to study the pathways that are conserved between the two species. In order to infer the functions of CVD-related orthologous genes in C. elegans, a PPI network was constructed. Orthologous gene PPI network analysis results reveal the hubs and important KRs: pmk-1, daf-21, gpb-1, crh-1, enpl-1, eef-1G, acdh-8, hif-1, pmk-2, and aha-1 in C. elegans. Modules were identified for determining the role of the orthologous genes at various levels in the created network. We also identified 9 commonly enriched pathways between humans and C. elegans linked with CVDs that include autophagy (animal), the ErbB signaling pathway, the FoxO signaling pathway, the MAPK signaling pathway, ABC transporters, the biosynthesis of unsaturated fatty acids, fatty acid metabolism, glutathione metabolism, and metabolic pathways. This study provides the first systematic genomic approach to explore the CVD-associated genes and pathways that are present in C. elegans, supporting the use of C. elegans as a prominent animal model organism for cardiovascular diseases.

Conserved Cardiovascular Network: Bioinformatics Insights into Genes and Pathways for Establishing Caenorhabditis elegans as an Animal Model for Cardiovascular Diseases

Cardiovascular disease (CVD) is a collective term for disorders of the heart and blood vessels. The molecular events and biochemical pathways associated with CVD are difficult to study in clinical settings on patients and in vitro conditions. Animal models play a pivotal and indispensable role in cardiovascular disease (CVD) research. Caenorhabditis elegans , a nematode species, has emerged as a prominent experimental organism widely utilised in various biomedical research fields. However, the specific number of CVD-related genes and pathways within the C. elegans genome remains undisclosed to date, limiting its in-depth utilisation for investigations. In the present study, we conducted a comprehensive analysis of genes and pathways related to CVD within the genomes of humans and C. elegans through a systematic bioinformatic approach. A total of 1113 genes in C. elegans orthologous to the most significant CVD-related genes in humans were identified, and the GO terms and pathways were compared to study the pathways that are conserved between the two species. In order to infer the functions of CVD-related orthologous genes in C. elegans, a PPI network was constructed. Orthologous gene PPI network analysis results reveal the hubs and important KRs: pmk-1, daf-21, gpb-1, crh-1, enpl-1, eef-1G, acdh-8, hif-1, pmk-2, and aha-1 in C. elegans. Modules were identified for determining the role of the orthologous genes at various levels in the created network. We also identified 9 commonly enriched pathways between humans and C. elegans linked with CVDs that include autophagy (animal), the ErbB signalling pathway, the FoxO signalling pathway, the MAPK signalling pathway, ABC transporters, the biosynthesis of unsaturated fatty acids, fatty acid metabolism, glutathione metabolism, and metabolic pathways. This study provides the first systematic genomic approach to explore the CVD-associated genes and pathways that are present in C. elegans, supporting the use of C. elegans as a prominent animal model organism for cardiovascular diseases.

Chemogenomic library design strategies for precision oncology, applied to phenotypic profiling of glioblastoma patient cells

Designing a targeted screening library of bioactive small molecules is a challenging task since most compounds modulate their effects through multiple protein targets with varying degrees of potency and selectivity. We implemented analytic procedures for designing anticancer compound libraries adjusted for library size, cellular activity, chemical diversity and availability, and target selectivity. The resulting compound collections cover a wide range of protein targets and biological pathways implicated in various cancers, making them widely applicable to precision oncology. We characterized the compound and target spaces of the virtual libraries, in comparison with a minimal screening library of 1,211 compounds for targeting 1,386 anticancer proteins. In a pilot screening study, we identified patient-specific vulnerabilities by imaging glioma stem cells from patients with glioblastoma (GBM), using a physical library of 789 compounds that cover 1,320 of the anticancer targets. The cell survival profiling revealed highly heterogeneous phenotypic responses across the patients and GBM subtypes.

A systems genomics approach to uncover patient-specific pathogenic pathways and proteins in ulcerative colitis

We describe a precision medicine workflow, the integrated single nucleotide polymorphism network platform (iSNP), designed to determine the mechanisms by which SNPs affect cellular regulatory networks, and how SNP co-occurrences contribute to disease pathogenesis in ulcerative colitis (UC). Using SNP profiles of 378 UC patients we map the regulatory effects of the SNPs to a human signalling network containing protein-protein, miRNA-mRNA and transcription factor binding interactions. With unsupervised clustering algorithms we group these patient-specific networks into four distinct clusters driven by PRKCB, HLA, SNAI1/CEBPB/PTPN1 and VEGFA/XPO5/POLH hubs. The pathway analysis identifies calcium homeostasis, wound healing and cell motility as key processes in UC pathogenesis. Using transcriptomic data from an independent patient cohort, with three complementary validation approaches focusing on the SNP-affected genes, the patient specific modules and affected functions, we confirm the regulatory impact of non-coding SNPs. iSNP identified regulatory effects for disease-associated non-coding SNPs, and by predicting the patient-specific pathogenic processes, we propose a systems-level way to stratify patients.

Construction and contextualization approaches for protein-protein interaction networks

Protein-protein interaction network (PPIN) analysis is a widely used method to study the contextual role of proteins of interest, to predict novel disease genes, disease or functional modules, and to identify novel drug targets. PPIN-based analysis uses both generic and context-specific networks. Multiple contextualization methodologies have been described, such as shortest-path algorithms, neighborhood-based methods, and diffusion/propagation algorithms. This review discusses these methods, provides intuitive representations of PPIN contextualization, and also examines how the quality of such context-specific networks could be improved by considering additional sources of evidence. As a heuristic, we observe that tasks such as identifying disease genes, drug targets, and protein complexes should consider local neighborhoods, while uncovering disease mechanisms and discovering disease-pathways would gain from diffusion-based construction.

Network Biology Approaches to Achieve Precision Medicine in Inflammatory Bowel Disease

Inflammatory bowel disease (IBD) is a chronic immune-mediated condition arising due to complex interactions between multiple genetic and environmental factors. Despite recent advances, the pathogenesis of the condition is not fully understood and patients still experience suboptimal clinical outcomes. Over the past few years, investigators are increasingly capturing multi-omics data from patient cohorts to better characterise the disease. However, reaching clinically translatable endpoints from these complex multi-omics datasets is an arduous task. Network biology, a branch of systems biology that utilises mathematical graph theory to represent, integrate and analyse biological data through networks, will be key to addressing this challenge. In this narrative review, we provide an overview of various types of network biology approaches that have been utilised in IBD including protein-protein interaction networks, metabolic networks, gene regulatory networks and gene co-expression networks. We also include examples of multi-layered networks that have combined various network types to gain deeper insights into IBD pathogenesis. Finally, we discuss the need to incorporate other data sources including metabolomic, histopathological, and high-quality clinical meta-data. Together with more robust network data integration and analysis frameworks, such efforts have the potential to realise the key goal of precision medicine in IBD.

An Evaluation of the Effectiveness of Image-based Texture Features Extracted from Static B-mode Ultrasound Images in Distinguishing between Benign and Malignant Ovarian Masses

Significant successes in machine learning approaches to image analysis for various applications have energized strong interest in automated diagnostic support systems for medical images. The evolving in-depth understanding of the way carcinogenesis changes the texture of cellular networks of a mass/tumor has been informing such diagnostics systems with use of more suitable image texture features and their extraction methods. Several texture features have been recently applied in discriminating malignant and benign ovarian masses by analysing B-mode images from ultrasound scan of the ovary with different levels of performance. However, comparative performance evaluation of these reported features using common sets of clinically approved images is lacking. This paper presents an empirical evaluation of seven commonly used texture features (histograms, moments of histogram, local binary patterns [256-bin and 59-bin], histograms of oriented gradients, fractal dimensions, and Gabor filter), using a collection of 242 ultrasound scan images of ovarian masses of various pathological characteristics. The evaluation examines not only the effectiveness of classification schemes based on the individual texture features but also the effectiveness of various combinations of these schemes using the simple majority-rule decision level fusion. Trained support vector machine classifiers on the individual texture features without any specific pre-processing, achieve levels of accuracy between 75% and 85% where the seven moments and the 256-bin LBP are at the lower end while the Gabor filter is at the upper end. Combining the classification results of the top k (k = 3, 5, 7) best performing features further improve the overall accuracy to a level between 86% and 90%. These evaluation results demonstrate that each of the investigated image-based texture features provides informative support in distinguishing benign or malignant ovarian masses.

Identification of the key target profiles underlying the drugs of narrow therapeutic index for treating cancer and cardiovascular disease

An appropriate therapeutic index is crucial for drug discovery and development since narrow therapeutic index (NTI) drugs with slight dosage variation may induce severe adverse drug reactions or potential treatment failure. To date, the shared characteristics underlying the targets of NTI drugs have been explored by several studies, which have been applied to identify potential drug targets. However, the association between the drug therapeutic index and the related disease has not been dissected, which is important for revealing the NTI drug mechanism and optimizing drug design. Therefore, in this study, two classes of disease (cancers and cardiovascular disorders) with the largest number of NTI drugs were selected, and the target property of the corresponding NTI drugs was analyzed. By calculating the biological system profiles and human protein–protein interaction (PPI) network properties of drug targets and adopting an AI-based algorithm, differentiated features between two diseases were discovered to reveal the distinct underlying mechanisms of NTI drugs in different diseases. Consequently, ten shared features and four unique features were identified for both diseases to distinguish NTI from NNTI drug targets. These computational discoveries, as well as the newly found features, suggest that in the clinical study of avoiding narrow therapeutic index in those diseases, the ability of target to be a hub and the efficiency of target signaling in the human PPI network should be considered, and it could thus provide novel guidance in the drug discovery and clinical research process and help to estimate the drug safety of cancer and cardiovascular disease.

Regulatory network analysis of Paneth cell and goblet cell enriched gut organoids using transcriptomics approaches

The epithelial lining of the small intestine consists of multiple cell types, including Paneth cells and goblet cells, that work in cohort to maintain gut health. 3D in vitro cultures of human primary epithelial cells, called organoids, have become a key model to study the functions of Paneth cells and goblet cells in normal and diseased conditions. Advances in these models include the ability to skew differentiation to particular lineages, providing a useful tool to study cell type specific function/dysfunction in the context of the epithelium. Here, we use comprehensive profiling of mRNA, microRNA and long non-coding RNA expression to confirm that Paneth cell and goblet cell enrichment of murine small intestinal organoids (enteroids) establishes a physiologically accurate model. We employ network analysis to infer the regulatory landscape altered by skewing differentiation, and using knowledge of cell type specific markers, we predict key regulators of cell type specific functions: Cebpa, Jun, Nr1d1 and Rxra specific to Paneth cells, Gfi1b and Myc specific for goblet cells and Ets1, Nr3c1 and Vdr shared between them. Links identified between these regulators and cellular phenotypes of inflammatory bowel disease (IBD) suggest that global regulatory rewiring during or after differentiation of Paneth cells and goblet cells could contribute to IBD aetiology. Future application of cell type enriched enteroids combined with the presented computational workflow can be used to disentangle multifactorial mechanisms of these cell types and propose regulators whose pharmacological targeting could be advantageous in treating IBD patients with Crohn's disease or ulcerative colitis.

QSAR-derived affinity fingerprints (part 2): modeling performance for potency prediction

Affinity fingerprints report the activity of small molecules across a set of assays, and thus permit to gather information about the bioactivities of structurally dissimilar compounds, where models based on chemical structure alone are often limited, and model complex biological endpoints, such as human toxicity and in vitro cancer cell line sensitivity. Here, we propose to model in vitro compound activity using computationally predicted bioactivity profiles as compound descriptors. To this aim, we apply and validate a framework for the calculation of QSAR-derived affinity fingerprints (QAFFP) using a set of 1360 QSAR models generated using K_i, K_d, IC₅₀ and EC₅₀ data from ChEMBL database. QAFFP thus represent a method to encode and relate compounds on the basis of their similarity in bioactivity space. To benchmark the predictive power of QAFFP we assembled IC₅₀ data from ChEMBL database for 18 diverse cancer cell lines widely used in preclinical drug discovery, and 25 diverse protein target data sets. This study complements part 1 where the performance of QAFFP in similarity searching, scaffold hopping, and bioactivity classification is evaluated. Despite being inherently noisy, we show that using QAFFP as descriptors leads to errors in prediction on the test set in the ~ 0.65-0.95 pIC₅₀ units range, which are comparable to the estimated uncertainty of bioactivity data in ChEMBL (0.76-1.00 pIC₅₀ units). We find that the predictive power of QAFFP is slightly worse than that of Morgan2 fingerprints and 1D and 2D physicochemical descriptors, with an effect size in the 0.02-0.08 pIC₅₀ units range. Including QSAR models with low predictive power in the generation of QAFFP does not lead to improved predictive power. Given that the QSAR models we used to compute the QAFFP were selected on the basis of data availability alone, we anticipate better modeling results for QAFFP generated using more diverse and biologically meaningful targets. Data sets and Python code are publicly available at https://github.com/isidroc/QAFFP_regression .

EMDIP: An Entropy Measure to Discover Important Proteins in PPI networks

Discovering important proteins in Protein-Protein Interaction (PPI) networks has attracted a lot of attention in recent years. Most of the previous work applies different network centrality measures such as Closeness, Betweenness, PageRank and many others to discover the most influential proteins in PPI networks. Although entropy is a well-known graph-based method in computer science, according to our knowledge, it is not used in the biology domain for this purpose. In this paper, first, we annotate the human PPI network with available annotation data. Second, we introduce a new concept called annotation-context that describes each protein according to annotation data of its neighbors. Third, we apply an entropy measure to discover proteins with varied annotation-context. Empirical results indicate that our proposed method succeeded in (1) differentiating essential and non-essential proteins in PPI networks with annotation data; (2) outperforming centrality measures in the task of discovering essential nodes; (3) predicting new annotated proteins based on existing annotation data.

Gene relevance based on multiple evidences in complex networks

MOTIVATION

Multi-omics approaches offer the opportunity to reconstruct a more complete picture of the molecular events associated with human diseases, but pose challenges in data analysis. Network-based methods for the analysis of multi-omics leverage the complex web of macromolecular interactions occurring within cells to extract significant patterns of molecular alterations. Existing network-based approaches typically address specific combinations of omics and are limited in terms of the number of layers that can be jointly analysed. In this study, we investigate the application of network diffusion to quantify gene relevance on the basis of multiple evidences (layers).

RESULTS

We introduce a gene score (mND) that quantifies the relevance of a gene in a biological process taking into account the network proximity of the gene and its first neighbours to other altered genes. We show that mND has a better performance over existing methods in finding altered genes in network proximity in one or more layers. We also report good performances in recovering known cancer genes. The pipeline described in this article is broadly applicable, because it can handle different types of inputs: in addition to multi-omics datasets, datasets that are stratified in many classes (e.g., cell clusters emerging from single cell analyses) or a combination of the two scenarios.

AVAILABILITY AND IMPLEMENTATION

The R package 'mND' is available at URL: https://www.itb.cnr.it/mnd.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Topological integration of RPPA proteomic data with multi-omics data for survival prediction in breast cancer via pathway activity inference

BACKGROUND

The analysis of integrated multi-omics data enables the identification of disease-related biomarkers that cannot be identified from a single omics profile. Although protein-level data reflects the cellular status of cancer tissue more directly than gene-level data, past studies have mainly focused on multi-omics integration using gene-level data as opposed to protein-level data. However, the use of protein-level data (such as mass spectrometry) in multi-omics integration has some limitations. For example, the correlation between the characteristics of gene-level data (such as mRNA) and protein-level data is weak, and it is difficult to detect low-abundance signaling proteins that are used to target cancer. The reverse phase protein array (RPPA) is a highly sensitive antibody-based quantification method for signaling proteins. However, the number of protein features in RPPA data is extremely low compared to the number of gene features in gene-level data. In this study, we present a new method for integrating RPPA profiles with RNA-Seq and DNA methylation profiles for survival prediction based on the integrative directed random walk (iDRW) framework proposed in our previous study. In the iDRW framework, each omics profile is merged into a single pathway profile that reflects the topological information of the pathway. In order to address the sparsity of RPPA profiles, we employ the random walk with restart (RWR) approach on the pathway network.

RESULTS

Our model was validated using survival prediction analysis for a breast cancer dataset from The Cancer Genome Atlas. Our proposed model exhibited improved performance compared with other methods that utilize pathway information and also out-performed models that did not include the RPPA data utilized in our study. The risk pathways identified for breast cancer in this study were closely related to well-known breast cancer risk pathways.

CONCLUSIONS

Our results indicated that RPPA data is useful for survival prediction for breast cancer patients under our framework. We also observed that iDRW effectively integrates RNA-Seq, DNA methylation, and RPPA profiles, while variation in the composition of the omics data can affect both prediction performance and risk pathway identification. These results suggest that omics data composition is a critical parameter for iDRW.

Integrating gene regulatory pathways into differential network analysis of gene expression data

The advent of next-generation sequencing has introduced new opportunities in analyzing gene expression data. Research in systems biology has taken advantage of these opportunities by gleaning insights into gene regulatory networks through the analysis of gene association networks. Contrasting networks from different populations can reveal the many different roles genes fill, which can lead to new discoveries in gene function. Pathologies can also arise from aberrations in these gene-gene interactions. Exposing these network irregularities provides a new avenue for understanding and treating diseases. A general framework for integrating known gene regulatory pathways into a differential network analysis between two populations is proposed. The framework importantly allows for any gene-gene association measure to be used, and inference is carried out through permutation testing. A simulation study investigates the performance in identifying differentially connected genes when incorporating known pathways, even if the pathway knowledge is partially inaccurate. Another simulation study compares the general framework with four state-of-the-art methods. Two RNA-seq datasets are analyzed to illustrate the use of this framework in practice. In both examples, the analysis reveals genes and pathways that are known to be biologically significant along with potentially novel findings that may be used to motivate future research.

Discrete dynamic network modeling of oncogenic signaling: Mechanistic insights for personalized treatment of cancer

Targeted drugs disrupting proteins that are dysregulated in cancer have emerged as promising treatments because of their specificity to cancer cell aberrations and thus their improved side effect profile. However, their success remains limited, largely due to existing or emergent therapy resistance. We suggest that this is due to limited understanding of the entire relevant cellular landscape. A class of mathematical models called discrete dynamic network models can be used to understand the integrated effect of an individual tumor's aberrations. We review the recent literature on discrete dynamic models of cancer and highlight their predicted therapeutic strategies. We believe dynamic network modeling can be used to drive treatment decision-making in a personalized manner to direct improved treatments in cancer.

•

Cancer is rooted in incorrect cellular decisions caused by genetic alterations.

•

Dynamic models of signaling networks can map the relevant repertoire of alterations.

•

Discrete dynamic network models can predict therapeutic interventions.

•

Progress in personalized medicine needs integration of multiple data and model types.

A side-effect free method for identifying cancer drug targets

Identifying effective drug targets, with little or no side effects, remains an ever challenging task. A potential pitfall of failing to uncover the correct drug targets, due to side effect of pleiotropic genes, might lead the potential drugs to be illicit and withdrawn. Simplifying disease complexity, for the investigation of the mechanistic aspects and identification of effective drug targets, have been done through several approaches of protein interactome analysis. Of these, centrality measures have always gained importance in identifying candidate drug targets. Here, we put forward an integrated method of analysing a complex network of cancer and depict the importance of k-core, functional connectivity and centrality (KFC) for identifying effective drug targets. Essentially, we have extracted the proteins involved in the pathways leading to cancer from the pathway databases which enlist real experimental datasets. The interactions between these proteins were mapped to build an interactome. Integrative analyses of the interactome enabled us to unearth plausible reasons for drugs being rendered withdrawn, thereby giving future scope to pharmaceutical industries to potentially avoid them (e.g. ESR1, HDAC2, F2, PLG, PPARA, RXRA, etc). Based upon our KFC criteria, we have shortlisted ten proteins (GRB2, FYN, PIK3R1, CBL, JAK2, LCK, LYN, SYK, JAK1 and SOCS3) as effective candidates for drug development.