Abstract:The germinal center response or reaction (GCR) is a hallmark event of adaptive humoral immunity. Unfolding in the B cell follicles of the secondary lymph organs, a GC culminates in the production of high-affinity antibody-secreting plasma cells along with memory B cells. By interacting with follicular dendritic cells (FDC) and T follicular helper (Tfh) cells, GC B cells exhibit complex spatiotemporal dynamics. Driving the B cell dynamics are the intracellular signal transduction and gene regulatory network that responds to cell surface signaling molecules, cytokines, and chemokines. As our knowledge of the GC continues to expand in depth and in scope, mathematical modeling has become an important tool to help disentangle the intricacy of the GCR and inform novel mechanistic and clinical insights. While the GC has been modeled at different granularities, a multiscale spatial simulation framework – integrating molecular, cellular, and tissue-level responses – is still rare. Here, we report our recent progress toward this end with a hybrid stochastic GC framework developed on the Cellular Potts Model-based CompuCell3D platform. Tellurium is used to simulate the B cell intracellular molecular network comprising NF-κB, FOXO1, MYC, AP4, CXCR4, and BLIMP1 that responds to B cell receptor (BCR) and CD40-mediated signaling. The molecular outputs of the network drive the spatiotemporal behaviors of B cells, including cyclic migration between the dark zone (DZ) and light zone (LZ) via chemotaxis; clonal proliferative bursts, somatic hypermutation, and DNA damage-induced apoptosis in the DZ; and positive selection, apoptosis via a death timer, and emergence of plasma cells in the LZ. Our simulations are able to recapitulate key molecular, cellular, and morphological GC events including B cell population growth, affinity maturation, and clonal dominance. This novel modeling framework provides an open-source, customizable, and multiscale virtual GC simulation platform that enables qualitative and quantitative investigations of a range of mechanic and applied research questions in future.

Mathematical modeling of lymphocytes selection in the germinal center

Calculating Germinal Centre Reactions

NetNeg: A Connectionist-Agent Integrated System for Representing Musical Knowledge

An analysis of B cell selection mechanisms in germinal centers

A mathematical model of tumor-immune interactions incorporated with danger model

Optimality of Mutation and Selection in Germinal Centers

Quantitative modeling of the terminal differentiation of B cells and mechanisms of lymphomagenesis

On Possible Indicators of Negative Selection in Germinal Centers

A Method to Analyze Models of the Dynamics of Germinal Centers

A Stochastic Model of the Germinal Center Integrating Local Antigen Competition, Individualistic T-B Interactions, and B Cell Receptor Signaling.

A Mathematical Model on Germinal Center Kinetics and Termination

Antibody Dynamics Simulation-A Mathematical Exploration of Clonal Deletion and Somatic Hypermutation.

Stochastic modeling of a gene regulatory network driving B cell development in germinal centers

Random walks on binary strings applied to the somatic hypermutation of B-cells

A structured population model of clonal selection in acute leukemias with multiple maturation stages

Generation of a Model to Predict Differentiation and Migration of Lymphocyte Subsets under Homeostatic and CNS Autoinflammatory Conditions

A Multiscale Spatial Modeling Framework for the Germinal Center Response

Rare event simulation in immune biology: Models of negative selection in T-cell maturation

Mathematical modeling of antigenicity for HIV dynamics

Recycling probability and dynamical properties of germinal center reactions

A branching stochastic evolutionary model of the B-cell repertoire