Information theory in systems biology. Part II: protein-protein interaction and signaling networks

The crosstalking immune cells network creates a collective function beyond the function of each cellular constituent during the progression of hepatocellular carcinoma

Abundance of data on the role of inflammatory immune responses in the progression or inhibition of hepatocellular carcinoma (HCC) has failed to offer a curative immunotherapy for HCC. This is largely because of focusing on detailed specific cell types and missing the collective function of the hepatic immune system. To discover the collective immune function, we take systems immunology approach by performing high-throughput analysis of snRNAseq data collected from the liver of DIAMOND mice during the progression of nonalcoholic fatty liver disease (NAFLD) to HCC. We report that mutual signaling interactions of the hepatic immune cells in a dominant-subdominant manner, as well as their interaction with structural cells shape the immunological pattern manifesting a collective function beyond the function of the cellular constituents. Such pattern discovery approach recognized direct role of the innate immune cells in the progression of NASH and HCC. These data suggest that discovery of the immune pattern not only detects the immunological mechanism of HCC in spite of dynamic changes in immune cells during the course of disease but also offers immune modulatory interventions for the treatment of NAFLD and HCC.

Rewiring of miRNA-mRNA bipartite co-expression network as a novel way to understand the prostate cancer related players

The differential expression and direct targeting of mRNA by miRNA are two main logics of the traditional approach to constructing the miRNA-mRNA network. This approach, could be led to the loss of considerable information and some challenges of direct targeting. To avoid these problems, we analyzed the rewiring network and constructed two miRNA-mRNA expression bipartite networks for both normal and primary prostate cancer tissue obtained from PRAD-TCGA. We then calculated beta-coefficient of the regression-model when miR was dependent and mRNA independent for each miR and mRNA and separately in both networks. We defined the rewired edges as a significant change in the regression coefficient between normal and cancer states. The rewired nodes through multinomial distribution were defined and network from rewired edges and nodes was analyzed and enriched. Of the 306 rewired edges, 112(37%) were new, 123(40%) were lost, 44(14%) were strengthened, and 27(9%) weakened connections were discovered. The highest centrality of 106 rewired mRNAs belonged to PGM5, BOD1L1, C1S, SEPG, TMEFF2, and CSNK2A1. The highest centrality of 68 rewired miRs belonged to miR-181d, miR-4677, miR-4662a, miR-9.3, and miR-1301. SMAD and beta-catenin binding were enriched as molecular functions. The regulation was a frequently repeated concept in the biological process. Our rewiring analysis highlighted the impact of β-catenin and SMAD signaling as also some transcript factors like TGFB1I1 in prostate cancer progression. Altogether, we developed a miRNA-mRNA co-expression bipartite network to identify the hidden aspects of the prostate cancer mechanism, which traditional analysis -like differential expression- was not detect it.

Impact of 5HydroxyMethylCytosine (5hmC) on reverse/direct association of cell-cycle, apoptosis, and extracellular matrix pathways in gastrointestinal cancers

BACKGROUND

Aberrant levels of 5-hydroxymethylcytosine (5-hmC) can lead to cancer progression. Identification of 5-hmC-related biological pathways in cancer studies can produce better understanding of gastrointestinal (GI) cancers. We conducted a network-based analysis on 5-hmC levels extracted from circulating free DNAs (cfDNA) in GI cancers including colon, gastric, and pancreatic cancers, and from healthy donors. The co-5-hmC network was reconstructed using the weighted-gene co-expression network method. The cancer-related modules/subnetworks were detected. Preservation of three detected 5-hmC-related modules was assessed in an external dataset. The 5-hmC-related modules were functionally enriched, and biological pathways were identified. The relationship between modules was assessed using the Pearson correlation coefficient (p-value < 0.05). An elastic network classifier was used to assess the potential of the 5-hmC modules in distinguishing cancer patients from healthy individuals. To assess the efficiency of the model, the Area Under the Curve (AUC) was computed using five-fold cross-validation in an external dataset.

RESULTS

The main biological pathways were the cell cycle, apoptosis, and extracellular matrix (ECM) organization. Direct association between the cell cycle and apoptosis, inverse association between apoptosis and ECM organization, and inverse association between the cell cycle and ECM organization were detected for the 5-hmC modules in GI cancers. An AUC of 92% (0.73-1.00) was observed for the predictive model including 11 genes.

CONCLUSION

The intricate association between biological pathways of identified modules may reveal the hidden significance of 5-hmC in GI cancers. The identified predictive model and new biomarkers may be beneficial in cancer detection and precision medicine using liquid biopsy in the early stages.

Large-scale k-mer-based analysis of the informational properties of genomes, comparative genomics and taxonomy

Information theoretic approaches are ubiquitous and effective in a wide variety of bioinformatics applications. In comparative genomics, alignment-free methods, based on short DNA words, or k-mers, are particularly powerful. We evaluated the utility of varying k-mer lengths for genome comparisons by analyzing their sequence space coverage of 5805 genomes in the KEGG GENOME database. In subsequent analyses on four k-mer lengths spanning the relevant range (11, 21, 31, 41), hierarchical clustering of 1634 genus-level representative genomes using pairwise 21- and 31-mer Jaccard similarities best recapitulated a phylogenetic/taxonomic tree of life with clear boundaries for superkingdom domains and high subtree similarity for named taxons at lower levels (family through phylum). By analyzing ~14.2M prokaryotic genome comparisons by their lowest-common-ancestor taxon levels, we detected many potential misclassification errors in a curated database, further demonstrating the need for wide-scale adoption of quantitative taxonomic classifications based on whole-genome similarity.

A rotamer relay information system in the epidermal growth factor receptor-drug complexes reveals clues to new paradigm in protein conformational change

Cancer cells can escape the effects of chemotherapy through mutations and upregulation of a tyrosine kinase protein called the epidermal growth factor receptor (EGFR). In the past two decades, four generations of tyrosine kinase inhibitors targeting EGFR have been developed. Using comparative structure analysis of 116 EGFR-drug complex crystal structures, cluster analysis produces two clans of 73 and 43 structures, respectively. The first clan of 73 structures is larger and is comprised mostly of the C-helix-IN conformation while the second clan of 43 structures correlates with the C-helix-OUT conformation. A deep rotamer analysis identifies 43 residues (18%) of the total of 237 residues spanning the kinase structures under investigation with significant rotamer variations between the C-helix-IN and C-helix-OUT clans. The locations of these rotamer variations take on the appearance of side chain conformational relays extending out from points of EGFR mutation to different regions of the EGFR kinase. Accordingly, we propose that key EGFR mutations act singly or together to induce drug resistant conformational changes in EGFR that are communicated via these side chain conformational relays. Accordingly, these side chain conformational relays appear to play a significant role in the development of tumour resistance. This phenomenon also suggests a new paradigm in protein conformational change that is mediated by supportive relays of rotamers on the protein surface, rather than through conventional backbone movements.

A propagation-based seed-centric local community detection for multilayer environment: The case study of colon adenocarcinoma

Regardless of all efforts on community discovery algorithms, it is still an open and challenging subject in network science. Recognizing communities in a multilayer network, where there are several layers (types) of connections, is even more complicated. Here, we concentrated on a specific type of communities called seed-centric local communities in the multilayer environment and developed a novel method based on the information cascade concept, called PLCDM. Our simulations on three datasets (real and artificial) signify that the suggested method outstrips two known earlier seed-centric local methods. Additionally, we compared it with other global multilayer and single-layer methods. Eventually, we applied our method on a biological two-layer network of Colon Adenocarcinoma (COAD), reconstructed from transcriptomic and post-transcriptomic datasets, and assessed the output modules. The functional enrichment consequences infer that the modules of interest hold biomolecules involved in the pathways associated with the carcinogenesis.

Identification of anticancer drug target genes using an outside competitive dynamics model on cancer signaling networks

Each cancer type has its own molecular signaling network. Analyzing the dynamics of molecular signaling networks can provide useful information for identifying drug target genes. In the present study, we consider an on-network dynamics model—the outside competitive dynamics model—wherein an inside leader and an opponent competitor outside the system have fixed and different states, and each normal agent adjusts its state according to a distributed consensus protocol. If any normal agent links to the external competitor, the state of each normal agent will converge to a stable value, indicating support to the leader against the impact of the competitor. We determined the total support of normal agents to each leader in various networks and observed that the total support correlates with hierarchical closeness, which identifies biomarker genes in a cancer signaling network. Of note, by experimenting on 17 cancer signaling networks from the KEGG database, we observed that 82% of the genes among the top 3 agents with the highest total support are anticancer drug target genes. This result outperforms those of four previous prediction methods of common cancer drug targets. Our study indicates that driver agents with high support from the other agents against the impact of the external opponent agent are most likely to be anticancer drug target genes.

Nodes with the highest control power play an important role at the final level of cooperation in directed networks

Controllability of complex networks aims to seek the lowest number of nodes (the driver nodes) that can control all the nodes by receiving the input signals. The concept of control centrality is used to determine the power of each node to control the network. The more a node controls the nodes through connections in the network, the more it has the power to control. Although the cooperative and free-rider strategies and the final level of cooperation in a population are considered and studied in the public goods game. However, it is yet to determine a solution to indicate the effectiveness of each member in changing the strategies of the other members. In a network, the choice of nodes effective in changing the other nodes' strategies, as free-riders, will lead to lower cooperation and vice versa. This paper uses simulated and real networks to investigate that the nodes with the highest control power are more effective than the hubs, local, and random nodes in changing the strategies of the other nodes and the final level of cooperation. Results indicate that the nodes with the highest control power as free-riders, compared to the other sets being under consideration, can lead to a lower level of cooperation and are, therefore, more effective in changing the strategies of the other nodes. The obtained results can be considered in the treatment of cancer. So that, destroying the tumoral cells with the highest control power should be a priority as these cells have a higher capability to change the strategies of the other cells from cooperators to free-riders (healthy to tumoral).

Cell-to-cell diversification in ERBB-RAS-MAPK signal transduction that produces cell-type specific growth factor responses

Growth factors regulate cell fates, including their proliferation, differentiation, survival, and death, according to the cell type. Even when the response to a specific growth factor is deterministic for collective cell behavior, significant levels of fluctuation are often observed between single cells. Statistical analyses of single-cell responses provide insights into the mechanism of cell fate decisions but very little is known about the distributions of the internal states of cells responding to growth factors. Using multi-color immunofluorescent staining, we have here detected the phosphorylation of seven elements in the early response of the ERBB-RAS-MAPK system to two growth factors. Among these seven elements, five were analyzed simultaneously in distinct combinations in the same single cells. Although principle component analysis suggested cell-type and input specific phosphorylation patterns, cell-to-cell fluctuation was large. Mutual information analysis suggested that each cell type uses multitrack (bush-like) signal transduction pathways under conditions in which clear fate changes have been reported. The clustering of single-cell response patterns indicated that the fate change in a cell population correlates with the large entropy of the response, suggesting a bet-hedging strategy is used in decision making. A comparison of true and randomized datasets further indicated that this large variation is not produced by simple reaction noise, but is defined by the properties of the signal-processing network.

Computation in bacterial communities

Bacteria across many scales are involved in a dynamic process of information exchange to coordinate activity and community structure within large and diverse populations. The molecular components bacteria use to communicate have been discovered and characterized, and recent efforts have begun to understand the potential for bacterial signal exchange to gather information from the environment and coordinate collective behaviors. Such computations made by bacteria to coordinate the action of a population of cells in response to information gathered by a multitude of inputs is a form of collective intelligence. These computations must be robust to fluctuations in both biological, chemical, and physical parameters as well as to operate with energetic efficiency. Given these constraints, what are the limits of computation by bacterial populations and what strategies have evolved to ensure bacterial communities efficiently work together? Here the current understanding of information exchange and collective decision making that occur in microbial populations will be reviewed. Looking toward the future, we consider how a deeper understanding of bacterial computation will inform future direction in microbiology, biotechnology, and biophysics.

Organ System Network Disruption Is Associated With Poor Prognosis in Patients With Chronic Liver Failure

BACKGROUND

A healthy individual has a high degree of functional connectivity between organ systems, which can be represented graphically in a network map. Disruption of this system connectivity is associated with mortality in life-threatening acute illnesses, demonstrated by a network approach. However, this approach has not been applied to chronic multisystem diseases and may be more reliable than conventional individual organ prognostic scoring methods. Cirrhosis is a chronic disease of the liver with multisystem involvement. Development of an efficient model for prediction of mortality in cirrhosis requires a profound understanding of the pathophysiologic processes that lead to poor prognosis. In the present study, we use a network approach to evaluate the differences in organ system connectivity between survivors and non-survivors in a group of well-characterized patients with cirrhosis.

METHODS

201 patients with cirrhosis originally referred to the Clinic five at the University Hospital of Padova, with 13 clinical variables available representing hepatic, metabolic, haematopoietic, immune, neural and renal organ systems were retrospectively enrolled and followed up for 3, 6, and 12 months. Software was designed to compute the correlation network maps of organ system interaction in survivors and non-survivors using analysis indices: A. Bonferroni corrected Pearson's correlation coefficient and B. Mutual Information. Network structure was quantitatively evaluated using the measures of edges, average degree of connectivity and closeness, and qualitatively using clinical significance. Pair-matching was also implemented as a further step after initial general analysis to control for sample size and Model for End-Stage Liver Disease (MELD-Na) score between the groups.

RESULTS

There was a higher number of significant correlations in survivors, as indicated by both the analysis indices of Bonferroni corrected Pearson's correlation coefficient and the Mutual Information analysis. The number of edges, average degree of connectivity and closeness were significantly higher in survivors compared to non-survivors group. Pair-matching for MELD-Na was also associated with increased connectivity in survivors compared to non-survivors over 3 and 6 months follow up.

CONCLUSION

This study demonstrates the application of a network approach in evaluating functional connectivity of organ systems in liver cirrhosis, demonstrating a significant degree of network disruption in organ systems in non-survivors. Network analysis of organ systems may provide insight in developing novel prognostic models for organ allocation in patients with cirrhosis.

Information Theory in Computational Biology: Where We Stand Today

"A Mathematical Theory of Communication" was published in 1948 by Claude Shannon to address the problems in the field of data compression and communication over (noisy) communication channels. Since then, the concepts and ideas developed in Shannon's work have formed the basis of information theory, a cornerstone of statistical learning and inference, and has been playing a key role in disciplines such as physics and thermodynamics, probability and statistics, computational sciences and biological sciences. In this article we review the basic information theory based concepts and describe their key applications in multiple major areas of research in computational biology-gene expression and transcriptomics, alignment-free sequence comparison, sequencing and error correction, genome-wide disease-gene association mapping, metabolic networks and metabolomics, and protein sequence, structure and interaction analysis.

Biological Network Approaches and Applications in Rare Disease Studies

Network biology has the capability to integrate, represent, interpret, and model complex biological systems by collectively accommodating biological omics data, biological interactions and associations, graph theory, statistical measures, and visualizations. Biological networks have recently been shown to be very useful for studies that decipher biological mechanisms and disease etiologies and for studies that predict therapeutic responses, at both the molecular and system levels. In this review, we briefly summarize the general framework of biological network studies, including data resources, network construction methods, statistical measures, network topological properties, and visualization tools. We also introduce several recent biological network applications and methods for the studies of rare diseases.

On signalling and estimation limits for molecular birth-processes

Understanding and uncovering the mechanisms or motifs that molecular networks employ to regulate noise is a key problem in cell biology. As it is often difficult to obtain direct and detailed insight into these mechanisms, many studies instead focus on assessing the best precision attainable on the signalling pathways that compose these networks. Molecules signal one another over such pathways to solve noise regulating estimation and control problems. Quantifying the maximum precision of these solutions delimits what is achievable and allows hypotheses about underlying motifs to be tested without requiring detailed biological knowledge. The pathway capacity, which defines the maximum rate of transmitting information along it, is a widely used proxy for precision. Here it is shown, for estimation problems involving elementary yet biologically relevant birth-process networks, that capacity can be surprisingly misleading. A time-optimal signalling motif, called birth-following, is derived and proven to better the precision expected from the capacity, provided the maximum signalling rate constraint is large and the mean one above a certain threshold. When the maximum constraint is relaxed, perfect estimation is predicted by the capacity. However, the true achievable precision is found highly variable and sensitive to the mean constraint. Since the same capacity can map to different combinations of rate constraints, it can only equivocally measure precision. Deciphering the rate constraints on a signalling pathway may therefore be more important than computing its capacity.

Exploring the information transmission properties of noise-induced dynamics: application to glioma differentiation

Background

Cells operate in an uncertain environment, where critical cell decisions must be enacted in the presence of biochemical noise. Information theory can measure the extent to which such noise perturbs normal cellular function, in which cells must perceive environmental cues and relay signals accurately to make timely and informed decisions. Using multivariate response data can greatly improve estimates of the latent information content underlying important cell fates, like differentiation.

Results

We undertake an information theoretic analysis of two stochastic models concerning glioma differentiation therapy, an alternative cancer treatment modality whose underlying intracellular mechanisms remain poorly understood. Discernible changes in response dynamics, as captured by summary measures, were observed at low noise levels. Mitigating certain feedback mechanisms present in the signaling network improved information transmission overall, as did targeted subsampling and clustering of response dynamics.

Conclusion

Computing the channel capacity of noisy signaling pathways present great probative value in uncovering the prevalent trends in noise-induced dynamics. Areas of high dynamical variation can provide concise snapshots of informative system behavior that may otherwise be overlooked. Through this approach, we can examine the delicate interplay between noise and information, from signal to response, through the observed behavior of relevant system components.

Electronic supplementary material

The online version of this article (10.1186/s12859-019-2970-7) contains supplementary material, which is available to authorized users.

Do cells use passwords in cell-state transitions? Is cell signaling sometimes encrypted?

Organisms must maintain proper regulation including defense and healing. Life-threatening problems may be caused by pathogens or by a multicellular organism's own cells through cancer or autoimmune disorders. Life evolved solutions to these problems that can be conceptualized through the lens of information security, which is a well-developed field in computer science. Here I argue that taking an information security view of cells is not merely semantics, but useful to explain features of signaling, regulation, and defense. An information security perspective also offers a conduit for cross-fertilization of advanced ideas from computer science and the potential for biology to inform computer science. First, I consider whether cells use passwords, i.e., initiation sequences that are required for subsequent signals to have effects, by analyzing the concept of pioneer transcription factors in chromatin regulation and cellular reprogramming. Second, I consider whether cells may encrypt signal transduction cascades. Encryption could benefit cells by making it more difficult for pathogens or oncogenes to hijack cell networks. By using numerous molecules, cells may gain a security advantage in particular against viruses, whose genome sizes are typically under selection pressure. I provide a simple conceptual argument for how cells may perform encryption through posttranslational modifications, complex formation, and chromatin accessibility. I invoke information theory to provide a criterion of an entropy spike to assess whether a signaling cascade has encryption-like features. I discuss how the frequently invoked concept of context dependency may oversimplify more advanced features of cell signaling networks, such as encryption. Therefore, by considering that biochemical networks may be even more complex than commonly realized we may be better able to understand defenses against pathogens and pathologies.

Robustness and Information Transfer within IL-6-induced JAK/STAT Signalling

Cellular communication via intracellular signalling pathways is crucial. Expression and activation of signalling proteins is heterogenous between isogenic cells of the same cell-type. However, mechanisms evolved to enable sufficient communication and to ensure cellular functions. We use information theory to clarify mechanisms facilitating IL-6-induced JAK/STAT signalling despite cell-to-cell variability. We show that different mechanisms enabling robustness against variability complement each other. Early STAT3 activation is robust as long as cytokine concentrations are low. Robustness at high cytokine concentrations is ensured by high STAT3 expression or serine phosphorylation. Later the feedback-inhibitor SOCS3 increases robustness. Channel Capacity of JAK/STAT signalling is limited by cell-to-cell variability in STAT3 expression and is affected by the same mechanisms governing robustness. Increasing STAT3 amount increases Channel Capacity and robustness, whereas increasing STAT3 tyrosine phosphorylation reduces robustness but increases Channel Capacity. In summary, we elucidate mechanisms preventing dysregulated signalling by enabling reliable JAK/STAT signalling despite cell-to-cell heterogeneity.

Symmetry and symmetry breaking in cancer: a foundational approach to the cancer problem

Symmetry and symmetry breaking concepts from physics and biology are applied to the problem of cancer. Three categories of symmetry breaking in cancer are examined: combinatorial, geometric, and functional. Within these categories, symmetry breaking is examined for relevant cancer features, including epithelial-mesenchymal transition (EMT); tumor heterogeneity; tensegrity; fractal geometric and information structure; functional interaction networks; and network stabilizability and attack tolerance. The new cancer symmetry concepts are relevant to homeostasis loss in cancer and to its origin, spread, treatment and resistance. Symmetry and symmetry breaking could provide a new way of thinking and a pathway to a solution of the cancer problem.

From protein-protein interactions to protein co-expression networks: a new perspective to evaluate large-scale proteomic data

The reductionist approach of dissecting biological systems into their constituents has been successful in the first stage of the molecular biology to elucidate the chemical basis of several biological processes. This knowledge helped biologists to understand the complexity of the biological systems evidencing that most biological functions do not arise from individual molecules; thus, realizing that the emergent properties of the biological systems cannot be explained or be predicted by investigating individual molecules without taking into consideration their relations. Thanks to the improvement of the current -omics technologies and the increasing understanding of the molecular relationships, even more studies are evaluating the biological systems through approaches based on graph theory. Genomic and proteomic data are often combined with protein-protein interaction (PPI) networks whose structure is routinely analyzed by algorithms and tools to characterize hubs/bottlenecks and topological, functional, and disease modules. On the other hand, co-expression networks represent a complementary procedure that give the opportunity to evaluate at system level including organisms that lack information on PPIs. Based on these premises, we introduce the reader to the PPI and to the co-expression networks, including aspects of reconstruction and analysis. In particular, the new idea to evaluate large-scale proteomic data by means of co-expression networks will be discussed presenting some examples of application. Their use to infer biological knowledge will be shown, and a special attention will be devoted to the topological and module analysis.