Monte Carlo simulation of a simple gene network yields new evolutionary insights

Modeling the evolution of molecular systems from a mechanistic perspective

Systems biology-inspired genotype-phenotype mapping models are increasingly being used to study the evolutionary properties of molecular biological systems, in particular the general emergent properties of evolving systems, such as modularity, robustness, and evolvability. However, the level of abstraction at which many of these models operate might not be sufficient to capture all relevant intricacies of biological evolution in sufficient detail. Here, we argue that in particular gene and genome duplications, both evolutionary mechanisms of potentially major importance for the evolution of molecular systems and of special relevance to plant evolution, are not adequately accounted for in most GPM modeling frameworks, and that more fine-grained mechanistic models may significantly advance understanding of how gen(om)e duplication impacts molecular systems evolution.

Hybrid T-helper cells: stabilizing the moderate center in a polarized system

Polarization of cell phenotypes, a common strategy to achieve cell type diversity in metazoa, results from binary cell-fate decisions in the branching pedigree of development. Such "either-or" fate decisions are controlled by two opposing cell fate-determining transcription factors. Each of the two distinct "master regulators" promotes differentiation of its respective sister lineage. But they also suppress one other, leading to their mutually exclusive expression in the two ensuing lineages. Thus, promiscuous coexistence of the antagonist regulators in the same cell, the hallmark of the common "undecided" progenitor of two sister lineages, is considered unstable. This antagonism ensures robust polarization into two discretely distinct cell types. But now the immune system's T-helper (Th) cells and their two canonical subtypes, Th1 and Th2 cells, tell a different story, as revealed in three papers recently published in PLOS Biology. The intermediate state that co-expresses the two opposing master regulators of the Th1 and Th2 subtypes, T-bet and Gata3, is highly stable and is not necessarily an undecided precursor. Instead, the Th1/Th2 hybrid cell is a robust new type with properties of both Th1 and Th2 cells. These hybrid cells are functionally active and possess the benefit of moderation: self-limitation of effector T cell function to prevent excessive inflammation, a permanent risk in host defense that can cause tissue damage or autoimmunity. Gene regulatory network analysis suggests that stabilization of the intermediate center in a polarizing system can be achieved by minor tweaking of the architecture of the mutual suppression gene circuit, and thus is a design option readily available to evolution.

A general model for binary cell fate decision gene circuits with degeneracy: indeterminacy and switch behavior in the absence of cooperativity

Background

The gene regulatory circuit motif in which two opposing fate-determining transcription factors inhibit each other but activate themselves has been used in mathematical models of binary cell fate decisions in multipotent stem or progenitor cells. This simple circuit can generate multistability and explains the symmetric “poised” precursor state in which both factors are present in the cell at equal amounts as well as the resolution of this indeterminate state as the cell commits to either cell fate characterized by an asymmetric expression pattern of the two factors. This establishes the two alternative stable attractors that represent the two fate options. It has been debated whether cooperativity of molecular interactions is necessary to produce such multistability.

Principal Findings

Here we take a general modeling approach and argue that this question is not relevant. We show that non-linearity can arise in two distinct models in which no explicit interaction between the two factors is assumed and that distinct chemical reaction kinetic formalisms can lead to the same (generic) dynamical system form. Moreover, we describe a novel type of bifurcation that produces a degenerate steady state that can explain the metastable state of indeterminacy prior to cell fate decision-making and is consistent with biological observations.

Conclusion

The general model presented here thus offers a novel principle for linking regulatory circuits with the state of indeterminacy characteristic of multipotent (stem) cells.

Multistable decision switches for flexible control of epigenetic differentiation

It is now recognized that molecular circuits with positive feedback can induce two different gene expression states (bistability) under the very same cellular conditions. Whether, and how, cells make use of the coexistence of a larger number of stable states (multistability) is however largely unknown. Here, we first examine how autoregulation, a common attribute of genetic master regulators, facilitates multistability in two-component circuits. A systematic exploration of these modules' parameter space reveals two classes of molecular switches, involving transitions in bistable (progression switches) or multistable (decision switches) regimes. We demonstrate the potential of decision switches for multifaceted stimulus processing, including strength, duration, and flexible discrimination. These tasks enhance response specificity, help to store short-term memories of recent signaling events, stabilize transient gene expression, and enable stochastic fate commitment. The relevance of these circuits is further supported by biological data, because we find them in numerous developmental scenarios. Indeed, many of the presented information-processing features of decision switches could ultimately demonstrate a more flexible control of epigenetic differentiation.

Extinction and resurrection in gene networks

When gene regulatory networks operate in regimes where the number of protein molecules is so small that the molecular species are on the verge of extinction, the death and resurrection of the species greatly modifies the attractor landscape. Deterministic models and the diffusion approximation to the master equation break down at the limits of protein populations in a way very analogous to the breakdown of geometrical optics that occurs at distances <1 wavelength of light from edges. Stable stochastic attractors arise from extinction and resurrection events that are not predicted by the deterministic description. With this view, we explore the attractors of the regular toggle switch and the exclusive switch, focusing on the effects of cooperative binding and the production of protein in bursts. Our arguments suggest that the stability of lysogeny in the lambda-phage may be influenced by such extinction phenomena.