Buffered Qualitative Stability explains the robustness and evolvability of transcriptional networks

The gene regulatory network (GRN) is the central decision‐making module of the cell. We have developed a theory called Buffered Qualitative Stability (BQS) based on the hypothesis that GRNs are organised so that they remain robust in the face of unpredictable environmental and evolutionary changes. BQS makes strong and diverse predictions about the network features that allow stable responses under arbitrary perturbations, including the random addition of new connections. We show that the GRNs of E. coli, M. tuberculosis, P. aeruginosa, yeast, mouse, and human all verify the predictions of BQS. BQS explains many of the small- and large‐scale properties of GRNs, provides conditions for evolvable robustness, and highlights general features of transcriptional response. BQS is severely compromised in a human cancer cell line, suggesting that loss of BQS might underlie the phenotypic plasticity of cancer cells, and highlighting a possible sequence of GRN alterations concomitant with cancer initiation. DOI: http://dx.doi.org/10.7554/eLife.02863.00

Cancer modelling: Getting to the heart of the problem

Paradoxically, improvements in healthcare that have enhanced the life expectancy of humans in the Western world have, indirectly, increased the prevalence of certain types of cancer such as prostate and breast. It remains unclear whether this phenomenon should be attributed to the ageing process itself or the cumulative effect of prolonged exposure to harmful environmental stimuli such as ultraviolet light, radiation and carcinogens (Franks and Teich, 1988). Equally, there is also compelling evidence that certain genetic abnormalities can predispose individuals to specific cancers (Ilyas et al., 1999). The variety of factors that have been implicated in the development of solid tumours stems, to a large extent, from the fact that ‘cancer’ is a generic term, often used to characterize a series of disorders that share common features. At this generic level of description, cancer may be viewed as a cellular disease in which controls that usually regulate growth and maintain homeostasis are disrupted. Cancer is typically initiated by genetic mutations that lead to enhanced mitosis of a cell lineage and the formation of an avascular tumour. Since it receives nutrients by diffusion from the surrounding tissue, the size of an avascular tumour is limited to several millimeters in diameter. Further growth relies on the tumour acquiring the ability to stimulate the ingrowth of a new, circulating blood supply from the host vasculature via a process termed angiogenesis (Folkman, 1974). Once vascularised, the tumour has access to a vast nutrient source and rapid growth ensues. Further, tumour fragments that break away from the primary tumour, on entering the vasculature, may be transported to other organs in which they may establish secondary tumours or metastases that further compromise the host. Invasion is another key feature of solid tumours whereby contact with the tissue stimulates the production of enzymes that digest the tissue, liberating space into which the tumour cells migrate. Thus, cancer is a complex, multiscale process. The spatial scales of interest range from the subcellular level, to the cellular and macroscopic (or tissue) levels while the timescales may vary from seconds (or less) for signal transduction pathways to months for tumour doubling times The variety of phenomena involved, the range of spatial and temporal scales over which they act and the complex way in which they are inter-related mean that the development of realistic theoretical models of solid tumour growth is extremely challenging. While there is now a large literature focused on modelling solid tumour growth (for a review, see, for example, Preziosi, 2003), existing models typically focus on a single spatial scale and, as a result, are unable to address the fundamental problem of how phenomena at different scales are coupled or to combine, in a systematic manner, data from the various scales. In this article, a theoretical framework will be presented that is capable of integrating a hierarchy of processes occurring at different scales into a detailed model of solid tumour growth (Alarcon et al., 2004). The model is formulated as a hybrid cellular automaton and contains interlinked elements that describe processes at each spatial scale: progress through the cell cycle and the production of proteins that stimulate angiogenesis are accounted for at the subcellular level; cell-cell interactions are treated at the cellular level; and, at the tissue scale, attention focuses on the vascular network whose structure adapts in response to blood flow and angiogenic factors produced at the subcellular level. Further coupling between the different spatial scales arises from the transport of blood-borne oxygen into the tissue and its uptake at the cellular level. Model simulations will be presented to illustrate the effect that spatial heterogeneity induced by blood flow through the vascular network has on the tumour’s growth dynamics and explain how the model may be used to compare the efficacy of different anti-cancer treatment protocols

Systems level circuit model of C. elegans undulatory locomotion: mathematical modeling and molecular genetics

To establish the relationship between locomotory behavior and dynamics of neural circuits in the nematode C. elegans we combined molecular and theoretical approaches. In particular, we quantitatively analyzed the motion of C. elegans with defective synaptic GABA and acetylcholine transmission, defective muscle calcium signaling, and defective muscles and cuticle structures, and compared the data with our systems level circuit model. The major experimental findings are: (i) anterior-to-posterior gradients of body bending flex for almost all strains both for forward and backward motion, and for neuronal mutants, also analogous weak gradients of undulatory frequency, (ii) existence of some form of neuromuscular (stretch receptor) feedback, (iii) invariance of neuromuscular wavelength, (iv) biphasic dependence of frequency on synaptic signaling, and (v) decrease of frequency with increase of the muscle time constant. Based on (i) we hypothesize that the Central Pattern Generator (CPG) is located in the head both for forward and backward motion. Points (i) and (ii) are the starting assumptions for our theoretical model, whose dynamical patterns are qualitatively insensitive to the details of the CPG design if stretch receptor feedback is sufficiently strong and slow. The model reveals that stretch receptor coupling in the body wall is critical for generation of the neuromuscular wave. Our model agrees with our behavioral data(iii), (iv), and (v), and with other pertinent published data, e.g., that frequency is an increasing function of muscle gap-junction coupling.Comment: Neural control of C. elegans motion with genetic perturbation

An excitable gene regulatory circuit induces transient cellular differentiation

Certain types of cellular differentiation are probabilistic and transient. In such systems individual cells can switch to an alternative state and, after some time, switch back again. In Bacillus subtilis, competence is an example of such a transiently differentiated state associated with the capability for DNA uptake from the environment. Individual genes and proteins underlying differentiation into the competent state have been identified, but it has been unclear how these genes interact dynamically in individual cells to control both spontaneous entry into competence and return to vegetative growth. Here we show that this behaviour can be understood in terms of excitability in the underlying genetic circuit. Using quantitative fluorescence time-lapse microscopy, we directly observed the activities of multiple circuit components simultaneously in individual cells, and analysed the resulting data in terms of a mathematical model. We find that an excitable core module containing positive and negative feedback loops can explain both entry into, and exit from, the competent state. We further tested this model by analysing initiation in sister cells, and by re-engineering the gene circuit to specifically block exit. Excitable dynamics driven by noise naturally generate stochastic and transient responses, thereby providing an ideal mechanism for competence regulation

The Molecular Clockwork of a Protein-Based Circadian Oscillator

The circadian clock of the cyanobacterium Synechococcuselongatus PCC 7942 is governed by a core oscillator consisting of the proteins KaiA, KaiB, and KaiC. Remarkably, circadian oscillations in the phosphorylation state of KaiC can be reconstituted in a test tube by mixing the three Kai proteins and adenosine triphosphate. The in vitro oscillator provides a well-defined system in which experiments can be combined with mathematical analysis to understand the mechanism of a highly robust biological oscillator. In this Review, we summarize the biochemistry of the Kai proteins and examine models that have been proposed to explain how oscillations emerge from the properties of the oscillator’s constituents.Molecular and Cellular Biolog

Modeling Non-Linear Dynamic Phenomena in Biochemical Networks

Facilitated by the development of high-throughput techniques, the focus of biological research has changed in the last decades from the investigation of single cell components to a system-level approach, which aims at an understanding of interactions between these cell components. This objective requires modeling and analysis methods for these regulatory networks. In this thesis, we investigate mechanisms causing qualitative dynamic behaviors of regulatory subsystems. For this purpose, we introduce a differential equation model based on underlying molecular binding reactions, whose parameters are estimated using time series concentration data. In the first part, the model is applied to subsystems with qualitatively different dynamic behaviors: The response of the Mycobacterium tuberculosis to DNA damages is described as the relaxation of a system to its steady state after external perturbation. Specific repression of genes in Escherichia coli by the global regulator protein H-NS is explained by the interrelation of feedback mechanisms. In order to prevent overfitting, a typical problem in network inference from experimental data, we introduce an approach based on Bayesian statistics, which includes prior knowledge about the system in terms of prior probability distributions. This approach is applied to simulated data and to the regulatory network of the Saccharomyces cerevisiae cell cycle. Motivated by results on the yeast cell cycle, the second part of this thesis investigates the robustness of periodic behavior in regulatory networks. The model presented belongs to a class of differential equations whose solutions tend to converge to a steady state. Accordingly, periodic behavior is not robust with respect to parameter variations. We explain this phenomenon by applying a bifurcation analysis and investigating the stability of steady states. It is shown that large time scale differences and an inclusion of time-delays can stabilize sustained oscillations, and we postulate that they are important to maintain oscillations in biological systems