Attractor landscape analysis of colorectal tumorigenesis and its reversion

Background

Colorectal cancer arises from the accumulation of genetic mutations that induce dysfunction of intracellular signaling. However, the underlying mechanism of colorectal tumorigenesis driven by genetic mutations remains yet to be elucidated.

Results

To investigate colorectal tumorigenesis at a system-level, we have reconstructed a large-scale Boolean network model of the human signaling network by integrating previous experimental results on canonical signaling pathways related to proliferation, metastasis, and apoptosis. Throughout an extensive simulation analysis of the attractor landscape of the signaling network model, we found that the attractor landscape changes its shape by expanding the basin of attractors for abnormal proliferation and metastasis along with the accumulation of driver mutations. A further hypothetical study shows that restoration of a normal phenotype might be possible by reversely controlling the attractor landscape. Interestingly, the targets of approved anti-cancer drugs were highly enriched in the identified molecular targets for the reverse control.

Conclusions

Our results show that the dynamical analysis of a signaling network based on attractor landscape is useful in acquiring a system-level understanding of tumorigenesis and developing a new therapeutic strategy.

     

Flickering as an early warning signal

Most work on generic early warning signals for critical transitions focuses on indicators of the phenomenon of critical slowing down that precedes a range of catastrophic bifurcation points. However, in highly stochastic environments, systems will tend to shift to alternative basins of attraction already far from such bifurcation points. In fact, strong perturbations (noise) may cause the system to “flicker” between the basins of attraction of the system’s alternative states. As a result, under such noisy conditions, critical slowing down is not relevant, and one would expect its related generic leading indicators to fail, signaling an impending transition. Here, we systematically explore how flickering may be detected and interpreted as a signal of an emerging alternative attractor. We show that—although the two mechanisms differ—flickering may often be reflected in rising variance, lag-1 autocorrelation and skewness in ways that resemble the effects of critical slowing down. In particular, we demonstrate how the probability distribution of a flickering system can be used to map potential alternative attractors and their resilience. Thus, while flickering systems differ in many ways from the classical image of critical transitions, changes in their dynamics may carry valuable information about upcoming major changes.     

Basins of attraction: population dynamics with two stable 4-cycles

We use the concepts of composite maps, basins of attraction, basin switching, and saddle fly-by's to make the ecological hypothesis of the existence of multiple attractors more accessible to experimental scrutiny. Specifically, in a periodically forced insect population growth model we identify multiple attractors, namely, two locally stable 4-cycles. Using the model-predicted basins of attraction, we examine data time series from a Tribolium experiment for evidence of the multiple attractors. We conclude that the multiple attractor hypothesis together with demographic stochasticity accounts for the experimental observations.     

Tumor progression: Chance and necessity in Darwinian and Lamarckian somatic (mutationless) evolution

Current investigation of cancer progression towards increasing malignancy focuses on the molecular pathways that produce the various cancerous traits of cells. Their acquisition is explained by the somatic mutation theory: tumor progression is the result of a neo-Darwinian evolution in the tissue. Herein cells are the units of selection. Random genetic mutations permanently affecting these pathways create malignant cell phenotypes that are selected for in the disturbed tissue. However, could it be that the capacity of the genome and its gene regulatory network to generate the vast diversity of cell types during development, i.e., to produce inheritable phenotypic changes without mutations, is harnessed by tumorigenesis to propel a directional change towards malignancy? Here we take an encompassing perspective, transcending the orthodoxy of molecular carcinogenesis and review mechanisms of somatic evolution beyond the Neo-Darwinian scheme. We discuss the central concept of "cancer attractors" - the hidden stable states of gene regulatory networks normally not occupied by cells. Noise-induced transitions into such attractors provide a source for randomness (chance) and regulatory constraints (necessity) in the acquisition of novel expression profiles that can be inherited across cell divisions, and hence, can be selected for. But attractors can also be reached in response to environmental signals - thus offering the possibility for inheriting acquired traits that can also be selected for. Therefore, we face the possibility of non-genetic (mutation-independent) equivalents to both Darwinian and Lamarckian evolution which may jointly explain the arrow of change pointing toward increasing malignancy.     

Noise-driven stem cell and progenitor population dynamics

Background

The balance between maintenance of the stem cell state and terminal differentiation is influenced by the cellular environment. The switching between these states has long been understood as a transition between attractor states of a molecular network. Herein, stochastic fluctuations are either suppressed or can trigger the transition, but they do not actually determine the attractor states.

Methodology/Principal Findings

We present a novel mathematical concept in which stem cell and progenitor population dynamics are described as a probabilistic process that arises from cell proliferation and small fluctuations in the state of differentiation. These state fluctuations reflect random transitions between different activation patterns of the underlying regulatory network. Importantly, the associated noise amplitudes are state-dependent and set by the environment. Their variability determines the attractor states, and thus actually governs population dynamics. This model quantitatively reproduces the observed dynamics of differentiation and dedifferentiation in promyelocytic precursor cells.

Conclusions/Significance

Consequently, state-specific noise modulation by external signals can be instrumental in controlling stem cell and progenitor population dynamics. We propose follow-up experiments for quantifying the imprinting influence of the environment on cellular noise regulation.     

A dynamical systems hypothesis of schizophrenia

We propose a top-down approach to the symptoms of schizophrenia based on a statistical dynamical framework. We show that a reduced depth in the basins of attraction of cortical attractor states destabilizes the activity at the network level due to the constant statistical fluctuations caused by the stochastic spiking of neurons. In integrate-and-fire network simulations, a decrease in the NMDA receptor conductances, which reduces the depth of the attractor basins, decreases the stability of short-term memory states and increases distractibility. The cognitive symptoms of schizophrenia such as distractibility, working memory deficits, or poor attention could be caused by this instability of attractor states in prefrontal cortical networks. Lower firing rates are also produced, and in the orbitofrontal and anterior cingulate cortex could account for the negative symptoms, including a reduction of emotions. Decreasing the GABA as well as the NMDA conductances produces not only switches between the attractor states, but also jumps from spontaneous activity into one of the attractors. We relate this to the positive symptoms of schizophrenia, including delusions, paranoia, and hallucinations, which may arise because the basins of attraction are shallow and there is instability in temporal lobe semantic memory networks, leading thoughts to move too freely round the attractor energy landscape.     

A comparison of two cell regulatory models entailing high dimensional attractors representing phenotype

Two models for mammalian cell regulation that invoke the concept of cellular phenotype represented by high dimensional dynamic attractor states are compared. In one model the attractors are derived from an experimentally determined genetic regulatory network (GRN) for the cell type. As the state space architecture within which the attractors are embedded is determined by the binding sites on proteins and the recognition sites on DNA the attractors can be described as “hard-wired” in the genome through the genomic DNA sequence. In the second model attractors arising from the interactions between active gene products (mainly proteins) and independent of the genomic sequence, are descended from a pre-cellular state from which life originated. As this model is based on the cell as an open system the attractor acts as the interface between the cell and its environment. Environmental sources of stress can serve to trigger attractor and therefore phenotypic, transitions without entailing genotypic sequence changes.It is asserted that the evidence from cell and molecular biological research and logic, favours the second model. If correct there are important implications for understanding how environmental factors impact on evolution and may be implicated in hereditary and somatic disease.     

Generating Boolean networks with a prescribed attractor structure

MOTIVATION: Dynamical modeling of gene regulation via network models constitutes a key problem for genomics. The long-run characteristics of a dynamical system are critical and their determination is a primary aspect of system analysis. In the other direction, system synthesis involves constructing a network possessing a given set of properties. This constitutes the inverse problem. Generally, the inverse problem is ill-posed, meaning there will be many networks, or perhaps none, possessing the desired properties. Relative to long-run behavior, we may wish to construct networks possessing a desirable steady-state distribution. This paper addresses the long-run inverse problem pertaining to Boolean networks (BNs). RESULTS: The long-run behavior of a BN is characterized by its attractors. The rest of the state transition diagram is partitioned into level sets, the j-th level set being composed of all states that transition to one of the attractor states in exactly j transitions. We present two algorithms for the attractor inverse problem. The attractors are specified, and the sizes of the predictor sets and the number of levels are constrained. Algorithm complexity and performance are analyzed. The algorithmic solutions have immediate application. Under the assumption that sampling is from the steady state, a basic criterion for checking the validity of a designed network is that there should be concordance between the attractor states of the model and the data states. This criterion can be used to test a design algorithm: randomly select a set of states to be used as data states; generate a BN possessing the selected states as attractors, perhaps with some added requirements such as constraints on the number of predictors and the level structure; apply the design algorithm; and check the concordance between the attractor states of the designed network and the data states. AVAILABILITY: The software and supplementary material is available at http://gsp.tamu.edu/Publications/BNs/bn.htm     

Robustness and state-space structure of Boolean gene regulatory models

Robustness to perturbation is an important characteristic of genetic regulatory systems, but the relationship between robustness and model dynamics has not been clearly quantified. We propose a method for quantifying both robustness and dynamics in terms of state-space structures, for Boolean models of genetic regulatory systems. By investigating existing models of the Drosophila melanogaster segment polarity network and the Saccharomyces cerevisiae cell-cycle network, we show that the structure of attractor basins can yield insight into the underlying decision making required of the system, and also the way in which the system maximises its robustness. In particular, gene networks implementing decisions based on a few genes have simple state-space structures, and their attractors are robust by virtue of their simplicity. Gene networks with decisions that involve many interacting genes have correspondingly more complicated state-space structures, and robustness cannot be achieved through the structure of the attractor basins, but is achieved by larger attractor basins that dominate the state space. These different types of robustness are demonstrated by the two models: the D. melanogaster segment polarity network is robust due to simple attractor basins that implement decisions based on spatial signals; the S. cerevisiae cell-cycle network has a complicated state-space structure, and is robust only due to a giant attractor basin that dominates the state space.     

Mesoscopic biochemical basis of isogenetic inheritance and canalization: stochasticity, nonlinearity, and emergent landscape

Biochemical reaction systems in mesoscopic volume, under sustained environmental chemical gradient(s), can have multiple stochastic attractors. Two distinct mechanisms are known for their origins: (a) Stochastic single-molecule events, such as gene expression, with slow gene on-off dynamics; and (b) nonlinear networks with feedbacks. These two mechanisms yield different volume dependence for the sojourn time of an attractor. As in the classic Arrhenius theory for temperature dependent transition rates, a landscape perspective provides a natural framework for the system's behavior. However, due to the nonequilibrium nature of the open chemical systems, the landscape, and the attractors it represents, are all themselves emergent properties of complex, mesoscopic dynamics. In terms of the landscape, we show a generalization of Kramers' approach is possible to provide a rate theory. The emergence of attractors is a form of self-organization in the mesoscopic system; stochastic attractors in biochemical systems such as gene regulation and cellular signaling are naturally inheritable via cell division. Delbrück-Gillespie's mesoscopic reaction system theory, therefore, provides a biochemical basis for spontaneous isogenetic switching and canalization.     

Expansion or extinction: deterministic and stochastic two-patch models with Allee effects

We investigate the impact of Allee effect and dispersal on the long-term evolution of a population in a patchy environment. Our main focus is on whether a population already established in one patch either successfully invades an adjacent empty patch or undergoes a global extinction. Our study is based on the combination of analytical and numerical results for both a deterministic two-patch model and a stochastic counterpart. The deterministic model has either two, three or four attractors. The existence of a regime with exactly three attractors only appears when patches have distinct Allee thresholds. In the presence of weak dispersal, the analysis of the deterministic model shows that a high-density and a low-density populations can coexist at equilibrium in nearby patches, whereas the analysis of the stochastic model indicates that this equilibrium is metastable, thus leading after a large random time to either a global expansion or a global extinction. Up to some critical dispersal, increasing the intensity of the interactions leads to an increase of both the basin of attraction of the global extinction and the basin of attraction of the global expansion. Above this threshold, for both the deterministic and the stochastic models, the patches tend to synchronize as the intensity of the dispersal increases. This results in either a global expansion or a global extinction. For the deterministic model, there are only two attractors, while the stochastic model no longer exhibits a metastable behavior. In the presence of strong dispersal, the limiting behavior is entirely determined by the value of the Allee thresholds as the global population size in the deterministic and the stochastic models evolves as dictated by their single-patch counterparts. For all values of the dispersal parameter, Allee effects promote global extinction in terms of an expansion of the basin of attraction of the extinction equilibrium for the deterministic model and an increase of the probability of extinction for the stochastic model.     

Synergistic impacts of global warming on the resilience of coral reefs

Recent epizootics have removed important functional species from Caribbean coral reefs and left communities vulnerable to alternative attractors. Global warming will impact reefs further through two mechanisms. A chronic mechanism reduces coral calcification, which can result in depressed somatic growth. An acute mechanism, coral bleaching, causes extreme mortality when sea temperatures become anomalously high. We ask how these two mechanisms interact in driving future reef state (coral cover) and resilience (the probability of a reef remaining within a coral attractor). We find that acute mechanisms have the greatest impact overall, but the nature of the interaction with chronic stress depends on the metric considered. Chronic and acute stress act additively on reef state but form a strong synergy when influencing resilience by intensifying a regime shift. Chronic stress increases the size of the algal basin of attraction (at the expense of the coral basin), whereas coral bleaching pushes the system closer to the algal attractor. Resilience can change faster—and earlier—than a change in reef state. Therefore, we caution against basing management solely on measures of reef state because a loss of resilience can go unnoticed for many years and then become disproportionately more difficult to restore.     

Attractor analysis of asynchronous Boolean models of signal transduction networks

Prior work on the dynamics of Boolean networks, including analysis of the state space attractors and the basin of attraction of each attractor, has mainly focused on synchronous update of the nodes’ states. Although the simplicity of synchronous updating makes it very attractive, it fails to take into account the variety of time scales associated with different types of biological processes. Several different asynchronous update methods have been proposed to overcome this limitation, but there have not been any systematic comparisons of the dynamic behaviors displayed by the same system under different update methods. Here we fill this gap by combining theoretical analysis such as solution of scalar equations and Markov chain techniques, as well as numerical simulations to carry out a thorough comparative study on the dynamic behavior of a previously proposed Boolean model of a signal transduction network in plants. Prior evidence suggests that this network admits oscillations, but it is not known whether these oscillations are sustained. We perform an attractor analysis of this system using synchronous and three different asynchronous updating schemes both in the case of the unperturbed (wild-type) and perturbed (node-disrupted) systems. This analysis reveals that while the wild-type system possesses an update-independent fixed point, any oscillations eventually disappear unless strict constraints regarding the timing of certain processes and the initial state of the system are satisfied. Interestingly, in the case of disruption of a particular node all models lead to an extended attractor. Overall, our work provides a roadmap on how Boolean network modeling can be used as a predictive tool to uncover the dynamic patterns of a biological system under various internal and environmental perturbations.     

The Problem of Colliding Networks and its Relation to Cell Fusion and?Cancer

Cell fusion, a process that merges two or more cells into one, is required for normal development and has been explored as a tool for stem cell therapy. It has also been proposed that cell fusion causes cancer and contributes to its progression. These functions rely on a poorly understood ability of cell fusion to create new cell types. We suggest that this ability can be understood by considering cells as attractor networks whose basic property is to adopt a set of distinct, stable, self-maintaining states called attractors. According to this view, fusion of two cell types is a collision of two networks that have adopted distinct attractors. To learn how these networks reach a consensus, we model cell fusion computationally. To do so, we simulate patterns of gene activities using a formalism developed to simulate patterns of memory in neural networks. We find that the hybrid networks can assume attractors that are unrelated to parental attractors, implying that cell fusion can create new cell types by nearly instantaneously moving cells between attractors. We also show that hybrid networks are prone to assume spurious attractors, which are emergent and sporadic network states. This finding means that cell fusion can produce abnormal cell types, including cancerous types, by placing cells into normally inaccessible spurious states. Finally, we suggest that the problem of colliding networks has general significance in many processes represented by attractor networks, including biological, social, and political phenomena.     

An olfactory recognition model based on spatio-temporal encoding of odor quality in the olfactory bulb

In order to study the problem how the olfactory neural system processes the odorant molecular information for constructing the olfactory image of each object, we present a dynamic model of the olfactory bulb constructed on the basis of well-established experimental and theoretical results. The information relevant to a single odor, i.e. its constituent odorant molecules and their mixing ratios, are encoded into a spatio-temporal pattern of neural activity in the olfactory bulb, where the activity pattern corresponds to a limit cycle attractor in the mitral cell network. The spatio-temporal pattern consists of a temporal sequence of spatial firing patterns: each constituent molecule is encoded into a single spatial pattern, and the order of magnitude of the mixing ratio is encoded into the temporal sequence. The formation of a limit cycle attractor under the application of a novel odor is carried out based on the intensity-to-time-delay encoding scheme. The dynamic state of the olfactory bulb, which has learned many odors, becomes a randomly itinerant state in which the current firing state of the bulb itinerates randomly among limit cycle attractors corresponding to the learned odors. The recognition of an odor is generated by the dynamic transition in the network from the randomly itinerant state to a limit cycle attractor state relevant to the odor, where the transition is induced by the short-term synaptic changes made according to the Hebbian rule under the application of the odor stimulus. Received: 28 July 1997 / Accepted in revised form: 6 May 1998     

Stochastic Boolean networks: An efficient approach to modeling gene regulatory networks

ABSTRACT: BACKGROUND: Various computational models have been of interest due to their use in the modelling of gene regulatory networks (GRNs). As a logical model, probabilistic Boolean networks (PBNs) consider molecular and genetic noise, so the study of PBNs provides significant insights into the understanding of the dynamics of GRNs. This will ultimately lead to advances in developing therapeutic methods that intervene in the process of disease development and progression. The applications of PBNs, however, are hindered by the complexities involved in the computation of the state transition matrix and the steady-state distribution of a PBN. For a PBN with n genes and N Boolean networks, the complexity to compute the state transition matrix is O(nN22n) or O(nN2n) for a sparse matrix. RESULTS: This paper presents a novel implementation of PBNs based on the notions of stochastic logic and stochastic computation. This stochastic implementation of a PBN is referred to as a stochastic Boolean network (SBN). An SBN provides an accurate and efficient simulation of a PBN without and with random gene perturbation. The state transition matrix is computed in an SBN with a complexity of O(nL2n), where L is a factor related to the stochastic sequence length. Since the minimum sequence length required for obtaining an evaluation accuracy approximately increases in a polynomial order with the number of genes, n, and the number of Boolean networks, N, usually increases exponentially with n, L is typically smaller than N, especially in a network with a large number of genes. Hence, the computational complexity of an SBN is primarily limited by the number of genes, but not directly by the total possible number of Boolean networks. Furthermore, a time-frame expanded SBN enables an efficient analysis of the steady-state distribution of a PBN. These findings are supported by the simulation results of a simplified p53 network, several randomly generated networks and a network inferred from a T cell immune response dataset. An SBN can also implement the function of an asynchronous PBN and is potentially useful in a hybrid approach in combination with a continuous or single-molecule level stochastic model. CONCLUSIONS: Stochastic Boolean networks (SBNs) are proposed as an efficient approach to modelling gene regulatory networks (GRNs). The SBN approach is able to recover biologically-proven regulatory behaviours, such as the oscillatory dynamics of the p53-Mdm2 network and the dynamic attractors in a T cell immune response network. The proposed approach can further predict the network dynamics when the genes are under perturbation, thus providing biologically meaningful insights for a better understanding of the dynamics of GRNs. The algorithms and methods described in this paper have been implemented in Matlab packages, which are attached as Additional files.     

Combined effects of LTP/LTD and synaptic scaling in formation of discrete and line attractors with persistent activity from non-trivial baseline

In this article, we analyze combined effects of LTP/LTD and synaptic scaling and study the creation of persistent activity from a periodic or chaotic baseline attractor. The bifurcations leading to the creation of new attractors have been detailed; this was achieved using a mean field approximation. Attractors encoding persistent activity can notably appear via generalized period-doubling bifurcations, tangent bifurcations of the second iterates or boundary crises, after which the basins of attraction become irregular. Synaptic scaling is shown to maintain the coexistence of a state of persistent activity and the baseline. According to the rate of change of the external inputs, different types of attractors can be formed: line attractors for rapidly changing external inputs and discrete attractors for constant external inputs.     

Uncovering the profile of somatic mtDNA mutations in Chinese colorectal cancer patients

In the past decade, a high incidence of somatic mitochondrial DNA (mtDNA) mutations has been observed, mostly based on a fraction of the molecule, in various cancerous tissues; nevertheless, some of them were queried due to problems in data quality. Obviously, without a comprehensive understanding of mtDNA mutational profile in the cancerous tissue of a specific patient, it is unlikely to disclose the genuine relationship between somatic mtDNA mutations and tumorigenesis. To achieve this objective, the most straightforward way is to directly compare the whole mtDNA genome variation among three tissues (namely, cancerous tissue, para-cancerous tissue, and distant normal tissue) from the same patient. Considering the fact that most of the previous studies on the role of mtDNA in colorectal tumor focused merely on the D-loop or partial segment of the molecule, in the current study we have collected three tissues (cancerous, para-cancerous and normal tissues) respectively recruited from 20 patients with colorectal tumor and completely sequenced the mitochondrial genome of each tissue. Our results reveal a relatively lower incidence of somatic mutations in these patients; intriguingly, all somatic mutations are in heteroplasmic status. Surprisingly, the observed somatic mutations are not restricted to cancer tissues, for the para-cancer tissues and distant normal tissues also harbor somatic mtDNA mutations with a lower frequency than cancerous tissues but higher than that observed in the general population. Our results suggest that somatic mtDNA mutations in cancerous tissues could not be simply explained as a consequence of tumorigenesis; meanwhile, the somatic mtDNA mutations in normal tissues might reflect an altered physiological environment in cancer patients.     

An efficient algorithm for identifying primary phenotype attractors of a large-scale Boolean network

Background

Boolean network modeling has been widely used to model large-scale biomolecular regulatory networks as it can describe the essential dynamical characteristics of complicated networks in a relatively simple way. When we analyze such Boolean network models, we often need to find out attractor states to investigate the converging state features that represent particular cell phenotypes. This is, however, very difficult (often impossible) for a large network due to computational complexity.

Results

There have been some attempts to resolve this problem by partitioning the original network into smaller subnetworks and reconstructing the attractor states by integrating the local attractors obtained from each subnetwork. But, in many cases, the partitioned subnetworks are still too large and such an approach is no longer useful. So, we have investigated the fundamental reason underlying this problem and proposed a novel efficient way of hierarchically partitioning a given large network into smaller subnetworks by focusing on some attractors corresponding to a particular phenotype of interest instead of considering all attractors at the same time. Using the definition of attractors, we can have a simplified update rule with fixed state values for some nodes. The resulting subnetworks were small enough to find out the corresponding local attractors which can be integrated for reconstruction of the global attractor states of the original large network.

Conclusions

The proposed approach can substantially extend the current limit of Boolean network modeling for converging state analysis of biological networks.