Optimal evolutionary decision-making to store immune memory

  1. Oskar H Schnaack
  2. Armita Nourmohammad  Is a corresponding author
  1. Max Planck Institute for Dynamics and Self-organization, Germany
  2. Department of Physics, University of Washington, United States
  3. Fred Hutchinson Cancer Research Center, United States
4 figures and 1 additional file

Figures

Immune memory or naive response upon infection.

(A) Schematic shows affinity maturation in germinal centers(right), where B-cell receptors acquire mutations and undergo selection, resulting in an increase in their affinity to an antigen (from light to dark receptors), indicated by the sharpening of receptors’ affinity profiles (on left). (B) Upon infection, the immune system can initiate a novel response (top) or a memory response (bottom). A novel B-cell response could involve affinity maturation to generate memory or high-affinity plasma cells (pink) that can secrete antibodies to battle the pathogen. A novel response can take 1–2 weeks, during which pathogen can replicate within a host and a patient can show symptoms from the disease (top, left). During this time, the proliferation of pathogens within a host incurs a cost associated with a naive response Ωτ, which is a monotonic function of the deliberation time τ (top, right). If the host carries memory from a previous infection or vaccination (bottom), the immune system can robustly and rapidly activate a memory response to battle the infection. The probability to mount such memory response Qmem. depends non-linearly on the relative utilities of memory versus naïve responses against a given infection ΔU=Umem.-Unaive (bottom, right). (C) Affinity profile Eα,θ(rm,υ)αexp[-(αd)θ] of a memory receptor rm is shown in orange as a function of the distance d=υr*-υ in the antigenic shape space, between the receptor’s cognate antigen υr* (orange) and an evolved novel target υi (red). The affinity of a receptor decays with increasing distance between targets and its cognate antigen. The antigenic range over which a receptor is reactive inversely depends on its specificity α. The shape of the binding profile is tuned by the factor θ, here shown for θ=2. The expected binding profile Eα,θ(i)(υ)¯ and the expected utility U for an immune response are weighted averages of these quantities over memory and naïve responses. The Kullback-Leibler distance between the expected profile Eα,θ(i)(υ)¯ and the profile centered around the infecting antigen Eα,θ(rυi,υ), in units of the deliberation factor β, defines the sub-optimality of a response, that is,, dissipation Kdiss (Equation 1). The net utility Unet measures the goodness of a decision to mount a memory vs. naive response against an infection (Equation 2). (D) Antigenic evolution of the H3N2 influenza virus is shown over 40 years along its first (most variable) antigenic dimension (data from Bedford et al., 2014). The decision of an immune system to utilize memory or to mount a novel response (B,C) is determined by the specificity α of receptors and the deliberation factor β. We characterize the optimal immune strategies (α*,β*) by maximizing the total net utility of immune responses against pathogens with different antigenic divergences, experienced over the lifetime of an organisms (Equation 3).

Figure 2 with 3 supplements
Optimal memory strategies against evolving pathogens.

(A) and (B) show the optimal specificity α^*α*/αmax and deliberation factor β^*β*/βmax, scaled by their respective upper bounds, as a function of the antigenic divergence per infection, scaled by the cross-reactive range (or inverse of maximum specificity) δ^δ/(αmax-1). Colors/markers indicate different naïve cost functions for deliberation, including no-cost Ω^Ω/Emax=0, linear cost Ω^=Ω^0β^, and quadratic cost Ω^=Ω^0β^2, with varying amplitudes Ω0. (C) The heat map shows the expected rescaled net utility U^net=Unet/Emax (Equation 2) per round of infection for an immune system with an optimal specificity α^*, as a function of rescaled antigenic divergence δ^ and deliberation factor β^. Rescaling by Emax sets the magnitude of net utility to one, for a response to conserved antigens (with δ^=0) and in the limit of zero deliberation cost Ω0. Boundaries indicate different levels of dissipation, with orange and blue encompassing regions of 40% and 70% of the maximum dissipation Kmax, respectively. The three modes of immune response are indicate based on the magnitude of dissipation and net utility in each reagion: (i) equilibrium memory, (ii) non-equilibrium memory, and (iii) equilibrium naive. Simulation parameters, (A–C): αmax=4, βmax=10, and θ=2, (C): linear deliberation cost function Ω^=Ω^0β^ with Ω^0=0.1. Results for other shape parameters θ and specificity thresholds αmax are shown in Figure 2—figure supplement 2, Figure 2—figure supplement 3, respectively.

Figure 2—figure supplement 1
Utility, dissipation, and usage of optimal memory.

Utility, dissipation, and usage of optimal memory. (A) and (B) show the scaled net utility U^netUnet/Emax (Equation 2) and dissipation K^dissKdiss/Emax (Equation 1) per round of infection as a function of the antigenic divergence δ^. Rescaling by Emax sets the magnitude of net utility for a response to conserved antigens (with δ^=0), and in the limit of zero deliberation cost Ω^0, to 1; see Figure 2 in the main text for comparison. (C) The expected number of rounds that a memory receptor is utilized prior to a novel response in an optimal system is shown to decay as the antigenic divergence δ^ increases. The results are evaluated for immune systems with optimized strategies (α^*, β^*) against pathogens with a given scaled antigenic divergence δ^; the corresponding strategies are shown in Figure 2. Colors/markers indicate different naïve cost functions for deliberation, including no-cost Ω^Ω/Emax=0, linear cost Ω^=Ω^0β^, and quadratic cost Ω^=Ω^0β^2, with varying amplitudes Ω^0. Simulation parameters: αmax=4, βmax=10, and θ=2.

Figure 2—figure supplement 2
Optimal memory strategies for different specificity shape factors θ.

Optimal memory strategies for different specificity shape factors θ. (A) Scaled specificity α^*α*/αmax, (B) scaled deliberation factor β^*β*/βmax, (C) scaled net utility U^netUnet/Emax, and (D) scaled dissipation are shown as a function of the scaled antigenic divergence per infection δ^=δ/(αmax-1) (similar to Figure 2). Colors indicate different shape factors θ of the specificity function, ranging from a double-sided exponential (θ=1), to Gaussian for θ=2 (as in Figure 2), and top-hat functions θ>2. The dependence of memory strategies on antigenic divergence is qualitatively insensitive to the shape factor of the specificity function. Simulation parameters: linear deliberation cost function Ω=Ω^0β^ with Ω^0=0.1, αmax=4, and βmax=10.

Figure 2—figure supplement 3
Optimal memory strategies for different specificity thresholds αmax.

Optimal memory strategies for different specificity thresholds α𝐦𝐚𝐱. (A) Scaled specificity α^*α*/αmax, (B) scaled deliberation factor β^*β*/βmax, (C) scaled net utility U^netUnet/Emax, and (D) scaled dissipation are shown as a function of the scaled antigenic divergence per infection δ^=δ/(αmax-1) (similar to Figure 2). Colors indicate different specificity thresholds αmax. Memory strategies are qualitatively insensitive to the specificity threshold. Simulation parameters: linear deliberation cost function Ω=Ω^0β^ with Ω^0=0.1 and βmax=10.

Figure 3 with 1 supplement
Mixed memory strategy against a mixture of pathogens with a broad range of evolutionary rates.

Distribution of scaled optimized specificities α^* for functional memory (purple) is shown for an immune system with a fixed deliberation factor β^=0.2. A mixture strategy with a bimodal distribution of specificities P(α^) is established to counter pathogens with a broad range of antigenic divergences. The dashed bar indicates stored memory with specificity α=0, which is not further used in response to infections. The solid line indicates the probability Pusage that a stored memory with a given specificity is utilized in future infections (Materials and methods). Optimization is done by maximizing the net utility of immune response averaged over encounters with 1000 independently evolving antigens with (scaled) antigenic divergences drawn uniformly from a range δ^(0,1.6) (Materials and methods). The distribution shows the ensemble statistics of functional memory accumulated from 200 independent optimizations, each starting from a flat prior for specificities (orange). The insert shows the optimized mixture strategy for one optimization with 3000 steps. Simulation parameters: αmax=4, βmax=10, and θ=2.

Figure 3—figure supplement 1
Mixed memory strategy against pathogens for different deliberation factors β^.

Mixed memory strategy against pathogens for different deliberation factors β^. Distribution of scaled optimized specificities α^* of functional memories is shown for an immune system with a fixed deliberation factor β^=0.2, in which a mixture strategy with a bimodal distribution of specificities P(α^) is established to counter pathogens with a broad range of antigenic divergences, drawn uniformly from a range δ^(0,1.6) (similar to Figure 3). The dashed bars indicate stored memory with specificity α=0, which is not further used in response to infections. Colors indicate different deliberation factors. Simulation parameters: αmax=4, and βmax=10.

Figure 4 with 1 supplement
Life expectancy influences the specificity of optimal memory.

(A,B) Memory strategies, that is, optimal rescaled specificity α^* (green) and deliberation factor β^* (orange) are shown as a function of the organism’s life expectancy (bottom axis) and the corresponding expected antigenic divergence over the organism’s life-time δ^lifetime (top axis). Antigenic divergence (per encounter) of the infecting pathogen is δ^=0.35 in (A) and δ^=0.5 in (B). Memory is highly specific in organisms with very short lifetimes, during which re-infections with evolved forms of a pathogen are unlikely (i.e. when life-expected antigenic divergence is smaller than 1, indicated by a dotted pink line). Memory becomes more cross-reactive with a smaller deliberation in organisms with (realistic) short lifetimes, up to a transition point c* (indicated by dotted purple line), after which specificity increases again. (C) Scaled net utility U^net is shown as a function of organism’s life span, whose immune strategies (α^*, β^*) are optimized for a specified life expectancy (colors as indicated in the legend). Net utility for memory optimized against pathogens with antigenic divergence δ^=0.35 (panel A) and δ^=0.5 (panel B) are shown by full and dashed lines, respectively. Life span and life expectancy are measured in units of the number of pathogenic encounters during lifetime. Simulation parameters: linear deliberation cost function Ω=Ω0β^ with an amplitude Ω^0=0.1, αmax=4, βmax=10, and θ=2.

Figure 4—figure supplement 1
Pathogen encounter threshold to transition between cross-reactive and specific memory.

Pathogen encounter threshold to transition between cross-reactive and specific memory. (A) The encounter threshold c*, shown in Figure 4A,B, decays as a function of the antigenic divergence (per encounter) δ^ and the amplitude of the naive cost Ω^0 (colors). (B) The expected antigenic divergence for the duration of c* (threshold) encounters δ^c* is shown as a function of antigenic divergence (per encounter) δ^. Simulation parameters: linear deliberation cost function Ω=Ω0β^, αmax=4, βmax=10, and θ=2.

Additional files

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Oskar H Schnaack
  2. Armita Nourmohammad
(2021)
Optimal evolutionary decision-making to store immune memory
eLife 10:e61346.
https://doi.org/10.7554/eLife.61346