As proposed earlier, the immune cells, cytokines and TD discussed above can be treated as network nodes. The aforementioned interactions between these nodes are then integrated into the network model shown in Fig 1 in the next section. The constructed network bears a multiple timescale feature. Cytokine regulation of cell function through signal transduction usually occurs on a sub-second timescale, for example, whereas cell production of cytokines takes minutes to hours . The cytokine regulation activity can thus be considered to be at steady state in the equations that describes the slow timescale activities of the cells. In this way a positive or a negative input can be modeled using an increasing or decreasing Hill function, respectively . The dynamics of the cytokines, the immune cells and TD can thus be described using a set of ordinary differential equations (ODEs). In this study, the system of ODEs is solved numerically using MATLAB (version R2013a Mathworks) with a variable order and multistep solver, ode15s. MATLAB is also used to plot the simulation data to generate the figures presented below.