PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of interfaceThe Royal Society PublishingInterfaceAboutBrowse by SubjectAlertsFree Trial
 
J R Soc Interface. 2017 July; 14(132): 20170095.
Published online 2017 July 12. doi:  10.1098/rsif.2017.0095
PMCID: PMC5550964

Modelling the host–pathogen interactions of macrophages and Candida albicans using Game Theory and dynamic optimization

Abstract

The release of fungal cells following macrophage phagocytosis, called non-lytic expulsion, is reported for several fungal pathogens. On one hand, non-lytic expulsion may benefit the fungus in escaping the microbicidal environment of the phagosome. On the other hand, the macrophage could profit in terms of avoiding its own lysis and being able to undergo proliferation. To analyse the causes of non-lytic expulsion and the relevance of macrophage proliferation in the macrophage–Candida albicans interaction, we employ Evolutionary Game Theory and dynamic optimization in a sequential manner. We establish a game-theoretical model describing the different strategies of the two players after phagocytosis. Depending on the parameter values, we find four different Nash equilibria and determine the influence of the systems state of the host upon the game. As our Nash equilibria are a direct consequence of the model parameterization, we can depict several biological scenarios. A parameter region, where the host response is robust against the fungal infection, is determined. We further apply dynamic optimization to analyse whether macrophage mitosis is relevant in the host–pathogen interaction of macrophages and C. albicans. For this, we study the population dynamics of the macrophage–C. albicans interactions and the corresponding optimal controls for the macrophages, indicating the best macrophage strategy of switching from proliferation to attacking fungal cells.

Keywords: game theory, Candida albicans, dynamic optimization, non-lytic expulsion, macrophage replication, pathogenic fungi

1. Introduction

Candida albicans is one of the most frequent opportunistic pathogens of humans, and can infect many body sites [13]. In response to environmental challenges, C. albicans can change its morphology from yeast to hyphal growth and back [35]. Hyphae formation and the invasion of host tissue are seen as determinants of the shift from a commensal to an invasive pathogen [6]. In immunocompetent hosts, the defence against C. albicans mainly relies on the innate immune system, especially neutrophils and macrophages [79]. Candida albicans, on the other hand, has developed a plethora of response and evasion mechanisms to circumvent recognition by the immune system and to control, evade and interfere with the host immune reactions [3,7]. An overview of the immunological cross-talk between C. albicans and the innate immune system can be found in Dühring et al. [10]. In this article, we focus on the interaction of C. albicans and macrophages.

Macrophages are dynamic immune cells that are disseminated in several tissues of the human host [11]. The number of macrophages at an infection site is mainly influenced by invading monocytes that differentiate into macrophages. The proliferation of mature macrophages is rather slow and rare but is observed repeatedly for different types of macrophages, e.g. human alveolar macrophages [1214], murine hepatic macrophages (Kupffer cells) [15,16], murine peritoneal macrophages [17,18] and the J774.1 macrophage cell line [19]. It is also of importance in chronic inflammation [13,14,18]. Being a special type of phagocyte, macrophages can phagocytose and eliminate the round yeast cells and relatively short hyphae of C. albicans and thereby limit the fungal burden early upon infection [8,11,20]. Inside the macrophage phagosome, some C. albicans cells are killed through oxidative and non-oxidative microbicidal mechanisms such as by diminution of the pH [2,8,11]. However, C. albicans hyphae and yeast cells can manipulate the phagosomes from the inside by raising the pH [9,21]. This leads to germination of C. albicans yeast cells and hyphae formation within the macrophage [7,9,22]. While the C. albicans yeast cells fail to damage the macrophage, the elongating C. albicans filaments can lead to macrophage death by lysis, allowing the fungus to escape [5,22,23]. Hence, macrophages and C. albicans can have a significant cytotoxic effect on each other [8]. Yet, the response of C. albicans to the phagocytosis by macrophages depends on the genomic background of the fungus. Distinct karyotypes (b and c) of C. albicans isolates differ in their ability to resist intracellular killing. The c karyotype isolates are reported to be more resistant to intracellular killing and to escape from within macrophages by killing the immune cell. This c karyotype is thus considered to behave in a more virulent (or aggressive) way than the b karyotype [2,24,25].

There also exists an alternative escape of hyphae-producing C. albicans cells from within macrophages, called non-lytic expulsion (also termed exocytosis, extrusion or vomocytosis). Shortly after a phagocytosed C. albicans cell starts to form hyphae, it is released into the medium by the macrophage. After non-lytic expulsion both C. albicans and the macrophage are capable of surviving. The host macrophage and the expelled C. albicans cell appear morphologically normal. The macrophage continues to undergo mitosis and the C. albicans hypha extends at normal rates [2629]. Non-lytic expulsion of fungal cells following phagocytosis by macrophages was first described in Cryptococcus neoformans [2931]. This process is also observed in Candida krusei [29,32,33]. Even though the frequency of non-lytic expulsion can be low (in C. albicans less than 1%), the event is observed under various experimental conditions [26,31,34,35]. Furthermore, Nicola et al. [34] showed in C. neoformans that non-lytic expulsion is not an in vitro artefact, but occurred in vivo in murine models with a frequency that is assumed to be much higher than that observed in vitro [29,33,34].

This gives rise to several questions: under which conditions is non-lytic expulsion beneficial to the macrophage and/or the fungus? On one hand, the pathogen is able to escape the hostile environment of the macrophage phagosome. On the other hand, the macrophage is able to avoid potential lysis. Another potential advantage of non-lytic expulsion to the macrophage is the possibility of undergoing mitosis, which can fail in macrophages with phagocytosed C. albicans cells [19,26]. This in turn leads to the question of under which conditions and to what proportions of the macrophage population is it beneficial to the macrophages to undergo mitosis? And further, is it better for macrophages to try to phagocytose and kill C. albicans cells straight away or to proliferate first and switch to phagocytosis at a later time point, when (locally) present in larger numbers? Both dynamic regimes have their advantages. Phagocytosing straight away reduces the number of C. albicans cells earlier upon infection. Proliferating first increases the number of macrophages, which provides an advantage later. To analyse the causes of non-lytic expulsion and the relevance of macrophage proliferation in the macrophage–C. albicans interaction, we employ two complementary approaches of mathematical modelling in a sequential manner: Evolutionary Game Theory and dynamic optimization.

Optimality principles are often used to study and explain biological processes [36,37]. Originating in engineering, dynamic optimization has been used successfully to find optimal regimes in several biological systems including infection processes [3842], protein assembly [43], metabolic pathways [44] and to optimize medical applications such as the treatment of cancer [45] or diabetes [46]. In the theoretical description and modelling of host–parasite interactions, Evolutionary Game Theory has turned out to be a very useful tool [4751]. Game Theory is used, in particular, to describe diverse C. albicans interactions [49,52,53]. Hummert et al. [49] studied the optimal survival strategy of C. albicans cells after phagocytosis by a macrophage. In their setting, C. albicans cells play against each other, while the macrophage is considered as a constant environment. Tyc et al. [52] analysed the coexistence of yeast and hyphal forms in a C. albicans population. In a later study by Tyc et al. [53], the colonization dynamics of C. albicans cells expressing different levels of EFG1 in response to the host immune status are presented. In contrast with these previous works, we consider the host (more specifically, macrophages) as an active player and not solely as an environment.

In §2.1, we establish a game-theoretical model describing the different strategies of macrophages and C. albicans after phagocytosis. Depending on the parameter values, we determine the Nash equilibria (solutions of the game) and analyse the influence of the systems state of the host upon the game (see §3.1). As our Nash equilibria are a direct consequence of the model parameterization, we can depict several biological scenarios. We further determine a parameter region, where the host response is robust against this fungal infection. In §2.2, we apply dynamic optimization to analyse whether macrophage mitosis is relevant in the host–pathogen interaction of macrophages and C. albicans. In this way, we attempt to clarify what regime (phagocytosing or proliferating first) the macrophage population should apply. For this, we determine the population dynamics and corresponding optimal controls indicating the best macrophage strategy of switching from proliferation to attacking fungal cells by macrophages in §3.2.

2. Material and methods

2.1. Characterization of the game

In Evolutionary Game Theory, both the host and pathogen are considered as evolutionary antagonistic players which can show different strategies to maximize their fitness. These strategies can be, for example, cellular traits like the up- and down-regulation of metabolic pathways, expression of virulence factors or generation of different splice variants [50]. A change in strategy can occur by mutation, epigenetic modifications, stochastic gene expression or due to the immunological cross-talk of the players. Each player's fitness is quantified as the net payoff of costs and benefits. Although it is often difficult to quantify the payoff, it is taken into account that the payoff for each player not only depends on its own strategy but also on that of the antagonistic counterpart(s). The solutions of the game are called Nash equilibria. Intuitively, the classical and frequently used concept of the Nash equilibrium is a situation in which neither of the players has an incentive to change strategy unilaterally [47,54].

Game Theory can even be applied if no changes between strategies occur, and should then be interpreted in view of population games. A population may consist of subpopulations determined, for example, by different alleles or karyotypes. These subpopulations can be characterized by different strategies. The final strategy observed in the pure Nash equilibrium arises because one subpopulation outcompetes the other. In the case of a mixed Nash equilibrium, all subpopulations coexist in oscillatory or stationary ways.

In this study, we followed the game-theoretical framework proposed by Renaud & De Meeüs [48] and adapted it to the specific situation after a macrophage (player I) has phagocytosed a C. albicans cell (player II). During this confrontation, each of the two players has two strategies (figure 1a). The macrophage can either ‘release’ the C. albicans cell (non-lytic expulsion) or ‘attack’ and try to kill the pathogen.

Figure 1.
Characterization of the game after phagocytosis: each player's strategies and possible outcomes (a); the payoff functions in terms of fitness F for the C. albicans cell (orange) and the macrophage (blue) (b).

The C. albicans cell can either be ‘less aggressive’ (e.g. karyotype b as described in §1) or ‘aggressive’ (e.g. karyotype c as described in §1). Hereafter, we assume that a macrophage is always able to win a confrontation with a ‘less aggressive’ C. albicans cell. This is not necessarily the case when a macrophage plays against an ‘aggressive’ C. albicans cell.

Pairing the different strategies leads to four distinct types of confrontation: ‘release’/‘less aggressive’, ‘release’/‘aggressive’, ‘attack’/‘less aggressive’ and ‘attack’/‘aggressive’. Each player's payoff for the different confrontations is recorded in terms of fitness F of individuals (figure 1b). The maximum fitness an individual can obtain in a population is denoted by Fmax. The macrophage's payoffs are referred to as FM, while the C. albicans' payoffs are referred to as FC. For our considered confrontations this gives the following pattern.

‘attack’/‘aggressive’: The payoff of the macrophage is

equation image
2.1

The payoff of the C. albicans cell is

equation image
2.2

Only one of the two players can survive in this setting. The probability that the macrophage wins the confrontation is denoted by PM. Hence 0 < PM < 1, with (1 − PM) being the probability that the fungus wins the confrontation. Both players have investment costs to play their strategy, denoted by IM and IC. These could be (but are not limited to) the costs for morphological adaptation, e.g. hyphae development. Each player also has to invest in resisting its opponent's attack. These resistance costs are given by RM and RC.

‘attack’/‘less aggressive’: The macrophage's payoff is

equation image
2.3

In this confrontation, the fungus is eliminated. Its payoff is therefore FC = 0.

‘release’/‘aggressive’: Here, the payoff of the macrophage is

equation image
2.4

The payoff of C. albicans is

equation image
2.5

‘release’/‘less aggressive’: In this confrontation, the payoff of the macrophage is

equation image
2.6

while the payoff of the pathogen is

equation image
2.7

As the ‘aggressive’ C. albicans cell is much more virulent and resistant to the attacks of the macrophages, we assume that the investment costs IC,2 are higher than the investment costs IC,1. Hence, IC,1 < IC,2.

For macrophages to release a C. albicans cell can be more or less severe depending on the constitution of the host and the severity of the fungal attack. We therefore assume that the fitness consequences of the macrophage for ‘releasing’ a C. albicans cell are directly linked to the overall conditions of the host. This is implemented by the systemic costs S1 (‘releasing’ a ‘less aggressive’ fungal cell) and S2 (‘releasing’ an ‘aggressive’ fungal cell). These systemic costs represent the costs of the host (and not just the macrophage). This enables us to investigate the outcome of the game (i.e. the location of the Nash equilibria) in the light of the overall host conditions. Taking into account the overall systemic consequences also allows us to investigate the game-theoretical conflicts between the cell level and organism level, e.g. the seemingly erroneous expulsion of an ‘aggressive’ C. albicans cell by a macrophage to escape lysis. We assume the systemic costs S1 and S2 to be low when the host benefits from ‘releasing’ the fungus. Lewis et al. [19] postulated that non-lytic expulsion might enable macrophage mitosis to proceed normally. As macrophage replication can fail while a C. albicans cell is ingested inside a macrophage, the proliferating macrophage has an incentive to release the fungal cell. Hence, in a scenario where macrophage proliferation is beneficial to the host the systemic costs S1 and S2 are assumed to be low.

The resulting Nash equilibria of our game are given in §3.1. In §2.2 and §3.2, we investigate under which conditions proliferation is relevant for the macrophage–C. albicans interaction.

2.2. The dynamic optimization model

In contrast with Evolutionary Game Theory, dynamic optimization tries to identify a time optimal control in biological systems. To predict an optimal regime of phagocytosis and proliferation, we started by setting up a first-order ordinary differential equation model to simulate the complex dynamics of the host–pathogen interactions (figure 2). We further defined constraints and an objective function to perform dynamic optimization.

Figure 2.
Candida albicans versus mammalian macrophages and the parameters of the differential equation model. C. albicans (orange) can proliferate and switch between its yeast and hyphal forms. Macrophages (blue) can undergo mitosis or phagocytose the fungus. ...

As macrophages are unable to replicate while they are phagocytosing, we are interested in an optimal strategy of the macrophage population with respect to mitosis versus phagocytosis. Hence, we are looking for a control u(t) [set membership] [0, 1] that describes the optimal proportion of the macrophage population undergoing phagocytosis at a given point in time t.

In our model, the macrophage population is denoted by m. The symbol y stands for the C. albicans yeast cell population, h for the C. albicans hyphal cell population and p for the phagocytosed C. albicans cell population. For all populations, our model considers countable cells instead of biomass. In this way, we can use experimental data from the literature, where relevant rates are given as cells per hour (table 1). The differential equation system describing the considered macrophage–C. albicans interactions is formulated as follows:

equation image
2.8
equation image
2.9
equation image
2.10
equation image
2.11

The proportion of proliferating macrophages is 1 − u(t). The macrophage proliferation rate is given by α (table 1). We assume that the number of fungal cells a macrophage is able to ingest is limited to a (table 1). Phagocytosed C. albicans cells kill macrophages at a rate ν and are killed at a rate λ. They are released via non-lytic expulsion at a rate μ. As a consequence, ν(p/a) macrophages are killed in each time step. The yeast cells proliferate at a rate δ and switch to the hyphal growth form at a rate ε. They are phagocytosed by macrophages at a rate β. Candida albicans hyphal cells are phagocytosed at rate γ. We here assume that all phagocytosed C. albicans cells form hyphae inside macrophages so that there are no cells remaining in yeast form. Those C. albicans cells that escaped from the phagosome by killing the macrophage are therefore added to the hyphae population. As a reasonable reduction in complexity, our model does not consider hyphal elongation. However, the qualitative behaviour of our results should not be affected by this, as macrophages are only able to phagocytose relatively short hyphae.

Table 1.
Summary of derived mean values for the parameter distribution of our optimization problem.

As the fungal burden should be minimal at all times, we minimize the integral of all C. albicans cells outside macrophages

equation image
2.12

subject to the differential equation system (equations (2.8)–(2.11)) and the following inequality constraints:

equation image
2.13
equation image
2.14
equation image
2.15
equation image
2.16

3. Results

3.1. Solution of the game

In this section, we determine the Nash equilibria for our game of the interaction between a macrophage (player I) and a C. albicans cell (player II).

The payoff matrix for the macrophage is

equation image
3.1

The payoff matrix for the fungus is

equation image
3.2

If the macrophage plays its ‘attack’ strategy, the C. albicans cell has to adopt the ‘aggressive’ strategy, whatever the average expected value of (1 − PM)(FC,maxIC,2RC) (equation (2.2)), as the fungus is otherwise eliminated. If the macrophage plays its ‘release’ strategy instead, the fungus will play its ‘less aggressive’ strategy as IC,1 < IC,2.

Depending on the values of the parameters IM,1, IM,2, S1, S2 and the probability PM of winning against the fungus the macrophage chooses either ‘attack’ or ‘release’.

There are four different cases of solutions to the game (table 2). In case 1 and case 3, we find the strategy pair ‘attack’/‘aggressive’ as a pure Nash equilibrium, whereas in case 2 and case 3 non-lytic expulsion occurs as a pure Nash equilibrium. In case 4, we find a mixed Nash equilibrium, where each strategy is played with a certain probability so that ‘attacking’/‘releasing’ and ‘aggressive’/‘less aggressive’ cells on both sides coexist with certain frequencies in the population.

Figure 3.
Nash equilibria of the game. Depicted are the four described cases with their corresponding Nash equilibria: ‘attack’/‘aggressive’ (a), ‘release’/‘less aggressive’ (b), ‘attack’/‘aggressive’ ...

Table 2.
Solutions of the game.

Referring to the payoff matrix A and B in equation (3.1) and equation (3.2), we can determine the mean expected payoff EM for the macrophage as a player and EC for the C. albicans player. The mean expected payoff for C. albicans is

equation image
3.3

The mean expected payoff EC depends on the macrophage's tendency to ‘release’ (k) or to ‘attack’ (1 − k) the C. albicans cell. It is, however, independent of the systemic costs of the host (figure 4). The mean expected payoff for the macrophage is

equation image
3.4

Here, d stands for the probability that an interacting C. albicans cell plays the ‘less aggressive’ strategy, whereas (1 − d) stands for the probability of playing the ‘aggressive’ strategy. The systemic costs S1 and S2 exclusively affect EM but not EC.

Figure 4.
The mean expected payoff for C. albicans EC is independent of the systemic costs of the host (orange line). The mean expected payoff of macrophages EM depends on the temptation of C. albicans cells to play the ‘less aggressive’ strategy ...

For EM > EC, the mean expected payoff of a macrophage is superior to that of a C. albicans cell (figure 4). If C. albicans cells exclusively play the ‘less aggressive’ strategy (d = 1), EM,1 falls monotonically in S1/S2. However, if interacting C. albicans cells exclusively play the ‘aggressive’ strategy, EM > EC is only true if S1[dbl greater-than sign]S2 (see EM,2 in figure 4). Let Q1 be defined as the point of intersection of EM,1 and EM,2. Further, let Q2 be the point of intersection of EM,1 with EC and Q3 be the point of intersection of EM,2 with EC. As long as the value of the mean expected payoff of the macrophage in point Q1 is bigger than the value of EC a parameter region can be determined, where the immune response is robust and independent of C. albicans’ temptation to play the ‘aggressive’ strategy (0 ≤ d ≤ 1). Simultaneously, Q2 and Q3 indicate the maximum deflection of the host's state, where the immune defence can be considered as robust. Since k determines the value of EC, a maximization of the area of robust solutions can be achieved if the macrophage is preferably playing the (pure) ‘attack’ strategy.

3.2. Solution of the dynamic optimization problem

Our optimization problem of the interaction of macrophages and C. albicans in §2.2 consists of a time-dependent, continuous control u(t) and state variables (cell populations). To solve this dynamic optimization problem we used a quasi-sequential approach with an extension to handle approximation errors and moving finite elements [71], as in previous works [43,72]. In this gradient-based approach the optimization is repeated 100 times for each parameter set with random initializations of the control variable to avoid ending in local optima.

To determine physiologically relevant model parameters, we used experimental data taken from the literature (table 1). Based on these parameter values, we determined the optimal replication strategy for macrophages with varying multiplicity of infection (MOI) and macrophage replication rates α.

Our initial scenario assumes a standard macrophage replication rate α = 0.0176 and an MOI of 1 : 1 (C. albicans cells : macrophages). For this parameter setting, the optimal strategy of macrophages involves only phagocytosis and no replication (see control u(t) in figure 5a). The population dynamics reveal that macrophages are sufficient to control the number of C. albicans cells. During a simulated period of 24 h, the macrophage population drops only a little due to lysis. Since the replication rate of macrophages is rather low (α = 0.01757) in this setting, we repeated the optimization for the maximal macrophage replication rate reported in the literature (α = 0.2660; figure 5b). This scenario leads to a different optimal strategy u(t) starting with a phase of exclusive phagocytosis. This is followed by a phase of replication and phagocytosis and ends with a phase of exclusive phagocytosis. Comparing both scenarios, the number of C. albicans cells over the whole simulated time span is nearly identical. This indicates that, for the macrophages, replication (even with high rates) is of small advantage when the MOI is balanced. To see whether the amount of C. albicans cells rather than the replication rate influences the macrophages' strategy, we changed the MOI to 7 : 1, assuming again a standard replication rate (α = 0.0176; figure 5c). The resulting dynamic shows that macrophages are not able to control the invading C. albicans cells. The optimal strategy of macrophages switches immediately from exclusive phagocytosis to almost complete replication to escape lysis by C. albicans.

Figure 5.
Dynamics of the population sizes of C. albicans cells and macrophages with standard and maximal macrophage replication rates and the corresponding optimal controls for different MOIs.

To quantify the influence of parameters on the optimal strategy of macrophages to combat C. albicans cells, we performed a parameter sensitivity analysis with varying MOIs (1 : 1 and 3 : 1) and macrophage replication rates (standard and high). As expected, the replication rate α has the highest impact on the amount of replication for the tested MOIs of 1 : 1 and 3 : 1 and replication rates α. This is followed by the C. albicans filamentation rate, ε, as well as the macrophage lysis rate by C. albicans, ν (electronic supplementary material, figures S1 and S2).

The corresponding optimal controls u(t) show for almost every parameter set no replication under the scenario MOI 1 : 1, where macrophages can efficiently control the C. albicans population (electronic supplementary material, figure S3). For the higher MOI of 3:1, more parameter sets lead to optimal strategies including phases of replication (electronic supplementary material, figure S4). This is because C. albicans cells exceed the manageable amount for macrophages in many parameter sets (electronic supplementary material, figures S1 and S2).

A closer inspection of both scenarios (MOI of 1 : 1, 3 : 1) reveals that replication is only effective if the replication rate exceeds a certain magnitude (electronic supplementary material, figures S1 and S2). The majority of optimal replication strategies have a similar time course to that seen in figure 5b. Therefore, we can conclude that an early phase of replication during the incubation of macrophages and C. albicans cells is only of advantage to the macrophage if the replication rate is very high. For more results on the sensitivity analysis see the electronic supplementary material.

4. Discussion

In this work, we used mathematical modelling to study non-lytic expulsion of C. albicans from within macrophages and analysed the relevance of macrophage replication in the macrophage–C. albicans interactions.

Applying dynamic optimization, we studied the population dynamics of the macrophage–C. albicans interactions and the corresponding optimal controls for the macrophages. To examine the conditions under which macrophage mitosis is beneficial to the macrophage population (and the host), we explored the dynamic regimes of phagocytosis and proliferation of macrophages in the light of an optimal response to fungal burden. For this, we looked into different infection scenarios represented by different MOIs (C. albicans : macrophages) and macrophage replication rates (see §2.2 and §3.2). Hereby, we considered cells instead of biomass for our populations.

Comparing the different scenarios, we found that replication (even with high rates) is of small advantage when the MOI is balanced (1 : 1). Our simulations indicate that, when macrophages are able to control the C. albicans invasion (balanced MOI of 1 : 1), the optimal strategy for macrophages includes almost no replication (figure 5a,b). Even with higher MOIs (3 : 1), replication is only effective if the replication rate exceeds a certain magnitude (electronic supplementary material, figures S2 and S4). When applying an MOI of 7 : 1, macrophages are no longer able to control the invading C. albicans cells regardless of the proliferation policy. In a healthy host with a standard macrophage replication rate (table 1), the great majority of macrophages solely phagocytose the fungal cells. Only with high macrophage replication rates, as, for example, seen in chronic inflammation, does the system's behaviour change and macrophage replication occurs in addition to phagocytosis.

We conclude that an early phase of macrophage replication during the incubation of macrophages and C. albicans cells is only of advantage to the macrophage if the replication rate is very high. In healthy individuals, macrophage replication is therefore of minor importance in this interaction.

We further analysed the conditions which render non-lytic expulsion beneficial for a given player (the macrophage and the fungus). For the fungus, non-lytic expulsion is always beneficial as the fungus escapes the hostile environment of the macrophage. For the macrophage, the situation is more nuanced. Using Evolutionary Game Theory, we found four different Nash equilibria depending on the investment costs of the macrophage, the systemic costs of the host and the probability that a macrophage is able to kill the fungus. From our game-theoretical model, we derived that non-lytic expulsion can occur as a pure Nash equilibrium (see cases 2 and 3 in table 2) as well as a mixed Nash equilibrium (see case 4 in table 2).

As our Nash equilibria are a direct consequence of the model parameterization, we can depict several biological scenarios. By adjusting the value of the systemic costs S1 for ‘releasing’ a ‘less aggressive’ fungal cell, it is possible to make a distinction between host conditions (i.e. healthy or weakened). In a healthy host, macrophage replication is of minor importance, as shown in our population dynamics in §3.2. Hence, the systemic costs S1 are high as the focus is on attacking the C. albicans cells and fungal clearance. The systemic costs S2 for ‘releasing’ an ‘aggressive’ fungal cell, on the other hand, are directly linked to the probability PM of the macrophage winning the confrontation with the fungus. With decreasing chances of the macrophage surviving a fungal ‘attack’, the systemic costs S2 also decrease. The systemic costs S2 can thus function as a measure of the severity of the fungal attack. Note that there can exist a lag phase in the linkage of PM and S2. During the advance of the C. albicans infection, the fungal cells become more difficult to handle for the macrophage because of hyphae formation. This can lead to a conflict between the host's cell and the organism level, especially when the macrophage's probability of losing to the fungus is increasing but the systemic costs S2 are still high. Our model predicts that, in this setting, the macrophage would still try to attack the fungus despite its decreasing chance of surviving such an attack. Only when the systemic costs S2 decrease to a low value does the strategy of the macrophage switch to non-lytic expulsion (see cases 2 and 4 in table 2) to avoid macrophage lysis. This loss in macrophage competence has to be compensated by other immune cells such as neutrophils (not part of the model presented here).

The scenario of an immune-competent host with a moderate fungal infection is best described by high systemic costs S1 while at the same time the systemic costs S2 and/or the probability PM of winning against the fungus are high. In this parameter setting, we find the pure Nash equilibrium ‘attack’/‘aggressive’ (see case 1 in table 2). The macrophages’ focus is on attacking the fungal cells. A deviation from this strategy is blocked by the prohibitively high systemic costs of this scenario. ‘Attacking’ the fungus also maximizes the robustness of the immune response to C. albicans, as shown in §3.1. This leads to the immune state becoming more tolerant of perturbations by C. albicans.

The model also provides a scenario with a healthy host suffering from an advanced infection (severe fungal attack). The systemic costs S1 are still high but either the systemic costs S2 and/or the chance of ‘winning’ an attack against the C. albicans player are low (see case 4 in table 2). In this scenario, we find the mixed Nash equilibrium, which includes the release of fungal cells. It is noteworthy that only in this model scenario is an ‘aggressive’ fungal cell released instead of being killed.

In a weakened host (low S1), independent of the severity of the fungal infection (S2), we find non-lytic expulsion as a pure Nash equilibrium (see cases 2 and 3 in table 2). This is a direct consequence of the macrophage's costs for ‘attacking’ a ‘less aggressive’ fungus (IM,1 and S1) being higher than the costs for ‘releasing’ this type of fungus (IM,2) (see cases 2 and 3 in table 2).

From this we conclude that, for non-lytic expulsion to be beneficial to the macrophage, either the fungal ‘attack’ must be severe with a low probability of the macrophages surviving (as in the case of the mixed Nash equilibrium) or the host needs to have a strong incentive to release the fungus, i.e. to undergo mitosis. But, as mitosis is of minor importance to the macrophage–C. albicans interaction in a healthy host, the host conditions need to be in a state where the macrophage replication rate is strongly elevated, e.g. by chronic inflammation or by using a fast replicating macrophage cell line such as J774.1. To simulate a severe fungal ‘attack’ in experiments, it is not sufficient to only have a high fungal burden but rather to use an aggressive karyotype of C. albicans. As our simulations indicate, macrophages are not able to control the fungal burden with MOIs of 7 : 1. Instead, one would expect to see fungal outgrowth with a high occurrence of macrophage lysis. It would be of interest to further study the situation of chronic inflammation and/or severe fungal infection in wet lab experiments, as non-lytic expulsion should appear more often under those conditions. Moreover, both the aggressive and less aggressive karyotypes of C. albicans could be used in co-infection experiments with slow and fast replicating macrophages to simulate the population game. Once the molecular mechanisms of non-lytic expulsion are better understood, overexpression and knockouts could be used to end up in specific Nash equilibria.

Our findings using Game Theory and dynamic optimization give a holistic perspective on fungal infection processes and are consistent with experimental observations [26]. They further explain why the frequency of non-lytic expulsion can be very low in experimental studies. Both our modelling approaches can be used for further studies of other fungal pathogens such as C. neoformans and C. krusei. For this, our models would need to be adapted to the specific situations of those infections and to the different lifestyles of those pathogens. An extension of our study could be the activation of macrophages through chemokines and cytokines and the recruitment of other immune cells like neutrophils. The question would then be whether it is more beneficial to invest in the attack of fungal cells or in recruitment. Also interesting for further extensions to our model are other macrophage strategies such as macrophage extracellular trap-like structures [11] and the competition within the C. albicans population.

Supplementary Material

Supplementary Material:

Acknowledgements

The authors thank Martin Bartl and Pu Li for stimulating discussions and providing the toolbox for solving the optimization problem.

Data accessibility

The datasets supporting this article are included in the main paper and/or have been uploaded as part of the electronic supplementary material.

Authors' contributions

Setting up the mathematical models: S.D., S.S. Implementing the optimization problem in the toolbox: J.E. Conception and design of the investigation and work: S.D., J.E., S.G., S.S. Drafting the manuscript: all. Revising it critically for important intellectual content and final approval of the version to be published: all. Agreement to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved: all.

Competing interests

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Funding

S.D., T.D. and S.S. acknowledge support by the Deutsche Forschungsgemeinschaft (DFG) CRC/Transregio 124 ‘Pathogenic fungi and their human host: Networks of interaction’ sub project B1 and J.E. by B2. C.K. acknowledges support by the Cluster of Excellence ‘Inflammation at Interfaces’ (DFG support code EXC306).

References

1. Wilson D, et al. 2009. Identifying infection-associated genes of Candida albicans in the postgenomic era. FEMS Yeast Res. 9, 688–700. (doi:10.1111/j.1567-1364.2009.00524.x) [PubMed]
2. Cheng S-C, Joosten LA, Kullberg B-J, Netea MG 2012. Interplay between Candida albicans and the mammalian innate host defense. Infect. Immun. 80, 1304–1313. (doi:10.1128/IAI.06146-11) [PMC free article] [PubMed]
3. Luo S, Skerka C, Kurzai O, Zipfel PF 2013. Complement and innate immune evasion strategies of the human pathogenic fungus Candida albicans. Mol. Immunol. 56, 161–169. (doi:10.1016/j.molimm.2013.05.218) [PubMed]
4. Yan L, Yang C, Tang J 2013. Disruption of the intestinal mucosal barrier in Candida albicans infections. Microbiol. Res. 168, 389–395. (doi:10.1016/j.micres.2013.02.008) [PubMed]
5. Uwamahoro N, et al. 2014. The pathogen Candida albicans hijacks pyroptosis for escape from macrophages. mBio 5, e00003–14. (doi:10.1128/mBio.00003-14) [PMC free article] [PubMed]
6. Mech F, Wilson D, Lehnert T, Hube B, Figge MT 2014. Epithelial invasion outcompetes hypha development during Candida albicans infection as revealed by an image-based systems biology approach. Cytometry A 85, 126–139. (doi:10.1002/cyto.a.22418) [PubMed]
7. Lopez CM. 2013. The roles of Candida albicans Gpm1p and Tef1p in immune evasion and tissue invasion of the human host. PhD thesis, Friedrich-Schiller-Universität Jena, Jena, Germany.
8. Krysan DJ, Sutterwala FS, Wellington M 2014. Catching fire: Candida albicans, macrophages, and pyroptosis. PLoS Pathog. 10, e1004139 (doi:10.1371/journal.ppat.1004139) [PMC free article] [PubMed]
9. Vylkova S, Lorenz MC 2014. Modulation of phagosomal pH by Candida albicans promotes hyphal morphogenesis and requires Stp2p, a regulator of amino acid transport. PLoS Pathog. 10, e1003995 (doi:10.1371/journal.ppat.1003995) [PMC free article] [PubMed]
10. Dühring S, Germerodt S, Skerka C, Zipfel PF, Dandekar T, Schuster S 2015. Host-pathogen interactions between the human innate immune system and Candida albicans—understanding and modeling defense and evasion strategies. Front. Microbiol. 6, 625 (doi:10.3389/fmicb.2015.00625) [PMC free article] [PubMed]
11. Liu P, et al. 2014. Escherichia coli and Candida albicans induced macrophage extracellular trap-like structures with limited microbicidal activity. PLoS ONE 9, e90042 (doi:10.1371/journal.pone.0090042) [PMC free article] [PubMed]
12. Evans MJ, Cabral LJ, Stephens RJ, Freeman G 1973. Cell division of alveolar macrophages in rat lung following exposure to NO2. Am. J. Pathol. 70, 199–208. [PubMed]
13. Bitterman PB, Saltzman LE, Adelberg S, Ferrans VJ, Crystal RG 1984. Alveolar macrophage replication. One mechanism for the expansion of the mononuclear phagocyte population in the chronically inflamed lung. J. Clin. Invest. 74, 460–469. (doi:10.1172/JCI111443) [PMC free article] [PubMed]
14. Barbers RG, Evans MJ, Gong H Jr, Tashkin DP 1991. Enhanced alveolar monocytic phagocyte (macrophage) proliferation in tobacco and marijuana smokers. Am. J. Respir. Crit. Care Med. 143, 1092–1095. (doi:10.1164/ajrccm/143.5_Pt_1.1092) [PubMed]
15. North RJ. 1969. The mitotic potential of fixed phagocytes in the liver as revealed during the development of cellular immunity. J. Exp. Med. 130, 315–326. (doi:10.1084/jem.130.2.315) [PMC free article] [PubMed]
16. More D, Nelson D 1972. Antigen-induced mitosis in liver macrophages of immunized mice. Experientia 28, 566–567. (doi:10.1007/BF01931880) [PubMed]
17. Khoo K, Mackaness G 1964. Macrophage proliferation in relation to acquired cellular resistance. Aust. J. Exp. Biol. Med. Sci. 42, 707–716. (doi:10.1038/icb.1964.67) [PubMed]
18. Murch A, Papadimitriou J 1981. The kinetics of murine peritoneal macrophage replication. J. Pathol. 133, 177–183. (doi:10.1002/path.1711330302) [PubMed]
19. Lewis LE, Bain JM, Lowes C, Gow NA, Erwig L-P 2012. Candida albicans infection inhibits macrophage cell division and proliferation. Fungal Genet. Biol. 49, 679–680. (doi:10.1016/j.fgb.2012.05.007) [PMC free article] [PubMed]
20. Jacobsen ID, Wilson D, Wächtler B, Brunke S, Naglik JR, Hube B 2012. Candida albicans dimorphism as a therapeutic target. Expert Rev. Anti Infect. Ther. 10, 85–93. (doi:10.1586/eri.11.152) [PubMed]
21. Vylkova S, Carman AJ, Danhof HA, Collette JR, Zhou H, Lorenz MC 2011. The fungal pathogen Candida albicans autoinduces hyphal morphogenesis by raising extracellular pH. MBio 2, e00055 (doi:10.1128/mBio.00055-11) [PMC free article] [PubMed]
22. Wellington M, Koselny K, Sutterwala FS, Krysan DJ 2014. Candida albicans triggers NLRP3-mediated pyroptosis in macrophages. Eukaryot. Cell 13, 329–340. (doi:10.1128/EC.00336-13) [PMC free article] [PubMed]
23. Dementhon K, El-Kirat-Chatel S, Noël T 2012. Development of an in vitro model for the multi-parametric quantification of the cellular interactions between Candida yeasts and phagocytes. PLoS ONE 7, e32621 (doi:10.1371/journal.pone.0032621) [PMC free article] [PubMed]
24. Tavanti A, Campa D, Bertozzi A, Pardini G, Naglik JR, Barale R, Senesi S 2006. Candida albicans isolates with different genomic backgrounds display a differential response to macrophage infection. Microbes Infect. 8, 791–800. (doi:10.1016/j.micinf.2005.09.016) [PubMed]
25. Neglia R, Colombari B, Peppoloni S, Orsi C, Tavanti A, Senesi S, Blasi E 2006. Adaptive response of microglial cells to in vitro infection by Candida albicans isolates with different genomic backgrounds. Microb. Pathog. 41, 251–256. (doi:10.1016/j.micpath.2006.08.001) [PubMed]
26. Bain JM, Lewis LE, Okai B, Quinn J, Gow NA, Erwig L-P 2012. Non-lytic expulsion/exocytosis of Candida albicans from macrophages. Fungal Genet. Biol. 49, 677–678. (doi:10.1016/j.fgb.2012.01.008) [PMC free article] [PubMed]
27. Brunke S, Hube B 2013. Two unlike cousins: Candida albicans and C. glabrata infection strategies. Cell Microbiol. 15, 701–708. (doi:10.1111/cmi.12091) [PMC free article] [PubMed]
28. Whittington A, Gow NA, Hube B 2014. 1 From commensal to pathogen: Candida albicans. In Human fungal pathogens (ed. Kurzai O, editor. ), 2nd edn The Mycota, vol. 12, pp. 3–18. Berlin, Germany: Springer.
29. Erwig LP, Gow NA 2016. Interactions of fungal pathogens with phagocytes. Nat. Rev. Microbiol. 14, 163–176. (doi:10.1038/nrmicro.2015.21) [PubMed]
30. Alvarez M, Casadevall A 2006. Phagosome extrusion and host-cell survival after Cryptococcus neoformans phagocytosis by macrophages. Curr. Biol. 16, 2161–2165. (doi:10.1016/j.cub.2006.09.061) [PubMed]
31. Ma H, Croudace JE, Lammas DA, May RC 2006. Expulsion of live pathogenic yeast by macrophages. Curr. Biol. 16, 2156–2160. (doi:10.1016/j.cub.2006.09.032) [PubMed]
32. García-Rodas R, González-Camacho F, Rodríguez-Tudela JL, Cuenca-Estrella M, Zaragoza O 2011. The interaction between Candida krusei and murine macrophages results in multiple outcomes, including intracellular survival and escape from killing. Infect. Immun. 79, 2136–2144. (doi:10.1128/IAI.00044-11) [PMC free article] [PubMed]
33. Smith LM, May RC 2013. Mechanisms of microbial escape from phagocyte killing. Biochem. Soc. Trans. 41, 475–490. (doi:10.1042/BST20130014) [PubMed]
34. Nicola AM, Robertson EJ, Albuquerque P, da Silveira Derengowski L, Casadevall A 2011. Nonlytic exocytosis of Cryptococcus neoformans from macrophages occurs in vivo and is influenced by phagosomal pH. mBio 2, e00167–11 (doi:10.1128/mBio.00167-11) [PMC free article] [PubMed]
35. Stukes SA, Cohen HW, Casadevall A 2014. Temporal kinetics and quantitative analysis of Cryptococcus neoformans nonlytic exocytosis. Infect. Immun. 82, 2059–2067. (doi:10.1128/IAI.01503-14) [PMC free article] [PubMed]
36. Lardon LA, Merkey BV, Martins S, Dötsch A, Picioreanu C, Kreft J-U, Smets BF 2011. iDynoMiCS: next-generation individual-based modelling of biofilms. Environ. Microbiol. 13, 2416–2434. (doi:10.1111/j.1462-2920.2011.02414.x) [PubMed]
37. Clegg RJ, Dyson RJ, Kreft J-U 2014. Repair rather than segregation of damage is the optimal unicellular aging strategy. BMC Biol. 12, 52 (doi:10.1186/s12915-014-0052-x) [PMC free article] [PubMed]
38. Segel LA. 1980. Mathematical models in molecular and cellular biology. Cambridge, UK: Cambridge University Press.
39. Lenhart S, Workman JT 2007. Optimal control applied to biological models. Boca Raton, FL: CRC Press.
40. Banga JR. 2008. Optimization in computational systems biology. BMC Syst. Biol. 2, 47 (doi:10.1186/1752-0509-2-47) [PMC free article] [PubMed]
41. Brauer F, Castillo-Chavez C, Castillo-Chavez C 2001. Mathematical models in population biology and epidemiology. Berlin, Germany: Springer.
42. Stengel RF, Ghigliazza RM, Kulkarni NV 2002. Optimal enhancement of immune response. Bioinformatics 18, 1227–1235. (doi:10.1093/bioinformatics/18.9.1227) [PubMed]
43. Ewald J, Kötzing M, Bartl M, Kaleta C 2015. Footprints of optimal protein assembly strategies in the operonic structure of prokaryotes. Metabolites 5, 252–269. (doi:10.3390/metabo5020252) [PMC free article] [PubMed]
44. Bartl M, Kötzing M, Schuster S, Li P, Kaleta C 2013. Dynamic optimization identifies optimal programmes for pathway regulation in prokaryotes. Nat. Commun. 4, 2243 (doi:10.1038/ncomms3243) [PubMed]
45. Swan GW. 1990. Role of optimal control theory in cancer chemotherapy. Math. Biosci. 101, 237–284. (doi:10.1016/0025-5564(90)90021-P) [PubMed]
46. Chen C-L, Tsai H-W 2009. Model-based insulin therapy scheduling: a mixed-integer nonlinear dynamic optimization approach. Ind. Eng. Chem. Res. 48, 8595–8604. (doi:10.1021/ie9005673)
47. Hofbauer J, Sigmund K 1998. Evolutionary games and population dynamics. Cambridge, UK: Cambridge University Press.
48. Renaud F, DeMeeüs T 1991. A simple model of host-parasite evolutionary relationships. Parasitism: compromise or conflict? J. Theor. Biol. 152, 319–327. (doi:10.1016/S0022-5193(05)80197-3) [PubMed]
49. Hummert S, Hummert C, Schröter A, Hube B, Schuster S 2010. Game theoretical modelling of survival strategies of Candida albicans inside macrophages. J. Theor. Biol. 264, 312–318. (doi:10.1016/j.jtbi.2010.01.022) [PubMed]
50. Hummert S, Bohl K, Basanta D, Deutsch A, Werner S, Theißen G, Schroeter A, Schuster S 2014. Evolutionary game theory: cells as players. Mol. Biosyst. 10, 3044–3065. (doi:10.1039/C3MB70602H) [PubMed]
51. Lindsay RJ, Kershaw MJ, Pawlowska BJ, Talbot NJ, Gudelj I 2016. Harbouring public good mutants within a pathogen population can increase both fitness and virulence. eLife 5, e18678 (doi:10.7554/eLife.18678) [PMC free article] [PubMed]
52. Tyc KM, Kühn C, Wilson D, Klipp E 2014. Assessing the advantage of morphological changes in Candida albicans: a game theoretical study. Front. Microbiol. 5, 41 (doi:10.3389/fmicb.2014.00041) [PMC free article] [PubMed]
53. Tyc KM, Herwald SE, Hogan JA, Pierce JV, Klipp E, Kumamoto CA 2016. The game theory of Candida albicans colonization dynamics reveals host status-responsive gene expression. BMC Syst. Biol. 10, 20 (doi:10.1186/s12918-016-0268-1) [PMC free article] [PubMed]
54. Barron EN. 2008. Game theory: an introduction, vol. 2 New York, NY: John Wiley & Sons.
55. Ishida K, de Mello JCP, Cortez DAG, Dias Filho BP, Ueda-Nakamura T, Nakamura CV 2006. Influence of tannins from Stryphnodendron adstringens on growth and virulence factors of Candida albicans. J. Antimicrob. Chemother. 58, 942–949. (doi:10.1093/jac/dkl377) [PubMed]
56. Gyanchandani A, Khan Z, Maitra S 2000. Arbortristosides modulate murine peritoneal macrophages for phagocytosis and intracellular killing of Candida albicans. Pharm. Biol. 38, 340–352. (doi:10.1076/phbi.38.5.340.5969)
57. Lewis LE, Bain JM, Lowes C, Gillespie C, Rudkin FM, Gow NA, Erwig L-P 2012. Stage specific assessment of Candida albicans phagocytosis by macrophages identifies cell wall composition and morphogenesis as key determinants. PLoS Pathog. 8, e1002578 (doi:10.1371/journal.ppat.1002578) [PMC free article] [PubMed]
58. Vera A, Olga S, Olga S, Natalia M 2013. Physiological parameters of clinical yeasts growth and isolation of specific antigens. Biochem. Physiol. Open Access 2013, 116 (doi:10.4172/2168-9652.1000116)
59. Baillie GS, Douglas LJ 1998. Iron-limited biofilms of Candida albicans and their susceptibility to amphotericin B. Antimicrob. Agents Chemother. 42, 2146–2149. [PMC free article] [PubMed]
60. Baillie GS, Douglas LJ 1999. Role of dimorphism in the development of Candida albicans biofilms. J. Med. Microbiol. 48, 671–679. (doi:10.1099/00222615-48-7-671) [PubMed]
61. Ozato K, Uesaka I 1974. The role of macrophages in Candida albicans infection in vitro. Jpn. J. Microbiol. 18, 29–35. (doi:10.1111/j.1348-0421.1974.tb00740.x) [PubMed]
62. McKenzie C, Koser U, Lewis L, Bain J, Mora-Montes H, Barker R, Gow N, Erwig L 2010. Contribution of Candida albicans cell wall components to recognition by and escape from murine macrophages. Infect. Immun. 78, 1650–1658. (doi:10.1128/IAI.00001-10) [PMC free article] [PubMed]
63. Ene IV, Cheng S-C, Netea MG, Brown AJ 2013. Growth of Candida albicans cells on the physiologically relevant carbon source lactate affects their recognition and phagocytosis by immune cells. Infect. Immun. 81, 238–248. (doi:10.1128/IAI.01092-12) [PMC free article] [PubMed]
64. Rudkin FM, Bain JM, Walls C, Lewis LE, Gow NA, Erwig LP 2013. Altered dynamics of Candida albicans phagocytosis by macrophages and PMNs when both phagocyte subsets are present. mBio 4, e00810 (doi:10.1128/mBio.00810-13) [PMC free article] [PubMed]
65. Loyola W, Gaziri DA, Gaziri LCJ, Felipe I 2002. Concanavalin A enhances phagocytosis and killing of Candida albicans by mice peritoneal neutrophils and macrophages. FEMS Immunol. Med. Microbiol. 33, 201–208. (doi:10.1111/j.1574-695X.2002.tb00591.x) [PubMed]
66. Prestes-Carneiro LE, Laraya RD, Silva PRC, Moliterno RA, Felipe I, Mathias PC 2006. Long-term effect of early protein malnutrition on growth curve, hematological parameters and macrophage function of rats. J. Nutr. Sci. Vitaminol. 52, 414–420. (doi:10.3177/jnsv.52.414) [PubMed]
67. Moresco T, Gaziri L, Yasumoto Y, Felipe I 2002. Phagocytic and candidacidal activities of macrophages from suckling and adult mice pretreated with concanavalin-A. Med. Mycol. 40, 393–397. (doi:10.1080/mmy.40.4.393.397) [PubMed]
68. Maródi L, Tournay C, Káposzta R, Johnston RB, Moguilevsky N 1998. Augmentation of human macrophage candidacidal capacity by recombinant human myeloperoxidase and granulocytemacrophage colony-stimulating factor. Infect. Immun. 66, 2750–2754. [PMC free article] [PubMed]
69. Redmond H, Hofmann K, Shou J, Leon P, Kelly C, Daly JM 1992. Effects of laparotomy on systemic macrophage function. Surgery 111, 647–655. [PubMed]
70. Brothers KM, Newman ZR, Wheeler RT 2011. Live imaging of disseminated candidiasis in zebrafish reveals role of phagocyte oxidase in limiting filamentous growth. Eukaryot. Cell. 10, 932–944. (doi:10.1128/EC.05005-11) [PMC free article] [PubMed]
71. Bartl M, Li P, Biegler LT 2011. Improvement of state profile accuracy in nonlinear dynamic optimization with the quasi-sequential approach. AIChE J. 57, 2185–2197. (doi:10.1002/aic.12437)
72. de Hijas-Liste GM, Balsa-Canto E, Ewald J, Bartl M, Li P, Banga JR, Kaleta C 2015. Optimal programs of pathway control: dissecting the influence of pathway topology and feedback inhibition on pathway regulation. BMC Bioinformatics 16, 163 (doi:10.1186/s12859-015-0587-z) [PMC free article] [PubMed]

Articles from Journal of the Royal Society Interface are provided here courtesy of The Royal Society