Synchronism vs asynchronism in Boolean automata networks Sylvain - - PowerPoint PPT Presentation

synchronism vs asynchronism in boolean automata networks
SMART_READER_LITE
LIVE PREVIEW

Synchronism vs asynchronism in Boolean automata networks Sylvain - - PowerPoint PPT Presentation

Synchronism vs asynchronism in Boolean automata networks Sylvain Sen MOVE seminar 18th January 2018 Outline Introduction 1 Main definitions 2 3 Deterministic periodic updates 4 Non-deterministic updates Sylvain Sen Synchronism vs


slide-1
SLIDE 1

Synchronism vs asynchronism in Boolean automata networks

Sylvain Sené

MOVE seminar 18th January 2018

slide-2
SLIDE 2

Outline

1

Introduction

2

Main definitions

3

Deterministic periodic updates

4

Non-deterministic updates

Sylvain Sené Synchronism vs asynchronism in BANs 2/27

slide-3
SLIDE 3

Introduction

Outline

1

Introduction

2

Main definitions

3

Deterministic periodic updates

4

Non-deterministic updates

Sylvain Sené Synchronism vs asynchronism in BANs 3/27

slide-4
SLIDE 4

Introduction

BANs, non formally

§ A discrete computational model of interaction systems. § From a theoretical standpoint:

§ Simple setting and representation. § Able to capture dynamically a lot of behavioural intricacies and

heterogeneities.

§ From a more practical/applied standpoint:

§ Originate from neural theoretical modelling (McCulloch, Pitts, 1943). § Developed in the context of genetics (Kauffman, 1969; Thomas, 1973). § The most used mathematical objects for genetic regulation qualitative

modelling.

Sylvain Sené Synchronism vs asynchronism in BANs 4/27

slide-5
SLIDE 5

Introduction

The (a-)synchronicity problematic(s)

§ The causality of events along time depends on the relation between

automata updates and “time” but...

§ How to define this relation? § How to study the causal perturbations due to changes of this relation?

§ Mathematical pertinence:

§ Neat problematic at the frontier of dynamical systems, combinatorics,

complexity and computability.

§ Biological pertinence:

§ Genetic expression and chromatin dynamics.

§ A remaining question: does model synchronicity stand for modelled

system simultaneity?

Sylvain Sené Synchronism vs asynchronism in BANs 5/27

slide-6
SLIDE 6

Main definitions

Outline

1

Introduction

2

Main definitions

3

Deterministic periodic updates

4

Non-deterministic updates

Sylvain Sené Synchronism vs asynchronism in BANs 6/27

slide-7
SLIDE 7

Main definitions

BANs and interaction graphs

A Boolean automata network (BAN) of size n is a function f : Bn Ñ Bn x “ px0,x1,...,xn´1q ÞÑ fpxq “ pf0pxq,f1pxq,...,fn´1pxqq , where @i P t0,...,n´1u, xi P B is the state of automaton i, and Bn is the set of configurations. The interaction graph of f is the signed digraph Gpfq : pV,E Ď V ˆVq where:

§ V “ t0,...,n´1u; § pi,jq P E is positive if Dx P Bn s.t.

fjpx0,...,xi´1,0,xi`1,...,xn´1q “ 0 and fjpx0,...,xi´1,1,xi`1,...,xn´1q “ 1;

§ pi,jq P E is negative if Dx P Bn s.t.

fjpx0,...,xi´1,0,xi`1,...,xn´1q “ 1 and fjpx0,...,xi´1,1,xi`1,...,xn´1q “ 0.

Sylvain Sené Synchronism vs asynchronism in BANs 7/27

slide-8
SLIDE 8

Main definitions

BANs and interaction graphs

A Boolean automata network (BAN) of size n is a function f : Bn Ñ Bn x “ px0,x1,...,xn´1q ÞÑ fpxq “ pf0pxq,f1pxq,...,fn´1pxqq , where @i P t0,...,n´1u, xi P B is the state of automaton i, and Bn is the set of configurations. f : B4 Ñ B4 f “ $ ’ ’ ’ & ’ ’ ’ % f0pxq “ x0 _x1 ^x3 f1pxq “ x0 ^px1 _x2q f2pxq “ x3 f3pxq “ x0 _x1 3 1 2

Sylvain Sené Synchronism vs asynchronism in BANs 7/27

slide-9
SLIDE 9

Main definitions

Automata updates

f0pxq “ x0 _x1 ^x3 f3pxq “ x0 _x1 f1pxq “ x0 ^px1 _x2q f2pxq “ x3

0101 1 2 3

Sylvain Sené Synchronism vs asynchronism in BANs 8/27

slide-10
SLIDE 10

Main definitions

Automata updates

t2u

f3pxq “ x0 _x1 f1pxq “ x0 ^px1 _x2q f2pxq “ x3

0101 1 2 3

f0pxq “ x0 _x1 ^x3

“ 0

Sylvain Sené Synchronism vs asynchronism in BANs 8/27

slide-11
SLIDE 11

Main definitions

Automata updates

0001

t3u t1u

f0pxq “ x0 _x1 ^x3

Asynchronous transitions

f3pxq “ x0 _x1 f1pxq “ x0 ^px1 _x2q f2pxq “ x3

0101 1 2 3 0000 0100 1000 1001 1100 1101

t2u t0u

Sylvain Sené Synchronism vs asynchronism in BANs 8/27

slide-12
SLIDE 12

Main definitions

Automata updates

1100 1101

t2u t0u

0001

t3u t1u t0,1u t0,1u

f0pxq “ x0 _x1 ^x3 f3pxq “ x0 _x1 f1pxq “ x0 ^px1 _x2q f2pxq “ x3

0101 1 3 0000 0100 2 1000 1001

“ 0 1 “

Sylvain Sené Synchronism vs asynchronism in BANs 8/27

slide-13
SLIDE 13

Main definitions

Automata updates

t1,2,3u

f3pxq “ x0 _x1 f1pxq “ x0 ^px1 _x2q f2pxq “ x3

0101 1 2 3 0000 0100 1000 1001 1100 1101

t2u t0u

0001

t3u t1u t1,3u t1,2u t2,3u t0,1,3u t0,2u t0,1,2,3u t0,1u t0,1,2u t0,3u t0,2,3u

f0pxq “ x0 _x1 ^x3

Synchronous transitions

Sylvain Sené Synchronism vs asynchronism in BANs 8/27

slide-14
SLIDE 14

Main definitions

Update modes and BAN behaviours

§ An update mode is a way of organising the automata updates along time. § It can be deterministic (periodic or not) or non-deterministic (stochastic

  • r not).

§ There exists an infinite number of update modes.

Sylvain Sené Synchronism vs asynchronism in BANs 9/27

slide-15
SLIDE 15

Main definitions

Update modes and BAN behaviours

§ An update mode is a way of organising the automata updates along time. § It can be deterministic (periodic or not) or non-deterministic (stochastic

  • r not).

§ There exists an infinite number of update modes. § The update mode defines the network behaviour. § The behaviour of a BAN f is described by a transition graph

G˛pfq “ pBn,T Ď Bn ˆpPpVqzHqˆt0,1unq, where ˛ represents a given “fair” update mode.

Sylvain Sené Synchronism vs asynchronism in BANs 9/27

slide-16
SLIDE 16

Main definitions

Some examples

f : B3 Ñ B3 f “ $ & % f0pxq “ x1 _x2 f1pxq “ x0 ^x2 f2pxq “ x2 ^px0 _x1q 1 2 Parallel evolution

000 001 010 011 100 101 110 111

toto

Sylvain Sené Synchronism vs asynchronism in BANs 10/27

slide-17
SLIDE 17

Main definitions

Some examples

f : B3 Ñ B3 f “ $ & % f0pxq “ x1 _x2 f1pxq “ x0 ^x2 f2pxq “ x2 ^px0 _x1q 1 2 Parallel evolution

000 001 010 011 100 101 110 111 § An attractor of pf,˛q is a terminal SCC of

G˛pfq.

§ A fixed point (stable configuration) is a

trivial attractor.

§ A limit cycle (stable oscillation) is a

non-trivial attractor.

Sylvain Sené Synchronism vs asynchronism in BANs 10/27

slide-18
SLIDE 18

Main definitions

Some examples

f : B3 Ñ B3 f “ $ & % f0pxq “ x1 _x2 f1pxq “ x0 ^x2 f2pxq “ x2 ^px0 _x1q 1 2 pt0u,t1u,t2uq-sequential evolution

000 000 000 001 010 011 100 100 101 101 101 110 110 111 111

Sylvain Sené Synchronism vs asynchronism in BANs 10/27

slide-19
SLIDE 19

Main definitions

Some examples

f : B3 Ñ B3 f “ $ & % f0pxq “ x1 _x2 f1pxq “ x0 ^x2 f2pxq “ x2 ^px0 _x1q 1 2 pt0u,t1u,t2uq-sequential evolution

000 001 010 011 100 101 110 111

Sylvain Sené Synchronism vs asynchronism in BANs 10/27

slide-20
SLIDE 20

Main definitions

Some examples

f : B3 Ñ B3 f “ $ & % f0pxq “ x1 _x2 f1pxq “ x0 ^x2 f2pxq “ x2 ^px0 _x1q 1 2 pt0,2u,t1uq-block-sequential evolution

000 000 001 001 010 011 100 100 101 110 110 111 111

Number

  • f
  • rdered

partitions:

Bord

n

n´1

ÿ

k“0

ˆn k ˙ Bord

k

,

with Bord “ 1.

Sylvain Sené Synchronism vs asynchronism in BANs 10/27

slide-21
SLIDE 21

Main definitions

Some examples

f : B3 Ñ B3 f “ $ & % f0pxq “ x1 _x2 f1pxq “ x0 ^x2 f2pxq “ x2 ^px0 _x1q 1 2 pt0,2u,t1uq-block-sequential evolution

000 001 010 011 100 101 110 111

Number

  • f
  • rdered

partitions:

Bord

n

n´1

ÿ

k“0

ˆn k ˙ Bord

k

,

with Bord “ 1.

Sylvain Sené Synchronism vs asynchronism in BANs 10/27

slide-22
SLIDE 22

Main definitions

Some examples

f : B3 Ñ B3 f “ $ & % f0pxq “ x1 _x2 f1pxq “ x0 ^x2 f2pxq “ x2 ^px0 _x1q 1 2 pt0,2u,t1uq-block-sequential evolution

000 001 010 011 100 101 110 111

Number

  • f
  • rdered

partitions:

Bord

n

n´1

ÿ

k“0

ˆn k ˙ Bord

k

,

with Bord “ 1.

Sylvain Sené Synchronism vs asynchronism in BANs 10/27

slide-23
SLIDE 23

Main definitions

Some examples

f : B3 Ñ B3 f “ $ & % f0pxq “ x1 _x2 f1pxq “ x0 ^x2 f2pxq “ x2 ^px0 _x1q 1 2 Asynchronous evolution ` t0,2u-synchronous transitions

000 001 010 011 100 101 110 111

1 2 1 1 1 1 1 1 1 2 2 2 2 Sylvain Sené Synchronism vs asynchronism in BANs 10/27

slide-24
SLIDE 24

Main definitions

Some examples

f : B3 Ñ B3 f “ $ & % f0pxq “ x1 _x2 f1pxq “ x0 ^x2 f2pxq “ x2 ^px0 _x1q 1 2 Asynchronous evolution ` t0,2u-synchronous transitions

000 001 010 011 100 101 110 111

1 2 1 1 1 1 1 1 1 2 2 2 2 Sylvain Sené Synchronism vs asynchronism in BANs 10/27

slide-25
SLIDE 25

Deterministic periodic updates

Outline

1

Introduction

2

Main definitions

3

Deterministic periodic updates

4

Non-deterministic updates

Sylvain Sené Synchronism vs asynchronism in BANs 11/27

slide-26
SLIDE 26

Deterministic periodic updates

Update graphs

Given an interaction graph G “ pV,Eq, a labelled graph is a graph pG,labq, with lab : E Ñ t‘,au. A labelled graph pG,labq is an update graph if there exist s : V Ñ t1,...,nu s.t. @pi,jq P E, labpi,jq “ # ‘ if spiq ě spjq a if spiq ă spjq .

2 1 3 1 2 3 3 2 1 2 3 1

‘ ‘ ‘ a ‘ ‘ ‘ ‘ a a ‘ a ‘ ‘ ‘ ‘ ‘ ‘

pt0,1,2,3uq pt0u,t1u,t2u,t3uq pt2,3u,t0,1uq

Sylvain Sené Synchronism vs asynchronism in BANs 12/27

slide-27
SLIDE 27

Deterministic periodic updates

Update graphs and dynamics

Let f be a BAN and Gpfq “ pV,Eq its interaction graph, let π be the parallel update mode, and let s ‰ s1 be two distinct block-sequential modes different from π. Theorem 1 (Aracena et al., 2009) If Gpf,labsq “ Gpf,labs1q then Gspfq “ Gs1pfq. Theorem 2 (Tchuente, 1988; Aracena et al., 2009) If s is defined as @j P t0,...,n´1u,@i s.t. pi,jq P E, spiq ě spjq then Gspfq “ Gπpfq. Theorem 3 (Aracena et al., 2009) Consider s and f s.t. all the loops in Gpfq are positive. Then there exists s1 such that Gspfq and Gs1pfq do not have any common limit cycle.

Sylvain Sené Synchronism vs asynchronism in BANs 13/27

slide-28
SLIDE 28

Deterministic periodic updates

Update graphs and dynamics

Theorem 1 (Aracena et al., 2009) If Gpf,labsq “ Gpf,labs1q then Gspfq “ Gs1pfq. s1 ” pt1u,t0u,t2u,t3uq s2 ” pt1u,t2u,t0u,t3uq s3 ” pt1u,t2u,t0,3uq 1 2 3 f “ $ ’ ’ ’ & ’ ’ ’ % f0pxq “ x1 ^x3 f1pxq “ x0 f2pxq “ x1 _x2 f3pxq “ x2 ^x3

Sylvain Sené Synchronism vs asynchronism in BANs 13/27

slide-29
SLIDE 29

Deterministic periodic updates

Update graphs and dynamics

Theorem 1 (Aracena et al., 2009) If Gpf,labsq “ Gpf,labs1q then Gspfq “ Gs1pfq.

‘ a ‘ a

s2 ” pt1u,t2u,t0u,t3uq s3 ” pt1u,t2u,t0,3uq 1 2 f “ $ ’ ’ ’ & ’ ’ ’ % f0pxq “ x1 ^x3 f1pxq “ x0 f2pxq “ x1 _x2 f3pxq “ x2 ^x3 3 1 2 3 s1 ” pt1u,t0u,t2u,t3uq

‘ ‘ a

1 4 5 8 10 12 6 14 2 7 3 9 11 13 15

Sylvain Sené Synchronism vs asynchronism in BANs 13/27

slide-30
SLIDE 30

Deterministic periodic updates

Interaction cycles

2 types of interaction cycles, the positive and the negative ones: an even number of negative arcs

5

C`

6 1 2 3 4 5

6 1 2 3 4

an odd number of negative arcs Seminal results: Theorem 4 (Robert, 1986) If Gpfq is acyclic, then f admits a unique attractor which is a fixed point. Theorem 5 (Thomas, 1981; Richard, Comet, 2007) If there are no positive cycles in Gpfq, f admits no more than one fixed point.

Sylvain Sené Synchronism vs asynchronism in BANs 14/27

slide-31
SLIDE 31

Deterministic periodic updates

Impact of update modes on cycles

sp1q “ 1 x1pt`1q “ f1px0ptqq x2pt`1q “ f2px1pt`1qq “ f2pf1px0ptqqq sp0q “ 1 Block-sequential mode s ” pt0,1u,t2uq 2 1 sp2q “ 2 x0pt`1q “ f0px2ptqq

Sylvain Sené Synchronism vs asynchronism in BANs 15/27

slide-32
SLIDE 32

Deterministic periodic updates

Impact of update modes on cycles

sp1q “ 1 x1pt`1q “ f1px0ptqq x2pt`1q “ f2px1pt`1qq “ f2pf1px0ptqqq sp0q “ 1 Block-sequential mode s ” pt0,1u,t2uq 2 1 sp2q “ 2 x0pt`1q “ f0px2ptqq

Sylvain Sené Synchronism vs asynchronism in BANs 15/27

slide-33
SLIDE 33

Deterministic periodic updates

Impact of update modes on cycles

sp1q “ 1 x1pt`1q “ f1px0ptqq x2pt`1q “ f2px1pt`1qq “ f2pf1px0ptqqq sp0q “ 1 Block-sequential mode s ” pt0,1u,t2uq 2 1 sp2q “ 2 x0pt`1q “ f0px2ptqq

Sylvain Sené Synchronism vs asynchronism in BANs 15/27

slide-34
SLIDE 34

Deterministic periodic updates

Impact of update modes on cycles

sp1q “ 1 x1pt`1q “ f1px0ptqq x2pt`1q “ f2px1pt`1qq “ f2pf1px0ptqqq Block-sequential mode sp0q “ 1 s ” pt0,1u,t2uq 2 1 sp2q “ 2 x0pt`1q “ f0px2ptqq 2 1 Interaction graph Gpf,sq “ pV,Epsqq

Each arc pi,jq P Epsq represents the dependence of xjpt`1q on xiptq.

Sylvain Sené Synchronism vs asynchronism in BANs 15/27

slide-35
SLIDE 35

Deterministic periodic updates

Impact of update modes on cycles

invpsq

“ tpi,i`1q | spiq ă spi`1qu

Sylvain Sené Synchronism vs asynchronism in BANs 15/27

slide-36
SLIDE 36

Deterministic periodic updates

Impact of update modes on cycles

invpsq

“ tpi,i`1q | spiq ă spi`1qu

Sylvain Sené Synchronism vs asynchronism in BANs 15/27

slide-37
SLIDE 37

Deterministic periodic updates

Impact of update modes on cycles

invpsq

“ tpi,i`1q | spiq ă spi`1qu

Theorems (Goles, Noual, 2010) ⊲ The dynamics induced by two update modes s and s1 are equal iff invpsq “ invps1q.

Given a cycle of size n, the total number of

distinct dynamics induced by block- sequential update modes is:

n´1

ÿ

k“0

ˆn k ˙ “ 2n ´1. ⊲ invpsq ‰ invps1q ù ñ no common limit cycles. ⊲ Iterating a cycle of size n with an update mode s with |invpsq| “ k corresponds to iterating a cycle of same sign and of size n´k in parallel.

Sylvain Sené Synchronism vs asynchronism in BANs 15/27

slide-38
SLIDE 38

Deterministic periodic updates

Impact of update modes on cycles

Theorem 6 (Goles, Noual, 2010) invpsq ‰ invps1q ù ñ no common limit cycles.

Proof First, let us note that @i,j P V, frj,is : # fj ˝fj´1 ˝¨¨¨˝fi if i ď j fj ˝fj´1 ˝¨¨¨˝f0 ˝fn´1 ˝¨¨¨˝fi if i ą j. Suppose that pi,i`1q P invpsqzinvps1q and that Dx “ xsptq “ xs1ptq s.t. xspt`1q “ xs1pt`1q. Then: xs

i`1pt`2q “ fi`1pxs ipt`2qq “ fri`1,i‹ `1spxs i‹pt`1qq,

and xs1

i`1pt`2q “ fi`1pxs1 i pt`1qq “ fi`1pxs ipt`1qq “ fri`1,i‹ `1spxs i‹ptqq,

where i‹ “ maxptk ă i | spkq ě spk`1quq. By the injectivity of fri`1,i‹ `1s, if xspt`2q “ xs1pt`2q then xi‹pt`1q “ xi‹ptq. Now, if x belongs to an attractor that is induced identically by both s and s1, then xsptq “ xs1ptq @t. As result, in this case, @t, xs

i‹pt`1q “ xs1 i‹ptq “ xs i‹ptq. In other terms, the state of

node i‹ is fixed in the attractor. Hence the states of all nodes are fixed in the attractor which therefore is a fixed point.

Sylvain Sené Synchronism vs asynchronism in BANs 15/27

slide-39
SLIDE 39

Deterministic periodic updates

Update graphs other related results

Q: Is a labelled graph an update graph?

3 2 3 1 2 3 1 Reversed labelled graph 4 pG,labq‘

R

Labelled graph pG,labq Reduced labelled graph pG,labq‘ 0,4 0,4 2 1

Sylvain Sené Synchronism vs asynchronism in BANs 16/27

slide-40
SLIDE 40

Deterministic periodic updates

Update graphs other related results

Q: Is a labelled graph an update graph?

3 Reversed labelled graph 2 pG,labq‘

R

3 1 2 3 1 4 Labelled graph pG,labq Reduced labelled graph pG,labq‘ 0,4 0,4 2 1

Theorem 7 (Aracena et al., 2011) A labelled digraph pG,labq is an update graph iff pG,labq‘

R does not contain

any forbidden cycle.

Idea sp1q ă sp2q sp1q ě psp0q “ sp4qq ^ sp3q ă psp0q “ sp4qq sp3q ě sp2q ù ñ sp1q ă sp3q sp3q ă sp1q

Sylvain Sené Synchronism vs asynchronism in BANs 16/27

slide-41
SLIDE 41

Deterministic periodic updates

Update graphs other related results

Q: How to find the most compact update mode on pG,labq?

3 2 3 1 2 pG1,labq‘

R

3

s ” pt0,4u,t1,3u,t2uq

1 4 0,4 0,4 pG1,labq pG1,labq‘ 2 1

Sylvain Sené Synchronism vs asynchronism in BANs 16/27

slide-42
SLIDE 42

Deterministic periodic updates

Update graphs other related results

Q: How to find the most compact update mode on pG,labq?

3 2 3 1 2 pG1,labq‘ 3 pG1,labq‘

R

1 4 0,4 0,4

s ” pt0,4u,t1,3u,t2uq

pG1,labq 2 1 Algorithm

  • Init. Take G1 :“ pG,labq‘

R and t :“ 1.

(1) Compute the paths Pa “ tP | #pa P Pq is max.u on G1. If Pa “ H, goto (4). (2) The targets T of the last negative arc of each P of Pa, and their successors SpTq are scheduled at time step t. t :“ t`1. (3) Remove T, SpTq and all their incoming arcs from G1, and go back to (1). (4) All the remaining nodes are scheduled all at once, at time step t.

Sylvain Sené Synchronism vs asynchronism in BANs 16/27

slide-43
SLIDE 43

Non-deterministic updates

Outline

1

Introduction

2

Main definitions

3

Deterministic periodic updates

4

Non-deterministic updates

Sylvain Sené Synchronism vs asynchronism in BANs 17/27

slide-44
SLIDE 44

Non-deterministic updates

Basic definitions and notations

@x “ px0,...,xn´1q P Bn,@i P V, xi “ px0,...,xi´1,xi,xi`1,...,xn´1q @x P Bn,@W “ W1 Ztiu Ď V, xW “ pxiq

W1

“ pxW1q

i

The sign of an influence of i on j in x is signxpi,jq “ fjpxq´fjpxiq xi ´xi

i

“ spxiq¨pfjpxq´fjpxiqq, where s : b P B ÞÑ b´b P t´1,1u. Given x,y P Bn, Dpx,yq “ ti P V | xi ‰ yiu and dpx,yq “ |Dpx,yq|. Epxq “ tpi,jq P V ˆV | signxpi,jq ‰ 0u represents the set of effective influences

  • f Gpfq in x, which formally means that

@i,j P V,Dx P Bn, fjpxq ‰ fjpxiq ð ñ pi,jq P E.

Sylvain Sené Synchronism vs asynchronism in BANs 18/27

slide-45
SLIDE 45

Non-deterministic updates

Monotonicity, unstabilities and frustrations

A local function fi is locally monotonic in j if either: @x, fipx0,...,xj´1,0,xj`1,...,xn´1q ď fipx0,...,xj´1,1,xj`1,...,xn´1q

  • r:

@x, fipx0,...,xj´1,0,xj`1,...,xn´1q ě fipx0,...,xj´1,1,xj`1,...,xn´1q . g “ $ ’ ’ ’ & ’ ’ ’ % g0pxq “ x1 ^x3 g1pxq “ x0 g2pxq “ x1 ‘ x2 g3pxq “ x2 _x3 1 2 3 1 2 3 is monotonic. f “ $ ’ ’ ’ & ’ ’ ’ % f0pxq “ x1 ^x3 f1pxq “ x0 f2pxq “ x1 _x2 f3pxq “ x2 _x3 is not.

Sylvain Sené Synchronism vs asynchronism in BANs 19/27

slide-46
SLIDE 46

Non-deterministic updates

Monotonicity, unstabilities and frustrations

A local function fi is locally monotonic in j if either: @x, fipx0,...,xj´1,0,xj`1,...,xn´1q ď fipx0,...,xj´1,1,xj`1,...,xn´1q

  • r:

@x, fipx0,...,xj´1,0,xj`1,...,xn´1q ě fipx0,...,xj´1,1,xj`1,...,xn´1q . An automaton i P V is unstable (resp. stable) in x P Bn if it belongs to the set

Upxq “ ti P V | fipxq ‰ xiu

(resp. Upxq “ VzUpxq). 1 f “ # f0pxq “ x1 f1pxq “ x0

x f0pxq f1pxq

Upxq

p0,0q 1 t0u p0,1q t1u p1,0q 1 1 t1u p1,1q 1 t0u

Sylvain Sené Synchronism vs asynchronism in BANs 19/27

slide-47
SLIDE 47

Non-deterministic updates

Monotonicity, unstabilities and frustrations

A local function fi is locally monotonic in j if either: @x, fipx0,...,xj´1,0,xj`1,...,xn´1q ď fipx0,...,xj´1,1,xj`1,...,xn´1q

  • r:

@x, fipx0,...,xj´1,0,xj`1,...,xn´1q ě fipx0,...,xj´1,1,xj`1,...,xn´1q . An automaton i P V is unstable (resp. stable) in x P Bn if it belongs to the set

Upxq “ ti P V | fipxq ‰ xiu

(resp. Upxq “ VzUpxq). An influence pi,jq P E is frustrated in x iff it belongs to

FRUSpxq “ tpi,jq P E | spxiq¨ spxjq “ ´signpi,jqu. 2 1 f “ $ ’ & ’ % f0pxq “ x2 f1pxq “ x0 _x1 f2pxq “ x0 ^x1

FRUSp000q “ tp0,2qu FRUSp001q “ tp1,2q,p2,0qu FRUSp010q “ tp0,1q,p0,2q,p1,2qu FRUSp011q “ tp0,1q,p2,0qu

Sylvain Sené Synchronism vs asynchronism in BANs 19/27

slide-48
SLIDE 48

Non-deterministic updates

Relations between unstabilities and frustrations

Remark (Noual, S., 2017) If j P Upxq then Di P V´pjq, pi,jq P FRUSpxq. 2 1 f “ $ ’ & ’ % f0pxq “ x0 f1pxq “ x0 _x2 f2pxq “ x1

FRUSp000q “ t (0,0) ,p2,1qu FRUSp001q “ t (0,0) ,p1,2qu FRUSp110q “ t (0,0) ,p1,2qu FRUSp111q “ t (0,0) ,p2,1qu

N.B: The reciprocal does not hold. 2 1 f “ $ ’ & ’ % f0pxq “ x2 f1pxq “ x0 _x1 f2pxq “ x0 ^x1

FRUSp000q “ t (0,2) u FRUSp001q “ tp1,2q,p2,0qu FRUSp010q “ tp0,1q,p0,2q,p1,2qu FRUSp011q “ tp0,1q,p2,0qu

Sylvain Sené Synchronism vs asynchronism in BANs 20/27

slide-49
SLIDE 49

Non-deterministic updates

Relations between unstabilities and frustrations

Lemma 1 (Noual, S., 2017) Adding frustrated influences incoming an unstable automaton cannot stabilise it. Formally, noting V´

FRUSpxqpjq “ V´pjqXti P V | pi,jq P FRUSpxqu, we

have: @x,y P Bn, j P Upxq^ ´ V´

FRUSpxqpjq Ď V´ FRUSpyqpjq

¯ ù ñ j P Upyq.

Proof Input provided by i to j: bj

ipxq “ bpsignpi,jq¨ spxiqq “

# xj if pi,jq R FRUSpxq xj

  • therwise

. By local monotonicity, fjpxq “ ľ

kďm

ckpxq “ ľ

kďm

´ ł

iPVj

k

bj

ipxq “

ł

iPVj

k

pi,jqPFRUSpxq

xj _ ł

iPVj

k

pi,jqRFRUSpxq

xj ¯ , where Vj

k is the set of in-neighbours of j involved in the kth clause.

Let x be unstable, admitting thus at least one frustrated incoming influence. Let y be such that it admits at least one more frustrated incoming influence than x. Since fj can be written as a conjunction of disjunctive clauses, the values of these clauses for y are necessarily the same as for x.

Sylvain Sené Synchronism vs asynchronism in BANs 20/27

slide-50
SLIDE 50

Non-deterministic updates

Critical cycles

Let f be a BAN, G “ pV,Eq its interaction graph, and x a configuration in Bn. A cycle C “ pVC,ECq of G is x-critical if EC Ď FRUSpxq. A cycle C is critical if it is x-critical for some x. Proposition 1 (Noual, S., 2017) A critical cycle is a NOPE-cycle, i.e. negative of odd length or positive of even length.

Proof Let x P Bn. By definition of frustrated influences, if C “ pVC,ECq is x-critical, has length ℓ and sign s then: ź

pi,jqPEC

´signpi,jq “ p´1qℓ ˆ s “ ź

pi,jqPEC

spxiq¨ spxjq “ 1.

1

x f0pxq f1pxq

FRUSpxq

p0,0q tp0,1q,p1,0qu p0,1q 1 H p1,0q 1 H p1,1q tp0,1q,p1,0qu

f “ # f0pxq “ x0 ^x1 f1pxq “ x0 ^x1

Sylvain Sené Synchronism vs asynchronism in BANs 21/27

slide-51
SLIDE 51

Non-deterministic updates

Transitions and trajectories

Name Notation Definition Asynchronous x y dpx,yq ď 1 Synchronous x y dpx,yq ą 1 Elementary x y x y P tx yuYtx yu Non-sequentialisable x y x y not decomposable into smaller elementary transitions

For all x,y P Bn s.t. x ‰ y, x is willing (resp. unwilling) towards y if Dpx,yq Ď Upxq (resp. Dpx,yqX Upxq “ H). A trajectory from x to y is a path x ... y in the transition graph. Let x “ xp0q xp1q ... xpm´1q y “ xpmq be a trajectory from x to y. If @t ă m, Dpxpt`1q,yq Ĺ Dpxptq,yq, this trajectory is direct. It performs no reversed changes, i.e. @t ă m, xptqi “ yi ù ñ @t ă t1 ď m, xpt1qi “ yi.

Sylvain Sené Synchronism vs asynchronism in BANs 22/27

slide-52
SLIDE 52

Non-deterministic updates

Results relating trajectories and critical cycles

Proposition 2 (Noual, S., 2017) Let x a willing configuration towards y.

  • 1. If there are no asynchronous trajectories from x to y, then Dpx,yq induces

a NOPE-cycle that is x-critical.

  • 2. If Dpx,yq does not induce an x-critical cycle, then there is a direct

asynchronous trajectory from x to y.

x “ 00 01 10 y “ 11

1 1 1 1

1 f “ # f0pxq “ x0 _x1 f1pxq “ x0 _x1

x f0pxq f1pxq

Upxq

p0,0q 1 1 Dpx,yq p0,1q 1 H p1,0q 1 H p1,1q 1 1 H

Sylvain Sené Synchronism vs asynchronism in BANs 23/27

slide-53
SLIDE 53

Non-deterministic updates

Results relating trajectories and critical cycles

Proposition 2 (Noual, S., 2017) Let x a willing configuration towards y.

  • 1. If there are no asynchronous trajectories from x to y, then Dpx,yq induces

a NOPE-cycle that is x-critical.

  • 2. If Dpx,yq does not induce an x-critical cycle, then there is a direct

asynchronous trajectory from x to y. 1

x f0pxq f1pxq

Upxq

p0,0q 1 1 t0,1u p0,1q t1u p1,0q 1 1 t1u p1,1q 1 1 H

f “ # f0pxq “ x0 _x1 f1pxq “ x0 _x1

x “ 00 01 10 y “ 11

1 1 1 1 1 1 Sylvain Sené Synchronism vs asynchronism in BANs 23/27

slide-54
SLIDE 54

Non-deterministic updates

Results relating trajectories and critical cycles

Proposition 2 (Noual, S., 2017) Let x a willing configuration towards y.

  • 1. If there are no asynchronous trajectories from x to y, then Dpx,yq induces

a NOPE-cycle that is x-critical.

  • 2. If Dpx,yq does not induce an x-critical cycle, then there is a direct

asynchronous trajectory from x to y. Implication When m local changes are possible in x, then, unless there is a NOPE-cycle of size m, these m changes can be made asynchronously without risking a deadlock, i.e. a situation in which some transitions would have transformed x into a configuration xptq from which y is not reachable anymore.

Sylvain Sené Synchronism vs asynchronism in BANs 23/27

slide-55
SLIDE 55

Non-deterministic updates

Results relating trajectories and critical cycles

Proposition 2 (Noual, S., 2017) Let x a willing configuration towards y.

  • 1. If there are no asynchronous trajectories from x to y, then Dpx,yq induces

a NOPE-cycle that is x-critical.

  • 2. If Dpx,yq does not induce an x-critical cycle, then there is a direct

asynchronous trajectory from x to y. Corollary 1 (Noual, S., 2017) If x y exists, then Dpx,yq induces a NOPE-cycle which is x-critical. Implication In a BAN with no NOPE-cycles of size smaller or equal than m P N, any synchronous change affecting no more than m automata states can be totally sequentialised.

Sylvain Sené Synchronism vs asynchronism in BANs 23/27

slide-56
SLIDE 56

Non-deterministic updates

Structural sensitivity: impact of synchronism

Class N Class F

“null” sensitivity “weak” sensitivity

1 2 3 4 5 6 7 8 9 1 2 3 4 5 6 7 8 9 Class G Class D

“medium” sensitivity “strong” sensitivity

1 2 3 4 5 6 7 8 9 1 2 3 4 5 6 7 8 9

Sylvain Sené Synchronism vs asynchronism in BANs 24/27

slide-57
SLIDE 57

Non-deterministic updates

Structural sensitivity: main result

Theorem 8 (Noual, S., 2017) 1) Synchronism-sensitivity requires the existence of a NOPE-cycle. 2) Significant sensitivity requires the existence of a NOPE-cycle of length strictly smaller than the BAN size as well as of a negative cycle. 3) In the absence of a Hamiltonian NOPE-cycle and positive loops on all automata, little sensitivity also requires a NOPE-cycle of length strictly smaller than the BAN size. A monotonic BAN belonging to sensitivity class D: 2 3 1 f “ $ ’ ’ ’ & ’ ’ ’ % f0pxq “ x2 _px0 ^x1q f1pxq “ x3 _px0 ^x1q f2pxq “ x0 ^x1 f3pxq “ x0 ^x1

Sylvain Sené Synchronism vs asynchronism in BANs 25/27

slide-58
SLIDE 58

Non-deterministic updates

Structural sensitivity: main result

Theorem 8 (Noual, S., 2017) 1) Synchronism-sensitivity requires the existence of a NOPE-cycle. 2) Significant sensitivity requires the existence of a NOPE-cycle of length strictly smaller than the BAN size as well as of a negative cycle. 3) In the absence of a Hamiltonian NOPE-cycle and positive loops on all automata, little sensitivity also requires a NOPE-cycle of length strictly smaller than the BAN size. A monotonic BAN belonging to sensitivity class D: tx P B4 | x0 _x1 “ 1u

0011 0000 0001 0010

asynchronous limit cycle fixed point

Sylvain Sené Synchronism vs asynchronism in BANs 25/27

slide-59
SLIDE 59

Non-deterministic updates

Class D and local (non-)monotonicity

Q: How are these two BANs related? f “ $ ’ ’ ’ & ’ ’ ’ % f0pxq “ x2 _px0 ^x1q f1pxq “ x3 _px0 ^x1q f2pxq “ x0 ^x1 f3pxq “ x0 ^x1 2 1 1 g “ # g0pxq “ x0 ‘x1 g1pxq “ x0 ‘x1 3 (S., 2012)

Sylvain Sené Synchronism vs asynchronism in BANs 26/27

slide-60
SLIDE 60

References

§ J. Aracena et al.. On the robustness of update schedules in Boolean

  • networks. BioSystems, 97:1–8, 2009.

§ J. Aracena et al.. Combinatorics on update digraphs in Boolean networks.

Discrete Applied Mathematics, 159:401–409, 2011.

§ E. Goles, M. Noual. Block-sequential update schedules and Boolean

automata circuits. Proceedings of AUTOMATA’2010, DMTCS, 41–50, 2010.

§ M. Noual, S.. Synchronism vs asynchronism in monotonic Boolean

automata networks. Natural Computing, doi:10.1007/s11047-016-9608-8, 2017.

§ S.. Sur la bio-informatique des réseaux d’automates, HDR, 2012. § M. Tchuente. Cycles generated by sequential iterations. Discrete Applied

Mathematics, 20:165–172, 1988.

Sylvain Sené Synchronism vs asynchronism in BANs 27/27