Figure 3a (Th17 Wild-type) has been reproduced here. The model was simulated using SBML odeSolver. The plots were obtained using Gnuplot.