A stochastic model of B cell affinity maturation and a network model of immune memory

Many events in the vertebrate immune system are influenced by some element of chance. The objective of the present work is to describe affinity maturation of B lymphocytes (in which random events are perhaps the most characteristic), and to study a possible network model of immune memory. In our model stochastic processes govern all events. A major novelty of this approach is that it permits studying random variations in the immune process. Four basic components are simulated in the model: non-immune self cells, nonself cells (pathogens), B lymphocytes, and bone marrow cells that produce naive B lymphocytes. A point in a generalized shape space plus the size of the corresponding population represents nonself and non-immune self cells. On the other hand, each individual B cell is represented by a disc that models its recognition region in the shape space. Infection is simulated by an “injection” of nonself cells into the system. Division of pathogens may instigate an attack of naive B cells, which in turn may induce clonal proliferation and hypermutation in the attacking B cells, and which eventually may slow down and stop the exponential growth of pathogens. Affinity maturation of newly produced B cells becomes expressed as a result of selection when the number of pathogens decreases. Under favorable conditions, the expanded primary B cell clones may stimulate the expansion of secondary B cell clones carrying complementary receptors to the stimulating B cells. Like in a hall of mirrors, the image of pathogens in the primary B cell clones then will be reflected in secondary B cell clones. This “ping-pong” game may survive for a long time even in the absence of the pathogen, creating a local network memory. This memory ensures that repeated infection by the same pathogen will be eliminated more efficiently.

💡 Research Summary

The paper presents a fully stochastic, agent‑based model of B‑cell affinity maturation and proposes a network‑based mechanism for long‑term immune memory. Four elementary entities are simulated: non‑immune self cells, non‑self (pathogen) cells, bone‑marrow derived naive B cells, and mature B cells. Both self and pathogen populations are represented as points in a high‑dimensional “shape space” together with their population size, while each B cell is modeled as a disc that defines its recognition region in the same space. An infection is introduced by injecting pathogen points; pathogen replication can trigger naive B cells to become activated. Activated B cells undergo clonal expansion and stochastic hyper‑mutation, producing offspring with altered recognition discs. The mutation rate is not fixed; it varies with pathogen load and the current B‑cell density, thereby mimicking antigen‑dependent mutation observed in vivo.

During clonal expansion, selection is driven by the distance between a B‑cell’s disc and pathogen points, i.e., its affinity. Cells with higher affinity proliferate faster, whereas low‑affinity clones either die or continue mutating in hopes of attaining better affinity. As the pathogen population declines, selection pressure intensifies, leading to the emergence of a highly affine primary clone.

A novel aspect of the model is the emergence of secondary clones that carry receptors complementary to those of the primary clone. The expanded primary clone creates a local “antigen image” in the shape space; B cells whose recognition regions are complementary to this image receive stimulation and begin to expand. This produces a “ping‑pong” interaction between primary and secondary clones that can persist for thousands of simulated time units even after the pathogen has been cleared. The authors interpret this sustained reciprocal activation as a local network memory, distinct from classical antibody‑antigen memory.

Simulation results demonstrate several key points: (1) stochastic hyper‑mutation coupled with affinity‑based selection efficiently drives the maturation of B‑cell repertoires; (2) the primary clone’s antigen image can seed complementary secondary clones, establishing a self‑reinforcing network that accelerates responses upon re‑infection; (3) this network can survive in the absence of antigen, providing a mechanistic basis for long‑lasting, cell‑intrinsic memory; and (4) a fully probabilistic framework captures both the random nature of somatic hyper‑mutation and the emergent collective dynamics that deterministic differential‑equation models overlook.

By integrating random mutation, selection, and inter‑clone network interactions, the model offers a quantitative tool for exploring how stochastic processes shape adaptive immunity and how durable memory may arise from the internal dynamics of B‑cell populations. The findings have implications for vaccine design, where inducing robust network memory could improve protective efficacy, and for autoimmune research, where dysregulated stochastic interactions might underlie pathological memory formation.