Double public goods games and acid-mediated tumor invasion
January 22, 2015 18 Comments
Although I’ve spent more time thinking about pairwise games, I’ve recently expanded my horizons to more serious considerations of public-goods games. They crop up frequently when we are modeling agents at the cellular level, since interacts are often indirect through production of some sort of common extra-cellular signal. Unlike the trivial to characterize two strategy pairwise games, two strategy public-goods have a more sophisticated range of possible dynamics. However, through a nice trick using the properties of Bernstein polynomials, Archetti (2013,2014) and Peña et al. (2014a) have greatly increased our understanding of the public good, and I will be borrowing heavily from the toolbag and extending it slightly in this post. I will discuss the obvious continuation of this work by considering more than two strategies and several public goods together. Unfortunately, the use of public goods games here — and of evolutionary game theory (EGT) more generally — is not without controversy. This extension is not meant to address the controversy of spatial structure (although for progress on this, see Peña et al., 2014b), but the rigorous qualitative analysis that I’ll use (mostly in a the next post on this project) will allow me to side-step much of the parameter-fitting issues.
Of course, having two public goods games is only interesting if we couple them to each other. In this case, we will have one public good from which everyone benefits, but the second good is anti-correlated in the sense that only those that don’t contribute to the first can benefit from the second. A more general analysis of all possible ways to correlate two public-goods game might be a fun future direction, but at this point it is not clear what other correlations would be useful for modeling; at least in mathematical oncology.
By the way, if you are curious what mathematical oncology research looks like, it is often just scribbles like this emailed back and forth:
I’ll use the rest of this post to guide you through the ideas behind the above sketch, and thus introduce you to the joint project that I am working on with Robert Vander Velde, David Basanta, and Jacob Scott. Treat this as a page from my open research notebook.
I’ll start with the second game. Let f(x) be the relative contribution of using a growth-factor — such as oxygen via aerobic respiration — to Alice’s fitness if x is the ratio of agents in Alice’s interaction neighborhood that are producing the growth-factor to those consuming it. The production of the growth-factor might be indirect — for instance producing VEGF to help vasculate nearby tissue and thus increase the abundance of oxygen — but it is not our goal to model all the potential complications, like time-delays, of this process. We will also assume that producing this growth-factor carries some constant cost c > 0.
Now moving onto the first game. Suppose that living without using the grown-factor produces some contaminant that is toxic to non-cancer cells, and maybe even to cancer cells in high enough quantities. In our aerobic vs. anaerobic respiration example, this might be the acid secreted by cells forfeiting aerobic respiration in favor of glycolysis. This is believed to help tumors achieve acid-mediated invasion (Gatenby & Gillies, 2004, 2007; Gatenby et al., 2006) and is one of the explanation for the Warburg effect in cancer progression. We’ll let m(x) be the relative contribution to Alice’s fitness of the contaminant if the proportion x of agents in Alice’s interaction is not using the growth-factor and thus producing the contaminant. There is no cost for producing the contaminant, except the implicit cost in not receiving the benefit f by functioning without using the growth-factor.
Based on the above considerations, we can write down the payoffs for the four strategies ‘‘ where the fist sign is ‘+‘ if the organism uses the growth-factor and the second sign is ‘+‘ if the organism produces — potentially indirectly — the growth factor. Note that in the case of the first sign, ‘–‘ means to cooperate, since by not using the growth-factor, the agent cooperates in producing the contaminant; while for the second sign, ‘+‘ means to cooperate by producing the growth-factor. Thus the double defector strategy is ‘-+‘ and the universal cooperator is ‘+-‘. Let be the proportion of agents with strategy ‘—‘, as strategy ‘++‘, as ‘+-‘, and as ‘-+‘.
In the above notation, the utilities (or finesses) for the four strategies are:
where the and are the multinomial coefficients. Note that c being positive means that we always have U(–) > U(-+) so the strategy ‘-+‘ will not be present in any stable equilibrium. For now, I am only interested in the equilibrium analysis of this system, so we can set and then factor our replicator equation with and . One of the first components we’ll need is the gradient for q:
where we go from the first to second line by applying the multinomial generalization of Bernstein’s approximation theorem (for example, see Duchoň, 2011; in a future post, I will provide a slightly cleaner probabilistic proof of this by generalizing Levasseur, 1984). Since the dynamic equation for q is given by , we know that — apart from the two trivial equilibria at q = 0 and q = 1 — we have an internal equilibrium when there exists a such that . Now, we can look for the gradient for p:
From the above we can characterize the dynamics of p. We know that — apart from the the two trivial quilibria at p = 0 and p = 1 — we have an internal equilibrium when there exists a such that . Note that if we are looking at q = 0 or q = q*, we have Q(q) simplifying to just be f(q).
Also, just to be pedantic I want to write down the simplified dynamic system of equations explicitly:
Note that for p = 0 we recover Archetti (2013) and for q = 0 or q = 1 we get Archetti (2014). In the other cases, if we want the equilibrium points (p*,q*) then we to solve the following system of equations:
At this point, it is important to do a sanity check by considering linear public goods and . This results in a stable internal equilibrium with and or in the simplex coordinates . From this, we can conclude that if then we will have all cooperators in the growth-factor public good. In our running example, this means that among the aerobic cells everybody will invest in VEGF. If the we will have all cooperators in the contamination public good. In our running example, this means that all the cells will revert to glycolosis. These cases are in agreement with more traditional approaches like Santos et al. (2008) hitchhikes.
The interesting case is that if both conditions are violated then unlike the single public-good (which would go toward all defectors), this anti-correlated pair of public goods will result in a coexistence of all 3 types of strategies. In our running example, we will have a tumor with some cells glycolotic, while others continue aerobic respiration with some of them producing VEGF. I am not sure if this co-existence is frequently observed in clinical cases, so I welcome any references. Note that this result would have been impossible to achieve with two decoupled linear public goods games, so the anti-correlated coupling of games is essential.
Of course, linear public goods are of limited appeal. In the next post on this project I will use the sign pattern technique of Gokhale & Traulsen (2010) and Peña et al. (2014a) to look at the general case of two anti-correlated non-linear public goods games. As far as I can tell, this post offers one of the first discussions of the EGT of multiple interacting public goods. Multiple public goods games have been briefly looked at in economics by Sable & Kling (2001), Ehlers (2001,2004), and Umezawa (2012) but as far as I can tell their interests in multiple public goods game is in the context of matching markets (which is fundamentally an equilibrium selection problem) while I am looking at them as social dilemmas (which is fundamentally playing with the tension of Nash Eq. that is not Pareto optimal). However, dear reader, if you know of existing work on the dynamics of multiple interacting public goods then please mention it in the comments!
Archetti, M. (2013). Evolutionary game theory of growth factor production: implications for tumour heterogeneity and resistance to therapies. British Journal of Cancer, 109(4): 1056-1062.
Archetti, M. (2014). Evolutionary dynamics of the Warburg effect: glycolysis as a collective action problem among cancer cells. Journal of Theoretical Biology, 341: 1-8
Duchoň, M. (2011). A generalized Bernstein approximation theorem. Tatra Mountains Mathematical Publications, 49(1): 99-109.
Ehlers, L. (2001). Independence axioms for the provision of multiple public goods as options. Mathematical Social Sciences, 41(2): 239-250.
Ehlers, L. (2002). Multiple public goods and lexicographic preferences: replacement principle. Journal of Mathematical Economics, 37(1): 1-15.
Gatenby, R. A., & Gillies, R. J. (2004). Why do cancers have high aerobic glycolysis? Nature Reviews Cancer, 4(11): 891-899.
Gatenby, R.A., Gawlinski, E.T., Gmitro, A.F., Kaylor, B., & Gillies, R.J. (2006). Acid-mediated tumor invasion: a multidisciplinary study. Cancer Research, 66(10): 5216-5223.
Gillies, R.J., & Gatenby, R.A. (2007). Adaptive landscapes and emergent phenotypes: why do cancers have high glycolysis?. Journal of Bioenergetics and Biomembranes, 39(3): 251-257.
Gokhale, C. S., & Traulsen, A. (2010). Evolutionary games in the multiverse. Proceedings of the National Academy of Sciences, 107(12): 5500-5504.
Levasseur, K.M. (1984). A probabilistic proof of the Weierstrass approximation theorem. Amer. Math. Monthly, 91: 249-250.
Peña, J., Lehmann, L., & Nöldeke, G. (2014a). Gains from switching and evolutionary stability in multi-player matrix games. Journal of Theoretical Biology, 346, 23-33 PMID: 24380778
Peña, J., Nöldeke, G., & Lehmann, L. (2014b) Relatedness and synergies of kind and scale in the evolution of helping. arXiv: 1412.1729
Sable, K. A., & Kling, R. W. (2001). The Double Public Good: A Conceptual Framework for“Shared Experience”Values Associated with Heritage Conservation. Journal of Cultural Economics, 25(2): 77-89.
Santos, F. C., Santos, M. D., & Pacheco, J. M. (2008). Social diversity promotes the emergence of cooperation in public goods games. Nature, 454(7201): 213-216.
Umezawa, M. (2012). The replacement principle for the provision of multiple public goods on tree networks. Social Choice and Welfare, 38(2): 211-235.