Abstract

The spatial range of cytokine communication is a crucial, yet poorly understood, functional property. Both containment of cytokine action in narrow junctions between immune cells (immunological synapses) and global signaling throughout entire lymph nodes have been proposed, but the conditions under which they might occur are not clear. Here we analyze spatially three-dimensional reaction-diffusion models for the dynamics of cytokine signaling at two successive scales: in immunological synapses and in dense multicellular environments. For realistic parameter values, we observe local spatial gradients, with the cytokine concentration around secreting cells decaying sharply across only a few cell diameters. Focusing on the well-characterized T-cell cytokine interleukin-2, we show how cytokine secretion and competitive uptake determine this signaling range. Uptake is shaped locally by the geometry of the immunological synapse. However, even for narrow synapses, which favor intrasynaptic cytokine consumption, escape fluxes into the extrasynaptic space are expected to be substantial (220% of secretion). Hence paracrine signaling will generally extend beyond the synapse but can be limited to cellular microenvi-ronments through uptake by target cells or strong competitors, such as regulatory T cells. By contrast, long-range cytokine signaling requires a high density of cytokine producers or weak consumption (9.9., by sparsely distributed target cells). Thus in a physiological setting, cytokine gradients between cells, and not bulk-phase concentrations, are crucial for cell-to-cell communication, emphasizing the need for spatially resolved data on cytokine

signaling.

Author Summary

The adaptive immune system fights pathogens through the activation of immune cell clones that specifically recognize a particular pathogen. Tight contacts, so-called immunological synapses, of immune cells With cells that present ‘digested’ pathogen molecules are pivotal for ensuring specificity. The discovery that immune responses are regulated by small diffusible proteins — the cytokines — has been surprising because cytokine diffusion to ‘bystander’ cells might compromise specificity. It has therefore been argued that cyto-kines are trapped in immunological synapses, Whereas other authors have found that cyto-kines act on a larger scale through entire lymph nodes. Measurements of cytokine concentrations with fine spatial resolution have not been achieved. Here, we study the spa-tio-temporal dynamics of cytokines through mathematical analysis and three-dimensional numerical simulation and identify key parameters that control signaling range. We predict that even tight immunological synapses leak a substantial portion of the secreted cyto-kines. Nevertheless, rapid cellular uptake will render cytokine signals short-range and thus incidental activation of bystander cells can be limited. Long-range signals will only occur with multiple secreting cells or/ and slow consumption by sparse target cells. Thus our study identifies key determinants of the spatial range of cytokine communication in realistic multicellular geometries.

Introduction

In particular, the release, sensing and uptake of cytokines, small signaling proteins, by cells is essential for the regulation of the mammalian immune system [1]. Prominent quantitative characteristics of cy-tokine signaling are high receptor specificity (with Kd % 10—10 nM) and low free cytokine concentrations in the picomolar range [2,3]. The physiological cytokine milieu regulates critical processes like the type and strength of the immune response. Quantitative understanding of such cytokine-driven cellular decisions is beginning to emerge [4—8] , yet the underlying spatio-temporal cytokine dynamics remain poorly understood. Cytokines act in a heterogeneous environment, typically with high cell-densities. It is not known how they diffuse under such conditions and, in turn, regulate immune responses. Specifically, how far cytokines can signal away from the producing cell is not clear. Perona-Wright et al. [9] have found that interleukin(IL)-4 is seen by most T cells in the lymph node upon parasite infection, including nonspecific ‘by-stander’ cells. In this case, many T cells throughout the lymph node could be IL-4 producers. By contrast, several observations suggest more localized cytokine communication [4,10—13].

Of note, 1 pM is about 1 molecule in 1700 um3, compared to ~500 um3 volume of a typical lymphocyte. Higher, systemically elevated cytokine levels arise only in certain immunopathologies, so-called ‘cytokine storms’, where they cause severe damage [14]. However, it has been demonstrated that cyto-kine concentrations are not always well mixed, and locally higher cytokine concentrations can occur also in eX vivo T cell cultures [12]. Therefore, we asked how and under which conditions such cytokine gradients arise, and if they are able to explain effective paracrine signals.

Immunological synapses are formed between immune cells by surface proteins after antigen recognition [15,16,19]. They have been observed between various cell types of the immune system, including immunological synapses between T cells and antigen presenting cells (APC, e.g. B cells [10] and dendritic cells [20]), and immunological synapses between T cells and T cells [12]. Many cytokines are secreted preferentially into the immunological synapse [10,21—23] , and a range of high-affinity cytokine receptors have been found to be specifically located in the immunological synapse, too [11]. Therefore, it is likely that the synapse has an important function for cytokine signaling, beyond its role for T cell receptor signaling on which theoretical studies have focused [24,25]. Cytokine signaling through immunological synapses might also explain the pleiotropic effects observed for most cytokines, as it would provide specificity of cytokine signaling by restriction of their action. Nevertheless, it is unlikely that paracrine cytokine signals are only possible between cells that are directly connected by an immunological synapse. For instance, Sanderson et al. [26] found that interferon-y can be seen by bystander cells other than the target cells to which the synapses are formed. To understand which parameters govern autocrine versus paracrine cytokine signaling, we analyzed in this study reaction-diffusion models of cytokine signaling at two scales: through the immunological synapse between two cells and in three-dimensional arrays of many (>100) cells.

IL-2 was first identified as a T cell growth factor [27] , but, paradoxically, is a critical mediator of immune tolerance [28—31]. It is secreted by T helper (Th) cells early after antigenic stimulation and taken up by high-affinity IL-2 receptor (IL-2R) on Th cells and regulatory T (Treg) cells [28—31]. Treg cells mediate immune tolerance and are critical for the prevention of autoimmune reactions [32,33]. IL-2 secretion is digital, i.e. upon receiving an antigen stimulus, only about one quarter of a Th cell population releases IL-2 molecules [34—36]. It is an open question if IL-2 and other cytokine signals act in an autocrine or paracrine manner [31,37]. In response to IL-2 uptake, Th cells and Treg cells upregulate CD25, the oc-subunit of the IL-2R. CD25 is often used as an activation marker of T cells, because it precedes proliferation of Th cells and subsequent recruitment of effector immune cells [30,38]. Although IL-2 secreting Th cells upregulate CD25, Long and Adler [37] reported that they lack phosphorylated STAT5, a key intermediate in the IL-2R signal transduction cascade. In the same experiment, other Th cells not secreting IL-2 also upregulate CD25 in response to IL-2, and in addition show fully functional signal transduction [5,37]. These data suggest that the dominant mode of IL-2 signaling is paracrine, in contrast to the presumed function of the immunological synapse in containing secreted cy-tokines [16,17]. However, unlike T cells, most APC do not express functional IL-2 receptor (IL-2R) [39]. Thus, both the study by Sanderson et al. [26] and the properties of IL-2 signaling suggest a role of the immunological synapse for cytokine signals that goes beyond signal amplification between the two cells associated by a synapse.

To this end, we considered the two key spatial scales, the sub-um scale of the immunological synapse and the supra-um scale of cell-to-cell communication. We investigated reaction-diffusion models on these two scales by analytical techniques and advanced finite-element computations in three spatial dimensions [40—44]. To be specific, we utilized a simple, experiment-based mathematical model for IL-2 signaling and gained more general insight through systematic variation of parameters. Our results show that paracrine cytokine signaling is possible in the presence of local concentration gradients combined with nonlinear signal amplification. The spatial range of cytokine signaling can be tuned from purely autocrine via intrasynaptic and short-range paracrine to long-range paracrine. For a wide array of parameters, we found that cytokine gradients in dense multicellular environments range over one to few cell diameters. These computational findings can inform novel experiments probing the spatio-temporal dynamics of cytokine signaling [45].

Results

Target cell density and receptor expression control paracrine cytokine signals

Thus, regulating the strength of cytokine signaling by cytokine receptor expression might also affect the extracellular cytokine concentration and hence, indirectly, signaling. To gain quantitative insight, we first studied a simple reaction-diffusion model, where a cyto-kine-secreting cell is surrounded by cells that can take up the cytokine. To allow for an analytical solution, we assume the surrounding cells to be placed on a spherical shell with the secreting cell in the center (Fig 1B, see Materials and Methods). For convenience, parameter values are summarized in Table 1. If the target cells are located far away (i.e., their density is low), the cytokine concentration experienced by the target cells is nearly independent of the

A O B O surrounding

secretion internalization

Paracrine cytokine signals depend on cell density and receptor number. (A) Cytokine secretion and uptake is followed by receptor internalization, leading to a reduced cytokine concentration in the medium. (B) Schematic of the high cell-density scenario. A cytokine secreting cell is surrounded by a layer of responder cells that provide a diffusion barrier for the cytokine. Signaling can be autocrine (Jauto), i.e. cytokine molecules are bound by receptors of the cytokine secreting cell, or paracrine (Jpara), i.e. bound by receptors on responder cells. (C) Cytokine concentration profile in the model with homogeneous secretion and uptake. (D) Cytokine concentration 1 pm away from the cytokine secreting cell, in the limits of high and low cell-density. (E) Paracrine signal Jpara and autocrine signal Jauto (inset) as fraction of secreted cytokine molecules. (F) Autocrine and paracrine uptake in the high cell-density scenario. Dotted lines are approximations in the limit of fast diffusion (see S1 Text). Parameter values: R = 4000 mol./cell, and see level of receptor expression (Fig 1C) because the dilution of the cytokine occurs primarily by diffusion in the three-dimensional tissue. On the other and, if the density of target cells is so high that they immediately surround the cytokine secreting cell, the cytokine concentration is practically homogeneous in the small intervening space, as the timescale of diffusion over such a short distance is fast compared to the timescale of cytokine uptake (Fig 1D and 81 Text). As a consequence, the cytokine concentration experienced by proximal target cells is set by the balance of secretion rate by the cytokine-producing cell and uptake rate. The autocrine and para-crine uptake rates Iauto and Ipara depend on the level of cytokine receptor expression on the target cells (Fig 1E), and are practically independent of the cell-to-cell distance even at high cell density (Fig IF; the low cell-density scenario is independent of the cell-to-cell distance by construction). Interestingly, cytokine concentration (Fig 1C) and uptake rates (Fig ID) are sensitive to receptor expression on proximal targets cells in the physiologic range of 100 to several 1000 receptor molecules per cell [5]. Thus, this simple model indicates that with a high density of target cells, cytokine receptor expression controls the amount of paracrine cytokine signal.

The immunological synapse controls type and strength of cytokine signals

However, T cells release IL-2 and other cytokines in a polarized fashion into the immunological synapse [10,21—23]. Therefore, we analyzed a model of cytokine secretion and uptake in the immunological synapse, represented by a small cylindrical region between a Th cell and an opposed APC or second T cell (Fig 2A, see Materials and Methods), extending previous work [46]. The distance between Th cell and opposed cell, in the following referred to as synaptic distance, is in the range of 10 to 40 nm [19,47]. This close contact between Th cell and opposed cell causes a cytokine concentration profile which is almost homogeneous between the two cells in the center of the synapse, and sharply falls off towards the outer boundaries through which cytokine molecules are lost practically irreversibly (Fig 2B, top). In the case of low receptor expression (Fig 2B, top left), the cytokine concentration reaches values in the nM

Thus, the synaptic cytokine secretion results in locally much higher concentrations than homogeneous secretion (see Fig 1C), in line with experimental data [12]. For comparison, consider cytokine secretion into a cylindrical region with length 2 pm, a typical value for nearby cells but much larger than the immunological synapse (Fig 2B, bottom). In this case, the cyto-kine concentration falls to less than 20 pM at the surface of the opposed cell. Hence, the very small synaptic distance in a fully formed immunological synapse is crucial for the establishment of high local cytokine concentrations.

Cytokine molecules may also induce autocrine signals by binding to receptors at the secretory cell (Iauto). Fig 2C shows the fractions of Iauto, Isynapse and Iescape, choosing IL-2 receptor densities that are characteristic of naive (IL-2R10W) or preactivated (IL-2Rhigh) IL-2 secreting T cells, and for opposed cells with different IL-2R expression. IL-2Rhigh cells recapture most of the secreted IL-2 molecules, irrespective of the type of opposed cell and the synaptic distance. Na’ive, IL-2Rlow cells show a strong dependence on the synaptic distance (Fig 2C, left). If the synaptic space is sufficiently narrow, Iescape is small; the escape flux could be even further reduced by adhesion molecules sealing off much of the synapse from the extracellular space. If the opposed cell is a second T cell, then the secretory cell and the opposed cell compete for the cytokine molecules. Treg cells outcompete Th cells due to their large receptor number. On the other hand, APCs do not express IL-2R and the IL-2 signals would be purely autocrine. For synapses with somewhat larger synaptic distances, a considerable amount of cytokine molecules can escape (enhanced Iescape) and provide paracrine signals to surrounding cells. Interestingly, the ratio of Iauto and Iescape is most sensitive to the synaptic distance in the physiologic range between 10 nm in the close-contact zone and 40 nm in the outer region of the immunological synapse [19,47]. In all cases, the fraction Jescape of cytokine molecules that escape for paracrine signaling is considerably smaller than in the case of homogeneously distributed cyto-kine secretion and receptor expression (Figs 2C and 1D). However, even if the cytokine secreting cell is pre-activated and hence autocrine IL-2 uptake is high (IL-2Rhigh), a sizeable fraction of cytokine molecules still diffuses out of the synapse (Jew,pe ~20%).

Both properties result from the high aspect ratio of the radius of the cell contact area and the synaptic distance (r and z in Fig 2A). The two properties have opposing effects on the strength of paracrine cytokine signals: While localized cytokine concentrations increase the likelihood of a local paracrine signal, the reduction in effective cy-tokine secretion reduces the potential of paracrine signals.

ln-silico Th cell culture exhibits localized paracrine lL-2 signaling

For example, the high aspect ratio of the immunological synapse evokes highly localized cytokine concentrations in the vicinity of cytokine secreting cells resembling secretion from a point source (see Fig 2B), and the high diffusion constant in relation to the receptor dynamics makes the system largely independent of the cell-to-cell distance (Fig 1E and IF). However, the simple models studied above cannot answer the question if effective paracrine signals are possible despite the low bulk cytokine concentrations. To illustrate this problem, consider a classical formula from Berg and Purcell for the timescale of ligand diffusion towards a receptor [48] (Materials and Methods). Measured cytokine concentrations in serum or in supernatants of ex-vivo T cell cultures are typically in the picomolar range [2,3]. Assuming a spatially uniform cytokine concentration of 10 pM and a receptor number of 100 per cell, as is typical for the high-affinity IL-2R on naive T cells, that calculation reveals that on average, every 7 min a receptor becomes bound by a cytokine molecule. Under these conditions it would take hours to induce a reliable signal, indicating that bulk cytokine concentrations might just be capable of, or even be too low for, stimulating signal transduction. However, it has been reported that IL-2 is subject to appreciable spatial gradients, with much higher concentrations at the surfaces of T cells [12,49].

As before, we focus on the cytokine IL-2, for which many parameters, including secretion and receptor expression rates, have been estimated from experiments [5,35,50], and experimentally tested models for the IL-2R dynamics are available [4,5,7]. To account for polarized secretion at the immunological synapse, we do not explicitly model synapse formation but consider the effect of discrete IL-2 sources from which IL-2 escapes into the extra-syn-aptic space (with rate qeff, corresponding to Iescape in the simplified model of the previous section). The position of the IL-2 source of a producing cell is a randomly chosen point at the cell surface. IL-2 secretion is all-or-nothing [34,36]: only about one quarter of antigen-stimu-lated T cells release IL-2 molecules, and among these cells, the IL-2 secretion rate is in the range of 10 molecules per second [35]. In accordance with experimental data, already activated IL-2 secreting cells have high IL-2R expression which, for simplicity, we take as constant [37]. Non-secreting cells are assumed to upregulate IL-2R expression in response to IL-2 homogeneously at their cell surface.

To focus on the role of IL-2 uptake by T cells, we do not consider the APC explicitly, but only its consequences for polarized secretion and uptake (see above) (Table 1). Despite this simplification, our simulations consider realistic extracellular volumes (as determined by the cell distances) between the cells as a basis for determining the extracellular concentrations of the secreted cytokines.

Stimulating IL-2 secretion in a fraction of T cells (Fig 3A and 3B), the IL-2 concentration increases rapidly and nearly homogeneously for several hours after stimulation (Fig 3C and 3D and SI Fig). Then, in response to the high IL-2 concentration resulting from paracrine signaling, IL-2R expression is upregulated in non-secreting cells (Fig 3E) and causes fast IL-2 uptake from the medium. As a result, concentration gradients occur: In large parts of the simulated region, the IL-2 concentration reaches a steady-state at around 10 pM while locally it is more than twice as large (Fig 3C, red regions at 30 h). This inhomoge-neity in IL-2 concentration corresponds to receptor upregulation (activation) of non-secreting Th cells: IL-2Rhigh cells are found near the regions with high IL-2 concentration. Analysis of the time course (Fig 3E) shows that all cells upregulate IL-2R levels in response to the increased IL-2 concentration in the first hours after antigenic stimulation and IL-2 secretion. However, as the high-affinity IL-2R is being upregulated, IL-2 becomes increasingly depleted in the medium. As a result, only a fraction of the cells receive a suff1cientIL-2 stimulus to sustain high IL-2R expression (IL-2Rhigh cells in Fig 3E), whereas the remaining cells downregulate IL-2R expression (IL-2R10W cells).

Later, at steady-state, the IL-2 concentration is somewhat higher in the microenvironment of IL-2R10W cells, because they do not consume as many IL-2 molecules. This form of local bistability, which occurs in the expression of IL-2R on Th cells, was observed already in Ref. [4]: Based on a quasi-stationary state assumption, Busse et al. showed that in the model without Treg cells, the IL-2R expression rate responds to the increase of the secretion rate in a digital way and the cells are activated only after a certain threshold is exceeded. A small bistable region around the threshold is observed. These findings were supported by experimental data from primary T cells cultured ex vivo [4]. Thus our present model with the immunological synapse and 3D diffusion matches the bistable system behavior seen in the simpler analytical model. Taken together, our simulations indicate that the amount of IL-2 escaping from the immunological synapse is sufficient to sustain paracrine signaling in at least a fraction of surrounding cells. However, competition for the cytokine can cause heterogeneity in the response of a cell population and result in bulk IL-2 levels that are much lower than local concentration peaks and in agreement With concentration levels measured by ELISA (see Discussion).

Competitive lL-2 uptake by regulatory T cells

To study the effect of Treg cells on the IL-2 dynamics after activation of conventional Th cells, we simulated a T-cell population consisting of antigen-stimulated IL-2 secreting and non-secreting Th cells as well as Treg cells (Fig 4A and 4B). Compared to the situation in the absence of Treg cells (cf. Fig 3), the IL-2 concentration attains a spatially inhomogeneous steady state more rapidly, with the overall IL-2 concentration being lower (Fig 4C and 4D and 82A Fig). Importantly, the non-secreting Th cells do not permanently upregulate IL-2R in the presence of Treg cells because the Treg cells suppress the paracrine IL-2 signal. The comparison with the simulations without Treg cells (Fig 3) imply that Th cells require for sustained IL-2 signaling both a transient strong and a stable weak IL-2 stimulus. The finding that Th cells can sustain IL-2 signaling at low cytokine concentration, but only after initial stimulation with high cytokine concentration, is a spatio-temporal phenomenon similar to hysteresis: Active cells express more cytokine receptors, which bind more cytokine molecules even at lower concentration and thus stabilize the active state once it is achieved. Treg cells can suppress prolonged IL-2 signaling in Th cells by inhibiting the strong initial IL-2 signal and the resulting upregulation of the high-affinity IL-2R.

Without Treg cells, the number of activated Th cells increases linearly with the effective IL-2 secretion rate qeff until, eventually, all cells in the simulated region become active (Fig 4E, left panel). By contrast, the presence of Treg cells creates a threshold at an effective secretion rate of qeff ~ 20000 molecules/h (about 5 molecules/s), below which there is no paracrine IL-2 signaling between Th cells. The same pattern is observed if we vary the fraction of cytokine secreting cells instead of the effective secretion rate (Fig 4E, middle panel), which reflects digital IL-2 secretion [34,36]. Hence the presence of Treg cells changes the paracrine IL-2 dynamics from a gradual to an all-or-none response: Either the paracrine signal is completely suppressed by competitive uptake, or suppression is overrun and all cells are activated. Note that qeff measures only the IL-2 molecules that escape from the immunological synapse; assuming a tight synapse, this would only be 20% of the total secretion (see Fig 2). However, measured IL-2 secretion rates are ~10 molecules/s [35,50], which is likely to be too small to titrate out the Treg cells in a physiological setting where IL-2 is secreted into the synapse.

This is because, as anticipated by the analytically treatable model (see Fig 1F), cytokine molecules can reach nearby cells rapidly by diffusion compared to the slower time scales of changes in IL-2R expression and IL-2 internalization. Thus, the exact cell-to-cell distance is unimportant in the physiological range.

lL-2 producers surrounded by lL-2-responsive cells produce short-range paracrine signals

Generally, we expect that the balance between cytokine secretion, dilution through diffusion in the three-dimensional extracellular space and cellular consumption will determine the signaling range. To understand the interplay of these three factors, we performed large-scale simulations of an area containing ~2000 cells, with a single IL-2-secreting Th cell surrounded by non-secreting Th cells which all are potential responders to the IL-2 (Fig 5A). Although we use the specific parameters for IL-2 here, this model is of more general interest and applies to other situations with few signaling cells and many responder cells (e.g., IL-4 secreting Th cells in a B cell population [9]), or can be thought of as representing a cluster of several cytokine secreting cells in a population with a small density of cytokine secreting cells elsewhere.

Remarkably, although secretion is assumed to be polarized through the synapse, the cytokine concentration is higher along the entire surface of the secreting cell, including the pole opposite to the synapse, than at nearby cells This is due to the absence of IL-2R on the surface of secreting cells (except for the synaptic space). For larger secretion rates (of the order to 106 molecules/h or 280 mole-cules/s), the IL-2 signal reaches hundreds of cells. However, with the experimental estimate for the IL-2 secretion rate (10 molecules/s [35,50]), of the order of a 100 secreting cells would be needed to realize such a high rate (assuming an effective secretion rate of 10-20% of the total rate, see Figure 2). Therefore, IL-2 from an individual producer will act locally whereas only large clusters of activated cells could cause long-range signals. The occurrence of two distinct spatial signaling regimes as a function of secretion rate is expected because the cellular uptake rate can be saturated by high cytokine concentrations (akin to an enzymatic Michaelis-Menten rate law where the cytokine receptors function as the enzyme). Below saturation the cytokine signal remains local. Interestingly, the spatial range scales linearly with the logarithm of the secretion rate (Fig 5C). Hence the signaling range exhibits a fold-change response to the effective secretion rate (see Discussion).

For this purpose, we simulated a pulsed, homogeneous stimulation (see Methods) in a region covering ~5000 cells. We found that despite the apparent short-range induction of effective paracrine signaling, the traveling distance has a broad distribution peaking around four cells away from the cytokine secreting cell in each direction (Fig 5D). Thus, in our reaction-diffusion system, the chemical reactions on the cell surface dominate the diffusion and determine the IL-2 gradient formation. We further compared the distribution of traveling distances with earlier analytical expressions obtained from a reaction-diffusion model of mor-phogen gradient formation [51]. Despite some differences in the model architecture (see Methods), and despite numerical limitations in simulating an ‘infinite domain’ as assumed by the analytical methods of Ref. [51], our simulations are in good agreement with those analytical results (S3 Fig). Taken together, our simulations of long-range cytokine diffusion and uptake show that long-range paracrine signals are possible in principle, but require exceptional circumstances (extremely high rates of cytokine production or large clusters of cytokine producing cells) that might not readily occur in vivo.

Discussion

Experimentally, however, cytokine signaling has not been probed directly at fine spatial resolution, although recent advances in synthetic biology could provide new tools in the near future [45]. Here, we used a computational approach to study cytokine signaling in realistic three-dimensional geometries. To this end, we considered two distinct spatial scales. First, we analyzed polarized signaling across narrow junctions — immunological synapses — between immune cells (nm scale). We find that synapses enhance autocrine signaling and signaling towards the cell connected by the synapse, but, importantly, cannot prevent substantial cytokine escape for paracrine communication by mere geometry. Second, we employed advanced simulation tools for partial differential equations to dissect the dynamics of this ‘spill-over’ paracrine signal in dense ensembles of hundreds of communicating cells (um scale). Using experimentally established parameters for the T-cell cytokine IL-2, we find that cytokine signals emanating from producing cells are short-range (one to few cell-to-cell distances) because of uptake by target cells or competitors. Long-range communication requires coherent secretion by tens to hundreds of producers or/ and sparse uptake. Thus we predict that gradients at the cellular length scale are a key property of cell-to-cell communication by cytokines. We note that the spatial range of diffusible signals is also of relevance for morphogen action [52,53]. In contrast to immune cell signaling with a typical time scale of many hours during which diffusive gradients reach steady state, the transient behavior on shorter time scales is of particular interest for morphogen gradients [51,54].

As low cytokine concentrations would imply long signaling times (see Eq 1 below), we hypothesized that paracrine cytokine signals rely on much higher cytokine concentrations in the microenvironment of target cells, which have indeed been detected by live cell imaging [12]. However, it is generally believed that cytokine signaling occurs in the regime of fast diffusion, which is reflected by our parameter values—the typical spatial range of fast diffusion, \ /D / k d (see e.g. [55]), spans 40 cells away from the cytokine secreting cell for our values (see Table 1). Therefore, we analyzed the spatiotemporal dynamics of cytokine signals by more detailed mathematical modeling and simulations. We found that spatial gradients do occur due to nonlinear receptor dynamics and polarized IL-2 secretion at the immunological synapse, despite fast diffusion. This was also quantified, e.g. in terms of the traveling distance of cytokine molecules (Fig 5D). We showed by extensive simulations in three spatial dimensions that such cytokine gradients can mediate paracrine signals targeting cells other than those connected by the immunological synapse, as previously suggested for interfer-on-y [26]. Moreover, we analyzed the parameters that control paracrine signaling on the different spatial scales.

IL-2 has initially been thought of as a prototypical autocrine signal facilitating self-activation of Th cells [27,30,56]. More recently, paracrine IL-2 signaling towards Treg cells was identified as essential to prevent autoimmune diseases [29,31], possibly due to competition with autocrine self-activation [4,5,28,57]. Recent experimental observations suggest that also paracrine IL-2 signals towards other Th cells are important for regulation of immune responses, while true au-tocrine IL-2 signals are suppressed by the intracellular signal transduction pathway [5,37]. A plausible explanation would be that IL-2 secreting cells are constitutively activated, i.e. prone to proliferation and differentiation, due to a strong signal from the T cell receptor, and do not rely on signals via the IL-2R. Th cells not secreting IL-2 may have received a weaker T cell receptor signal, and are only fully activated if they receive additional stimulation from the IL-2R. In this theoretical study, we cannot address the question to what extent such a mechanism is responsible for the activation of T cell populations in vivo. However, our simulations show that the need for paracrine cytokine signals provides several checkpoints for the induction of immune responses downstream of the T cell receptor.

First, cytokine receptors have high affinity and are internalized after binding of cytokine molecules. That allows for control of paracrine cytokine signals by expression of cytokine receptors (Fig 1). Second, the effective rate of cytokine secretion, i.e. the paracrine cytokine signal escaping the immunological synapse, sensitively depends on the configuration of the immunological synapse, in terms of the exact synaptic distance (Fig 2). Therefore, we propose that regulation of cytokine signals is an important function of the immunological synapse (see also Refs. [15—18] ), along with regulation of the strength of T cell receptor activation [47] and the exchange of microvesicles between T cell and AFC [58]. Note that using the synaptic distance is an idealization; in reality the influence of the immunological synapse on cytokine diffusion is more complex, due to its structure consisting of several layers with different types of surface proteins [16,24]. Third, in our simulations, Treg cells efficiently suppress paracrine IL-2 signals, because they express high basal levels of IL-2R, preventing the strong transient cytokine signal. In line with earlier work from us and others [4,5,29,57], this suggests that suppression of IL-2 signals is an important mechanism contributing to immune tolerance mediated by Treg cells. Of note, Treg cells most likely interfere with T cell activation in several other ways, e.g. by release of antiinflammatory cytokines like IL-10, by forming immunological synapses with T cells, and by other mechanisms yet to be discovered [12,33]. Interestingly, a fourth system property one might expect to have a large influence on the dynamics of the system, the cell density or cell-to-cell distance, is unimportant for the results of our simulations (Fig 4E). This property results from the timescale separation between cytokine diffusion and cytokine uptake (see Fig 1E and SI Text), and explains recent experimental data [7].

Such a transient cytokine signal has been observed by single-cell IL-2 capture assays [49,59], and recently also by ELISA in cell supernatants [7], although with conflicting time-scales: IL-2 capture assays evoked a peak in the number of IL-2 secreting cells at 1—6 hr after antigen stimulation [49,59], while Tkach et al. report a peak in the IL-2 concentration measured in vitro after ~50 hr [7]. Our simulations point to an IL-2 peak in the first 10 hr after stimulation, and thus support the earlier suggestion [49] that ELISA studies have limitations in reflecting the time-course of in vivo cytokine signals, although the study of Tkach et al. provides valuable quantitative insight into the dose-response characteristics of IL-2 signals. A reason might be that in culture, cells form thin layers on the bottom of the well, and therefore cytokine molecules are detected by ELISA in the supernatant after a certain delay.

That means, the cell population recognizes relative rather than absolute increases in the stimulus strength (here, the amount of secreted cyto-kine molecules per time). Fold-changes in sensory biological systems are a classical phenomenon referred to as Weber’s law, and were recently observed in various intracellular signal transduction pathways [60—63]. As a consequence of the fold-change response, sensory systems can act over a broad range of stimulus intensities, from nearly detectable to very intense stimulations. Our computer simulations suggest a similar mechanism for paracrine cyto-kine signals: Moderate effective secretion by a small fraction of cells allows for short-range signals inside an immunological synapse, larger effective secretion rates may evoke paracrine signals that reach bystander cells in close vicinity but not connected by a synapse, and very high secretion rates or large clusters of secreting cells may evoke an organ-wide cytokine signal or ‘cytokine storm’ [14]. Adaptive immune responses must be rapid and effective in the case of strong infection, but also carefully controlled to avoid autoimmune diseases. In our simulations, the spatial distribution of cytokine secretion and uptake within a population of immune cells had a huge impact on the cellular response, generating multiple layers of plasticity that can be exploited for appropriate regulation of immune responses.

Materials and Methods

Software

The system was discretized in time by the damped Crank-Nicolson method. The intercellular area was discretized with an unstructured adaptive mesh, which describes each cell with at least 342 degrees of freedom (64 in the long range simulations in Fig 5) by a Galerkin approach using continuous finite elements (Q1). The discretized system was solved efficiently by controlling the error with adaptive space and time grids by means of the Dual Weighted Residual (DWR) method[42,44]. To allow for larger time steps, the equations were solved in a fully coupled fashion and not with the commonly applied iterative segregating approach. We linearized the nonlinear equations with Newton's method and applied Krylov-Space methods (GMRES) with a geometric multilevel preconditioner [40,43] to solve the resulting linear equations. In the simulations, the secreting Th cells and the Treg cells and the synapse on the cell surface of secreting cells were positioned randomly. We checked the influence of this cell positioning on the simulations with different randomly chosen positions and found that the variations between simulations were negligible.

For the graphical representation of the three-dimesional scalar data, here the IL-2 distribution in space, two methods were applied, the visualization of iso-surfaces using topological methods [64,65] and volume rendering [64]. With the first method specific isosurfaces are visualized by varying the transparency for different isovalues to get an impression of the 3D data set (Figs 5A, 81 and 82A). To choose these specific isosurfaces with important features, topological information (Morse complex, persistent homology classes and Betti numbers) is computed. The rendering was performed by using the Visualization Toolkit VTK (http://www.vtk.org) which allows rotation in real time. The second method, volume rendering, produces the image directly from the data without an intermediate geometrical representation. A play with transparency of the whole data set makes the inner structures visible (Figs 3D and 4D). With flexible mapping of the data on colors and opacity, different structures can be visualized efficiently and a realistic representation is obtained (Figs 3C and 3D and 4C and 4D). Difficulties in the data-representation were the wide range of the values over several orders of magnitude and the porous domain (extracellular domain).

The simulation is then run until the concentration reaches zero in the whole area. The fraction of the released IL-2 concentration bound by a certain responder cell is equivalent to the probability that the ‘secret-ed molecule’ was bound. This probability was calculated for the successive layers of responder cells surrounding the secretory cell, in order to obtain the distribution of the traveling distance. Analytical calculations were supported by Wolfram’s Mathematica. Matlab from Math-works was used to generate plots and to calculate the special functions applied in Fig 2.

Characteristic time of cytokine signaling

A classical formula derived by Berg and Purcell approximates the characteristic time T of a ligand diffusing towards a receptor [48]:

Note that in Eq 1 and in the following, cytokine concentrations (nM) are implicitly converted to molecules/um3 by Avogadro’s constant N A, wherever necessary, as follows: nM = 10'9 mol/l = 10'9NA molecules/(1015um3) = 6/ 10 molecules/um3. Note that the time to diffuse towards a T cell (first term in Eq 1) is less than a second, but the mean time to reach a receptor at the cell surface (second term in Eq 1) is in the order of minutes due to the small number of receptors on naive T cells.

Homogeneous cytokine secretion and uptake

The cytokine secreting cell has R cytokine receptors, and responder cells have Rresp cytokine receptors, both binding cyto-kine molecules in their immediate Vicinity with rate kon. We assume homogeneous cytokine secretion and uptake, so that the system has radial symmetry. As diffusion is fast (D = 10um2/s, see Table 1), it reaches a steady state after about LZ/D = 0.5 s, Where L is the cell-to-cell distance in the case of high cell-density. Thus, it is sufficient to consider the diffusion equation in steady state in the extracellular domain with flux boundary condition at the cell surface: c(r) is the cytokine concentration at distance r from the center of the cell, A is the Laplace operator in spherical coordinates, p is the cell radius, and q is the cytokine secretion rate. Note that cytokine concentrations are implicitly converted from unit nM to unit molecules/um3, as above. The boundary condition on the outer boundary is either (low cell-density limit) or (high cell-density limit) where N is the number of IL-2 consuming responder cells. In both cases, the problem can be solved analytically for the cytokine concentration c(r) and eventually for the uptake rates

Cytokine diffusion in the immunological synapse

We consider stationary cytokine diffusion in a cylindrical region between a cytokine secreting Th cell and a responder cell, both potentially expressing cytokine receptors (see Fig 2A). This leads to the following boundary conditions at the cytokine secreting cell (2 = 0) and the

At the outer boundary of the synapse, we assume c(a,z) = 0, which means that cytokine molecules Which escape the cylindrical region do not return to it. The cytokine concentration, and the uptake rates Iauto, Iescape and Isynapse resulting from this model, can be calculated analytically using Bessel functions (see 81 Text).

In silico T cell population

Moreover, due to recent experimental observations [5,35], we discard the previously assumed positive feedback from IL-2 uptake to IL-2 secretion, and we set the IL-2 secretion rate to 10 molecules/s and the fraction of IL-2 secreting cells to about 25% (see Table 1). In brief, the model [4] considers interactions of three kinds of cells: Secretory Th cells, responder Th cells and Treg cells. All three cell types express IL-2R molecules on the cell surface. Responder Th cells and Treg cells express IL-2R homogeneously at the cell surface, Treg cells at higher levels than responder Th cells. IL-2 signaling leads to the expression of the oc subunit of the IL-2 receptor that is required for high-affinity IL-2 binding in both responder Th cells and Treg cells. Hence both cell types enhance their rate of IL-2R expression (v) upon IL-2 uptake, which we model, following Busse et al. [4], by a Hill equation with a moderate Hill coefficient of 3: Here, v0 and v1 are the basal and the IL-2 induced rates of IL-2R expression, K is the half-saturation constant, and C(t) is the number of IL-2/IL-2R complexes, which is a dynamic variable of the model (Table 1). For details and the full model see 81 Text.

Supporting Information

Alternative Visualization of simulations shown in Fig 3D—3E, in terms of isosurfaces of the IL-2 concentration (see Section Materials and Methods).

(A) Alternative Visualization of simulations shown in Fig 4C and 4D, in terms of isosurfaces of the IL-2 concentration (see Section Materials and Methods). (BC) Time course of receptor number and IL-2 concentration at the cell surface, see Fig 3E and 3F. (PDF)

Comparison of cytokine traveling distances (Fig 5D) with earlier analytical results by Coppey et. al. [51]. According to their theory, the probability of a ligand binding in this limited domain amounts to 64%, therefore we normalized the results of the 3D-simulations to this amount for the comparison to the theoretical probability distribution. Our simulations comply best with Coppey et al. in the center of the distribution, i.e. in the Vicinity of the maximum. The larger somewhat number of ligand trapping points at the outer boundary of the simulated domain can be explained by the limited number of simulated cells. The smaller number close to the secretion source stems from the difference between an analytical model of spatially homogenized receptors and our discrete model with receptors only on the surface of individual cells.

Acknowledgments

We thank Luca Mariani, Dorothea Busse, Uwe Benary, Rolf Rannacher and Susanne Kromker for discussions, Thomas Carraro for his numerical support and expertise, Lisa Kolb for the rendering of isosurfaces Via topological methods, Marcus Schaber for the volume rendering and Alexander Skupin for critical reading of the manuscript.

Author Contributions

Performed the experiments: KT DG. Analyzed the data: KT DG. Contributed reagents/ materials/ analysis tools: KT DG. Wrote the paper: KT DG EF TH.