Boolean model of anchorage dependence and contact inhibition points to coordinated inhibition but semi-independent induction of proliferation and migration

Unlocking mitochondrial dysfunction-associated senescence (MiDAS) with NAD+ - A Boolean model of mitochondrial dynamics and cell cycle control

The steady accumulation of senescent cells with aging creates tissue environments that aid cancer evolution. Aging cell states are highly heterogeneous. 'Deep senescent' cells rely on healthy mitochondria to fuel a strong proinflammatory secretome, including cytokines, growth and transforming signals. Yet, the physiological triggers of senescence such as reactive oxygen species (ROS) can also trigger mitochondrial dysfunction, and sufficient energy deficit to alter their secretome and cause chronic oxidative stress - a state termed Mitochondrial Dysfunction-Associated Senescence (MiDAS). Here, we offer a mechanistic hypothesis for the molecular processes leading to MiDAS, along with testable predictions. To do this we have built a Boolean regulatory network model that qualitatively captures key aspects of mitochondrial dynamics during cell cycle progression (hyper-fusion at the G1/S boundary, fission in mitosis), apoptosis (fission and dysfunction) and glucose starvation (reversible hyper-fusion), as well as MiDAS in response to SIRT3 knockdown or oxidative stress. Our model reaffirms the protective role of NAD+ and external pyruvate. We offer testable predictions about the growth factor- and glucose-dependence of MiDAS and its reversibility at different stages of reactive oxygen species (ROS)-induced senescence. Our model provides mechanistic insights into the distinct stages of DNA-damage induced senescence, the relationship between senescence and epithelial-to-mesenchymal transition in cancer and offers a foundation for building multiscale models of tissue aging.

Unlocking Mitochondrial Dysfunction-Associated Senescence (MiDAS) with NAD + - a Boolean Model of Mitochondrial Dynamics and Cell Cycle Control

The steady accumulation of senescent cells with aging creates tissue environments that aid cancer evolution. Aging cell states are highly heterogeneous. 'Deep senescent' cells rely on healthy mitochondria to fuel a strong proinflammatory secretome, including cytokines, growth and transforming signals. Yet, the physiological triggers of senescence such as the reactive oxygen species (ROS) can also trigger mitochondrial dysfunction, and sufficient energy deficit to alter their secretome and cause chronic oxidative stress - a state termed Mitochondrial Dysfunction-Associated Senescence (MiDAS). Here, we offer a mechanistic hypothesis for the molecular processes leading to MiDAS, along with testable predictions. To do this we have built a Boolean regulatory network model that qualitatively captures key aspects of mitochondrial dynamics during cell cycle progression (hyper-fusion at the G1/S boundary, fission in mitosis), apoptosis (fission and dysfunction) and glucose starvation (reversible hyper-fusion), as well as MiDAS in response to SIRT3 knockdown or oxidative stress. Our model reaffirms the protective role of NAD + and external pyruvate. We offer testable predictions about the growth factor- and glucose-dependence of MiDAS and its reversibility at different stages of reactive oxygen species (ROS)-induced senescence. Our model provides mechanistic insights into the distinct stages of DNA-damage induced senescence, the relationship between senescence and epithelial-to-mesenchymal transition in cancer and offers a foundation for building multiscale models of tissue aging.

Highlights

Boolean regulatory network model reproduces mitochondrial dynamics during cell cycle progression, apoptosis, and glucose starvation. Model offers a mechanistic explanation for the positive feedback loop that locks in Mitochondrial Dysfunction-Associated Senescence (MiDAS), involving autophagy-resistant, hyperfused, dysfunctional mitochondria. Model reproduces ROS-mediated mitochondrial dysfunction and suggests that MiDAS is part of the early phase of damage-induced senescence. Model predicts that cancer-driving mutations that bypass the G1/S checkpoint generally increase the incidence of MiDAS, except for p53 loss.

Network analysis reveals that the tumor suppressor lncRNA GAS5 acts as a double-edged sword in response to DNA damage in gastric cancer

The lncRNA GAS5 acts as a tumor suppressor and is downregulated in gastric cancer (GC). In contrast, E2F1, an important transcription factor and tumor promoter, directly inhibits miR-34c expression in GC cell lines. Furthermore, in the corresponding GC cell lines, lncRNA GAS5 directly targets E2F1. However, lncRNA GAS5 and miR-34c remain to be studied in conjunction with GC. Here, we present a dynamic Boolean network to classify gene regulation between these two non-coding RNAs (ncRNAs) in GC. This is the first study to show that lncRNA GAS5 can positively regulate miR-34c in GC through a previously unknown molecular pathway coupling lncRNA/miRNA. We compared our network to several in-vivo/in-vitro experiments and obtained an excellent agreement. We revealed that lncRNA GAS5 regulates miR-34c by targeting E2F1. Additionally, we found that lncRNA GAS5, independently of p53, inhibits GC proliferation through the ATM/p38 MAPK signaling pathway. Accordingly, our results support that E2F1 is an engaging target of drug development in tumor growth and aggressive proliferation of GC, and favorable results can be achieved through tumor suppressor lncRNA GAS5/miR-34c axis in GC. Thus, our findings unlock a new avenue for GC treatment in response to DNA damage by these ncRNAs.

Relative importance of composition structures and biologically meaningful logics in bipartite Boolean models of gene regulation

Boolean networks have been widely used to model gene networks. However, such models are coarse-grained to an extent that they abstract away molecular specificities of gene regulation. Alternatively, bipartite Boolean network models of gene regulation explicitly distinguish genes from transcription factors (TFs). In such bipartite models, multiple TFs may simultaneously contribute to gene regulation by forming heteromeric complexes, thus giving rise to composition structures. Since bipartite Boolean models are relatively recent, an empirical investigation of their biological plausibility is lacking. Here, we estimate the prevalence of composition structures arising through heteromeric complexes. Moreover, we present an additional mechanism where composition structures may arise as a result of multiple TFs binding to cis-regulatory regions and provide empirical support for this mechanism. Next, we compare the restriction in BFs imposed by composition structures and by biologically meaningful properties. We find that though composition structures can severely restrict the number of Boolean functions (BFs) driving a gene, the two types of minimally complex BFs, namely nested canalyzing functions (NCFs) and read-once functions (RoFs), are comparatively more restrictive. Finally, we find that composition structures are highly enriched in real networks, but this enrichment most likely comes from NCFs and RoFs.

Probabilistic edge weights fine-tune Boolean network dynamics

Biological systems are noisy by nature. This aspect is reflected in our experimental measurements and should be reflected in the models we build to better understand these systems. Noise can be especially consequential when trying to interpret specific regulatory interactions, i.e. regulatory network edges. In this paper, we propose a method to explicitly encode edge-noise in Boolean dynamical systems by probabilistic edge-weight (PEW) operators. PEW operators have two important features: first, they introduce a form of edge-weight into Boolean models through the noise, second, the noise is dependent on the dynamical state of the system, which enables more biologically meaningful modeling choices. Moreover, we offer a simple-to-use implementation in the already well-established BooleanNet framework. In two application cases, we show how the introduction of just a few PEW operators in Boolean models can fine-tune the emergent dynamics and increase the accuracy of qualitative predictions. This includes fine-tuning interactions which cause non-biological behaviors when switching between asynchronous and synchronous update schemes in dynamical simulations. Moreover, PEW operators also open the way to encode more exotic cellular dynamics, such as cellular learning, and to implementing edge-weights for regulatory networks inferred from omics data.

In Silico Pleiotropy Analysis in KEGG Signaling Networks Using a Boolean Network Model

Pleiotropy, which refers to the ability of different mutations on the same gene to cause different pathological effects in human genetic diseases, is important in understanding system-level biological diseases. Although some biological experiments have been proposed, still little is known about pleiotropy on gene–gene dynamics, since most previous studies have been based on correlation analysis. Therefore, a new perspective is needed to investigate pleiotropy in terms of gene–gene dynamical characteristics. To quantify pleiotropy in terms of network dynamics, we propose a measure called in silico Pleiotropic Scores (sPS), which represents how much a gene is affected against a pair of different types of mutations on a Boolean network model. We found that our model can identify more candidate pleiotropic genes that are not known to be pleiotropic than the experimental database. In addition, we found that many types of functionally important genes tend to have higher sPS values than other genes; in other words, they are more pleiotropic. We investigated the relations of sPS with the structural properties in the signaling network and found that there are highly positive relations to degree, feedback loops, and centrality measures. This implies that the structural characteristics are principles to identify new pleiotropic genes. Finally, we found some biological evidence showing that sPS analysis is relevant to the real pleiotropic data and can be considered a novel candidate for pleiotropic gene research. Taken together, our results can be used to understand the dynamics pleiotropic characteristics in complex biological systems in terms of gene–phenotype relations.

Dynamical Analysis of a Boolean Network Model of the Oncogene Role of lncRNA ANRIL and lncRNA UFC1 in Non-Small Cell Lung Cancer

Long non-coding RNA (lncRNA) such as ANRIL and UFC1 have been verified as oncogenic genes in non-small cell lung cancer (NSCLC). It is well known that the tumor suppressor microRNA-34a (miR-34a) is downregulated in NSCLC. Furthermore, miR-34a induces senescence and apoptosis in breast, glioma, cervical cancer including NSCLC by targeting Myc. Recent evidence suggests that these two lncRNAs act as a miR-34a sponge in corresponding cancers. However, the biological functions between these two non-coding RNAs (ncRNAs) have not yet been studied in NSCLC. Therefore, we present a Boolean model to analyze the gene regulation between these two ncRNAs in NSCLC. We compared our model to several experimental studies involving gain- or loss-of-function genes in NSCLC cells and achieved an excellent agreement. Additionally, we predict three positive circuits involving miR-34a/E2F1/ANRIL, miR-34a/E2F1/UFC1, and miR-34a/Myc/ANRIL. Our circuit- perturbation analysis shows that these circuits are important for regulating cell-fate decisions such as senescence and apoptosis. Thus, our Boolean network permits an explicit cell-fate mechanism associated with NSCLC. Therefore, our results support that ANRIL and/or UFC1 is an attractive target for drug development in tumor growth and aggressive proliferation of NSCLC, and that a valuable outcome can be achieved through the miRNA-34a/Myc pathway.

Hepatectomy-Induced Alterations in Hepatic Perfusion and Function - Toward Multi-Scale Computational Modeling for a Better Prediction of Post-hepatectomy Liver Function

Liver resection causes marked perfusion alterations in the liver remnant both on the organ scale (vascular anatomy) and on the microscale (sinusoidal blood flow on tissue level). These changes in perfusion affect hepatic functions via direct alterations in blood supply and drainage, followed by indirect changes of biomechanical tissue properties and cellular function. Changes in blood flow impose compression, tension and shear forces on the liver tissue. These forces are perceived by mechanosensors on parenchymal and non-parenchymal cells of the liver and regulate cell-cell and cell-matrix interactions as well as cellular signaling and metabolism. These interactions are key players in tissue growth and remodeling, a prerequisite to restore tissue function after PHx. Their dysregulation is associated with metabolic impairment of the liver eventually leading to liver failure, a serious post-hepatectomy complication with high morbidity and mortality. Though certain links are known, the overall functional change after liver surgery is not understood due to complex feedback loops, non-linearities, spatial heterogeneities and different time-scales of events. Computational modeling is a unique approach to gain a better understanding of complex biomedical systems. This approach allows (i) integration of heterogeneous data and knowledge on multiple scales into a consistent view of how perfusion is related to hepatic function; (ii) testing and generating hypotheses based on predictive models, which must be validated experimentally and clinically. In the long term, computational modeling will (iii) support surgical planning by predicting surgery-induced perfusion perturbations and their functional (metabolic) consequences; and thereby (iv) allow minimizing surgical risks for the individual patient. Here, we review the alterations of hepatic perfusion, biomechanical properties and function associated with hepatectomy. Specifically, we provide an overview over the clinical problem, preoperative diagnostics, functional imaging approaches, experimental approaches in animal models, mechanoperception in the liver and impact on cellular metabolism, omics approaches with a focus on transcriptomics, data integration and uncertainty analysis, and computational modeling on multiple scales. Finally, we provide a perspective on how multi-scale computational models, which couple perfusion changes to hepatic function, could become part of clinical workflows to predict and optimize patient outcome after complex liver surgery.