A 3-D agent based model (ABM) (where the agents are dendritic cells (DCs)) is presented that simulates the stochastic response of DCs to viral infection up to 15 hours post infection. Simulation results are in good agreement with the experimentally measured DDX58 transcript and IFN-β transcript distributions at 6 and 11 hours post infection.
The results show that both the spatially random distribution of infected DCs (only infected DCs can express IFNB1) and the broad distribution of times the infected DCs begin to secrete IFN-β lead to an initially heterogeneous distribution in IFN-β in the extracellular medium. During the first 6 hours post infection, the spatial heterogeneity of the IFN-β distribution is amplified in the spatial distribution of the bound IFNAR receptors across the DC population by extracellular communication (via autocrine signaling), see Figure 1. By 11 hours post infection the spatial distribution of IFNAR receptors begins to homogenize due to paracrine signaling between DCs. Bound IFNAR receptors activate the JAK/STAT pathway that induces the interferon stimulated genes (ISGs) such as DDX58 (by modulating the rate DDX58 transitions from a state that results in a basal rate of transcript production to a state that results in an enhanced rate of transcript production). The results show that the spatial heterogeneity of the bound IFNAR distribution leads to a broad distribution of DDX58 transcripts from 7 to 15 hours post infection. The ABM predicts the other ISGs’ transcript distributions will also be broad during these times, which may result in the functional consequence of a broad range of antiviral activity throughout these times. 
The stochastic intracellular processes of both infected and uninfected DCs are simulated using the standard Gillespie algorithm . The reactions are depicted diagrammatically in Figure 2 (see the supplementary section of  for detailed tables of reactions, rate constants and parameters, and the initial conditions). The method used to simulate the extracellular diffusion of IFN-β is described in detail in . The code is written in C++ and requires the user to download the GSL library in order to compile it. An annotated version of the source code for the 3-D ABM can be downloaded below. See Figure 1 for an example of the spatial distributions that can be produced using this program. For a more comprehensive description of the model see the methods section of .