In the driver model, each of the two daughter cells acquires *k*_{d} ∼ Pois(*m*_{d}/2)
driver mutations in each cell division, and can accumulate *N*_{d} driver mutations at maximum.
Each driver mutation increases the cell division probability by *f* fold: *g* = *g*_{0}*f*^{nd},
where *d*_{0} is the base cell division probability and *n*_{d} is the number of driver mutations accumulated in a cell.

Here, we incorporate punctuated evolution into the driver model to build the “punctuated” model. For the models considered so far, we assumed that cells can infinitely grows without a decrease of growth speed. However, it is more natural to assume that there exists a limit of population because of the resource limitation and the growth speed gradually slows down as the population size approaches the limit. The limit of population sizes called the carrying capacity is employed in the logistic equation. By mimicking the logistic equation, we modified the division probabilit as *g* = *g*_{0}*f*^{nd}(1 − *p*/*p*_{c}), where *p*_{c} is the carrying capacity. To reproduce punctuated evolution, we additionally employ an “explosive" driver mutation, which negates the effect of the carrying capacity. After a cell accumulates driver mutations up to the maximum *N*_{d}, the explosive driver mutation is introduced at a probability *m*_{e}. For a cell that has the explosive driver mutation, the carrying capacity *p*_{c} is set to infinite; that is, it is assumed that the explosive driver mutation rapidly evolves the cell so that it can conquer the growth limit and attain infinite proliferation ability.

Information of variables and parameters are listed in Tables 1 and 2.
In MASSIVE, we converted *m*_{d},*m*_{e}, *p*_{c} and *P* to log scale, i.e., *m*_{d}' = log_{10} *m*_{d},
*m*_{e}' = log_{10} *m*_{e}, *p*_{c}' = log_{10} *p*_{c}, and *P*' = log_{10} *P*, and then tested every combination of
*m*_{d}' ∈ {-2, -1.5, -1},
*m*_{e}' ∈ {-6, -5.5, -5,...,-3},
*p*_{e}' ∈ {3, 3.5, 4} and
*P*' ∈ {3, 4, 5, 6}.
All results are explorable in the focused and comparative view modes of the MASSIVE viewer.

Table 1. a list of the variabes

symbol | description |
---|---|

k_{d} |
number of driver mutations obtained in a cell division |

n_{d} |
number of driver mutations accumulated in a cell |

p |
population size |

t |
number of time steps |

g |
cell division probability |

Table 1. a list of the parameters

symbol | description | value |
---|---|---|

m_{d} |
expected number of driver mutations generated per cell division | {10^{-2}, 10^{-1.5}, 10^{-1}} |

m_{e} |
probability of acquiring an explosive mutation | {10^{-6}, 10^{-5.5}, 10^{-5},...,10^{-3}} |

N_{d} |
maximum number of driver mutations accumulated in a cell | 3 |

f |
increase of the cell division probability per driver mutation | 10^{0.3} |

g_{0} |
base cell division probability | 10^{-4} |

p_{c} |
carrying capacity | {10^{3}, 10^{3.5}, 10^{4}} |

P |
maximum population size | {10^{3}, 10^{4}, 10^{5}, 10^{6}} |

T |
maximum number of time steps | 10^{6} |