Active regulation in gene networks poses mathematical challenges that have led to conflicting approaches to analysis. Competing regulation that keeps concentrations of some transcription factors at or near threshold values leads to so-called singular dynamics when steeply sigmoidal interactions are approximated by step functions. An extension, due to Artstein and coauthors, of the classical singular perturbation approach was suggested as an appropriate way to handle the complex situation where non-trivial dynamics, such as a limit cycle, of fast variables occur in switching domains. This non-trivial behaviour can occur when a gene regulates multiple other genes at the same threshold. Here, it is shown that it is possible for nonuniqueness to arise in such a system in the case of limiting step-function interactions. This nonuniqueness is reminiscent of but not identical to the nonuniqueness of Filippov solutions. More realistic gene network models have sigmoidal interactions, however, and in the example considered here, it is shown numerically that the corresponding phenomenon in smooth systems is a sensitivity to initial conditions that leads in the limit to densely interwoven basins of attraction of different fixed point attractors.

1.
Z.
Artstein
, “
On singularly perturbed ordinary differential equations with measure-valued limits
,”
Math. Bohemica
127
(
2
),
139
152
(
2002
).
2.
Z.
Artstein
,
J.
Linshiz
, and
E. S.
Titi
, “
Young measure approach to computing slowly advancing fast oscillations
,”
Multiscale Model. Simul.
6
,
1085
1097
(
2007
).
3.
Z.
Artstein
and
A.
Vigodner
, “
Singularly perturbed ordinary differential equations with dynamic limits
,”
Proc. R. Soc. Edinburgh
126
,
541
569
(
1996
).
4.
R.
Casey
,
H.
de Jong
, and
J.-L.
Gouzé
, “
Piecewise-linear models of genetic regulatory networks: Equilibria and their stability
,”
J. Math. Biol.
52
(
1
),
27
56
(
2006
).
5.
H.
de Jong
,
J.
Geiselmann
,
C.
Hernandez
, and
M.
Page
, “
Genetic network analyzer: Qualitative simulation of genetic regulatory networks
,”
Bioinformatics
19
(
3
),
336
344
(
2003
).
6.
H.
de Jong
,
J.-L.
Gouzé
,
C.
Hernandez
,
M.
Page
,
T.
Sari
, and
J.
Geiselmann
, “
Qualitative simulations of genetic regulatory networks using piecewise linear models
,”
Bull. Math. Biol.
66
(
2
),
301
340
(
2004
).
7.
R.
Edwards
, “
Analysis of continuous-time switching networks
,”
Physica D
146
(
1
),
165
199
(
2000
).
8.
R.
Edwards
and
L.
Glass
, “
A calculus for relating the dynamics and structure of complex biological networks
,” in
Adventures in Chemical Physics
, edited by
R. S.
Berry
and
J.
Jortner
(
Wiley
,
Hoboken, New Jersey
,
2006
), Vol.
132
, pp.
151
178
.
9.
E.
Farcot
and
J.-L.
Gouzé
, “
Limit cycles in piecewise-affine gene network models with multiple interaction loops
,”
Int. J. Syst. Sci.
41
,
119
130
(
2010
).
10.
A. F.
Filippov
,
Differential Equations with Discontinuous Right-Hand Sides
(
Nauka
,
Moscow
,
1985
) (Russian); English transl. (Kluwer, Dordrecht, 1998).
11.
L.
Glass
, “
Global analysis of nonlinear chemical kinetics
,”
Statistical Mechanics, Part B
(
Plenum Press
,
New York
,
1977
), pp.
311
349
.
12.
L.
Glass
and
S.
Kauffman
, “
The logical analysis of continuous non-linear biochemical control networks
,”
J. Theor. Biol.
39
(
1
),
103
129
(
1973
).
13.
J.-L.
Gouzé
and
T.
Sari
, “
A class of piecewise linear differential equations arising in biological models
,”
Dyn. Syst.
17
(
4
),
299
316
(
2002
).
14.
L.
Ironi
,
L.
Panzeri
,
E.
Plahte
, and
V.
Simoncini
, “
Dynamics of actively regulated gene networks
,”
Physica D
240
(
8
),
779
794
(
2011
).
15.
A.
Machina
,
R.
Edwards
, and
P.
van den Driessche
, “
Singular dynamics in gene network models
,”
SIAM J. Appl. Dyn. Syst.
12
(
1
),
95
125
(
2013
).
16.
A.
Machina
and
A.
Ponosov
, “
Stability of stationary solutions of piecewise affine differential equations describing gene regulatory networks
,”
J. Math. Anal. Appl.
380
(
2
),
736
749
(
2011
).
17.
A.
Machina
and
A.
Ponosov
, “
Filippov solutions in the analysis of piecewise linear models describing gene regulatory networks
,”
Nonlinear Anal.
74
,
882
900
(
2011
).
18.
J.
Monod
and
F.
Jacob
, “
Teleonomic mechanisms in cellular metabolism, growth, and differentiation
,”
Cold Spring Harbor Symp. Quant. Biol.
26
,
389
401
(
1961
).
19.
E.
Plahte
and
S.
Kjøglum
, “
Analysis and generic properties of gene regulatory networks with graded response functions
,”
Physica D
201
,
150
176
(
2005
).
20.
E.
Plahte
,
T.
Mestl
, and
S. W.
Omholt
, “
A methodological basis for the description and analysis of systems with complex switch-like interactions
,”
J. Math. Biol.
36
,
321
348
(
1998
).
21.
S.
Veflingstad
and
E.
Plahte
, “
Analysis of gene regulatory network models with graded and binary transcriptional responses
,”
BioSystems
90
,
323
339
(
2007
).
You do not currently have access to this content.