Amount of antigen, T follicular helper cells and affinity of founder cells shape the diversity of germinal center B cells: A computational study

Soluble ACE2 correlates with severe COVID-19 and can impair antibody responses

Identifying immune modulators that impact neutralizing antibody responses against severe acute respiratory syndrome coronavirus type 2 (SARS-CoV-2) is of great relevance. We postulated that high serum concentrations of soluble angiotensin-converting enzyme 2 (sACE2) might mask the spike and interfere with antibody maturation toward the SARS-CoV-2-receptor-binding motif (RBM). We tested 717 longitudinal samples from 295 COVID-19 patients and showed a 2- to 10-fold increase of enzymatically active sACE2 (a-sACE2), with up to 1 μg/mL total sACE2 in moderate and severe patients. Fifty percent of COVID-19 sera inhibited ACE2 activity, in contrast to 1.3% of healthy donors and 4% of non-COVID-19 pneumonia patients. A mild inverse correlation of a-sACE2 with RBM-directed serum antibodies was observed. In silico, we show that sACE2 concentrations measured in COVID-19 sera can disrupt germinal center formation and inhibit timely production of high-affinity antibodies. We suggest that sACE2 is a biomarker for COVID-19 and that soluble receptors may contribute to immune suppression informing vaccine design.

A Multiscale Spatial Modeling Framework for the Germinal Center Response

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 in silico investigations of a range of mechanic and applied research questions in future.

Germinal centers are permissive to subdominant antibody responses

Introduction

A protective humoral response to pathogens requires the development of high affinity antibodies in germinal centers (GC). The combination of antigens available during immunization has a strong impact on the strength and breadth of the antibody response. Antigens can display various levels of immunogenicity, and a hierarchy of immunodominance arises when the GC response to an antigen dampens the response to other antigens. Immunodominance is a challenge for the development of vaccines to mutating viruses, and for the development of broadly neutralizing antibodies. The extent by which antigens with different levels of immunogenicity compete for the induction of high affinity antibodies and therefore contribute to immunodominance is not known.

Methods

Here, we perform in silico simulations of the GC response, using a structural representation of antigens with complex surface amino acid composition and topology. We generate antigens with complex domains of different levels of immunogenicity and perform simulations with combinations of these domains.

Results

We found that GC dynamics were driven by the most immunogenic domain and immunodominance arose as affinity maturation to less immunogenic domain was inhibited. However, this inhibition was moderate since the less immunogenic domain exhibited a weak GC response in the absence of the most immunogenic domain. Less immunogenic domains reduced the dominance of GC responses to more immunogenic domains, albeit at a later time point.

Discussion

The simulations suggest that increased vaccine valency may decrease immunodominance of the GC response to strongly immunogenic domains and therefore, act as a potential strategy for the natural induction of broadly neutralizing antibodies in GC reactions.