Topological estimation of signal flow in complex signaling networks

Control of Cellular Differentiation Trajectories for Cancer Reversion

Cellular differentiation is controlled by intricate layers of gene regulation, involving the modulation of gene expression by various transcriptional regulators. Due to the complexity of gene regulation, identifying master regulators across the differentiation trajectory has been a longstanding challenge. To tackle this problem, a computational framework, single-cell Boolean network inference and control (BENEIN), is presented. Applying BENEIN to human large intestinal single-cell transcriptome data, MYB, HDAC2, and FOXA2 are identified as the master regulators whose inhibition induces enterocyte differentiation. It is found that simultaneous knockdown of these master regulators can revert colorectal cancer cells into normal-like enterocytes by synergistically inducing differentiation and suppressing malignancy, which is validated by in vitro and in vivo experiments.

A gray box framework that optimizes a white box logical model using a black box optimizer for simulating cellular responses to perturbations

Predicting cellular responses to perturbations requires interpretable insights into molecular regulatory dynamics to perform reliable cell fate control, despite the confounding non-linearity of the underlying interactions. There is a growing interest in developing machine learning-based perturbation response prediction models to handle the non-linearity of perturbation data, but their interpretation in terms of molecular regulatory dynamics remains a challenge. Alternatively, for meaningful biological interpretation, logical network models such as Boolean networks are widely used in systems biology to represent intracellular molecular regulation. However, determining the appropriate regulatory logic of large-scale networks remains an obstacle due to the high-dimensional and discontinuous search space. To tackle these challenges, we present a scalable derivative-free optimizer trained by meta-reinforcement learning for Boolean network models. The logical network model optimized by the trained optimizer successfully predicts anti-cancer drug responses of cancer cell lines, while simultaneously providing insight into their underlying molecular regulatory mechanisms.

Reducing State Conflicts between Network Motifs Synergistically Enhances Cancer Drug Effects and Overcomes Adaptive Resistance

Inducing apoptosis in cancer cells is a primary goal in anti-cancer therapy, but curing cancer with a single drug is unattainable due to drug resistance. The complex molecular network in cancer cells causes heterogeneous responses to single-target drugs, thereby inducing an adaptive drug response. Here, we showed that targeted drug perturbations can trigger state conflicts between multi-stable motifs within a molecular regulatory network, resulting in heterogeneous drug responses. However, we revealed that properly regulating an interconnecting molecule between these motifs can synergistically minimize the heterogeneous responses and overcome drug resistance. We extracted the essential cellular response dynamics of the Boolean network driven by the target node perturbation and developed an algorithm to identify a synergistic combinatorial target that can reduce heterogeneous drug responses. We validated the proposed approach using exemplary network models and a gastric cancer model from a previous study by showing that the targets identified with our algorithm can better drive the networks to desired states than those with other control theories. Of note, our approach suggests a new synergistic pair of control targets that can increase cancer drug efficacy to overcome adaptive drug resistance.

From time-series transcriptomics to gene regulatory networks: A review on inference methods

Inference of gene regulatory networks has been an active area of research for around 20 years, leading to the development of sophisticated inference algorithms based on a variety of assumptions and approaches. With the ever increasing demand for more accurate and powerful models, the inference problem remains of broad scientific interest. The abstract representation of biological systems through gene regulatory networks represents a powerful method to study such systems, encoding different amounts and types of information. In this review, we summarize the different types of inference algorithms specifically based on time-series transcriptomics, giving an overview of the main applications of gene regulatory networks in computational biology. This review is intended to give an updated reference of regulatory networks inference tools to biologists and researchers new to the topic and guide them in selecting the appropriate inference method that best fits their questions, aims, and experimental data.

Cell-Type-Specific Gene Regulatory Networks of Pro-Inflammatory and Pro-Resolving Lipid Mediator Biosynthesis in the Immune System

Lipid mediators are important regulators in inflammatory responses, and their biosynthetic pathways are targeted by commonly used anti-inflammatory drugs. Switching from pro-inflammatory lipid mediators (PIMs) to specialized pro-resolving (SPMs) is a critical step toward acute inflammation resolution and preventing chronic inflammation. Although the biosynthetic pathways and enzymes for PIMs and SPMs have now been largely identified, the actual transcriptional profiles underlying the immune cell type-specific transcriptional profiles of these mediators are still unknown. Using the Atlas of Inflammation Resolution, we created a large network of gene regulatory interactions linked to the biosynthesis of SPMs and PIMs. By mapping single-cell sequencing data, we identified cell type-specific gene regulatory networks of the lipid mediator biosynthesis. Using machine learning approaches combined with network features, we identified cell clusters of similar transcriptional regulation and demonstrated how specific immune cell activation affects PIM and SPM profiles. We found substantial differences in regulatory networks in related cells, accounting for network-based preprocessing in functional single-cell analyses. Our results not only provide further insight into the gene regulation of lipid mediators in the immune response but also shed light on the contribution of selected cell types in their biosynthesis.

NETISCE: a network-based tool for cell fate reprogramming

The search for effective therapeutic targets in fields like regenerative medicine and cancer research has generated interest in cell fate reprogramming. This cellular reprogramming paradigm can drive cells to a desired target state from any initial state. However, methods for identifying reprogramming targets remain limited for biological systems that lack large sets of experimental data or a dynamical characterization. We present NETISCE, a novel computational tool for identifying cell fate reprogramming targets in static networks. In combination with machine learning algorithms, NETISCE estimates the attractor landscape and predicts reprogramming targets using signal flow analysis and feedback vertex set control, respectively. Through validations in studies of cell fate reprogramming from developmental, stem cell, and cancer biology, we show that NETISCE can predict previously identified cell fate reprogramming targets and identify potentially novel combinations of targets. NETISCE extends cell fate reprogramming studies to larger-scale biological networks without the need for full model parameterization and can be implemented by experimental and computational biologists to identify parts of a biological system relevant to the desired reprogramming task.

Identifying molecular targets for reverse aging using integrated network analysis of transcriptomic and epigenomic changes during aging

Aging is associated with widespread physiological changes, including skeletal muscle weakening, neuron system degeneration, hair loss, and skin wrinkling. Previous studies have identified numerous molecular biomarkers involved in these changes, but their regulatory mechanisms and functional repercussions remain elusive. In this study, we conducted next-generation sequencing of DNA methylation and RNA sequencing of blood samples from 51 healthy adults between 20 and 74 years of age and identified aging-related epigenetic and transcriptomic biomarkers. We also identified candidate molecular targets that can reversely regulate the transcriptomic biomarkers of aging by reconstructing a gene regulatory network model and performing signal flow analysis. For validation, we screened public experimental data including gene expression profiles in response to thousands of chemical perturbagens. Despite insufficient data on the binding targets of perturbagens and their modes of action, curcumin, which reversely regulated the biomarkers in the experimental dataset, was found to bind and inhibit JUN, which was identified as a candidate target via signal flow analysis. Collectively, our results demonstrate the utility of a network model for integrative analysis of omics data, which can help elucidate inter-omics regulatory mechanisms and develop therapeutic strategies against aging.

OCSANA+: optimal control and simulation of signaling networks from network analysis

SUMMARY

OCSANA+ is a Cytoscape app for identifying nodes to drive the system toward a desired long-term behavior, prioritizing combinations of interventions in large-scale complex networks, and estimating the effects of node perturbations in signaling networks, all based on the analysis of the network's structure. OCSANA+ includes an update to optimal combinations of interventions from network analysis software tool with cutting-edge and rigorously tested algorithms, together with recently developed structure-based control algorithms for non-linear systems and an algorithm for estimating signal flow. All these algorithms are based on the network's topology. OCSANA+ is implemented as a Cytoscape app to enable a user interface for running analyses and visualizing results.

AVAILABILITY AND IMPLEMENTATION

OCSANA+ app and its tutorial can be downloaded from the Cytoscape App Store or https://veraliconaresearchgroup.github.io/OCSANA-Plus/. The source code and computations are available in https://github.com/VeraLiconaResearchGroup/OCSANA-Plus_SourceCode.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Signal flow control of complex signaling networks

Complex disease such as cancer is often caused by genetic mutations that eventually alter the signal flow in the intra-cellular signaling network and result in different cell fate. Therefore, it is crucial to identify control targets that can most effectively block such unwanted signal flow. For this purpose, systems biological analysis provides a useful framework, but mathematical modeling of complicated signaling networks requires massive time-series measurements of signaling protein activity levels for accurate estimation of kinetic parameter values or regulatory logics. Here, we present a novel method, called SFC (Signal Flow Control), for identifying control targets without the information of kinetic parameter values or regulatory logics. Our method requires only the structural information of a signaling network and is based on the topological estimation of signal flow through the network. SFC will be particularly useful for a large-scale signaling network to which parameter estimation or inference of regulatory logics is no longer applicable in practice. The identified control targets have significant implication in drug development as they can be putative drug targets.