The model should reproduce the figure 1C of the article (successfully reproduced in MathSBML). If your software does not support the variable "time", you can replace the assignmentRule:
n = n0 * [ exp(-kbN*time) + kappa * (1 - exp(-kbN*time))]
by
n = n...
