Robustness and fragility of Boolean models for genetic regulatory networks

Exploring gene regulation and biological processes in insects: Insights from omics data using gene regulatory network models

Gene regulatory network (GRN) comprises complicated yet intertwined gene-regulator relationships. Understanding the GRN dynamics will unravel the complexity behind the observed gene expressions. Insect gene regulation is often complicated due to their complex life cycles and diverse ecological adaptations. The main interest of this review is to have an update on the current mathematical modelling methods of GRNs to explain insect science. Several popular GRN architecture models are discussed, together with examples of applications in insect science. In the last part of this review, each model is compared from different aspects, including network scalability, computation complexity, robustness to noise and biological relevancy.

Spiking image processing unit based on neural analog of Boolean logic operations

McCulloch and Pitts hypothesized in 1943 that the brain is entirely composed of logic gates, akin to current computers' IP cores, which led to several neural analogs of Boolean logic. The current study proposes a spiking image processing unit (SIPU) based on spiking frequency gates and coordinate logic operations, as a dynamical model of synapses and spiking neurons. SIPU can imitate DSP functions like edge recognition, picture magnification, noise reduction, etc. but can be extended to cater for more advanced computing tasks. The proposed spiking Boolean logic platform can be used to develop advanced applications without relying on learning or specialized datasets. It could aid in gaining a deeper understanding of complex brain functions and spur new forms of neural analogs.

Predicting attractors from spectral properties of stylized gene regulatory networks

How the architecture of gene regulatory networks shapes gene expression patterns is an open question, which has been approached from a multitude of angles. The dominant strategy has been to identify nonrandom features in these networks and then argue for the function of these features using mechanistic modeling. Here we establish the foundation of an alternative approach by studying the correlation of network eigenvectors with synthetic gene expression data simulated with a basic and popular model of gene expression dynamics: Boolean threshold dynamics in signed directed graphs. We show that eigenvectors of the network adjacency matrix can predict collective states (attractors). However, the overall predictive power depends on details of the network architecture, namely the fraction of positive 3-cycles, in a predictable fashion. Our results are a set of statistical observations, providing a systematic step towards a further theoretical understanding of the role of network eigenvectors in dynamics on graphs.

Robust KALMAN Filter State Estimation for Gene Regulatory Networks

This paper proposes a revised version of the robust generalized maximum likelihood (GM)-type unscented KALMAN filter (GM-UKF) for the state estimation of gene regulatory networks (GRNs) in the presence of different types of deviations from assumptions. As known, the parameters and the power of the assumed noises within the GRN model may change abruptly as a result of jump behavior and bursting process in transcription and translation phases. Moreover, there may be outlying samples among genomic measurement data. Some other outliers may also occur in the model dynamics. The outliers may be misinterpreted by the filtering method if not detected and downweighted. To deal with all such deviations, a robust GM-UKF is designed that includes some modifications to address the challenges in calculating the projection statistics in GRNs such as the nonlinear behavior and the natural distance of the states. The proposed filter is compared to four Bayesian filters, i.e., the conventional UKF, the H _∞-UKF, the downweighting UKF (DW-UKF), and a modified version of the GM-UKF, the so-called maximum-likelihood UKF(M-UKF). The outcome results demonstrate that the GM-UKF outperforms other methods for all outlier types while the H _∞-UKF is appropriate for the changes in noise powers.

Functional Resilience of Mutually Repressing Motifs Embedded in Larger Networks

Elucidating the design principles of regulatory networks driving cellular decision-making has important implications for understanding cell differentiation and guiding the design of synthetic circuits. Mutually repressing feedback loops between 'master regulators' of cell fates can exhibit multistable dynamics enabling "single-positive" phenotypes: (high A, low B) and (low A, high B) for a toggle switch, and (high A, low B, low C), (low A, high B, low C) and (low A, low B, high C) for a toggle triad. However, the dynamics of these two motifs have been interrogated in isolation in silico, but in vitro and in vivo, they often operate while embedded in larger regulatory networks. Here, we embed these motifs in complex larger networks of varying sizes and connectivity to identify hallmarks under which these motifs maintain their canonical dynamical behavior. We show that an increased number of incoming edges onto a motif leads to a decay in their canonical stand-alone behaviors. We also show that this decay can be exacerbated by adding self-inhibition but not self-activation loops on the 'master regulators'. These observations offer insights into the design principles of biological networks containing these motifs and can help devise optimal strategies for the integration of these motifs into larger synthetic networks.

Attractor-Specific and Common Expression Values in Random Boolean Network Models (with a Preliminary Look at Single-Cell Data)

Random Boolean Networks (RBNs for short) are strongly simplified models of gene regulatory networks (GRNs), which have also been widely studied as abstract models of complex systems and have been used to simulate different phenomena. We define the “common sea” (CS) as the set of nodes that take the same value in all the attractors of a given network realization, and the “specific part” (SP) as the set of all the other nodes, and we study their properties in different ensembles, generated with different parameter values. Both the CS and of the SP can be composed of one or more weakly connected components, which are emergent intermediate-level structures. We show that the study of these sets provides very important information about the behavior of the model. The distribution of distances between attractors is also examined. Moreover, we show how the notion of a “common sea” of genes can be used to analyze data from single-cell experiments.

Effective connectivity determines the critical dynamics of biochemical networks

Living systems comprise interacting biochemical components in very large networks. Given their high connectivity, biochemical dynamics are surprisingly not chaotic but quite robust to perturbations-a feature C.H. Waddington named canalization. Because organisms are also flexible enough to evolve, they arguably operate in a critical dynamical regime between order and chaos. The established theory of criticality is based on networks of interacting automata where Boolean truth values model presence/absence of biochemical molecules. The dynamical regime is predicted using network connectivity and node bias (to be on/off) as tuning parameters. Revising this to account for canalization leads to a significant improvement in dynamical regime prediction. The revision is based on effective connectivity, a measure of dynamical redundancy that buffers automata response to some inputs. In both random and experimentally validated systems biology networks, reducing effective connectivity makes living systems operate in stable or critical regimes even though the structure of their biochemical interaction networks predicts them to be chaotic. This suggests that dynamical redundancy may be naturally selected to maintain living systems near critical dynamics, providing both robustness and evolvability. By identifying how dynamics propagates preferably via effective pathways, our approach helps to identify precise ways to design and control network models of biochemical regulation and signalling.

The Multiple Dimensions of Networks in Cancer: A Perspective

This perspective article gathers the latest developments in mathematical and computational oncology tools that exploit network approaches for the mathematical modelling, analysis, and simulation of cancer development and therapy design. It instigates the community to explore new paths and synergies under the umbrella of the Special Issue “Networks in Cancer: From Symmetry Breaking to Targeted Therapy”. The focus of the perspective is to demonstrate how networks can model the physics, analyse the interactions, and predict the evolution of the multiple processes behind tumour-host encounters across multiple scales. From agent-based modelling and mechano-biology to machine learning and predictive modelling, the perspective motivates a methodology well suited to mathematical and computational oncology and suggests approaches that mark a viable path towards adoption in the clinic.

A comprehensive survey of developmental programs reveals a dearth of tree-like lineage graphs and ubiquitous regeneration

Background

Multicellular organisms are characterized by a wide diversity of forms and complexity despite a restricted set of key molecules and mechanisms at the base of organismal development. Development combines three basic processes—asymmetric cell division, signaling, and gene regulation—in a multitude of ways to create this overwhelming diversity of multicellular life forms. Here, we use a generative model to test the limits to which such processes can be combined to generate multiple differentiation paths during development, and attempt to chart the diversity of multicellular organisms generated.

Results

We sample millions of biologically feasible developmental schemes, allowing us to comment on the statistical properties of cell differentiation trajectories they produce. We characterize model-generated “organisms” using the graph topology of their cell type lineage maps. Remarkably, tree-type lineage differentiation maps are the rarest in our data. Additionally, a majority of the “organisms” generated by our model appear to be endowed with the ability to regenerate using pluripotent cells.

Conclusions

Our results indicate that, in contrast to common views, cell type lineage graphs are unlikely to be tree-like. Instead, they are more likely to be directed acyclic graphs, with multiple lineages converging on the same terminal cell type. Furthermore, the high incidence of pluripotent cells in model-generated organisms stands in line with the long-standing hypothesis that whole body regeneration is an epiphenomenon of development. We discuss experimentally testable predictions of our model and some ways to adapt the generative framework to test additional hypotheses about general features of development.

Supplementary Information

The online version contains supplementary material available at (10.1186/s12915-021-01013-4).

Qualitative Modeling, Analysis and Control of Synthetic Regulatory Circuits

Qualitative modeling approaches are promising and still underexploited tools for the analysis and design of synthetic circuits. They can make predictions of circuit behavior in the absence of precise, quantitative information. Moreover, they provide direct insight into the relation between the feedback structure and the dynamical properties of a network. We review qualitative modeling approaches by focusing on two specific formalisms, Boolean networks and piecewise-linear differential equations, and illustrate their application by means of three well-known synthetic circuits. We describe various methods for the analysis of state transition graphs, discrete representations of the network dynamics that are generated in both modeling frameworks. We also briefly present the problem of controlling synthetic circuits, an emerging topic that could profit from the capacity of qualitative modeling approaches to rapidly scan a space of design alternatives.

Canalization and Robustness in Human Genetics and Disease

Canalization refers to the evolution of populations such that the number of individuals who deviate from the optimum trait, or experience disease, is minimized. In the presence of rapid cultural, environmental, or genetic change, the reverse process of decanalization may contribute to observed increases in disease prevalence. This review starts by defining relevant concepts, drawing distinctions between the canalization of populations and robustness of individuals. It then considers evidence pertaining to three continuous traits and six domains of disease. In each case, existing genetic evidence for genotype-by-environment interactions is insufficient to support a strong inference of decanalization, but we argue that the advent of genome-wide polygenic risk assessment now makes an empirical evaluation of the role of canalization in preventing disease possible. Finally, the contributions of both rare and common variants to congenital abnormality and adult onset disease are considered in light of a new kerplunk model of genetic effects.

General method to find the attractors of discrete dynamic models of biological systems

Analyzing the long-term behaviors (attractors) of dynamic models of biological networks can provide valuable insight. We propose a general method that can find the attractors of multilevel discrete dynamical systems by extending a method that finds the attractors of a Boolean network model. The previous method is based on finding stable motifs, subgraphs whose nodes' states can stabilize on their own. We extend the framework from binary states to any finite discrete levels by creating a virtual node for each level of a multilevel node, and describing each virtual node with a quasi-Boolean function. We then create an expanded representation of the multilevel network, find multilevel stable motifs and oscillating motifs, and identify attractors by successive network reduction. In this way, we find both fixed point attractors and complex attractors. We implemented an algorithm, which we test and validate on representative synthetic networks and on published multilevel models of biological networks. Despite its primary motivation to analyze biological networks, our motif-based method is general and can be applied to any finite discrete dynamical system.

Dynamics of Bacterial Gene Regulatory Networks

The ability of bacterial cells to adjust their gene expression program in response to environmental perturbation is often critical for their survival. Recent experimental advances allowing us to quantitatively record gene expression dynamics in single cells and in populations coupled with mathematical modeling enable mechanistic understanding on how these responses are shaped by the underlying regulatory networks. Here, we review how the combination of local and global factors affect dynamical responses of gene regulatory networks. Our goal is to discuss the general principles that allow extrapolation from a few model bacteria to less understood microbes. We emphasize that, in addition to well-studied effects of network architecture, network dynamics are shaped by global pleiotropic effects and cell physiology.

Global dynamics for switching systems and their extensions by linear differential equations

Switching systems use piecewise constant nonlinearities to model gene regulatory networks. This choice provides advantages in the analysis of behavior and allows the global description of dynamics in terms of Morse graphs associated to nodes of a parameter graph. The parameter graph captures spatial characteristics of a decomposition of parameter space into domains with identical Morse graphs. However, there are many cellular processes that do not exhibit threshold-like behavior and thus are not well described by a switching system. We consider a class of extensions of switching systems formed by a mixture of switching interactions and chains of variables governed by linear differential equations. We show that the parameter graphs associated to the switching system and any of its extensions are identical. For each parameter graph node, there is an order-preserving map from the Morse graph of the switching system to the Morse graph of any of its extensions. We provide counterexamples that show why possible stronger relationships between the Morse graphs are not valid.

Method for identification of sensitive nodes in Boolean models of biological networks

Biological systems are often represented as Boolean networks and analysed to identify sensitive nodes which on perturbation disproportionately change a predefined output. There exist different kinds of perturbation methods: perturbation of function, perturbation of state and perturbation in update scheme. Nodes may have defects in interpretation of the inputs from other nodes and calculation of the node output. To simulate these defects and systematically assess their effect on the system output, two new function perturbations, referred to as ‘not of function’ and ‘function of not’, are introduced. In the former, the inputs are assumed to be correctly interpreted but the output of the update rule is perturbed; and in the latter, each input is perturbed but the correct update rule is applied. These and previously used perturbation methods were applied to two existing Boolean models, namely the human melanogenesis signalling network and the fly segment polarity network. Through mathematical simulations, it was found that these methods successfully identified nodes earlier found to be sensitive using other methods, and were also able to identify sensitive nodes which were previously unreported.

Boolean network modeling in systems pharmacology

Quantitative systems pharmacology (QSP) is an emerging discipline that aims to discover how drugs modulate the dynamics of biological components in molecular and cellular networks and the impact of those perturbations on human pathophysiology. The integration of systems-based experimental and computational approaches is required to facilitate the advancement of this field. QSP models typically consist of a series of ordinary differential equations (ODE). However, this mathematical framework requires extensive knowledge of parameters pertaining to biological processes, which is often unavailable. An alternative framework that does not require knowledge of system-specific parameters, such as Boolean network modeling, could serve as an initial foundation prior to the development of an ODE-based model. Boolean network models have been shown to efficiently describe, in a qualitative manner, the complex behavior of signal transduction and gene/protein regulatory processes. In addition to providing a starting point prior to quantitative modeling, Boolean network models can also be utilized to discover novel therapeutic targets and combinatorial treatment strategies. Identifying drug targets using a network-based approach could supplement current drug discovery methodologies and help to fill the innovation gap across the pharmaceutical industry. In this review, we discuss the process of developing Boolean network models and the various analyses that can be performed to identify novel drug targets and combinatorial approaches. An example for each of these analyses is provided using a previously developed Boolean network of signaling pathways in multiple myeloma. Selected examples of Boolean network models of human (patho-)physiological systems are also reviewed in brief.

A toolbox for discrete modelling of cell signalling dynamics

We present a library of network motifs for the development of complex and realistic biological network models using the BioModelAnalyzer, and demonstrate their wider value by using them to construct a model of the cell cycle.

Understanding the mTOR signaling pathway via mathematical modeling

The mechanistic target of rapamycin (mTOR) is a central regulatory pathway that integrates a variety of environmental cues to control cellular growth and homeostasis by intricate molecular feedbacks. In spite of extensive knowledge about its components, the molecular understanding of how these function together in space and time remains poor and there is a need for Systems Biology approaches to perform systematic analyses. In this work, we review the recent progress how the combined efforts of mathematical models and quantitative experiments shed new light on our understanding of the mTOR signaling pathway. In particular, we discuss the modeling concepts applied in mTOR signaling, the role of multiple feedbacks and the crosstalk mechanisms of mTOR with other signaling pathways. We also discuss the contribution of principles from information and network theory that have been successfully applied in dissecting design principles of the mTOR signaling network. We finally propose to classify the mTOR models in terms of the time scale and network complexity, and outline the importance of the classification toward the development of highly comprehensive and predictive models. WIREs Syst Biol Med 2017, 9:e1379. doi: 10.1002/wsbm.1379 For further resources related to this article, please visit the WIREs website.

Inferring Intracellular Signal Transduction Circuitry from Molecular Perturbation Experiments

The development of network inference methodologies that accurately predict connectivity in dysregulated pathways may enable the rational selection of patient therapies. Accurately inferring an intracellular network from data remains a very challenging problem in molecular systems biology. Living cells integrate extremely robust circuits that exhibit significant heterogeneity, but still respond to external stimuli in predictable ways. This phenomenon allows us to introduce a network inference methodology that integrates measurements of protein activation from perturbation experiments. The methodology relies on logic-based networks to provide a predictive approximation of the transfer of signals in a network. The approach presented was validated in silico with a set of test networks and applied to investigate the epidermal growth factor receptor signaling of a breast epithelial cell line, MFC10A. In our analysis, we predict the potential signaling circuitry most likely responsible for the experimental readouts of several proteins in the mitogen-activated protein kinase and phosphatidylinositol-3 kinase pathways. The approach can also be used to identify additional necessary perturbation experiments to distinguish between a set of possible candidate networks.

Mathematical and Computational Modeling in Complex Biological Systems

The biological process and molecular functions involved in the cancer progression remain difficult to understand for biologists and clinical doctors. Recent developments in high-throughput technologies urge the systems biology to achieve more precise models for complex diseases. Computational and mathematical models are gradually being used to help us understand the omics data produced by high-throughput experimental techniques. The use of computational models in systems biology allows us to explore the pathogenesis of complex diseases, improve our understanding of the latent molecular mechanisms, and promote treatment strategy optimization and new drug discovery. Currently, it is urgent to bridge the gap between the developments of high-throughput technologies and systemic modeling of the biological process in cancer research. In this review, we firstly studied several typical mathematical modeling approaches of biological systems in different scales and deeply analyzed their characteristics, advantages, applications, and limitations. Next, three potential research directions in systems modeling were summarized. To conclude, this review provides an update of important solutions using computational modeling approaches in systems biology.

Facing the challenges of multiscale modelling of bacterial and fungal pathogen-host interactions

Recent and rapidly evolving progress on high-throughput measurement techniques and computational performance has led to the emergence of new disciplines, such as systems medicine and translational systems biology. At the core of these disciplines lies the desire to produce multiscale models: mathematical models that integrate multiple scales of biological organization, ranging from molecular, cellular and tissue models to organ, whole-organism and population scale models. Using such models, hypotheses can systematically be tested. In this review, we present state-of-the-art multiscale modelling of bacterial and fungal infections, considering both the pathogen and host as well as their interaction. Multiscale modelling of the interactions of bacteria, especially Mycobacterium tuberculosis, with the human host is quite advanced. In contrast, models for fungal infections are still in their infancy, in particular regarding infections with the most important human pathogenic fungi, Candida albicans and Aspergillus fumigatus. We reflect on the current availability of computational approaches for multiscale modelling of host-pathogen interactions and point out current challenges. Finally, we provide an outlook for future requirements of multiscale modelling.

Causality analysis detects the regulatory role of maternal effect genes in the early Drosophila embryo

In developmental studies, inferring regulatory interactions of segmentation genetic network play a vital role in unveiling the mechanism of pattern formation. As such, there exists an opportune demand for theoretical developments and new mathematical models which can result in a more accurate illustration of this genetic network. Accordingly, this paper seeks to extract the meaningful regulatory role of the maternal effect genes using a variety of causality detection techniques and to explore whether these methods can suggest a new analytical view to the gene regulatory networks. We evaluate the use of three different powerful and widely-used models representing time and frequency domain Granger causality and convergent cross mapping technique with the results being thoroughly evaluated for statistical significance. Our findings show that the regulatory role of maternal effect genes is detectable in different time classes and thereby the method is applicable to infer the possible regulatory interactions present among the other genes of this network.

Transient dynamics and their control in time-delay autonomous Boolean ring networks

Biochemical systems with switch-like interactions, such as gene regulatory networks, are well modeled by autonomous Boolean networks. Specifically, the topology and logic of gene interactions can be described by systems of continuous piecewise-linear differential equations, enabling analytical predictions of the dynamics of specific networks. However, most models do not account for time delays along links associated with spatial transport, mRNA transcription, and translation. To address this issue, we have developed an experimental test bed to realize a time-delay autonomous Boolean network with three inhibitory nodes, known as a repressilator, and use it to study the dynamics that arise as time delays along the links vary. We observe various nearly periodic oscillatory transient patterns with extremely long lifetime, which emerge in small network motifs due to the delay, and which are distinct from the eventual asymptotically stable periodic attractors. For repeated experiments with a given network, we find that stochastic processes give rise to a broad distribution of transient times with an exponential tail. In some cases, the transients are so long that it is doubtful the attractors will ever be approached in a biological system that has a finite lifetime. To counteract the long transients, we show experimentally that small, occasional perturbations applied to the time delays can force the trajectories to rapidly approach the attractors.

Compensatory interactions to stabilize multiple steady states or mitigate the effects of multiple deregulations in biological networks

Complex diseases can be modeled as damage to intracellular networks that results in abnormal cell behaviors. Network-based dynamic models such as Boolean models have been employed to model a variety of biological systems including those corresponding to disease. Previous work designed compensatory interactions to stabilize an attractor of a Boolean network after single node damage. We generalize this method to a multinode damage scenario and to the simultaneous stabilization of multiple steady state attractors. We classify the emergent situations, with a special focus on combinatorial effects, and characterize each class through simulation. We explore how the structural and functional properties of the network affect its resilience and its possible repair scenarios. We demonstrate the method's applicability to two intracellular network models relevant to cancer. This work has implications in designing prevention strategies for complex disease.

COMPUTATIONAL METHODS FOR ASYNCHRONOUS BASINS

For a Boolean network we consider asynchronous updates and define the exclusive asynchronous basin of attraction for any steady state or cyclic attractor. An algorithm based on commutative algebra is presented to compute the exclusive basin. Finally its use for targeting desirable attractors by selective intervention on network nodes is illustrated with two examples, one cell signalling network and one sensor network measuring human mobility.

Estimation of delays in generalized asynchronous Boolean networks

A new generalized asynchronous Boolean network (GABN) model has been proposed in this paper. This continuous-time discrete-state model captures the biological reality of cellular dynamics without compromising the computational efficiency of the Boolean framework. The GABN synthesis procedure is based on the prior knowledge of the logical structure of the regulatory network, and the experimental transcriptional parameters. The novelty of the proposed methodology lies in considering different delays associated with the activation and deactivation of a particular protein (especially the transcription factors). A few illustrative examples of some well-studied network motifs have been provided to explore the scope of using the GABN model for larger networks. The GABN model of the p53-signaling pathway in response to γ-irradiation has also been simulated in the current paper to provide an indirect validation of the proposed schema.

Stacking the odds for Golgi cisternal maturation

What is the minimal set of cell-biological ingredients needed to generate a Golgi apparatus? The compositions of eukaryotic organelles arise through a process of molecular exchange via vesicle traffic. Here we statistically sample tens of thousands of homeostatic vesicle traffic networks generated by realistic molecular rules governing vesicle budding and fusion. Remarkably, the plurality of these networks contain chains of compartments that undergo creation, compositional maturation, and dissipation, coupled by molecular recycling along retrograde vesicles. This motif precisely matches the cisternal maturation model of the Golgi, which was developed to explain many observed aspects of the eukaryotic secretory pathway. In our analysis cisternal maturation is a robust consequence of vesicle traffic homeostasis, independent of the underlying details of molecular interactions or spatial stacking. This architecture may have been exapted rather than selected for its role in the secretion of large cargo.

DOI:http://dx.doi.org/10.7554/eLife.16231.001

Logical Modeling and Dynamical Analysis of Cellular Networks

The logical (or logic) formalism is increasingly used to model regulatory and signaling networks. Complementing these applications, several groups contributed various methods and tools to support the definition and analysis of logical models. After an introduction to the logical modeling framework and to several of its variants, we review here a number of recent methodological advances to ease the analysis of large and intricate networks. In particular, we survey approaches to determine model attractors and their reachability properties, to assess the dynamical impact of variations of external signals, and to consistently reduce large models. To illustrate these developments, we further consider several published logical models for two important biological processes, namely the differentiation of T helper cells and the control of mammalian cell cycle.

Coherent organization in gene regulation: a study on six networks

Structural and dynamical fingerprints of evolutionary optimization in biological networks are still unclear. Here we analyze the dynamics of genetic regulatory networks responsible for the regulation of cell cycle and cell differentiation in three organisms or cell types each, and show that they follow a version of Hebb's rule which we have termed coherence. More precisely, we find that simultaneously expressed genes with a common target are less likely to act antagonistically at the attractors of the regulatory dynamics. We then investigate the dependence of coherence on structural parameters, such as the mean number of inputs per node and the activatory/repressory interaction ratio, as well as on dynamically determined quantities, such as the basin size and the number of expressed genes.

NACE: A web-based tool for prediction of intercompartmental efficiency of human molecular genetic networks

Molecular genetic processes generally involve proteins from distinct intracellular localisations. Reactions that follow the same process are distributed among various compartments within the cell. In this regard, the reaction rate and the efficiency of biological processes can depend on the subcellular localisation of proteins. Previously, the authors proposed a method of evaluating the efficiency of biological processes based on the analysis of the distribution of protein subcellular localisation (Popik et al., 2014). Here, NACE is presented, which is an open access web-oriented program that implements this method and allows the user to evaluate the intercompartmental efficiency of human molecular genetic networks. The method has been extended by a new feature that provides the evaluation of the tissue-specific efficiency of networks for more than 2800 anatomical structures. Such assessments are important in cases when molecular genetic pathways in different tissues proceed with the participation of various proteins with a number of intracellular localisations. For example, an analysis of KEGG pathways, conducted using the developed program, showed that the efficiencies of many KEGG pathways are tissue-specific. Analysis of efficiencies of regulatory pathways in the liver, linking proteins of the hepatitis C virus with human proteins involved in the KEGG apoptosis pathway, showed that intercompartmental efficiency might play an important role in host-pathogen interactions. Thus, the developed tool can be useful in the study of the effectiveness of functioning of various molecular genetic networks, including metabolic, regulatory, host-pathogen interactions and others taking into account tissue-specific gene expression. The tool is available via the following link: http://www-bionet.sscc.ru/nace/.

Wine glasses and hourglasses: Non-adaptive complexity of vesicle traffic in microbial eukaryotes

•

We are motivated by the diversity of vesicle traffic systems in microbial parasites.

•

We present a mathematical model of vesicle traffic in a manner accessible to a broad audience.

•

We show that many complex features of vesicle traffic systems arise spontaneously due to molecular interactions.

•

Traffic features such as compartmental maturation might arise non-adaptively and later be selected for function.

Microbial eukaryotes present a stunning diversity of endomembrane organization. From specialized secretory organelles such as the rhoptries and micronemes of apicomplexans, to peroxisome-derived metabolic compartments such as the glycosomes of kinetoplastids, different microbial taxa have explored different solutions to the compartmentalization and processing of cargo. The basic secretory and endocytic system, comprising the ER, Golgi, endosomes, and plasma membrane, as well as diverse taxon-specific specialized endomembrane organelles, are coupled by a complex network of cargo transport via vesicle traffic. It is tempting to connect form to function, ascribing biochemical roles to each compartment and vesicle of such a system. Here we argue that traffic systems of high complexity could arise through non-adaptive mechanisms via purely physical constraints, and subsequently be exapted for various taxon-specific functions. Our argument is based on a Boolean mathematical model of vesicle traffic: we specify rules of how compartments exchange vesicles; these rules then generate hypothetical cells with different types of endomembrane organization. Though one could imagine a large number of hypothetical vesicle traffic systems, very few of these are consistent with molecular interactions. Such molecular constraints are the bottleneck of a metaphorical hourglass, and the rules that make it through the bottleneck are expected to generate cells with many special properties. Sampling at random from among such rules represents an evolutionary null hypothesis: any properties of the resulting cells must be non-adaptive. We show by example that vesicle traffic systems generated in this random manner are reminiscent of the complex trafficking apparatus of real cells.

Combinatorial representation of parameter space for switching networks

We describe the theoretical and computational framework for the Dynamic Signatures for Genetic Regulatory Network ( DSGRN) database. The motivation stems from urgent need to understand the global dynamics of biologically relevant signal transduction/gene regulatory networks that have at least 5 to 10 nodes, involve multiple interactions, and decades of parameters. The input to the database computations is a regulatory network, i.e. a directed graph with edges indicating up or down regulation. A computational model based on switching networks is generated from the regulatory network. The phase space dimension of this model equals the number of nodes and the associated parameter space consists of one parameter for each node (a decay rate), and three parameters for each edge (low level of expression, high level of expression, and threshold at which expression levels change). Since the nonlinearities of switching systems are piece-wise constant, there is a natural decomposition of phase space into cells from which the dynamics can be described combinatorially in terms of a state transition graph. This in turn leads to a compact representation of the global dynamics called an annotated Morse graph that identifies recurrent and nonrecurrent dynamics. The focus of this paper is on the construction of a natural computable finite decomposition of parameter space into domains where the annotated Morse graph description of dynamics is constant. We use this decomposition to construct an SQL database that can be effectively searched for dynamical signatures such as bistability, stable or unstable oscillations, and stable equilibria. We include two simple 3-node networks to provide small explicit examples of the type of information stored in the DSGRN database. To demonstrate the computational capabilities of this system we consider a simple network associated with p53 that involves 5 nodes and a 29-dimensional parameter space.

Network modelling reveals the mechanism underlying colitis-associated colon cancer and identifies novel combinatorial anti-cancer targets

The connection between inflammation and tumourigenesis has been well established. However, the detailed molecular mechanism underlying inflammation-associated tumourigenesis remains unknown because this process involves a complex interplay between immune microenvironments and epithelial cells. To obtain a more systematic understanding of inflammation-associated tumourigenesis as well as to identify novel therapeutic approaches, we constructed a knowledge-based network describing the development of colitis-associated colon cancer (CAC) by integrating the extracellular microenvironment and intracellular signalling pathways. Dynamic simulations of the CAC network revealed a core network module, including P53, MDM2, and AKT, that may govern the malignant transformation of colon epithelial cells in a pro-tumor inflammatory microenvironment. Furthermore, in silico mutation studies and experimental validations led to a novel finding that concurrently targeting ceramide and PI3K/AKT pathway by chemical probes or marketed drugs achieves synergistic anti-cancer effects. Overall, our network model can guide further mechanistic studies on CAC and provide new insights into the design of combinatorial cancer therapies in a rational manner.

Fanconi anemia cells with unrepaired DNA damage activate components of the checkpoint recovery process

Background

The FA/BRCA pathway repairs DNA interstrand crosslinks. Mutations in this pathway cause Fanconi anemia (FA), a chromosome instability syndrome with bone marrow failure and cancer predisposition. Upon DNA damage, normal and FA cells inhibit the cell cycle progression, until the G2/M checkpoint is turned off by the checkpoint recovery, which becomes activated when the DNA damage has been repaired. Interestingly, highly damaged FA cells seem to override the G2/M checkpoint. In this study we explored with a Boolean network model and key experiments whether checkpoint recovery activation occurs in FA cells with extensive unrepaired DNA damage.

Methods

We performed synchronous/asynchronous simulations of the FA/BRCA pathway Boolean network model. FA-A and normal lymphoblastoid cell lines were used to study checkpoint and checkpoint recovery activation after DNA damage induction. The experimental approach included flow cytometry cell cycle analysis, cell division tracking, chromosome aberration analysis and gene expression analysis through qRT-PCR and western blot.

Results

Computational simulations suggested that in FA mutants checkpoint recovery activity inhibits the checkpoint components despite unrepaired DNA damage, a behavior that we did not observed in wild-type simulations. This result implies that FA cells would eventually reenter the cell cycle after a DNA damage induced G2/M checkpoint arrest, but before the damage has been fixed. We observed that FA-A cells activate the G2/M checkpoint and arrest in G2 phase, but eventually reach mitosis and divide with unrepaired DNA damage, thus resolving the initial checkpoint arrest. Based on our model result we look for ectopic activity of checkpoint recovery components. We found that checkpoint recovery components, such as PLK1, are expressed to a similar extent as normal undamaged cells do, even though FA-A cells harbor highly damaged DNA.

Conclusions

Our results show that FA cells, despite extensive DNA damage, do not loss the capacity to express the transcriptional and protein components of checkpoint recovery that might eventually allow their division with unrepaired DNA damage. This might allow cell survival but increases the genomic instability inherent to FA individuals and promotes cancer.

Systemic modeling myeloma-osteoclast interactions under normoxic/hypoxic condition using a novel computational approach

Interaction of myeloma cells with osteoclasts (OC) can enhance tumor cell expansion through activation of complex signaling transduction networks. Both cells reside in the bone marrow, a hypoxic niche. How OC-myeloma interaction in a hypoxic environment affects myeloma cell growth and their response to drug treatment is poorly understood. In this study, we i) cultured myeloma cells in the presence/absence of OCs under normoxia and hypoxia conditions and did protein profiling analysis using reverse phase protein array; ii) computationally developed an Integer Linear Programming approach to infer OC-mediated myeloma cell-specific signaling pathways under normoxic and hypoxic conditions. Our modeling analysis indicated that in the presence OCs, (1) cell growth-associated signaling pathways, PI3K/AKT and MEK/ERK, were activated and apoptotic regulatory proteins, BAX and BIM, down-regulated under normoxic condition; (2) β1 Integrin/FAK signaling pathway was activated in myeloma cells under hypoxic condition. Simulation of drug treatment effects by perturbing the inferred cell-specific pathways showed that targeting myeloma cells with the combination of PI3K and integrin inhibitors potentially (1) inhibited cell proliferation by reducing the expression/activation of NF-κB, S6, c-Myc, and c-Jun under normoxic condition; (2) blocked myeloma cell migration and invasion by reducing the expression of FAK and PKC under hypoxic condition.

An Extended, Boolean Model of the Septation Initiation Network in S.Pombe Provides Insights into Its Regulation

Cytokinesis in fission yeast is controlled by the Septation Initiation Network (SIN), a protein kinase signaling network using the spindle pole body as scaffold. In order to describe the qualitative behavior of the system and predict unknown mutant behaviors we decided to adopt a Boolean modeling approach. In this paper, we report the construction of an extended, Boolean model of the SIN, comprising most SIN components and regulators as individual, experimentally testable nodes. The model uses CDK activity levels as control nodes for the simulation of SIN related events in different stages of the cell cycle. The model was optimized using single knock-out experiments of known phenotypic effect as a training set, and was able to correctly predict a double knock-out test set. Moreover, the model has made in silico predictions that have been validated in vivo, providing new insights into the regulation and hierarchical organization of the SIN.

Inference of Network Dynamics and Metabolic Interactions in the Gut Microbiome

We present a novel methodology to construct a Boolean dynamic model from time series metagenomic information and integrate this modeling with genome-scale metabolic network reconstructions to identify metabolic underpinnings for microbial interactions. We apply this in the context of a critical health issue: clindamycin antibiotic treatment and opportunistic Clostridium difficile infection. Our model recapitulates known dynamics of clindamycin antibiotic treatment and C. difficile infection and predicts therapeutic probiotic interventions to suppress C. difficile infection. Genome-scale metabolic network reconstructions reveal metabolic differences between community members and are used to explore the role of metabolism in the observed microbial interactions. In vitro experimental data validate a key result of our computational model, that B. intestinihominis can in fact slow C. difficile growth.

The community of bacteria that live in our intestines (called the “gut microbiome”) is important to normal intestinal function, and destruction of this community has a causative role in diseases including obesity, diabetes, and even neurological disorders. Clostridum difficile is an opportunistic pathogenic bacterium that causes potentially life-threatening intestinal inflammation and diarrhea and frequently occurs after antibiotic treatment, which wipes out the normal intestinal bacterial community. We use a mathematical model to identify how the normal bacterial community interacts and how this community changes with antibiotic treatment and C. difficile infection. We use this model to identify bacteria that may inhibit C. difficile growth. Our model and subsequent experiments indicate that Barnesiella intestinihominis inhibits C. difficile growth. This result suggests that B. intestinihominis could potentially be used as a probiotic to treat or prevent C. difficile infection.

Boolean modeling: a logic-based dynamic approach for understanding signaling and regulatory networks and for making useful predictions

The biomolecules inside or near cells form a complex interacting system. Cellular phenotypes and behaviors arise from the totality of interactions among the components of this system. A fruitful way of modeling interacting biomolecular systems is by network-based dynamic models that characterize each component by a state variable, and describe the change in the state variables due to the interactions in the system. Dynamic models can capture the stable state patterns of this interacting system and can connect them to different cell fates or behaviors. A Boolean or logic model characterizes each biomolecule by a binary state variable that relates the abundance of that molecule to a threshold abundance necessary for downstream processes. The regulation of this state variable is described in a parameter free manner, making Boolean modeling a practical choice for systems whose kinetic parameters have not been determined. Boolean models integrate the body of knowledge regarding the components and interactions of biomolecular systems, and capture the system's dynamic repertoire, for example the existence of multiple cell fates. These models were used for a variety of systems and led to important insights and predictions. Boolean models serve as an efficient exploratory model, a guide for follow-up experiments, and as a foundation for more quantitative models.

Controllability and observability of Boolean networks arising from biology

Boolean networks are currently receiving considerable attention as a computational scheme for system level analysis and modeling of biological systems. Studying control-related problems in Boolean networks may reveal new insights into the intrinsic control in complex biological systems and enable us to develop strategies for manipulating biological systems using exogenous inputs. This paper considers controllability and observability of Boolean biological networks. We propose a new approach, which draws from the rich theory of symbolic computation, to solve the problems. Consequently, simple necessary and sufficient conditions for reachability, controllability, and observability are obtained, and algorithmic tests for controllability and observability which are based on the Gröbner basis method are presented. As practical applications, we apply the proposed approach to several different biological systems, namely, the mammalian cell-cycle network, the T-cell activation network, the large granular lymphocyte survival signaling network, and the Drosophila segment polarity network, gaining novel insights into the control and/or monitoring of the specific biological systems.

MONOMIALS AND BASIN CYLINDERS FOR NETWORK DYNAMICS

We describe methods to identify cylinder sets inside a basin of attraction for Boolean dynamics of biological networks. Such sets are used for designing regulatory interventions that make the system evolve towards a chosen attractor, for example initiating apoptosis in a cancer cell. We describe two algebraic methods for identifying cylinders inside a basin of attraction, one based on the Groebner fan that finds monomials that define cylinders and the other on primary decomposition. Both methods are applied to current examples of gene networks.

Analysis of signaling networks distributed over intracellular compartments based on protein-protein interactions

Background

Biological processes are usually distributed over various intracellular compartments. Proteins from diverse cellular compartments are often involved in similar signaling networks. However, the difference in the reaction rates between similar proteins among different compartments is usually quite high. We suggest that the estimation of frequency of intracompartmental as well as intercompartmental protein-protein interactions is an appropriate approach to predict the efficiency of a pathway.

Results

Using data from the databases STRING, ANDSystem, IntAct and UniProt, a PPI frequency matrix of intra/inter-compartmental interactions efficiencies was constructed. This matrix included 15 human-specific cellular compartments. An approach for estimating pathway efficiency using the matrix of intra/inter-compartmental PPI frequency, based on analysis of reactions efficiencies distribution was suggested. An investigation of KEGG pathway efficiencies was conducted using the developed method. The clusterization and the ranking of KEGG pathways based on their efficiency were performed. "Amino acid metabolism" and "Genetic information processing" revealed the highest efficiencies among other functional classes of KEGG pathways. "Nervous system" and "Signaling molecules interaction" contained the most inefficient pathways. Statistically significant differences were found between efficiencies of KEGG and randomly-generated pathways. Based on these observations, the validity of this approach was discussed.

Conclusion

The estimation of efficiency of signaling networks is a complicated task because of the need for the data on the kinetic reactions. However, the proposed method does not require such data and can be used for preliminary analysis of different protein networks.

Systems approaches in integrative cardiac biology: illustrations from cardiac heterocellular signalling studies

Understanding the complexity of cardiac physiology requires system-level studies of multiple cardiac cell types. Frequently, however, the end result of published research lacks the detail of the collaborative and integrative experimental design process, and the underlying conceptual framework. We review the recent progress in systems modelling and omics analysis of the heterocellular heart environment through complementary forward and inverse approaches, illustrating these conceptual and experimental frameworks with case studies from our own research program. The forward approach begins by collecting curated information from the niche cardiac biology literature, and connecting the dots to form mechanistic network models that generate testable system-level predictions. The inverse approach starts from the vast pool of public omics data in recent cardiac biological research, and applies bioinformatics analysis to produce novel candidates for further investigation. We also discuss the possibility of combining these two approaches into a hybrid framework, together with the benefits and challenges. These interdisciplinary research frameworks illustrate the interplay between computational models, omics analysis, and wet lab experiments, which holds the key to making real progress in improving human cardiac wellbeing.

Boolean modelling reveals new regulatory connections between transcription factors orchestrating the development of the ventral spinal cord

We have assembled a network of cell-fate determining transcription factors that play a key role in the specification of the ventral neuronal subtypes of the spinal cord on the basis of published transcriptional interactions. Asynchronous Boolean modelling of the network was used to compare simulation results with reported experimental observations. Such comparison highlighted the need to include additional regulatory connections in order to obtain the fixed point attractors of the model associated with the five known progenitor cell types located in the ventral spinal cord. The revised gene regulatory network reproduced previously observed cell state switches between progenitor cells observed in knock-out animal models or in experiments where the transcription factors were overexpressed. Furthermore the network predicted the inhibition of Irx3 by Nkx2.2 and this prediction was tested experimentally. Our results provide evidence for the existence of an as yet undescribed inhibitory connection which could potentially have significance beyond the ventral spinal cord. The work presented in this paper demonstrates the strength of Boolean modelling for identifying gene regulatory networks.

Multi-level modeling of light-induced stomatal opening offers new insights into its regulation by drought

Plant guard cells gate CO2 uptake and transpirational water loss through stomatal pores. As a result of decades of experimental investigation, there is an abundance of information on the involvement of specific proteins and secondary messengers in the regulation of stomatal movements and on the pairwise relationships between guard cell components. We constructed a multi-level dynamic model of guard cell signal transduction during light-induced stomatal opening and of the effect of the plant hormone abscisic acid (ABA) on this process. The model integrates into a coherent network the direct and indirect biological evidence regarding the regulation of seventy components implicated in stomatal opening. Analysis of this signal transduction network identified robust cross-talk between blue light and ABA, in which [Ca2+]c plays a key role, and indicated an absence of cross-talk between red light and ABA. The dynamic model captured more than 10(31) distinct states for the system and yielded outcomes that were in qualitative agreement with a wide variety of previous experimental results. We obtained novel model predictions by simulating single component knockout phenotypes. We found that under white light or blue light, over 60%, and under red light, over 90% of all simulated knockouts had similar opening responses as wild type, showing that the system is robust against single node loss. The model revealed an open question concerning the effect of ABA on red light-induced stomatal opening. We experimentally showed that ABA is able to inhibit red light-induced stomatal opening, and our model offers possible hypotheses for the underlying mechanism, which point to potential future experiments. Our modelling methodology combines simplicity and flexibility with dynamic richness, making it well suited for a wide class of biological regulatory systems.

Inferring the Effects of Honokiol on the Notch Signaling Pathway in SW480 Colon Cancer Cells

In a tumor cell, the development of acquired therapeutic resistance and the ability to survive in extracellular environments that differ from the primary site are the result of molecular adaptations in potentially highly plastic molecular networks. The accurate prediction of intracellular networks in a tumor remains a difficult problem in cancer informatics. In order to make truly rational patient-driven therapeutic decisions, it will be critical to develop methodologies that can accurately infer the molecular circuitry in the cells of a specific tumor. Despite enormous heterogeneity, cellular networks elicit deterministic digital-like responses. We discuss the use and limitations of methodologies that model molecular networks in cancer cells as a digital circuit. We also develop a network model of Notch signaling in colon cancer using a novel reverse engineering logic-based method and published western blot data to elucidate the interactions likely present in the circuits of the SW480 colon cancer cell line. Within this framework, we make predictions related to the role that honokiol may be playing as an anti-cancer drug.

A computational paradigm for dynamic logic-gates in neuronal activity

In 1943 McCulloch and Pitts suggested that the brain is composed of reliable logic-gates similar to the logic at the core of today's computers. This framework had a limited impact on neuroscience, since neurons exhibit far richer dynamics. Here we propose a new experimentally corroborated paradigm in which the truth tables of the brain's logic-gates are time dependent, i.e., dynamic logic-gates (DLGs). The truth tables of the DLGs depend on the history of their activity and the stimulation frequencies of their input neurons. Our experimental results are based on a procedure where conditioned stimulations were enforced on circuits of neurons embedded within a large-scale network of cortical cells in-vitro. We demonstrate that the underlying biological mechanism is the unavoidable increase of neuronal response latencies to ongoing stimulations, which imposes a non-uniform gradual stretching of network delays. The limited experimental results are confirmed and extended by simulations and theoretical arguments based on identical neurons with a fixed increase of the neuronal response latency per evoked spike. We anticipate our results to lead to better understanding of the suitability of this computational paradigm to account for the brain's functionalities and will require the development of new systematic mathematical methods beyond the methods developed for traditional Boolean algebra.