Characterizing multistationarity regimes in biochemical reaction networks

Switch like responses appear as common strategies in the regulation of cellular systems. Here we present a method to characterize bistable regimes in biochemical reaction networks that can be of use to both direct and reverse engineering of biological switches. In the design of a synthetic biological switch, it is important to study the capability for bistability of the underlying biochemical network structure. Chemical Reaction Network Theory (CRNT) may help at this level to decide whether a given network has the capacity for multiple positive equilibria, based on their structural properties. However, in order to build a working switch, we also need to ensure that the bistability property is robust, by studying the conditions leading to the existence of two different steady states. In the reverse engineering of biological switches, knowledge collected about the bistable regimes of the underlying potential model structures can contribute at the model identification stage to a drastic reduction of the feasible region in the parameter space of search. In this work, we make use and extend previous results of the CRNT, aiming not only to discriminate whether a biochemical reaction network can exhibit multiple steady states, but also to determine the regions within the whole space of parameters capable of producing multistationarity. To that purpose we present and justify a condition on the parameters of biochemical networks for the appearance of multistationarity, and propose an efficient and reliable computational method to check its satisfaction through the parameter space.     

Modular response analysis of cellular regulatory networks

The sheer complexity of intracellular regulatory networks, which involve signal transducing, metabolic, and genetic circuits, hampers our ability to carry out a quantitative analysis of their functions. Here, we describe an approach that greatly simplifies this type of analysis by capitalizing on the modular organization of such networks. Steady-state responses of the network as a whole are accounted for in terms of intermodular interactions between the modules alone; processes operating solely within modules need not be considered when analysing signal transfer through the entire network. The intermodular interactions are quantified through (local) response coefficients which populate an interaction map (matrix). This matrix can be derived from a biochemical or molecular biological analysis of (macro) molecular interactions that constitute the regulatory network. The approach is illustrated by two examples: (i) mitogenic signalling through the mitogen-activated protein kinase cascade in the epidermal growth factor receptor network and (ii) regulation of ammonium assimilation in Escherichia coli.     

From Boolean Network Model to Continuous Model Helps in Design of Functional Circuits

Computational circuit design with desired functions in a living cell is a challenging task in synthetic biology. To achieve this task, numerous methods that either focus on small scale networks or use evolutionary algorithms have been developed. Here, we propose a two-step approach to facilitate the design of functional circuits. In the first step, the search space of possible topologies for target functions is reduced by reverse engineering using a Boolean network model. In the second step, continuous simulation is applied to evaluate the performance of these topologies. We demonstrate the usefulness of this method by designing an example biological function: the SOS response of E. coli. Our numerical results show that the desired function can be faithfully reproduced by candidate networks with different parameters and initial conditions. Possible circuits are ranked according to their robustness against perturbations in parameter and gene expressions. The biological network is among the candidate networks, yet novel designs can be generated. Our method provides a scalable way to design robust circuits that can achieve complex functions, and makes it possible to uncover design principles of biological networks.     

Rational design of complex phenotype via network models

We demonstrate a modeling and computational framework that allows for rapid screening of thousands of potential network designs for particular dynamic behavior. To illustrate this capability we consider the problem of hysteresis, a prerequisite for construction of robust bistable switches and hence a cornerstone for construction of more complex synthetic circuits. We evaluate and rank most three node networks according to their ability to robustly exhibit hysteresis where robustness is measured with respect to parameters over multiple dynamic phenotypes. Focusing on the highest ranked networks, we demonstrate how additional robustness and design constraints can be applied. We compare our results to more traditional methods based on specific parameterization of ordinary differential equation models and demonstrate a strong qualitative match at a small fraction of the computational cost.     

Efficient parametric analysis of the chemical master equation through model order reduction

ABSTRACT: BACKGROUND: Stochastic biochemical reaction networks are commonly modelled by the chemical master equation, and can be simulated as first order linear differential equations through a finite state projection. Due to the very high state space dimension of these equations, numerical simulations are computationally expensive. This is a particular problem for analysis tasks requiring repeated simulations for different parameter values. Such tasks are computationally expensive to the point of infeasibility with the chemical master equation. RESULTS: In this article, we apply parametric model order reduction techniques in order to construct accurate low-dimensional parametric models of the chemical master equation. These surrogate models can be used in various parametric analysis task such as identifiability analysis, parameter estimation, or sensitivity analysis. As biological examples, we consider two models for gene regulation networks, a bistable switch and a network displaying stochastic oscillations. CONCLUSIONS: The results show that the parametric model reduction yields efficient models of stochastic biochemical reaction networks, and that these models can be useful for systems biology applications involving parametric analysis problems such as parameter exploration, optimization, estimation or sensitivity analysis.     

Coupled feedback loops form dynamic motifs of cellular networks

Cellular networks are composed of complicated interconnections among components, and some subnetworks of particular functioning are often identified as network motifs. Among such network motifs, feedback loops are thought to play important dynamical roles. Intriguingly, such feedback loops are very often found as a coupled structure in cellular circuits. Therefore, we integrated all the scattered information regarding the coupled feedbacks in various cellular circuits and investigated the dynamical role of each coupled structure. Finally, we discovered that coupled positive feedbacks enhance signal amplification and bistable characteristics; coupled negative feedbacks realize enhanced homeostasis; coupled positive and negative feedbacks enable reliable decision-making by properly modulating signal responses and effectively dealing with noise.     

Stochasticity,Bistability and the Wisdom of Crowds: A Model for Associative Learning in Genetic Regulatory Networks

It is generally believed that associative memory in the brain depends on multistable synaptic dynamics, which enable the synapses to maintain their value for extended periods of time. However, multistable dynamics are not restricted to synapses. In particular, the dynamics of some genetic regulatory networks are multistable, raising the possibility that even single cells, in the absence of a nervous system, are capable of learning associations. Here we study a standard genetic regulatory network model with bistable elements and stochastic dynamics. We demonstrate that such a genetic regulatory network model is capable of learning multiple, general, overlapping associations. The capacity of the network, defined as the number of associations that can be simultaneously stored and retrieved, is proportional to the square root of the number of bistable elements in the genetic regulatory network. Moreover, we compute the capacity of a clonal population of cells, such as in a colony of bacteria or a tissue, to store associations. We show that even if the cells do not interact, the capacity of the population to store associations substantially exceeds that of a single cell and is proportional to the number of bistable elements. Thus, we show that even single cells are endowed with the computational power to learn associations, a power that is substantially enhanced when these cells form a population.     

Ingeneue: a versatile tool for reconstituting genetic networks,with examples from the segment polarity network

Here we describe a software tool for synthesizing molecular genetic data into models of genetic networks. Our software program Ingeneue, written in Java, lets the user quickly turn a map of a genetic network into a dynamical model consisting of a set of ordinary differential equations. We developed Ingeneue as part of an ongoing effort to explore the design and evolvability of genetic networks. Ingeneue has three principal advantages over other available mathematical software: it automates instantiation of the same network model in each cell in a 2-D sheet of cells; it constructs model equations from pre-made building blocks corresponding to common biochemical processes; and it automates searches through parameter space, sensitivity analyses, and other common tasks. Here we discuss the structure of the software and some of the issues we have dealt with. We conclude with some examples of results we have achieved with Ingeneue for the Drosophila segment polarity network.     

Dual selection of a genetic switch by a single selection marker

Forward engineering of synthetic genetic circuits in living cells is expected to deliver various applications in biotechnology and medicine and to provide valuable insights into the design principles of natural gene networks. However, lack of biochemical data and complexity of biological environment complicate rational design of such circuits based on quantitative simulation. Previously, we have shown that directed evolution can complement our weakness in designing genetic circuits by screening or selecting functional circuits from a large pool of nonfunctional ones. Here we describe a dual selection strategy that allows selection of both ON and OFF states of genetic circuits using tetA as a single selection marker. We also describe a successful demonstration of a genetic switch selection from a 2000-fold excess background of nonfunctional switches in three rounds of iterative selection. The dual selection system is more robust than the previously reported selection system employing three genes, with no observed false positive mutants during the simulated selections.     

Network design meets in silico evolutionary biology

Cell fate is programmed through gene regulatory networks that perform several calculations to take the appropriate decision. In silico evolutionary optimization mimics the way Nature has designed such gene regulatory networks. In this review we discuss the basic principles of these evolutionary approaches and how they can be applied to engineer synthetic networks. We summarize the basic guidelines to implement an in silico evolutionary design method, the operators for mutation and selection that iteratively drive the network architecture towards a specified dynamical behavior. Interestingly, as it happens in natural evolution, we show the existence of patterns of punctuated evolution. In addition, we highlight several examples of models that have been designed using automated procedures, together with different objective functions to select for the proper behavior. Finally, we briefly discuss the modular designability of gene regulatory networks and its potential application in biotechnology.     

Methods for and results from the study of design principles in molecular systems

Most aspects of molecular biology can be understood in terms of biological design principles. These principles can be loosely defined as qualitative and quantitative features that emerge in evolution and recur more frequently than one would expect by chance alone in biological systems that perform a given type of process or function. Furthermore, such recurrence can be rationalized in terms of the functional advantage that the design provides to the system when compared with possible alternatives. This paper focuses on those design features that can be related to improved functional effectiveness of molecular and regulatory networks. We begin by reviewing assumptions and methods that underlie the study of such principles in molecular networks. We follow by discussing many of the design principles that have been found in genetic, metabolic, and signal transduction circuits. We concentrate mainly on results in the context of Biochemical Systems Theory, although we also briefly discuss other work. We conclude by discussing the importance of these principles for both, understanding the natural evolution of complex networks at the molecular level and for creating artificial biological systems with specific features.     

Proteomics of calcium-signaling components in plants

Calcium functions as a versatile messenger in mediating responses to hormones, biotic/abiotic stress signals and a variety of developmental cues in plants. The Ca(2+)-signaling circuit consists of three major "nodes"--generation of a Ca(2+)-signature in response to a signal, recognition of the signature by Ca2+ sensors and transduction of the signature message to targets that participate in producing signal-specific responses. Molecular genetic and protein-protein interaction approaches together with bioinformatic analysis of the Arabidopsis genome have resulted in identification of a large number of proteins at each "node"--approximately 80 at Ca2+ signature, approximately 400 sensors and approximately 200 targets--that form a myriad of Ca2+ signaling networks in a "mix and match" fashion. In parallel, biochemical, cell biological, genetic and transgenic approaches have unraveled functions and regulatory mechanisms of a few of these components. The emerging paradigm from these studies is that plants have many unique Ca2+ signaling proteins. The presence of a large number of proteins, including several families, at each "node" and potential interaction of several targets by a sensor or vice versa are likely to generate highly complex networks that regulate Ca(2+)-mediated processes. Therefore, there is a great demand for high-throughput technologies for identification of signaling networks in the "Ca(2+)-signaling-grid" and their roles in cellular processes. Here we discuss the current status of Ca2+ signaling components, their known functions and potential of emerging high-throughput genomic and proteomic technologies in unraveling complex Ca2+ circuitry.     

Competing to coordinate cell fate decisions: the MST2-Raf-1 signaling device

How do biochemical signaling pathways generate biological specificity? This question is fundamental to modern biology, and its enigma has been accentuated by the discovery that most proteins in signaling networks serve multifunctional roles. An answer to this question may lie in analyzing network properties rather than individual traits of proteins in order to elucidate design principles of biochemical networks that enable biological decision-making. We discuss how this is achieved in the MST2/Hippo-Raf-1 signaling network with the help of mathematical modeling and model-based analysis, which showed that competing protein interactions with affinities controlled by dynamic protein modifications can function as Boolean computing devices that determine cell fate decisions. In addition, we discuss areas of interest for future research and highlight how systems approaches would be of benefit.     

Emergence of bimodal cell population responses from the interplay between analog single-cell signaling and protein expression noise

ABSTRACT: BACKGROUND: Cell-to-cell variability in protein expression can be large, and its propagation through signaling networks affects biological outcomes. Here, we apply deterministic and probabilistic models and biochemical measurements to study how network topologies and cell-to-cell protein abundance variations interact to shape signaling responses. RESULTS: We observe bimodal distributions of extracellular signal-regulated kinase (ERK) responses to epidermal growth factor (EGF) stimulation, which are generally thought to indicate bistable or ultrasensitive signaling behavior in single cells. Surprisingly, we find that a simple MAPK/ERK-cascade model with negative feedback that displays graded, analog ERK responses at a single cell level can explain the experimentally observed bimodality at the cell population level. Model analysis suggests that a conversion of graded input--output responses in single cells to digital responses at the population level is caused by a broad distribution of ERK pathway activation thresholds brought about by cell-to-cell variability in protein expression. CONCLUSIONS: Our results show that bimodal signaling response distributions do not necessarily imply digital (ultrasensitive or bistable) single cell signaling, and the interplay between protein expression noise and network topologies can bring about digital population responses from analog single cell dose responses. Thus, cells can retain the benefits of robustness arising from negative feedback, while simultaneously generating population-level on/off responses that are thought to be critical for regulating cell fate decisions.     

In vitro regulatory models for systems biology

The reductionist approach has revolutionized biology in the past 50 years. Yet its limits are being felt as the complexity of cellular interactions is gradually revealed by high-throughput technology. In order to make sense of the deluge of “omic data”, a hypothesis-driven view is needed to understand how biomolecular interactions shape cellular networks. We review recent efforts aimed at building in vitro biochemical networks that reproduce the flow of genetic regulation. We highlight how those efforts have culminated in the rational construction of biochemical oscillators and bistable memories in test tubes. We also recapitulate the lessons learned about in vivo biochemical circuits such as the importance of delays and competition, the links between topology and kinetics, as well as the intriguing resemblance between cellular reaction networks and ecosystems.