ODe model describing the role of E2F1 in IL-6/IL-6R/STAT3 feedback loop system

We utilized previously published models (Vera et al. 2008, DOI: 10.1186/1752-0509-2-38, and Vera et al. 2013 https://doi.org/10.1158/0008-5472.CAN-12-4095 ) to derive a mathematical model in ordinary differential equations describing the activation of the signalling pathway IL-6/IL-6R/STAT3/IL-10. We created two instances of the model, one reflecting the activation of the pathway in melanoma cells, and another the activation in CD4 T cells, which differ in the expression levels attributed to the components of the IL-6 signalling pathway as indicated by transcriptomics data from monoculture experiments. 

For the case of simulating melanoma cells, we set up the model to a switch off configuration, modified iteratively the values of E2F1 expression in the model (from 0.1 to 100 a.u.), performed simulations, and computed the values of IL-6 at time 100 hr. We repeated the procedure but  modifying iteratively the expression levels of E2F1 and STAT3 (from 0.1 to 10 a.u.),  and computed the values of IL-6 at time 100 hr.

For the case of simulating CD4 T cells, we set up the model to a switch off configuration, utilized as input variable the vector of IL-6 abundance values generated in panel B left, performed simulations, and computed the values of IL-10 at time 100 hr.