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

A Boolean network model of hypoxia, mechanosensing and TGF-β signaling captures the role of phenotypic plasticity and mutations in tumor metastasis

The tumor microenvironment aids cancer progression by promoting several cancer hallmarks, independent of cancer-related mutations. Biophysical properties of this environment, such as the stiffness of the matrix cells adhere to and local cell density, impact proliferation, apoptosis, and the epithelial to mesenchymal transition (EMT). The latter is a rate-limiting step for invasion and metastasis, enhanced in hypoxic tumor environments but hindered by soft matrices and/or high cell densities. As these influences are often studied in isolation, the crosstalk between hypoxia, biomechanical signals, and the classic EMT driver TGF-β is not well mapped, limiting our ability to predict and anticipate cancer cell behaviors in changing tumor environments. To address this, we built a Boolean regulatory network model that integrates hypoxic signaling with a mechanosensitive model of EMT, which includes the EMT-promoting crosstalk of mitogens and biomechanical signals, cell cycle control, and apoptosis. Our model reproduces the requirement of Hif-1α for proliferation, the anti-proliferative effects of strong Hif-1α stabilization during hypoxia, hypoxic protection from anoikis, and hypoxia-driven mechanosensitive EMT. We offer experimentally testable predictions about the effect of VHL loss on cancer hallmarks, with or without secondary oncogene activation. Taken together, our model serves as a predictive framework to synthesize the signaling responses associated with tumor progression and metastasis in healthy vs. mutant cells. Our single-cell model is a key step towards more extensive regulatory network models that cover damage-response and senescence, integrating most cell-autonomous cancer hallmarks into a single model that can, in turn, control the behavior of in silico cells within a tissue model of epithelial homeostasis and carcinoma.

Detection of Oral Submucous Fibrosis at the Crossroads of Altering Nucleogeometry, Nuclear Spatial Distribution and Nuclear Molecular Abundance

The manuscript concentrates on spatial distribution and various nucleomorphometric parameters of epithelium to diagnose Oral Submucous Fibrosis (OSF). Histologically confirmed OSF and normal submucosa tissue samples were procured and stained with diamidino phenylindole (DAPI) to visualize nuclei. E-cadherins and p63 were also immunohistochemically stained. Microphotographs were analyzed to quantify the spatial distance among the nuclei of the stained tissue samples. Intensity of the DAPI stained nuclei and p63 was quantified. In addition, morphometric analysis of nuclei was done with the help of ImageJ software to quantify the geometric alterations in OSF tissue. Spatial distances among the nuclei of OSF tissue samples were found to be significantly higher than that of normal tissue. We also observed a significant decrease in the mean intensity of DAPI and p63 in OSF tissue samples. In addition, we have found statistically significant alterations of various morphometric quantifications in OSF tissue nuclei. There was a considerable change in the spatial distribution of nuclei, as well as some distinct changes in the nucleogeometry of OSF tissue, which corresponds to histological abnormalities. Decreased intensity of DAPI and p63 advocate disease progression. The biomarkers in this study are the accountable role-players for early detection of oral carcinoma.

Attractor-Transition Control of Complex Biological Networks: A Constant Control Approach

This article presents attractor-transition control of complex biological networks represented by Boolean networks (BNs) wherein the BN is steered from a prescribed initial attractor toward a desired one. The proposed approach leverages the similarity between attractors and Boolean algebraic properties embedded in the underlying state transition equations. To enhance the clarity of expression regarding stabilization toward the desired attractor, a simple coordinate transformation is performed on the considered BN. Based on the characteristics of transformed state equations, self-stabilizing state variables requiring no control efforts are derived in the first. Next, by applying the feedback vertex set (FVS) control scheme, control inputs stabilizing the remaining state variables are determined. The proposed control scheme exhibits versatility by accommodating both fixed-point and cyclic attractors. We validate the effectiveness of the proposed strategy through extensive numerical experiments conducted on random BNs as well as complex biological systems. In adherence to the reproducible research initiative, detailed results of numerical experiments and all the implementation codes are provided on the authors' website: https://github.com/choonlog/AttractorTransition.

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.