TD4: Toggle Switch

  1. Competing syntheses with cooperativity
  2. Stochastic semantics: toggle switch
  3. Continuous semantics: first-order approximation of the mean stochastic behavior

F. Fages 25 Jan 2019

Exclusive Switch (Loinger et al., Phy. Rev. E, 2007) Two genes with common promoter Slow or fast binding to the promoter with rate constant k Screen%20Shot%202019-01-24%20at%2022.42.07.png

Let us first consider a model with irreversible binding to the promotor

In [1]:
MA(s1) for DNA => DNA + X1.
In [2]:
MA(s2) for DNA => DNA + X2.
In [3]:
MA(d1) for X1 => _.
In [4]:
MA(d2) for X2 => _.
In [5]:
MA(b1) for DNA + X1 => DNAX1.
In [6]:
MA(b2) for DNA + X2 => DNAX2.
In [7]:
MA(sb1) for DNAX1 => DNAX1 + X1.
In [8]:
MA(sb2) for DNAX2 => DNAX2 + X2.
In [9]:
parameter(s1=1, d1=1, b1=1, sb1=10).
parameter(s2=1, d2=1, b2=1, sb2=10).
In [10]:
present(DNA, d). parameter(d=0.01).
In [11]:
draw_reactions.
In [13]:
draw_influences.
In [14]:
search_conservations.
Out[14]:
DNA+DNAX1+DNAX2
1 complex invariant(s)

Stochastic semantics

In [15]:
option(time:200, method:ssa).
In [16]:
numerical_simulation. plot.

Slow vs fast binding rates of the promotors X1 X2

  • Fast binding rates create cooperativity and a toggle switch stochastic behavior
  • Slow binding rates delete cooperativity and create oscillating behavior
In [17]:
%slider b1 b2 sb1 sb2

Boolean semantics

  • possible stability of DNAX1 or DNAX2 or no free DNA
  • possible boolean oscillations for X1 and X2
  • (oscillations are not expressible in Computation Tree Logic CTL, an over-approximation formula is used)
In [18]:
generate_ctl_not.
Out[18]:
reachable(stable('DNAX1'))
reachable(stable('DNAX2'))
reachable(stable(not'DNA'))
reachable(steady('DNA'))
reachable(steady('X1'))
reachable(steady('X2'))
reachable(steady(not'X1'))
reachable(steady(not'X2'))
oscil('X1')
oscil('X2')
In [19]:
expand_ctl(oscil(f)).
Out[19]:
EU(not f,f/\EU(f,not f/\EU(not f,f/\EU(f,not f/\EU(not f,f)))))/\EG(EF(f)/\EF(not f))
In [20]:
list_stable_states.
Out[20]:
[DNA-0,DNAX1-0,DNAX2-0,X1-0,X2-0]

Continuous semantics

  • Deterministic semantics
  • First-order approximation of the mean stochastic behavior
  • Same picture for slow and fast binding rates
In [21]:
parameter(b1=1, b2=1). option(method:bsimp).
In [22]:
numerical_simulation. plot.
In [23]:
plot(show:X2, against:X1).
In [24]:
%slider b1 b2
In [25]:
plot(show:X2, against:X1).

Questions

1) add unbinding reactions from the promotor with kinetic parameters u1, u2

In [26]:
MA(ub1) for DNAX1 => DNA + X1.
In [27]:
MA(ub2) for DNAX2 => DNA + X2.
In [31]:
parameter(ub1=2, ub2=2).
present(DNA, d). parameter(d=0.01).

2) Do stochastic simulations and explain the differences with the previous stochastic behavior

use the cells below for stochastic simulations, possibly Boolean model-checking, and

explain your answer here

...

In [32]:
option(time:200, method:ssa).
In [33]:
numerical_simulation. plot.

As expected, the bistability phenomenon is far less pronounced, since the promoters can unbind, so having one of the $X_i$ binding the DNA is not sufficient anymore to keep producing $X_i$ forever (in other words, the DNAXi states are not stable anymore).

3) Give parameter values to b1, b2, sb1, sb2 to recover similar stochastic behavior as before

use the cells below for stochastic simulations with sliders and

give your answer here

...

In [38]:
%slider b1 b2 sb1 sb2 ub1 ub2

The bistability phenomenon reappaears when $b_1, b_2 >> ub_1, ub_2$ (as expected: when binding is way more likely than unbinding, resulting in a situation similar to the without-unbinding one).

For example, with:

Parameter Value
$b_1, b_2$ 100
$ub_1, ub_2$ 1
$sb_1, sb_2$ 1000