Using the topology of the basic apoptosis network (Figure 1), constructed from the literature, we derived a computational model by generating rate equations for the network. The complete set of differential equations for the basic network N1 is given in Table 2. Concentration of each molecule is represented by enclosing the molecule name with square brackets (e.g. [Casp8]). Reaction rate constants are represented as rc1, rc2 etc. For example, based on the basic network N1 in Figure 1, the differential equations for the change of concentration of AKT and p53 were derived as: