Briefings in Bioinformatics Advance Access published online on July 11, 2007
Briefings in Bioinformatics, doi:10.1093/bib/bbm028
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Dynamical roles of biological regulatory circuits1
Corresponding author. Denis Thieffry, INSERM ERM 206 - TAGC & Université de la Méditerranée, Campus Scientifique de Luminy, case 928, 13288 Marseille, Cedex 09 France. E-mail: thieffry{at}tagc.univ-mrs.fr
| ABSTRACT |
|---|
|
|
|---|
Regulatory circuits are found at the basis of all non-trivial dynamical properties of biological networks. More specifically, positive circuits are involved in the generation of multiple differentiated states, whereas negative circuits can generate cyclic or homeostatic behaviours. These notions are briefly reviewed, from initial biological formulations to mathematical formalisations, further encompassing their application to the design of synthetic regulatory systems. Finally, current challenges for the analysis of increasingly complex regulatory networks are indicated, as well as prospects for our understanding of development and evolution.
Keywords: Regulatory circuits, multistability, homeostasis, oscillations, synthetic biology, systems biology, dynamics
| BIOLOGICAL FEEDBACKS |
|---|
|
|
|---|
Most remarkable properties of biological regulatory networks can be related to the occurrence of (generally non-linear) regulatory feedbacks. In biochemistry, the role of feedback inhibition of catalytic enzymes by biosynthetic products has been underlined already 50 years ago [1]. In parallel, positive feedback was postulated to account for enzymatic induction phenomena in bacteria [2]. Soon, the distinction between regulatory and structural genes, together with the specification of a concrete mechanism for gene regulation (operon model), led Monod and Jacob [3] to imagine different regulatory schemes accounting for the generation of multiple cell differentiation states by a single genotype (Figure 1).
|
During the 1970s, molecular analyses and modelling studies progressively enabled general characterisations of simple biochemical and genetic networks [46].
Notably, proper definition and classification of regulatory circuits (or feedback loops) were proposed by Thomas [7], together with general rules about circuit requirement to generate non-trivial dynamical behaviour. Indeed, defined as simple circular chains of oriented interactions, feedback circuits can be classified into positive or negative circuits, depending on the parity of the number of negative interactions (i.e. the sign of a circuit is given by the product of the signs of its constitutive interactions). Positive circuits are necessary to generate alternative cellular states (multiple attractors, using the terminology of Dynamical Systems), whereas negative circuits are needed to generate homeostasis or sustained oscillatory behaviour (Figure 2; cf. [8] for an extensive review).
|
| FORMAL DEFINITIONS OF FEEDBACK CIRCUITS |
|---|
|
|
|---|
The biological regulatory schemes found in the literature encompass arrows representing different types of molecular processes, e.g. metabolic reactions, transcriptional regulation, degradation, transport, etc. In order to investigate the dynamical roles of regulatory circuits, such regulatory schemes have to be translated into an homogenous mathematical formalism. In the present context, modellers often refer to a qualitative formalisation in terms of graphs or logical equations, or yet, to a quantitative representation in terms of ordinary differential equations (ODEs).
In the case of a graph-based representation, the notion of regulatory circuit directly derives from the graph-theoretic circuit concept (interaction signs are then defined as arc labels). For given logical rules specifying the behaviour of each regulatory components, it is possible to induce the regulatory interactions, including their signs (positive, negative or yet dual, i.e. context-dependent sign) (cf. [9]; for another qualitative approach using piece-wise differential equations, see [10]).
In the differential case, regulatory circuits can be rigorously defined on the basis of the Jacobian matrix of the ODE system (
, stating that the evolution or the derivative of a component i depends on the levels of other components of the network, represented by a state vector
). Indeed, the terms (aij =
Fi /
xj) of this matrix provides information about pair-wise influences (from regulatory component j onto component i). Consequently, any set of non-zero elements aij of the Jacobian matrix of the differential system, such that the i (row) and j (column) indices form a circular permutation, defines a regulatory circuit, whose sign is given by the product of the signs of these elements (cf. example in Figure 3). Note that, for most biological models, the consideration of passive, linear degradation of molecular components gives rise to negative diagonal terms and thus to auto-inhibitory circuits.
|
The Jacobian matrix is most often defined at the level of a specific steady state to determine its dynamical properties (linear stability analysis; see e.g. [11]). In this respect, it is worth noting that the mathematical objects used for steady state characterization (characteristic equation, eigenvalues) exclusively rely on terms involved in feedback circuits. However, this does not impede other terms to play important roles in the location or the number of steady states (for a more detailed discussion on the relationship between the Jacobian matrix and steady state properties, see [12]).
| FROM THOMAS RULES TO MATHEMATICAL THEOREMS AND DEMONSTRATIONS |
|---|
|
|
|---|
Thomas rules [7] have attracted the interest of a number of mathematicians, who translated them into proper theorems. In the differential context, a series of increasingly general theorems have been proposed, stating the necessity of a positive circuit to generate multistationarity [1317], or the necessity of a negative circuit to generate a cyclic attractor [1416]. More recently, the same rules have been demonstrated in a biologically relevant, discrete, multilevel framework, considering transitions between discrete states affecting at most one component at a time, switching it to a neighbouring value (i.e. the Hamming distance between two following states is exactly 1) [9, 18].
Lately, Thomas has proposed a third rule, stating that a negative circuit is necessary to have an attractor in the case of a differential system [12].
Altogether, these rules and theorems specify necessary conditions, which can be directly applied to specific networks. For example, whatever the complexity of the network considered, if it encompasses no positive circuit, one can directly conclude that the system can generate at most one attractor. Similarly, in the absence of negative circuit, one cannot expect a cyclic attractor.
However, the sole presence of a circuit in a regulatory graph (or in the Jacobian matrix of a differential system) does not necessarily imply the corresponding dynamical behaviour. Indeed, it is well known that at least some non-linearity is further necessary in continuous systems, as well as specific constraints on relevant parameters (kinetic parameters in the case of differential systems, logical rules directing component behaviours in the case of discrete systems). In the differential framework, these conditions can be established through the analysis of the sensitivity of specific state properties depending on crucial parameter values (parameter bifurcation analysis, cf. [19] and [20]). In the discrete case, it is possible to specify necessary conditions to obtain non-trivial dynamical behaviour with isolated (functional) feedback circuits, whatever the number and signs of involved interactions (see examples in Figure 4) [21].
|
| ENGINEERING OF SYNTHETIC REGULATORY CIRCUITS |
|---|
|
|
|---|
The appreciation of the crucial dynamical roles of regulatory circuits in real systems was at the root of the design of synthetic genetic circuits endowed with specific dynamical properties, including a cross-inhibitory positive circuit [22], a negative circuit made of the sequential inhibitions of three repressor genes [23], as well as a simple auto-inhibitory circuit [24] (see Figure 5; an example of circuit engineering in Yeast is found in [25]; for a review, see [26]).
|
A whole new field is presently emerging from these seminal studies, Synthetic Biology, which combines relatively standard molecular genetic protocols with mathematical modelling, under an engineering perspective.
In this respect, the MIT promotes a yearly international competition (iGEM) on biological circuit engineering, and organises the distribution of a growing collection of basic molecular genetic components or modules, which can be further combined with self-engineered regulatory modules in vivo to produce more sophisticated regulatory systems.
| DISENTANGLING COMPLEX NETWORKS |
|---|
|
|
|---|
As molecular data on regulatory components and interactions are quickly accumulating, biologists are facing increasingly complex networks, which are defying not only our intuition, but also standard modelling approaches.
In this context, an important challenge lies in the development of strategies to decompose large networks into functional cross-regulatory modules, thereby enabling stepwise dynamical modelling and analysis, while managing procedures to reassemble sub-models into a comprehensive model and to efficiently compute collective properties.
In this respect, a first instructive step would be to systematically identify the feedback circuits at the origin of remarkable dynamical properties of complex networks (functional circuits). Several groups are currently addressing this issue, either in terms of novel algorithmic developments, or through the formulation of mathematical theorems (sufficient conditions for circuit functionality). Preliminary results in the discrete framework suggest that only a relatively small fraction of (often short) regulatory circuits are functional in large networks, but that intertwined circuits might collectively contribute to generate specific properties (e.g. multistability, cf. [27, 28]).
This last point leads to the question of whether other, higher level feedback structures should be systematically defined, which could be endowed with more sophisticated dynamical properties. Indeed, combinations of feedback circuits have already been evocated in the context of steady-state linear stability analyses for biological systems [29]. More recently, working on abstract dynamical systems, Thomas and Kaufman [12] have introduced the concept of nucleus, defined as a combination of disjoint circuits, which together involve all the components of a regulatory system. In differential systems, these nuclei appear to play a crucial role in the specification of steady-state properties. However, the significance of these nuclei and their potential practical application to systems biology remain to be properly assessed.
Finally, as our knowledge progresses and diversifies, it becomes possible to compare the design of regulatory networks controlling similar dynamical processes in different organisms, or yet similar processes at different time or place in the same organism. Such comparative approach is already underway in the case of networks controlling cell cycle [30]. Tentatively, comparative model analyses focusing on network design in relation with dynamical behaviour should ultimately contribute to our understanding of the essential molecular processes at the basis of biological development and evolution.
Key Points
|
Funding
Supported by the European Commission (contract LSHG-CT-2004-512143) and the French Ministry of Research (ACI IMPbio and ANR BioSys).
| Acknowledgements |
|---|
|
|
|---|
I also wish to thank Claudine Chaouiya for her critical reading of a draft of this article.
| FOOTNOTES |
|---|
|
|
|---|
1This article is dedicated to René Thomas, for his inspiring pioneering contributions to Systems Biology
Denis Thieffry, with a background in Molecular Biology, is a Professor in Computational Biology at the Université de la Méditerranée and project leader at the INSERM ERM 206 research unit. He is widely recognised for his contributions to the field of qualitative dynamical modelling of genetic regulatory networks.
Submitted: March 31, 2007. Received (in revised form): June 4, 2007.
| References |
|---|
|
|
|---|
- Umbarger HE. Evidence for a negative-feedback mechanism in the biosynthesis of isoleucine. Science (1956) 123:848.
[Free Full Text] - Novick A, Weiner M. Enzyme induction as an all or none phenomenon. Proc Natl Acad Sci USA (1957) 43:55366.
[Free Full Text] - Monod J, Jacob F. Teleonomic mechanisms in cellular metabolism, growth, and differentiation. Cold Spring Harb Symp Quant Biol (1961) 26:389401.
[Abstract/Free Full Text] - Goldberger RF. Autogenous regulation of gene expression. Science (1974) 183:8106.
[Abstract/Free Full Text] - Meinhardt H, Gierer A. Applications of a theory of biological pattern formation based on lateral inhibition. J Cell Sci (1974) 15:32146.
[Abstract/Free Full Text] - Savageau MA. Design of molecular control mechanisms and the demand for gene expression. Proc Natl Acad Sci USA (1977) 74:564751.
[Abstract/Free Full Text] - Thomas R. On the relation between the logical structure of systems and their ability to generate multiple steady states of sustained oscillations. Springer Series Synergetics (1981) 9:18093.
- Thomas R, DAri R. Biological Feedback (1990) Boca Raton, Forida: CRC Press.
- Remy R, Ruet, Thieffry D. Positive or negative regulatory circuit inference from multilevel dynamics. Lect Notes Control Inf Sci (2006) 341:26370.
- de Jong H, Ropers D. Strategies for dealing with incomplete information in the modeling of molecular interaction networks. Brief Bioinform (2006) 7:35463.
[Abstract/Free Full Text] - Tyson J. Classification of instabilities in chemical reactions systems. Chem Phys (1975) 62:101015.[CrossRef]
- Thomas R, Kaufman M. Frontier diagrams: partition of the phase space according to the signs of the eigenvalues or to the sign patterns of the circuits. Int J Bifurcation Chaos (2005) 15:124.[CrossRef][Web of Science]
- Thomas R. Positive feedback circuits are a necessary conditions for positive eigenvalues of the Jacobian matrix. Ber Bungenses Phys Chem (1994) 98:114851.
- Plahte E, Mestl T, Omholt SW. Feedback loops, stability and multistationarity in dynamical systems. J Biol Syst (1995) 3:40913.[CrossRef]
- Snoussi EH. Necessary conditions for multistationarity and stable periodicity. J Biol Syst (1998) 6:39.[CrossRef]
- Gouzé J-L. Positive and negative circuits in dynamical systems. J Biol Syst (1998) 6:1115.[CrossRef]
- Soulé C. Graphic requirements for multistationarity. ComPlexUs (2003) 1:12333.[CrossRef]
- Richard A. Modèle formel pour les réseaux de régulation génétique et influence des circuits de rétroaction. (2006) France: University of Evry. Thesis.
- Goldbeter A. Computational approaches to cellular rhythms. Nature (2002) 420:23845.[CrossRef][Medline]
- Tyson JJ, Chen K, Novak B. Network dynamics and cel physiology. Nat Rev Mol Cell Biol (2001) 2:90816.[CrossRef][Web of Science][Medline]
- Remy R Mossé B, Chaouiya C, Thieffry D. Discrete dynamics of regulatory feedback circuits. Bioinformatics (2003) 10:ii1728.
- Gardner TS, Cantor CR, Collins JJ. Construction of a genetic toggle switch in Escherichia coli. Nature (2000) 403:33942.[CrossRef][Medline]
- Elowitz MB, Leibler S. A synthetic oscillatory network of transcriptional regulators. Nature (2000) 403:3358.[CrossRef][Medline]
- Becskei A, Serrano L. Engineering stability in gene networks by autoregulation. Nature (2000) 405:5903.[CrossRef][Medline]
- Becskei A, Seraphin B, Serrano L. Positive feedback in eukaryotic gene networks: cell differentiation by graded to binary response conversion. EMBO J (2001) 20:252835.[CrossRef][Web of Science][Medline]
- Wall ME, Hlavacek WS, Savageau MA. Design of gene circuits: lessons from bacteria. Nat Rev Genet (2004) 5:3442.[CrossRef][Web of Science][Medline]
- Cinquin O, Demongeot J. High-dimensional switches and the modelling of cellular differentiation. J Theor Biol (2005) 233:391411.[CrossRef][Web of Science][Medline]
- Sánchez L, Thieffry D. Segmenting the fly embryo: a logical analysis of the pair-rule cross-regulatory module. J theor Biol (2003) 224:51737.[CrossRef][Web of Science][Medline]
- Eisenfeld J, DeLisi C. On conditions for qualitative instability of regulatory circuits with application to immunological loops. In: Mathematics and Computers in Medical ApplicationsEisenfeld J, DeLisi C, eds. (1985) Elsevier Science Publishers. 3953.
- Csikasz-Nagy A, Battogtokh D, Chen KC, et al. Analysis of a generic model of eukaryotic cell-cycle regulation. Biophys J (2006) 90:436179.[CrossRef][Web of Science][Medline]
This article has been cited by other articles:
![]() |
A. Csikasz-Nagy Computational systems biology of the cell cycle Brief Bioinform, July 1, 2009; 10(4): 424 - 434. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. Remy and P. Ruet From minimal signed circuits to the dynamics of Boolean regulatory networks Bioinformatics, August 15, 2008; 24(16): i220 - i226. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||



ij and
i denote (constant and strictly positive) kinetic parameters. According to these equations, the first gene (whose product level is given by x1) activates the expression of the second gene (x2), which in turn inhibits the third gene (x3); finally, genes 2 and 3 both activate the first gene independently (sum). All regulatory products are further linearly degraded (i.e. proportionally to product concentration). In the bottom panel, the Jacobian matrix of this ODE system is displayed. Each term aij of this matrix gives the partial derivative of the evolution term i with respect to the variable j. As all these terms involve only positive parameters and variable squares, they have fixed signs. The system encompasses five circuits indicated by dotted circles and arrows on the matrix: three one-element negative circuits, one two-element positive circuit involving genes 1 and 2, and one three-element negative circuit.


