Fig. 4
Bifurcation analysis of key system parameters V max1, V max2 and the complex parameter c 6/k YciRde. The stable states are shown by continuous lines while unstable fixed points are indicated by dashed lines. Note that each stable level of active YdaM is characterised by a distinct stable amount of c-di-GMP and YciR II (not shown). The parameter regions inducing bistability are marked by grey areas. Inlay plots show stochastic simulations of the system within the bistable region, indicating the probability (relative amount of cells) to contain a high or low amount of active YdaM molecules. The initial condition for the stochastic simulations was set to x 1=100,x 2=0,x 3=768 and 103 simulations were conducted for each parameter combination of the bifurcation analysis. A sufficiently large simulation time of 100 seconds was chosen in order to ensure equilibration. In each experiment A-C all the other parameters were kept fixed corresponding to the parameter set 1 in Additional file 3 (first data row of the table)