Adaptation and learning of molecular networks as a description of cancer development at the systems-level: Potential use in anti-cancer therapies

Collective intelligence: A unifying concept for integrating biology across scales and substrates

A defining feature of biology is the use of a multiscale architecture, ranging from molecular networks to cells, tissues, organs, whole bodies, and swarms. Crucially however, biology is not only nested structurally, but also functionally: each level is able to solve problems in distinct problem spaces, such as physiological, morphological, and behavioral state space. Percolating adaptive functionality from one level of competent subunits to a higher functional level of organization requires collective dynamics: multiple components must work together to achieve specific outcomes. Here we overview a number of biological examples at different scales which highlight the ability of cellular material to make decisions that implement cooperation toward specific homeodynamic endpoints, and implement collective intelligence by solving problems at the cell, tissue, and whole-organism levels. We explore the hypothesis that collective intelligence is not only the province of groups of animals, and that an important symmetry exists between the behavioral science of swarms and the competencies of cells and other biological systems at different scales. We then briefly outline the implications of this approach, and the possible impact of tools from the field of diverse intelligence for regenerative medicine and synthetic bioengineering.

Cellular signaling pathways as plastic, proto-cognitive systems: Implications for biomedicine

Many aspects of health and disease are modeled using the abstraction of a "pathway"-a set of protein or other subcellular activities with specified functional linkages between them. This metaphor is a paradigmatic case of a deterministic, mechanistic framework that focuses biomedical intervention strategies on altering the members of this network or the up-/down-regulation links between them-rewiring the molecular hardware. However, protein pathways and transcriptional networks exhibit interesting and unexpected capabilities such as trainability (memory) and information processing in a context-sensitive manner. Specifically, they may be amenable to manipulation via their history of stimuli (equivalent to experiences in behavioral science). If true, this would enable a new class of biomedical interventions that target aspects of the dynamic physiological "software" implemented by pathways and gene-regulatory networks. Here, we briefly review clinical and laboratory data that show how high-level cognitive inputs and mechanistic pathway modulation interact to determine outcomes in vivo. Further, we propose an expanded view of pathways from the perspective of basal cognition and argue that a broader understanding of pathways and how they process contextual information across scales will catalyze progress in many areas of physiology and neurobiology. We argue that this fuller understanding of the functionality and tractability of pathways must go beyond a focus on the mechanistic details of protein and drug structure to encompass their physiological history as well as their embedding within higher levels of organization in the organism, with numerous implications for data science addressing health and disease. Exploiting tools and concepts from behavioral and cognitive sciences to explore a proto-cognitive metaphor for the pathways underlying health and disease is more than a philosophical stance on biochemical processes; at stake is a new roadmap for overcoming the limitations of today's pharmacological strategies and for inferring future therapeutic interventions for a wide range of disease states.

Spatial-Temporal Genome Regulation in Stress-Response and Cell-Fate Change

Complex functioning of the genome in the cell nucleus is controlled at different levels: (a) the DNA base sequence containing all relevant inherited information; (b) epigenetic pathways consisting of protein interactions and feedback loops; (c) the genome architecture and organization activating or suppressing genetic interactions between different parts of the genome. Most research so far has shed light on the puzzle pieces at these levels. This article, however, attempts an integrative approach to genome expression regulation incorporating these different layers. Under environmental stress or during cell development, differentiation towards specialized cell types, or to dysfunctional tumor, the cell nucleus seems to react as a whole through coordinated changes at all levels of control. This implies the need for a framework in which biological, chemical, and physical manifestations can serve as a basis for a coherent theory of gene self-organization. An international symposium held at the Biomedical Research and Study Center in Riga, Latvia, on 25 July 2022 addressed novel aspects of the abovementioned topic. The present article reviews the most recent results and conclusions of the state-of-the-art research in this multidisciplinary field of science, which were delivered and discussed by scholars at the Riga symposium.

The Transcriptome and Proteome Networks of Malignant Tumours Reveal Atavistic Attractors of Polyploidy-Related Asexual Reproduction

The expression of gametogenesis-related (GG) genes and proteins, as well as whole genome duplications (WGD), are the hallmarks of cancer related to poor prognosis. Currently, it is not clear if these hallmarks are random processes associated only with genome instability or are programmatically linked. Our goal was to elucidate this via a thorough bioinformatics analysis of 1474 GG genes in the context of WGD. We examined their association in protein-protein interaction and coexpression networks, and their phylostratigraphic profiles from publicly available patient tumour data. The results show that GG genes are upregulated in most WGD-enriched somatic cancers at the transcriptome level and reveal robust GG gene expression at the protein level, as well as the ability to associate into correlation networks and enrich the reproductive modules. GG gene phylostratigraphy displayed in WGD+ cancers an attractor of early eukaryotic origin for DNA recombination and meiosis, and one relative to oocyte maturation and embryogenesis from early multicellular organisms. The upregulation of cancer-testis genes emerging with mammalian placentation was also associated with WGD. In general, the results suggest the role of polyploidy for soma-germ transition accessing latent cancer attractors in the human genome network, which appear as pre-formed along the whole Evolution of Life.

An integrated network representation of multiple cancer-specific data for graph-based machine learning

Genomic profiles of cancer cells provide valuable information on genetic alterations in cancer. Several recent studies employed these data to predict the response of cancer cell lines to drug treatment. Nonetheless, due to the multifactorial phenotypes and intricate mechanisms of cancer, the accurate prediction of the effect of pharmacotherapy on a specific cell line based on the genetic information alone is problematic. Emphasizing on the system-level complexity of cancer, we devised a procedure to integrate multiple heterogeneous data, including biological networks, genomics, inhibitor profiling, and gene-disease associations, into a unified graph structure. In order to construct compact, yet information-rich cancer-specific networks, we developed a novel graph reduction algorithm. Driven by not only the topological information, but also the biological knowledge, the graph reduction increases the feature-only entropy while preserving the valuable graph-feature information. Subsequent comparative benchmarking simulations employing a tissue level cross-validation protocol demonstrate that the accuracy of a graph-based predictor of the drug efficacy is 0.68, which is notably higher than those measured for more traditional, matrix-based techniques on the same data. Overall, the non-Euclidean representation of the cancer-specific data improves the performance of machine learning to predict the response of cancer to pharmacotherapy. The generated data are freely available to the academic community at https://osf.io/dzx7b/.

Cancer progression as a learning process

Drug resistance and metastasis-the major complications in cancer-both entail adaptation of cancer cells to stress, whether a drug or a lethal new environment. Intriguingly, these adaptive processes share similar features that cannot be explained by a pure Darwinian scheme, including dormancy, increased heterogeneity, and stress-induced plasticity. Here, we propose that learning theory offers a framework to explain these features and may shed light on these two intricate processes. In this framework, learning is performed at the single-cell level, by stress-driven exploratory trial-and-error. Such a process is not contingent on pre-existing pathways but on a random search for a state that diminishes the stress. We review underlying mechanisms that may support this search, and show by using a learning model that such exploratory learning is feasible in a high-dimensional system as the cell. At the population level, we view the tissue as a network of exploring agents that communicate, restraining cancer formation in health. In this view, disease results from the breakdown of homeostasis between cellular exploratory drive and tissue homeostasis.

Overexpression of NREP Promotes Migration and Invasion in Gastric Cancer Through Facilitating Epithelial-Mesenchymal Transition

The identification of biomarkers and effective therapeutic targets for gastric cancer (GC), the most common cause of cancer-related deaths around the world, is currently a major focus area in research. Here, we examined the utility of Neuronal Regeneration Related Protein (NREP) as a prognostic biomarker and therapeutic target for GC. We assessed the clinical relevance, function, and molecular role of NREP in GC using bioinformatics analysis and experimental validation. Our results showed that in GC, NREP overexpression was significantly associated with a poor prognosis. Our findings also suggested that NREP may be involved in the activation of cancer-associated fibroblasts and the epithelial-mesenchymal transition (EMT), with transforming growth factor β1 mediating both processes. In addition, NREP expression showed a positive correlation with the abundance of M2 macrophages, which are potent immunosuppressors. Together, these results indicate that NREP is overexpressed in GC and affects GC prognosis. Thus, NREP could be a prognostic biomarker and therapeutic target for GC.

Cell Differentiation Trajectory-Associated Molecular Classification of Osteosarcoma

This study aims to investigate the differentiation trajectory of osteosarcoma cells and to construct molecular subtypes with their respective characteristics and generate a multi-gene signature for predicting prognosis. Integrated single-cell RNA-sequencing (scRNA-seq) data, bulk RNA-seq data and microarray data from osteosarcoma samples were used for analysis. Via scRNA-seq data, time-related as well as differentiation-related genes were recognized as osteosarcoma tumor stem cell-related genes (OSCGs). In Gene Expression Omnibus (GEO) cohort, osteosarcoma patients were classified into two subtypes based on prognostic OSCGs and it was found that molecular typing successfully predicted overall survival, tumor microenvironment and immune infiltration status. Further, available drugs for influencing osteosarcoma via prognostic OSCGs were revealed. A 3-OSCG-based prognostic risk score signature was generated and by combining other clinic-pathological independent prognostic factor, stage at diagnosis, a nomogram was established to predict individual survival probability. In external independent TARGET cohort, the molecular types, the 3-gene signature as well as nomogram were validated. In conclusion, osteosarcoma cell differentiation occupies a crucial position in many facets, such as tumor prognosis and microenvironment, suggesting promising therapeutic targets for this disease.

Using proteomic and transcriptomic data to assess activation of intracellular molecular pathways

Analysis of molecular pathway activation is the recent instrument that helps to quantize activities of various intracellular signaling, structural, DNA synthesis and repair, and biochemical processes. This may have a deep impact in fundamental research, bioindustry, and medicine. Unlike gene ontology analyses and numerous qualitative methods that can establish whether a pathway is affected in principle, the quantitative approach has the advantage of exactly measuring the extent of a pathway up/downregulation. This results in emergence of a new generation of molecular biomarkers-pathway activation levels, which reflect concentration changes of all measurable pathway components. The input data can be the high-throughput proteomic or transcriptomic profiles, and the output numbers take both positive and negative values and positively reflect overall pathway activation. Due to their nature, the pathway activation levels are more robust biomarkers compared to the individual gene products/protein levels. Here, we review the current knowledge of the quantitative gene expression interrogation methods and their applications for the molecular pathway quantization. We consider enclosed bioinformatic algorithms and their applications for solving real-world problems. Besides a plethora of applications in basic life sciences, the quantitative pathway analysis can improve molecular design and clinical investigations in pharmaceutical industry, can help finding new active biotechnological components and can significantly contribute to the progressive evolution of personalized medicine. In addition to the theoretical principles and concepts, we also propose publicly available software for the use of large-scale protein/RNA expression data to assess the human pathway activation levels.

GraphGR: A graph neural network to predict the effect of pharmacotherapy on the cancer cell growth.. [DOI: 10.1101/2020.05.20.107458]

Genomic profiles of cancer cells provide valuable information on genetic alterations in cancer. Several recent studies employed these data to predict the response of cancer cell lines to treatment with drugs. Nonetheless, due to the multifactorial phenotypes and intricate mechanisms of cancer, the accurate prediction of the effect of pharmacotherapy on a specific cell line based on the genetic information alone is problematic. High prediction accuracies reported in the literature likely result from significant overlaps among training, validation, and testing sets, making many predictors inapplicable to new data. To address these issues, we developed GraphGR, a graph neural network with sophisticated attention propagation mechanisms to predict the therapeutic effects of kinase inhibitors across various tumors. Emphasizing on the system-level complexity of cancer, GraphGR integrates multiple heterogeneous data, such as biological networks, genomics, inhibitor profiling, and genedisease associations, into a unified graph structure. In order to construct diverse and information-rich cancer-specific networks, we devised a novel graph reduction protocol based on not only the topological information, but also the biological knowledge. The performance of GraphGR, properly cross-validated at the tissue level, is 0.83 in terms of the area under the receiver operating characteristics, which is notably higher than those measured for other approaches on the same data. Finally, several new predictions are validated against the biomedical literature demonstrating that GraphGR generalizes well to unseen data, i.e. it can predict therapeutic effects across a variety of cancer cell lines and inhibitors. GraphGR is freely available to the academic community at https://github.com/pulimeng/GraphGR.

Role of miR-142 in the pathogenesis of osteosarcoma and its potential as therapeutic approach

Osteosarcoma (OS) is the most common primary malignant tumor of the bone with a strong tendency to early metastasis, and occurs in growing bones more commonly in children and adolescents. Considering the limited therapeutic methods and lack of 100% success of these methods, developing innovative therapies with high efficacy and lower side effects is needed. Meanwhile, miRNAs and the studies indicating the involvement of miRNAs in OS development have attracted attentions as a result of the frequent abnormalities in expression of miRNAs in cancer. miRNAs are noncoding short sequences with lengths ranging from 18 to 25 nucleotides that play a very important role in cellular processes, such as proliferation, differentiation, migration, and apoptosis. MiRNAs can have either oncogenic or tumor suppressive role based on cellular function and targets. This review aimed to have overview on miR-142 as a tumor suppressor in OS. Moreover, the genes involved in the disease, such as RAC1, HMAG1, MMP9, MMP2, and E-cadherin, which have irregularities as a result of change in miR-142 expression, and, thereby, result in increasing the proliferation, invasion, and metastasis of the cells in the tissues and OS cells will be discussed.

Network dynamics: quantitative analysis of complex behavior in metabolism, organelles, and cells, from experiments to models and back

Advancing from two core traits of biological systems: multilevel network organization and nonlinearity, we review a host of novel and readily available techniques to explore and analyze their complex dynamic behavior within the framework of experimental-computational synergy. In the context of concrete biological examples, analytical methods such as wavelet, power spectra, and metabolomics-fluxomics analyses, are presented, discussed, and their strengths and limitations highlighted. Further shown is how time series from stationary and nonstationary biological variables and signals, such as membrane potential, high-throughput metabolomics, O₂ and CO₂ levels, bird locomotion, at the molecular, (sub)cellular, tissue, and whole organ and animal levels, can reveal important information on the properties of the underlying biological networks. Systems biology-inspired computational methods start to pave the way for addressing the integrated functional dynamics of metabolic, organelle and organ networks. As our capacity to unravel the control and regulatory properties of these networks and their dynamics under normal or pathological conditions broadens, so is our ability to address endogenous rhythms and clocks to improve health-span in human aging, and to manage complex metabolic disorders, neurodegeneration, and cancer. WIREs Syst Biol Med 2017, 9:e1352. doi: 10.1002/wsbm.1352 For further resources related to this article, please visit the WIREs website.

On Having No Head: Cognition throughout Biological Systems

The central nervous system (CNS) underlies memory, perception, decision-making, and behavior in numerous organisms. However, neural networks have no monopoly on the signaling functions that implement these remarkable algorithms. It is often forgotten that neurons optimized cellular signaling modes that existed long before the CNS appeared during evolution, and were used by somatic cellular networks to orchestrate physiology, embryonic development, and behavior. Many of the key dynamics that enable information processing can, in fact, be implemented by different biological hardware. This is widely exploited by organisms throughout the tree of life. Here, we review data on memory, learning, and other aspects of cognition in a range of models, including single celled organisms, plants, and tissues in animal bodies. We discuss current knowledge of the molecular mechanisms at work in these systems, and suggest several hypotheses for future investigation. The study of cognitive processes implemented in aneural contexts is a fascinating, highly interdisciplinary topic that has many implications for evolution, cell biology, regenerative medicine, computer science, and synthetic bioengineering.

Evaluating the evidence for targeting FOXO3a in breast cancer: a systematic review

Background

Tumour cells show greater dependency on glycolysis so providing a sufficient and rapid energy supply for fast growth. In many breast cancers, estrogen, progesterone and epidermal growth factor receptor-positive cells proliferate in response to growth factors and growth factor antagonists are a mainstay of treatment. However, triple negative breast cancer (TNBC) cells lack receptor expression, are frequently more aggressive and are resistant to growth factor inhibition. Downstream of growth factor receptors, signal transduction proceeds via phosphatidylinositol 3-kinase (PI3k), Akt and FOXO3a inhibition, the latter being partly responsible for coordinated increases in glycolysis and apoptosis resistance. FOXO3a may be an attractive therapeutic target for TNBC. Therefore we have undertaken a systematic review of FOXO3a as a target for breast cancer therapeutics.

Methods

Articles from NCBI were retrieved systematically when reporting primary data about FOXO3a expression in breast cancer cells after cytotoxic drug treatment.

Results

Increased FOXO3a expression is common following cytotoxic drug treatment and is associated with apoptosis and cell cycle arrest. There is some evidence that metabolic enzyme expression is also altered and that this effect is also elicited in TNBC cells. FOXO3a expression serves as a positive prognostic marker, especially in estrogen (ER) receptor positive cells.

Discussion

FOXO3a is upregulated by a number of receptor-dependent and -independent anti-cancer drugs and associates with apoptosis. The identification of microRNA that regulate FOXO3a directly suggest that it offers a tangible therapeutic target that merits wider evaluation.

Cancer stem cells display extremely large evolvability: alternating plastic and rigid networks as a potential Mechanism: network models, novel therapeutic target strategies, and the contributions of hypoxia, inflammation and cellular senescence

Cancer is increasingly perceived as a systems-level, network phenomenon. The major trend of malignant transformation can be described as a two-phase process, where an initial increase of network plasticity is followed by a decrease of plasticity at late stages of tumor development. The fluctuating intensity of stress factors, like hypoxia, inflammation and the either cooperative or hostile interactions of tumor inter-cellular networks, all increase the adaptation potential of cancer cells. This may lead to the bypass of cellular senescence, and to the development of cancer stem cells. We propose that the central tenet of cancer stem cell definition lies exactly in the indefinability of cancer stem cells. Actual properties of cancer stem cells depend on the individual "stress-history" of the given tumor. Cancer stem cells are characterized by an extremely large evolvability (i.e. a capacity to generate heritable phenotypic variation), which corresponds well with the defining hallmarks of cancer stem cells: the possession of the capacity to self-renew and to repeatedly re-build the heterogeneous lineages of cancer cells that comprise a tumor in new environments. Cancer stem cells represent a cell population, which is adapted to adapt. We argue that the high evolvability of cancer stem cells is helped by their repeated transitions between plastic (proliferative, symmetrically dividing) and rigid (quiescent, asymmetrically dividing, often more invasive) phenotypes having plastic and rigid networks. Thus, cancer stem cells reverse and replay cancer development multiple times. We describe network models potentially explaining cancer stem cell-like behavior. Finally, we propose novel strategies including combination therapies and multi-target drugs to overcome the Nietzschean dilemma of cancer stem cell targeting: "what does not kill me makes me stronger".