Identification of genes and critical control proteins associated with inflammatory breast cancer using network controllability

A practically efficient algorithm for identifying critical control proteins in directed probabilistic biological networks

Network controllability is unifying the traditional control theory with the structural network information rooted in many large-scale biological systems of interest, from intracellular networks in molecular biology to brain neuronal networks. In controllability approaches, the set of minimum driver nodes is not unique, and critical nodes are the most important control elements because they appear in all possible solution sets. On the other hand, a common but largely unexplored feature in network control approaches is the probabilistic failure of edges or the uncertainty in the determination of interactions between molecules. This is particularly true when directed probabilistic interactions are considered. Until now, no efficient algorithm existed to determine critical nodes in probabilistic directed networks. Here we present a probabilistic control model based on a minimum dominating set framework that integrates the probabilistic nature of directed edges between molecules and determines the critical control nodes that drive the entire network functionality. The proposed algorithm, combined with the developed mathematical tools, offers practical efficiency in determining critical control nodes in large probabilistic networks. The method is then applied to the human intracellular signal transduction network revealing that critical control nodes are associated with important biological features and perturbed sets of genes in human diseases, including SARS-CoV-2 target proteins and rare disorders. We believe that the proposed methodology can be useful to investigate multiple biological systems in which directed edges are probabilistic in nature, both in natural systems or when determined with large uncertainties in-silico.

Measuring criticality in control of complex biological networks

Recent controllability analyses have demonstrated that driver nodes tend to be associated to genes related to important biological functions as well as human diseases. While researchers have focused on identifying critical nodes, intermittent nodes have received much less attention. Here, we propose a new efficient algorithm based on the Hamming distance for computing the importance of intermittent nodes using a Minimum Dominating Set (MDS)-based control model. We refer to this metric as criticality. The application of the proposed algorithm to compute criticality under the MDS control framework allows us to unveil the biological importance and roles of the intermittent nodes in different network systems, from cellular level such as signaling pathways and cell-cell interactions such as cytokine networks, to the complete nervous system of the nematode worm C. elegans. Taken together, the developed computational tools may open new avenues for investigating the role of intermittent nodes in many biological systems of interest in the context of network control.

CASP9 As a Prognostic Biomarker and Promising Drug Target Plays a Pivotal Role in Inflammatory Breast Cancer

Background

Inflammatory breast cancer (IBC) is one of the most rare and aggressive subtypes of primary breast cancer (BC). Our study aimed to explore hub genes related to the pathogenesis of IBC, which could be considered as novel molecular biomarkers for IBC diagnosis and prognosis. Material and Methods. Two datasets from gene expression omnibus database (GEO) were selected. Enrichment analysis and protein-protein interaction (PPI) network for the DEGs were performed. We analyzed the prognostic values of hub genes in the Kaplan-Meier Plotter. Connectivity Map (CMap) and Comparative Toxicogenomics Database (CTD) was used to find candidate small molecules capable to reverse the gene status of IBC.

Results

157 DEGs were selected in total. We constructed the PPI network with 154 nodes interconnected by 128 interactions. The KEGG pathway analysis indicated that the DEGs were enriched in apoptosis, pathways in cancer and insulin signaling pathway. PTEN, PSMF1, PSMC6, AURKB, FZR1, CASP9, CASP6, CASP8, BAD, AKR7A2, ZNF24, SSX2IP, SIGLEC1, MS4A4A, and VSIG4 were selected as hub genes based on the high degree of connectivity. Six hub genes (PSMC6, AURKB, CASP9, BAD, ZNF24, and SSX2IP) that were significantly associated with the prognosis of breast cancer. The expression of CASP9 protein was associated with prognosis and immune cells infiltration of breast cancer. CASP9- naringenin (NGE) is expected to be the most promising candidate gene-compound interaction for the treatment of IBC.

Conclusion

Taken together, CASP9 can be used as a prognostic biomarker and a novel therapeutic target in IBC.

Domination based classification algorithms for the controllability analysis of biological interaction networks

Deciding the size of a minimum dominating set is a classic NP-complete problem. It has found increasing utility as the basis for classifying vertices in networks derived from protein-protein, noncoding RNA, metabolic, and other biological interaction data. In this context it can be helpful, for example, to identify those vertices that must be present in any minimum solution. Current classification methods, however, can require solving as many instances as there are vertices, rendering them computationally prohibitive in many applications. In an effort to address this shortcoming, new classification algorithms are derived and tested for efficiency and effectiveness. Results of performance comparisons on real-world biological networks are reported.

Pharmacologically controlling protein-protein interactions through epichaperomes for therapeutic vulnerability in cancer

Cancer cell plasticity due to the dynamic architecture of interactome networks provides a vexing outlet for therapy evasion. Here, through chemical biology approaches for systems level exploration of protein connectivity changes applied to pancreatic cancer cell lines, patient biospecimens, and cell- and patient-derived xenografts in mice, we demonstrate interactomes can be re-engineered for vulnerability. By manipulating epichaperomes pharmacologically, we control and anticipate how thousands of proteins interact in real-time within tumours. Further, we can essentially force tumours into interactome hyperconnectivity and maximal protein-protein interaction capacity, a state whereby no rebound pathways can be deployed and where alternative signalling is supressed. This approach therefore primes interactomes to enhance vulnerability and improve treatment efficacy, enabling therapeutics with traditionally poor performance to become highly efficacious. These findings provide proof-of-principle for a paradigm to overcome drug resistance through pharmacologic manipulation of proteome-wide protein-protein interaction networks.

Highly connected, non-redundant microRNA functional control in breast cancer molecular subtypes

Breast cancer is a complex, heterogeneous disease at the phenotypic and molecular level. In particular, the transcriptional regulatory programs are known to be significantly affected and such transcriptional alterations are able to capture some of the heterogeneity of the disease, leading to the emergence of breast cancer molecular subtypes. Recently, it has been found that network biology approaches to decipher such abnormal gene regulation programs, for instance by means of gene co-expression networks, have been able to recapitulate the differences between breast cancer subtypes providing elements to further understand their functional origins and consequences. Network biology approaches may be extended to include other co-expression patterns, like those found between genes and non-coding transcripts such as microRNAs (miRs). As is known, miRs play relevant roles in the establishment of normal and anomalous transcription processes. Commodore miRs (cdre-miRs) have been defined as miRs that, based on their connectivity and redundancy in co-expression networks, are potential control elements of biological functions. In this work, we reconstructed miR–gene co-expression networks for each breast cancer molecular subtype, from high throughput data in 424 samples from the Cancer Genome Atlas consortium. We identified cdre-miRs in three out of four molecular subtypes. We found that in each subtype, each cdre-miR was linked to a different set of associated genes, as well as a different set of associated biological functions. We used a systematic literature validation strategy, and identified that the associated biological functions to these cdre-miRs are hallmarks of cancer such as angiogenesis, cell adhesion, cell cycle and regulation of apoptosis. The relevance of such cdre-miRs as actionable molecular targets in breast cancer is still to be determined from functional studies.

Constraint-based models for dominating protein interaction networks

The minimum dominating set (MDSet) comprises the smallest number of graph nodes, where other graph nodes are connected with at least one MDSet node. The MDSet has been successfully applied to extract proteins that control protein–protein interaction (PPI) networks and to reveal the correlation between structural analysis and biological functions. Although the PPI network contains many MDSets, the identification of multiple MDSets is an NP‐complete problem, and it is difficult to determine the best MDSets, enriched with biological functions. Therefore, the MDSet model needs to be further expanded and validated to find constrained solutions that differ from those generated by the traditional models. Moreover, by identifying the critical set of the network, the set of nodes common to all MDSets can be time‐consuming. Herein, the authors adopted the minimisation of metabolic adjustment (MOMA) algorithm to develop a new framework, called maximisation of interaction adjustment (MOIA). In MOIA, they provide three models; the first one generates two MDSets with a minimum number of shared proteins, the second model generates constrained multiple MDSets (k‐MDSets), and the third model generates user‐defined MDSets, containing the maximum number of essential genes and/or other important genes of the PPI network. In practice, these models significantly reduce the cost of finding the critical set and classifying the graph nodes. Herein, the authors termed the critical set as the k‐critical set, where k is the number of MDSets generated by the proposed model. Then, they defined a new set of proteins called the (k−1)‐critical set, where each node belongs to (k−1) MDSets. This set has been shown to be as important as the k‐critical set and contains many essential genes, transcription factors, and protein kinases as the k‐critical set. The (k−1)‐critical set can be used to extend the search for drug target proteins. Based on the performance of the MOIA models, the authors believe the proposed methods contribute to answering key questions about the MDSets of PPI networks, and their results and analysis can be extended to other network types.

Control Analysis of Protein-Protein Interaction Network Reveals Potential Regulatory Targets for MYCN

BACKGROUND

MYCN is an oncogenic transcription factor of the MYC family and plays an important role in the formation of tissues and organs during development before birth. Due to the difficulty in drugging MYCN directly, revealing the molecules in MYCN regulatory networks will help to identify effective therapeutic targets.

METHODS

We utilized network controllability theory, a recent developed powerful tool, to identify the potential drug target around MYCN based on Protein-Protein interaction network of MYCN. First, we constructed a Protein-Protein interaction network of MYCN based on public databases. Second, network control analysis was applied on network to identify driver genes and indispensable genes of the MYCN regulatory network. Finally, we developed a novel integrated approach to identify potential drug targets for regulating the function of the MYCN regulatory network.

RESULTS

We constructed an MYCN regulatory network that has 79 genes and 129 interactions. Based on network controllability theory, we analyzed driver genes which capable to fully control the network. We found 10 indispensable genes whose alternation will significantly change the regulatory pathways of the MYCN network. We evaluated the stability and correlation analysis of these genes and found EGFR may be the potential drug target which closely associated with MYCN.

CONCLUSION

Together, our findings indicate that EGFR plays an important role in the regulatory network and pathways of MYCN and therefore may represent an attractive therapeutic target for cancer treatment.

Target controllability with minimal mediators in complex biological networks

Controllability of a complex network system is related to finding a set of minimum number of nodes, known as drivers, controlling which allows having a full control on the dynamics of the network. For some applications, only a portion of the network is required to be controlled, for which target control has been proposed. Often, along the controlling route from driver nodes to target nodes, some mediators (intermediate nodes) are also unwillingly controlled, which might cause various side effects. In controlling cancerous cells, unwillingly controlling healthy cells, might result in weakening them, thus affecting the immune system against cancer. This manuscript proposes a suitable candidate solution to the problem of finding minimum number of driver nodes under minimal mediators. Although many others have attempted to develop algorithms to find minimum number of drivers for target control, the newly proposed algorithm is the first one that is capable of achieving this goal and at the same time, keeping the number of the mediators to a minimum. The proposed controllability condition, based on path lengths between node pairs, meets Kalman's controllability rank condition and can be applied on directed networks. Our results show that the path length is a major determinant of in properties of the target control under minimal mediators. As the average path length becomes larger, the ratio of drivers to target nodes decreases and the ratio of mediators to targets increases. The proposed methodology has potential applications in biological networks. The source code of the algorithm and the networks that have been used are available from the following link: https://github.com/LBBSoft/Target-Control-with-Minimal-Mediators.git.

The Hidden Control Architecture of Complex Brain Networks

The brain controls various cognitive functions in a robust and efficient way. What is the control architecture of brain networks that enables such robust and optimal control? Is this brain control architecture distinct from that of other complex networks? Here, we developed a framework to delineate a control architecture of a complex network that is compatible with the behavior of the network and applied the framework to structural brain networks and other complex networks. As a result, we revealed that the brain networks have a distributed and overlapping control architecture governed by a small number of control nodes, which may be responsible for the robust and efficient brain functions. Moreover, our artificial network evolution analysis showed that the distributed and overlapping control architecture of the brain network emerges when it evolves toward increasing both robustness and efficiency.

•

We develop a framework to delineate the control architecture of brain networks

•

The control architecture of brain networks is compared with other complex networks

•

Brain networks have a distributed and overlapping control architecture

•

Robust and efficient brain functions might be rooted in its control architecture