AbstractMicroglia function is orchestrated through highly-coupled signaling pathways that depend on calcium (Ca2+). In response to extracellular adenosine triphosphate (ATP), transient increases in intracellular Ca2+ driven through the activation of purinergic receptors, P2X and P2Y, are sufficient to promote cytokine synthesis and potentially their release. While steps comprising the pathways bridging purinergic receptor activation with transcriptional responses have been probed in great detail, a quantitative model for how these steps collectively control cytokine production has not been established. Here we developed a minimal computational model that quantitatively links extracellular stimulation of two prominent ionotropic puriner-gic receptors, P2X4 and P2X7, with the graded production of a gene product, namely the tumor necrosis factor α (TNFα) cytokine. In addition to Ca2+ handling mechanisms common to eukaryotic cells, our model includes microglia-specific processes including ATP-dependent P2X4 and P2X7 activation, activation of NFAT transcription factors, and TNFα production. Parameters for this model were optimized to reproduce published data for these processes, where available. With this model, we determined the propensity for TNFα production in microglia, subject to a wide range of ATP exposure amplitudes, frequencies and durations that the cells could encounter in vivo. Furthermore, we have investigated the extent to which modulation of the signal transduction pathways influence TNFα production. Our key findings are that TNFα production via P2X4 is maximized at low ATP when subject to high frequency ATP stimulation, whereas P2X7 contributes most significantly at millimolar ATPranges. Given that Ca2+ homeostasis in microglia is profoundly important to its function, this computational model provides a quantitative framework to explore hypotheses pertaining to microglial physiology.
Modelling human skull growth: a validated computational model
