Manipulating the selection forces during affinity maturation to generate cross-reactive HIV antibodies

Antibody Fc receptor binding and T cell responses to homologous and heterologous immunization with inactivated or mRNA vaccines against SARS-CoV-2

Whole virion inactivated vaccine CoronaVac (C) and Spike (S) mRNA BNT162b2 (B) vaccines differ greatly in their ability to elicit neutralizing antibodies but have somewhat comparable effectiveness in protecting from severe COVID-19. We conducted further analyses for a randomized trial (Cobovax study, NCT05057169) of third dose homologous and heterologous booster vaccination, i.e. four interventions CC-C, CC-B, BB-C and BB-B. Here, we assess vaccine immunogenicity beyond neutralizing function, including S and non-S antibodies with Fc receptor (FcR) binding, antibody avidity and T cell specificity to 6 months post-vaccination. Ancestral and Omicron S-specific IgG and FcR binding are significantly higher by BNT162b2 booster than CoronaVac, regardless of first doses. Nucleocapsid (N) antibodies are only increased in homologous boosted CoronaVac participants (CC-C). CoronaVac primed participants have lower baseline S-specific CD4+ IFNγ+ cells, but are significantly increased by either CoronaVac or BNT162b2 boosters. Priming vaccine content defined T cell peptide specificity preference, with S-specific T cells dominating B primed groups and non-S structural peptides contributing more in C primed groups, regardless of booster type. S-specific CD4+ T cell responses, N-specific antibodies, and antibody effector functions via Fc receptor binding may contribute to protection and compensate for less potent neutralizing responses in CoronaVac recipients.

Modelling HIV-1 control and remission

Remarkable advances are being made in developing interventions for eliciting long-term remission of HIV-1 infection. The success of these interventions will obviate the need for lifelong antiretroviral therapy, the current standard-of-care, and benefit the millions living today with HIV-1. Mathematical modelling has made significant contributions to these efforts. It has helped elucidate the possible mechanistic origins of natural and post-treatment control, deduced potential pathways of the loss of such control, quantified the effects of interventions, and developed frameworks for their rational optimization. Yet, several important questions remain, posing challenges to the translation of these promising interventions. Here, we survey the recent advances in the mathematical modelling of HIV-1 control and remission, highlight their contributions, and discuss potential avenues for future developments.

Parallel HIV-1 evolutionary dynamics in humans and rhesus macaques who develop broadly neutralizing antibodies

Human immunodeficiency virus (HIV)-1 exhibits remarkable genetic diversity. For this reason, an effective HIV-1 vaccine must elicit antibodies that can neutralize many variants of the virus. While broadly neutralizing antibodies (bnAbs) have been isolated from HIV-1 infected individuals, a general understanding of the virus-antibody coevolutionary processes that lead to their development remains incomplete. We performed a quantitative study of HIV-1 evolution in two individuals who developed bnAbs. We observed strong selection early in infection for mutations affecting HIV-1 envelope glycosylation and escape from autologous strain-specific antibodies, followed by weaker selection for bnAb resistance later in infection. To confirm our findings, we analyzed data from rhesus macaques infected with viruses derived from the same two individuals. We inferred remarkably similar fitness effects of HIV-1 mutations in humans and macaques. Moreover, we observed a striking pattern of rapid HIV-1 evolution, consistent in both humans and macaques, that precedes the development of bnAbs. Our work highlights strong parallels between infection in rhesus macaques and humans, and it reveals a quantitative evolutionary signature of bnAb development.

Universal characterization of epitope immunodominance from a multiscale model of clonal competition in germinal centers

We introduce a multiscale model for affinity maturation, which aims to capture the intraclonal, interclonal, and epitope-specific organization of the B-cell population in a germinal center. We describe the evolution of the B-cell population via a quasispecies dynamics, with species corresponding to unique B-cell receptors (BCRs), where the desired multiscale structure is reflected on the mutational connectivity of the accessible BCR space, and on the statistical properties of its fitness landscape. Within this mathematical framework, we study the competition among classes of BCRs targeting different antigen epitopes, and we construct an effective immunogenic space where epitope immunodominance relations can be universally characterized. We finally study how varying the relative composition of a mixture of antigens with variable and conserved domains allows for a parametric exploration of this space, and we identify general principles for the rational design of two-antigen cocktails.

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 lymphoid 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 mechanistic and applied research questions on the adaptive humoral immune response in the future.

Broad antibody and T cell responses to Ebola, Sudan, and Bundibugyo ebolaviruses using mono- and multi-valent adjuvanted glycoprotein vaccines

Currently, there are two approved vaccine regimens designed to prevent Ebola virus (EBOV) disease (EVD). Both are virus-vectored, and concerns about cold-chain storage and pre-existing immunity to the vectors warrant investigating additional vaccine strategies. Here, we have explored the utility of adjuvanted recombinant glycoproteins (GPs) from ebolaviruses Zaire (EBOV), Sudan (SUDV), and Bundibugyo (BDBV) for inducing antibody (Ab) and T cell cross-reactivity. Glycoproteins expressed in insect cells were administered to C57BL/6 mice as free protein or bound to the surface of liposomes, and formulated with toll-like receptor agonists CpG and MPLA (agonists for TLR 9 and 4, respectively), with or without the emulsions AddaVax or TiterMax. The magnitude of Ab cross-reactivity in binding and neutralization assays, and T cell cross-reactivity in antigen recall assays, correlated with phylogenetic relatedness. While most adjuvants screened induced IgG responses, a combination of CpG, MPLA and AddaVax emulsion ("IVAX-1") was the most potent and polarized in an IgG2c (Th1) direction. Breadth was also achieved by combining GPs into a trivalent (Tri-GP) cocktail with IVAX-1, which did not compromise antibody responses to individual components in binding and neutralizing assays. Th1 signature cytokines in T cell recall assays were undetectable after Tri-GP/IVAX-1 administration, despite a robust IgG2c response, although administration of Tri-GP on lipid nanoparticles in IVAX-1 elevated Th1 cytokines to detectable levels. Overall, the data indicate an adjuvanted trivalent recombinant GP approach may represent a path toward a broadly reactive, deployable vaccine against EVD.

Repeated vaccination with homologous influenza hemagglutinin broadens human antibody responses to unmatched flu viruses

The on-going diversification of influenza virus necessicates annual vaccine updating. The vaccine antigen, the viral spike protein hemagglutinin (HA), tends to elicit strain-specific neutralizing activity, predicting that sequential immunization with the same HA strain will boost antibodies with narrow coverage. However, repeated vaccination with homologous SARS-CoV-2 vaccine eventually elicits neutralizing activity against highly unmatched variants, questioning this immunological premise. We evaluated a longitudinal influenza vaccine cohort, where each year the subjects received the same, novel H1N1 2009 pandemic vaccine strain. Repeated vaccination gradually enhanced receptor-blocking antibodies (HAI) to highly unmatched H1N1 strains within individuals with no initial memory recall against these historical viruses. An in silico model of affinity maturation in germinal centers integrated with a model of differentiation and expansion of memory cells provides insight into the mechanisms underlying these results and shows how repeated exposure to the same immunogen can broaden the antibody response against diversified targets.

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.

A Quantitative Clinical Pharmacology-Based Framework For Model-Informed Vaccine Development

Historically, vaccine development and dose optimization have followed mostly empirical approaches without clinical pharmacology and model-informed approaches playing a major role, in contrast to conventional drug development. This is attributed to the complex cascade of immunobiological mechanisms associated with vaccines and a lack of quantitative frameworks for extracting dose-exposure-efficacy-toxicity relationships. However, the Covid-19 pandemic highlighted the lack of sufficient immunogenicity due to suboptimal vaccine dosing regimens and the need for well-designed, model-informed clinical trials which enhance the probability of selection of optimal vaccine dosing regimens. In this perspective, we attempt to develop a quantitative clinical pharmacology-based approach that integrates vaccine dose-efficacy-toxicity across various stages of vaccine development into a unified framework that we term as model-informed vaccine dose-optimization and development (MIVD). We highlight scenarios where the adoption of MIVD approaches may have a strategic advantage compared to conventional practices for vaccines.

Steering and controlling evolution - from bioengineering to fighting pathogens

Control interventions steer the evolution of molecules, viruses, microorganisms or other cells towards a desired outcome. Applications range from engineering biomolecules and synthetic organisms to drug, therapy and vaccine design against pathogens and cancer. In all these instances, a control system alters the eco-evolutionary trajectory of a target system, inducing new functions or suppressing escape evolution. Here, we synthesize the objectives, mechanisms and dynamics of eco-evolutionary control in different biological systems. We discuss how the control system learns and processes information about the target system by sensing or measuring, through adaptive evolution or computational prediction of future trajectories. This information flow distinguishes pre-emptive control strategies by humans from feedback control in biotic systems. We establish a cost-benefit calculus to gauge and optimize control protocols, highlighting the fundamental link between predictability of evolution and efficacy of pre-emptive control.

Engaging an HIV vaccine target through the acquisition of low B cell affinity

Low affinity is common for germline B cell receptors (BCR) seeding development of broadly neutralizing antibodies (bnAbs) that engage hypervariable viruses, including HIV. Antibody affinity selection is also non-homogenizing, insuring the survival of low affinity B cell clones. To explore whether this provides a natural window for expanding human B cell lineages against conserved vaccine targets, we deploy transgenic mice mimicking human antibody diversity and somatic hypermutation (SHM) and immunize with simple monomeric HIV glycoprotein envelope immunogens. We report an immunization regimen that focuses B cell memory upon the conserved CD4 binding site (CD4bs) through both conventional affinity maturation and reproducible expansion of low affinity BCR clones with public patterns in SHM. In the latter instance, SHM facilitates target acquisition by decreasing binding strength. This suggests that permissive B cell selection enables the discovery of antibody epitopes, in this case an HIV bnAb site.

'Persistent germinal center responses: slow-growing trees bear the best fruits'

Germinal centers (GCs) are key microanatomical sites in lymphoid organs where responding B cells mature and undergo affinity-based selection. The duration of the GC reaction has long been assumed to be relatively brief, but recent studies in humans, nonhuman primates, and mice indicate that GCs can last for weeks to months after initial antigen exposure. This review examines recent studies investigating the factors that influence GC duration, including antigen persistence, T-follicular helper cells, and mode of immunization. Potential mechanisms for how persistent GCs influence the B-cell repertoire are considered. Overall, these studies provide a blueprint for how to design better vaccines that elicit persistent GC responses.

Inspecting the interaction between human immunodeficiency virus and the immune system through genetic turnover

Chronic infections of the human immunodeficiency virus (HIV) create a very complex coevolutionary process, where the virus tries to escape the continuously adapting host immune system. Quantitative details of this process are largely unknown and could help in disease treatment and vaccine development. Here we study a longitudinal dataset of ten HIV-infected people, where both the B-cell receptors and the virus are deeply sequenced. We focus on simple measures of turnover, which quantify how much the composition of the viral strains and the immune repertoire change between time points. At the single-patient level, the viral-host turnover rates do not show any statistically significant correlation, however, they correlate if one increases the amount of statistics by aggregating the information across patients. We identify an anti-correlation: large changes in the viral pool composition come with small changes in the B-cell receptor repertoire. This result seems to contradict the naïve expectation that when the virus mutates quickly, the immune repertoire needs to change to keep up. However, a simple model of antagonistically evolving populations can explain this signal. If it is sampled at intervals comparable with the sweep time, one population has had time to sweep while the second cannot start a counter-sweep, leading to the observed anti-correlation. This article is part of the theme issue 'Interdisciplinary approaches to predicting evolutionary biology'.

Antigen presentation dynamics shape the antibody response to variants like SARS-CoV-2 Omicron after multiple vaccinations with the original strain

The Omicron variant of SARS-CoV-2 is not effectively neutralized by most antibodies elicited by two doses of mRNA vaccines, but a third dose increases anti-Omicron neutralizing antibodies. We reveal mechanisms underlying this observation by combining computational modeling with data from vaccinated humans. After the first dose, limited antigen availability in germinal centers (GCs) results in a response dominated by B cells that target immunodominant epitopes that are mutated in an Omicron-like variant. After the second dose, these memory cells expand and differentiate into plasma cells that secrete antibodies that are thus ineffective for such variants. However, these pre-existing antigen-specific antibodies transport antigen efficiently to secondary GCs. They also partially mask immunodominant epitopes. Enhanced antigen availability and epitope masking in secondary GCs together result in generation of memory B cells that target subdominant epitopes that are less mutated in Omicron. The third dose expands these cells and boosts anti-variant neutralizing antibodies.

Quantitatively Visualizing Bipartite Datasets

As experiments continue to increase in size and scope, a fundamental challenge of subsequent analyses is to recast the wealth of information into an intuitive and readily interpretable form. Often, each measurement conveys only the relationship between a pair of entries, and it is difficult to integrate these local interactions across a dataset to form a cohesive global picture. The classic localization problem tackles this question, transforming local measurements into a global map that reveals the underlying structure of a system. Here, we examine the more challenging bipartite localization problem, where pairwise distances are available only for bipartite data comprising two classes of entries (such as antibody-virus interactions, drug-cell potency, or user-rating profiles). We modify previous algorithms to solve bipartite localization and examine how each method behaves in the presence of noise, outliers, and partially observed data. As a proof of concept, we apply these algorithms to antibody-virus neutralization measurements to create a basis set of antibody behaviors, formalize how potently inhibiting some viruses necessitates weakly inhibiting other viruses, and quantify how often combinations of antibodies exhibit degenerate behavior.

Mechanisms that promote the evolution of cross-reactive antibodies upon vaccination with designed influenza immunogens

Immunogens that elicit broadly neutralizing antibodies targeting the conserved receptor-binding site (RBS) on influenza hemagglutinin may serve as candidates for a universal influenza vaccine. Here, we develop a computational model to interrogate antibody evolution by affinity maturation after immunization with two types of immunogens: a heterotrimeric "chimera" hemagglutinin that is enriched for the RBS epitope relative to other B cell epitopes and a cocktail composed of three non-epitope-enriched homotrimers of the monomers that comprise the chimera. Experiments in mice find that the chimera outperforms the cocktail for eliciting RBS-directed antibodies. We show that this result follows from an interplay between how B cells engage these antigens and interact with diverse helper T cells and requires T cell-mediated selection of germinal center B cells to be a stringent constraint. Our results shed light on antibody evolution and highlight how immunogen design and T cells modulate vaccination outcomes.

A Germline-Targeting Chimpanzee SIV Envelope Glycoprotein Elicits a New Class of V2-Apex Directed Cross-Neutralizing Antibodies

HIV-1 and its SIV precursors share a broadly neutralizing antibody (bNAb) epitope in variable loop 2 (V2) at the envelope glycoprotein (Env) trimer apex. Here, we tested the immunogenicity of germ line-targeting versions of a chimpanzee SIV (SIVcpz) Env in human V2-apex bNAb heavy-chain precursor-expressing knock-in mice and as chimeric simian-chimpanzee immunodeficiency viruses (SCIVs) in rhesus macaques (RMs). Trimer immunization of knock-in mice induced V2-directed NAbs, indicating activation of V2-apex bNAb precursor-expressing mouse B cells. SCIV infection of RMs elicited high-titer viremia, potent autologous tier 2 neutralizing antibodies, and rapid sequence escape in the canonical V2-apex epitope. Six of seven animals also developed low-titer heterologous plasma breadth that mapped to the V2-apex. Antibody cloning from two of these animals identified multiple expanded lineages with long heavy chain third complementarity determining regions that cross-neutralized as many as 7 of 19 primary HIV-1 strains, but with low potency. Negative stain electron microscopy (NSEM) of members of the two most cross-reactive lineages confirmed V2 targeting but identified an angle of approach distinct from prototypical V2-apex bNAbs, with antibody binding either requiring or inducing an occluded-open trimer. Probing with conformation-sensitive, nonneutralizing antibodies revealed that SCIV-expressed, but not wild-type SIVcpz Envs, as well as a subset of primary HIV-1 Envs, preferentially adopted a more open trimeric state. These results reveal the existence of a cryptic V2 epitope that is exposed in occluded-open SIVcpz and HIV-1 Env trimers and elicits cross-neutralizing responses of limited breadth and potency. IMPORTANCE An effective HIV-1 vaccination strategy will need to stimulate rare precursor B cells of multiple bNAb lineages and affinity mature them along desired pathways. Here, we searched for V2-apex germ line-targeting Envs among a large set of diverse primate lentiviruses and identified minimally modified versions of one chimpanzee SIV Env that bound several human V2-apex bNAb precursors and stimulated one of these in a V2-apex bNAb precursor-expressing knock-in mouse. We also generated chimeric simian-chimpanzee immunodeficiency viruses and showed that they elicit low-titer V2-directed heterologous plasma breadth in six of seven infected rhesus macaques. Characterization of this antibody response identified a new class of weakly cross-reactive neutralizing antibodies that target the V2-apex, but only in occluded-open Env trimers. The existence of this cryptic epitope, which in some Env backgrounds is immunodominant, needs to be considered in immunogen design.

Antibodies to repeat-containing antigens in Plasmodium falciparum are exposure-dependent and short-lived in children in natural malaria infections

Protection against Plasmodium falciparum, which is primarily antibody-mediated, requires recurrent exposure to develop. The study of both naturally acquired limited immunity and vaccine induced protection against malaria remains critical for ongoing eradication efforts. Towards this goal, we deployed a customized P. falciparum PhIP-seq T7 phage display library containing 238,068 tiled 62-amino acid peptides, covering all known coding regions, including antigenic variants, to systematically profile antibody targets in 198 Ugandan children and adults from high and moderate transmission settings. Repeat elements - short amino acid sequences repeated within a protein - were significantly enriched in antibody targets. While breadth of responses to repeat-containing peptides was twofold higher in children living in the high versus moderate exposure setting, no such differences were observed for peptides without repeats, suggesting that antibody responses to repeat-containing regions may be more exposure dependent and/or less durable in children than responses to regions without repeats. Additionally, short motifs associated with seroreactivity were extensively shared among hundreds of antigens, potentially representing cross-reactive epitopes. PfEMP1 shared motifs with the greatest number of other antigens, partly driven by the diversity of PfEMP1 sequences. These data suggest that the large number of repeat elements and potential cross-reactive epitopes found within antigenic regions of P. falciparum could contribute to the inefficient nature of malaria immunity.

Hierarchical sequence-affinity landscapes shape the evolution of breadth in an anti-influenza receptor binding site antibody

Broadly neutralizing antibodies (bnAbs) that neutralize diverse variants of a particular virus are of considerable therapeutic interest. Recent advances have enabled us to isolate and engineer these antibodies as therapeutics, but eliciting them through vaccination remains challenging, in part due to our limited understanding of how antibodies evolve breadth. Here, we analyze the landscape by which an anti-influenza receptor binding site (RBS) bnAb, CH65, evolved broad affinity to diverse H1 influenza strains. We do this by generating an antibody library of all possible evolutionary intermediates between the unmutated common ancestor (UCA) and the affinity-matured CH65 antibody and measure the affinity of each intermediate to three distinct H1 antigens. We find that affinity to each antigen requires a specific set of mutations - distributed across the variable light and heavy chains - that interact non-additively (i.e., epistatically). These sets of mutations form a hierarchical pattern across the antigens, with increasingly divergent antigens requiring additional epistatic mutations beyond those required to bind less divergent antigens. We investigate the underlying biochemical and structural basis for these hierarchical sets of epistatic mutations and find that epistasis between heavy chain mutations and a mutation in the light chain at the V_H-V_L interface is essential for binding a divergent H1. Collectively, this is the first work to comprehensively characterize epistasis between heavy and light chain mutations and shows that such interactions are both strong and widespread. Together with our previous study analyzing a different class of anti-influenza antibodies, our results implicate epistasis as a general feature of antibody sequence-affinity landscapes that can potentiate and constrain the evolution of breadth.

Long-Term Vaccine Delivery and Immunological Responses Using Biodegradable Polymer-Based Carriers

Biodegradable polymers are largely employed in the biomedical field, ranging from tissue regeneration to drug/vaccine delivery. The biodegradable polymers are highly biocompatible and possess negligible toxicity. In addition, biomaterial-based vaccines possess adjuvant properties, thereby enhancing immune responses. This Review introduces the use of different biodegradable polymers and their degradation mechanism. Different kinds of vaccines, as well as the interaction between the carriers with the immune system, then are highlighted. Natural and synthetic biodegradable micro-/nanoplatforms, hydrogels, and scaffolds for local or targeted and controlled vaccine release are subsequently discussed.

Moving the needle: Employing deep reinforcement learning to push the boundaries of coarse-grained vaccine models

Highly mutable infectious disease pathogens (hm-IDPs) such as HIV and influenza evolve faster than the human immune system can contain them, allowing them to circumvent traditional vaccination approaches and causing over one million deaths annually. Agent-based models can be used to simulate the complex interactions that occur between immune cells and hm-IDP-like proteins (antigens) during affinity maturation-the process by which antibodies evolve. Compared to existing experimental approaches, agent-based models offer a safe, low-cost, and rapid route to study the immune response to vaccines spanning a wide range of design variables. However, the highly stochastic nature of affinity maturation and vast sequence space of hm-IDPs render brute force searches intractable for exploring all pertinent vaccine design variables and the subset of immunization protocols encompassed therein. To address this challenge, we employed deep reinforcement learning to drive a recently developed agent-based model of affinity maturation to focus sampling on immunization protocols with greater potential to improve the chosen metrics of protection, namely the broadly neutralizing antibody (bnAb) titers or fraction of bnAbs produced. Using this approach, we were able to coarse-grain a wide range of vaccine design variables and explore the relevant design space. Our work offers new testable insights into how vaccines should be formulated to maximize protective immune responses to hm-IDPs and how they can be minimally tailored to account for major sources of heterogeneity in human immune responses and various socioeconomic factors. Our results indicate that the first 3 to 5 immunizations, depending on the metric of protection, should be specially tailored to achieve a robust protective immune response, but that beyond this point further immunizations require only subtle changes in formulation to sustain a durable bnAb response.

'Mix and Match' vaccination: Is dengue next?

The severity of the COVID-19 pandemic and the development of multiple SARS-CoV-2 vaccines expedited vaccine 'mix and match' trials in humans and demonstrated the benefits of mixing vaccines that vary in formulation, strength, and immunogenicity. Heterologous sequential vaccination may be an effective approach for protecting against dengue, as this strategy would mimic the natural route to broad dengue protection and may overcome the imbalances in efficacy of the individual leading live attenuated dengue vaccines. Here we review 'mix and match' vaccination trials against SARS-CoV-2, HIV, and dengue virus and discuss the possible advantages and concerns of future heterologous immunization with the leading dengue vaccines. COVID-19 trials suggest that priming with a vaccine that induces strong cellular responses, such as an adenoviral vectored product, followed by heterologous boost may optimize T cell immunity. Moreover, heterologous vaccination may induce superior humoral immunity compared to homologous vaccination when the priming vaccine induces a narrower response than the boost. The HIV trials reported that heterologous vaccination was associated with broadened antigen responses and that the sequence of the vaccines significantly impacts the regimen's immunogenicity and efficacy. In heterologous dengue immunization trials, where at least one dose was with a live attenuated vaccine, all reported equivalent or increased immunogenicity compared to homologous boost, although one study reported increased reactogenicity. The three leading dengue vaccines have been evaluated for safety and efficacy in thousands of study participants but not in combination in heterologous dengue vaccine trials. Various heterologous regimens including different combinations and sequences should be trialed to optimize cellular and humoral immunity and the breadth of the response while limiting reactogenicity. A blossoming field dedicated to more accurate correlates of protection and enhancement will help confirm the safety and efficacy of these strategies.

Vaccine models predict rules for updating vaccines against evolving pathogens such as SARS-CoV-2 and influenza in the context of pre-existing immunity

Currently, vaccines for SARS-CoV-2 and influenza viruses are updated if the new vaccine induces higher antibody-titers to circulating variants than current vaccines. This approach does not account for complex dynamics of how prior immunity skews recall responses to the updated vaccine. We: (i) use computational models to mechanistically dissect how prior immunity influences recall responses; (ii) explore how this affects the rules for evaluating and deploying updated vaccines; and (iii) apply this to SARS-CoV-2. Our analysis of existing data suggests that there is a strong benefit to updating the current SARS-CoV-2 vaccines to match the currently circulating variants. We propose a general two-dose strategy for determining if vaccines need updating as well as for vaccinating high-risk individuals. Finally, we directly validate our model by reanalysis of earlier human H5N1 influenza vaccine studies.

Harnessing Antibody Polyspecificity for Cancer Immunotherapy

Targeting the diverse glycan repertoire expressed on tumor cells is considered a viable therapeutic strategy to deal with tumor cell heterogeneity. Inherently polyspecific, natural, glycan-reactive antibodies are purported to be protective in thwarting infections and in cancer immunotherapy. Tumor-associated carbohydrate antigens (TACAs) are related to pathogen glycans, to which nascent or natural antibodies exist and IgM responses are elicited. To capture the polyspecific nature of anticarbohydrate responses, we have focused on the rational design of carbohydrate mimetic peptides (CMPs) cross-reactive with TACA reactive antibodies. In particular, we have focused on the development of CMPs that display reactivity to GD2 and Lewis Y (LeY) reactive monoclonal antibodies. They would serve as templates for pan-immunogens inducing biosimilar polyreactive antibodies. In the design, we relied on structural analyses of CMP's enhanced binding to the templates using molecular modeling. Glycan reactivity patterns of affinity CMP-purified human antibodies further refined specificity profiles in comparison with the immune response to the CMP in clinical trials. In this study, we further define the molecular characteristics for this mimicry by considering the polyspecificity of LeY and GD2 reactive antibodies binding to the lacto-ceramide core Galβ(1,4)Glcβ(1-1')Cer. Binding to this minimum building block can be capitalized on for cancer therapy and diagnostics and illustrates a new approach in designing cancer vaccines taking advantage of the latent polyspecificity of antibodies and the relevance of natural antibodies in antigen discovery and design.

Design of immunogens for eliciting antibody responses that may protect against SARS-CoV-2 variants

The rise of SARS-CoV-2 variants and the history of outbreaks caused by zoonotic coronaviruses point to the need for next-generation vaccines that confer protection against variant strains. Here, we combined analyses of diverse sequences and structures of coronavirus spikes with data from deep mutational scanning to design SARS-CoV-2 variant antigens containing the most significant mutations that may emerge. We trained a neural network to predict RBD expression and ACE2 binding from sequence, which allowed us to determine that these antigens are stable and bind to ACE2. Thus, they represent viable variants. We then used a computational model of affinity maturation (AM) to study the antibody response to immunization with different combinations of the designed antigens. The results suggest that immunization with a cocktail of the antigens is likely to promote evolution of higher titers of antibodies that target SARS-CoV-2 variants than immunization or infection with the wildtype virus alone. Finally, our analysis of 12 coronaviruses from different genera identified the S2’ cleavage site and fusion peptide as potential pan-coronavirus vaccine targets.

SARS-CoV-2 variants have already emerged and future variants may pose greater threats to the efficacy of current vaccines. Rather than using a reactive approach to vaccine development that would lag behind the evolution of the virus, such as updating the sequence in the vaccine with a current variant, we sought to use a proactive approach that predicts some of the mutations that could arise that could evade current immune responses. Then, by including these mutations in a new vaccine antigen, we might be able to protect against those potential variants before they appear. Toward this end, we used various computational methods including sequence analysis and machine learning to design such antigens. We then used simulations of antibody development, and the results suggest that immunization with our designed antigens is likely to result in an antibody response that is better able to target SARS-CoV-2 variants than current vaccines. We also leveraged our sequence analysis to suggest that a particular site on the spike protein could serve as a useful target for a pan-coronavirus vaccine.

Antigen presentation dynamics shape the response to emergent variants like SARS-CoV-2 Omicron strain after multiple vaccinations with wild type strain

The Omicron variant of SARS-CoV-2 evades neutralization by most serum antibodies elicited by two doses of mRNA vaccines, but a third dose of the same vaccine increases anti-Omicron neutralizing antibodies. By combining computational modeling with data from vaccinated humans we reveal mechanisms underlying this observation. After the first dose, limited antigen availability in germinal centers results in a response dominated by B cells with high germline affinities for immunodominant epitopes that are significantly mutated in an Omicron-like variant. After the second dose, expansion of these memory cells and differentiation into plasma cells shape antibody responses that are thus ineffective for such variants. However, in secondary germinal centers, pre-existing higher affinity antibodies mediate enhanced antigen presentation and they can also partially mask dominant epitopes. These effects generate memory B cells that target subdominant epitopes that are less mutated in Omicron. The third dose expands these cells and boosts anti-variant neutralizing antibodies.

The role of geometric features in a germinal center

The germinal center (GC) is a self-organizing structure produced in the lymphoid follicle during the T-dependent immune response and is an important component of the humoral immune system. However, the impact of the special structure of GC on antibody production is not clear. According to the latest biological experiments, we establish a spatiotemporal stochastic model to simulate the whole self-organization process of the GC including the appearance of two specific zones: the dark zone (DZ) and the light zone (LZ), the development of which serves to maintain an effective competition among different cells and promote affinity maturation. A phase transition is discovered in this process, which determines the critical GC volume for a successful growth in both the stochastic and the deterministic model. Further increase of the volume does not make much improvement on the performance. It is found that the critical volume is determined by the distance between the activated B cell receptor (BCR) and the target epitope of the antigen in the shape space. The observation is confirmed in both 2D and 3D simulations and explains partly the variability of the observed GC size.

Modeling suggests that multiple immunizations or infections will reveal the benefits of updating SARS-CoV-2 vaccines

When should vaccines to evolving pathogens such as SARS-CoV-2 be updated? Our computational models address this focusing on updating SARS-CoV-2 vaccines to the currently circulating Omicron variant. Current studies typically compare the antibody titers to the new variant following a single dose of the original-vaccine versus the updated-vaccine in previously immunized individuals. These studies find that the updated-vaccine does not induce higher titers to the vaccine-variant compared with the original-vaccine, suggesting that updating may not be needed. Our models recapitulate this observation but suggest that vaccination with the updated-vaccine generates qualitatively different humoral immunity, a small fraction of which is specific for unique epitopes to the new variant. Our simulations suggest that these new variant-specific responses could dominate following subsequent vaccination or infection with either the currently circulating or future variants. We suggest a two-dose strategy for determining if the vaccine needs updating and for vaccinating high-risk individuals.

From Antibody Repertoires to Cell-Cell Interactions to Molecular Networks: Bridging Scales in the Germinal Center

Antibody-mediated adaptive immunity must provide effective long-term protection with minimal adverse effects, against rapidly mutating pathogens, in a human population with diverse ages, genetics, and immune histories. In order to grasp and leverage the complexities of the antibody response, we advocate for a mechanistic understanding of the multiscale germinal center (GC) reaction - the process by which precursor B-cells evolve high-affinity antigen-specific antibodies, forming an effector repertoire of plasma and memory cells for decades-long protection. The regulatory dynamics of B-cells within the GC are complex, and unfold across multiple interacting spatial and temporal scales. At the organism scale, over weeks to years, the antibody sequence repertoire formed by various B-cell clonal lineages modulates antibody quantity and quality over time. At the tissue and cellular scale, over hours to weeks, B-cells undergo selection via spatially distributed interactions with local stroma, antigen, and helper T-cells. At the molecular scale, over seconds to days, intracellular signaling, transcriptional, and epigenetic networks modulate B-cell fates and shape their clonal lineages. We summarize our current understanding within each of these scales, and identify missing links in connecting them. We suggest that quantitative multi-scale mathematical models of B-cell and GC reaction dynamics provide predictive frameworks that can apply basic immunological knowledge to practical challenges such as rational vaccine design.

Multiscale affinity maturation simulations to elicit broadly neutralizing antibodies against HIV

The design of vaccines against highly mutable pathogens, such as HIV and influenza, requires a detailed understanding of how the adaptive immune system responds to encountering multiple variant antigens (Ags). Here, we describe a multiscale model of B cell receptor (BCR) affinity maturation that employs actual BCR nucleotide sequences and treats BCR/Ag interactions in atomistic detail. We apply the model to simulate the maturation of a broadly neutralizing Ab (bnAb) against HIV. Starting from a germline precursor sequence of the VRC01 anti-HIV Ab, we simulate BCR evolution in response to different vaccination protocols and different Ags, which were previously designed by us. The simulation results provide qualitative guidelines for future vaccine design and reveal unique insights into bnAb evolution against the CD4 binding site of HIV. Our model makes possible direct comparisons of simulated BCR populations with results of deep sequencing data, which will be explored in future applications.

Vaccination has saved more lives than any other medical procedure. But, we do not have robust ways to develop vaccines against highly mutable pathogens. For example, there is no effective vaccine against HIV, and a universal vaccine against diverse strains of influenza is also not available. The development of immunization strategies to elicit antibodies that can neutralize diverse strains of highly mutable pathogens (so-called ‘broadly neutralizing antibodies’, or bnAbs) would enable the design of universal vaccines against such pathogens, as well as other viruses that may emerge in the future. In this paper, we present an agent-based model of affinity maturation–the Darwinian process by which antibodies evolve against a pathogen–that, for the first time, enables the in silico investigation of real germline nucleotide sequences of antibodies known to evolve into potent bnAbs, evolving against real amino acid sequences of HIV-based vaccine-candidate proteins. Our results provide new insights into bnAb evolution against HIV, and can be used to qualitatively guide the future design of vaccines against highly mutable pathogens.

Optimizing clinical dosing of combination broadly neutralizing antibodies for HIV prevention

Broadly neutralizing antibodies (bNAbs) are promising agents to prevent HIV infection and achieve HIV remission without antiretroviral therapy (ART). As with ART, bNAb combinations are likely needed to cover HIV's extensive diversity. Not all bNAbs are identical in terms of their breadth, potency, and in vivo longevity (half-life). Given these differences, it is important to optimally select the composition, or dose ratio, of combination bNAb therapies for future clinical studies. We developed a model that synthesizes 1) pharmacokinetics, 2) potency against a wide HIV diversity, 3) interaction models for how drugs work together, and 4) correlates that translate in vitro potency to clinical protection. We found optimization requires drug-specific balances between potency, longevity, and interaction type. As an example, tradeoffs between longevity and potency are shown by comparing a combination therapy to a bi-specific antibody (a single protein merging both bNAbs) that takes the better potency but the worse longevity of the two components. Then, we illustrate a realistic dose ratio optimization of a triple combination of VRC07, 3BNC117, and 10-1074 bNAbs. We apply protection estimates derived from both a non-human primate (NHP) challenge study meta-analysis and the human antibody mediated prevention (AMP) trials. In both cases, we find a 2:1:1 dose emphasizing VRC07 is nearly optimal. Our approach can be immediately applied to optimize the next generation of combination antibody prevention and cure studies.

A Coarse-Grained Model of Affinity Maturation Indicates the Importance of B-Cell Receptor Avidity in Epitope Subdominance

The elicitation of broadly neutralizing antibodies (bnAbs) is a major goal in the design of vaccines against rapidly-mutating viruses. In the case of influenza, many bnAbs that target conserved epitopes on the stem of the hemagglutinin protein (HA) have been discovered. However, these antibodies are rare, are not boosted well upon reinfection, and often have low neutralization potency, compared to strain-specific antibodies directed to the HA head. Different hypotheses have been proposed to explain this phenomenon. We use a coarse-grained computational model of the germinal center reaction to investigate how B-cell receptor binding valency affects the growth and affinity maturation of competing B-cells. We find that receptors that are unable to bind antigen bivalently, and also those that do not bind antigen cooperatively, have significantly slower rates of growth, memory B-cell production, and, under certain conditions, rates of affinity maturation. The corresponding B-cells are predicted to be outcompeted by B-cells that bind bivalently and cooperatively. We use the model to explore strategies for a universal influenza vaccine, e.g., how to boost the concentrations of the slower growing cross-reactive antibodies directed to the stem. The results suggest that, upon natural reinfections subsequent to vaccination, the protectiveness of such vaccines would erode, possibly requiring regular boosts. Collectively, our results strongly support the importance of bivalent antibody binding in immunodominance, and suggest guidelines for developing a universal influenza vaccine.

Affinity maturation for an optimal balance between long-term immune coverage and short-term resource constraints

Humoral immunity relies on the mutation and selection of B cells to better recognize pathogens. This affinity maturation process produces cells with diverse recognition capabilities. Examining optimal immune strategies that maximize the long-term immune coverage at a minimal metabolic cost, we show when the immune system should mount a de novo response rather than rely on existing memory cells. Our theory recapitulates known modes of the B cell response, predicts the empirical form of the distribution of clone sizes, and rationalizes as a trade-off between metabolic and immune costs the antigenic imprinting effects that limit the efficacy of vaccines (original antigenic sin). Our predictions provide a framework to interpret experimental results that could be used to inform vaccination strategies.

In order to target threatening pathogens, the adaptive immune system performs a continuous reorganization of its lymphocyte repertoire. Following an immune challenge, the B cell repertoire can evolve cells of increased specificity for the encountered strain. This process of affinity maturation generates a memory pool whose diversity and size remain difficult to predict. We assume that the immune system follows a strategy that maximizes the long-term immune coverage and minimizes the short-term metabolic costs associated with affinity maturation. This strategy is defined as an optimal decision process on a finite dimensional phenotypic space, where a preexisting population of cells is sequentially challenged with a neutrally evolving strain. We show that the low specificity and high diversity of memory B cells—a key experimental result—can be explained as a strategy to protect against pathogens that evolve fast enough to escape highly potent but narrow memory. This plasticity of the repertoire drives the emergence of distinct regimes for the size and diversity of the memory pool, depending on the density of de novo responding cells and on the mutation rate of the strain. The model predicts power-law distributions of clonotype sizes observed in data and rationalizes antigenic imprinting as a strategy to minimize metabolic costs while keeping good immune protection against future strains.

Immunological Evidence of Variation in Exposure and Immune Response to Bacillus anthracis in Herbivores of Kruger and Etosha National Parks

Exposure and immunity to generalist pathogens differ among host species and vary across spatial scales. Anthrax, caused by a multi-host bacterial pathogen, Bacillus anthracis, is enzootic in Kruger National Park (KNP), South Africa and Etosha National Park (ENP), Namibia. These parks share many of the same potential host species, yet the main anthrax host in one (greater kudu (Tragelaphus strepsiceros) in KNP and plains zebra (Equus quagga) in ENP) is only a minor host in the other. We investigated species and spatial patterns in anthrax mortalities, B. anthracis exposure, and the ability to neutralize the anthrax lethal toxin to determine if observed host mortality differences between locations could be attributed to population-level variation in pathogen exposure and/or immune response. Using serum collected from zebra and kudu in high and low incidence areas of each park (18- 20 samples/species/area), we estimated pathogen exposure from anti-protective antigen (PA) antibody response using enzyme-linked immunosorbent assay (ELISA) and lethal toxin neutralization with a toxin neutralization assay (TNA). Serological evidence of pathogen exposure followed mortality patterns within each system (kudus: 95% positive in KNP versus 40% in ENP; zebras: 83% positive in ENP versus 63% in KNP). Animals in the high-incidence area of KNP had higher anti-PA responses than those in the low-incidence area, but there were no significant differences in exposure by area within ENP. Toxin neutralizing ability was higher for host populations with lower exposure prevalence, i.e., higher in ENP kudus and KNP zebras than their conspecifics in the other park. These results indicate that host species differ in their exposure to and adaptive immunity against B. anthracis in the two parks. These patterns may be due to environmental differences such as vegetation, rainfall patterns, landscape or forage availability between these systems and their interplay with host behavior (foraging or other risky behaviors), resulting in differences in exposure frequency and dose, and hence immune response.

Designing spatial and temporal control of vaccine responses

Vaccines are the key technology to combat existing and emerging infectious diseases. However, increasing the potency, quality and durability of the vaccine response remains a challenge. As our knowledge of the immune system deepens, it becomes clear that vaccine components must be in the right place at the right time to orchestrate a potent and durable response. Material platforms, such as nanoparticles, hydrogels and microneedles, can be engineered to spatially and temporally control the interactions of vaccine components with immune cells. Materials-based vaccination strategies can augment the immune response by improving innate immune cell activation, creating local inflammatory niches, targeting lymph node delivery and controlling the time frame of vaccine delivery, with the goal of inducing enhanced memory immunity to protect against future infections. In this Review, we highlight the biological mechanisms underlying strong humoral and cell-mediated immune responses and explore materials design strategies to manipulate and control these mechanisms.

Increased B Cell Selection Stringency In Germinal Centers Can Explain Improved COVID-19 Vaccine Efficacies With Low Dose Prime or Delayed Boost

The efficacy of COVID-19 vaccines appears to depend in complex ways on the vaccine dosage and the interval between the prime and boost doses. Unexpectedly, lower dose prime and longer prime-boost intervals have yielded higher efficacies in clinical trials. To elucidate the origins of these effects, we developed a stochastic simulation model of the germinal center (GC) reaction and predicted the antibody responses elicited by different vaccination protocols. The simulations predicted that a lower dose prime could increase the selection stringency in GCs due to reduced antigen availability, resulting in the selection of GC B cells with higher affinities for the target antigen. The boost could relax this selection stringency and allow the expansion of the higher affinity GC B cells selected, improving the overall response. With a longer dosing interval, the decay in the antigen with time following the prime could further increase the selection stringency, amplifying this effect. The effect remained in our simulations even when new GCs following the boost had to be seeded by memory B cells formed following the prime. These predictions offer a plausible explanation of the observed paradoxical effects of dosage and dosing interval on vaccine efficacy. Tuning the selection stringency in the GCs using prime-boost dosages and dosing intervals as handles may help improve vaccine efficacies.

Structure-guided envelope trimer design in HIV-1 vaccine development: a narrative review

INTRODUCTION

The development of a human immunodeficiency virus 1 (HIV-1) vaccine remains a formidable challenge. An effective vaccine likely requires the induction of broadly neutralizing antibodies (bNAbs), which likely involves the use of native-like HIV-1 envelope (Env) trimers at some or all stages of vaccination. Development of such trimers has been very difficult, but much progress has been made in the past decade, starting with the BG505 SOSIP trimer, elucidation of its atomic structure and implementing subsequent design iterations. This progress facilitated understanding the weaknesses of the Env trimer, fuelled structure-guided HIV-1 vaccine design and assisted in the development of new vaccine designs. This review summarizes the relevant literature focusing on studies using structural biology to reveal and define HIV-1 Env sites of vulnerability; to improve Env trimers, by creating more stable versions; understanding antibody responses in preclinical vaccination studies at the atomic level; understanding the glycan shield; and to improve "on-target" antibody responses versus "off-target" responses.

METHODS

The authors conducted a narrative review of recently published articles that made a major contribution to HIV-1 structural biology and vaccine design efforts between the years 2000 and 2021.

DISCUSSION

The field of structural biology is evolving at an unprecedented pace, where cryo-electron microscopy (cryo-EM) and X-ray crystallography provide complementary information. Resolving protein structures is necessary for defining which Env surfaces are accessible for the immune system and can be targeted by neutralizing antibodies. Recently developed techniques, such as electron microscopy-based polyclonal epitope mapping (EMPEM) are revolutionizing the way we are analysing immune responses and shed light on the immunodominant targets on new vaccine immunogens. Such information accelerates iterative vaccine design; for example, by reducing undesirable off-target responses, while improving immunogens to drive the more desirable on-target responses.

CONCLUSIONS

Resolving high-resolution structures of the HIV-1 Env trimer was instrumental in understanding and improving recombinant HIV-1 Env trimers that mimic the structure of viral HIV-1 Env spikes. Newly emerging techniques in structural biology are aiding vaccine design efforts and improving immunogens. The role of structural biology in HIV-1 vaccine design has indeed become very prominent and is unlikely to diminish any time soon.

System-Level Scenarios for the Elucidation of T Cell-Mediated Germinal Center B Cell Differentiation

Germinal center (GC) reactions are vital to the correct functioning of the adaptive immune system, through formation of high affinity, class switched antibodies. GCs are transient anatomical structures in secondary lymphoid organs where specific B cells, after recognition of antigen and with T cell help, undergo class switching. Subsequently, B cells cycle between zones of proliferation and somatic hypermutation and zones where renewed antigen acquisition and T cell help allows for selection of high affinity B cells (affinity maturation). Eventually GC B cells first differentiate into long-lived memory B cells (MBC) and finally into plasma cells (PC) that partially migrate to the bone marrow to encapsulate into long-lived survival niches. The regulation of GC reactions is a highly dynamically coordinated process that occurs between various cells and molecules that change in their signals. Here, we present a system-level perspective of T cell-mediated GC B cell differentiation, presenting and discussing the experimental and computational efforts on the regulation of the GCs. We aim to integrate Systems Biology with B cell biology, to advance elucidation of the regulation of high-affinity, class switched antibody formation, thus to shed light on the delicate functioning of the adaptive immune system. Specifically, we: i) review experimental findings of internal and external factors driving various GC dynamics, such as GC initiation, maturation and GCBC fate determination; ii) draw comparisons between experimental observations and mathematical modeling investigations; and iii) discuss and reflect on current strategies of modeling efforts, to elucidate B cell behavior during the GC tract. Finally, perspectives are specifically given on to the areas where a Systems Biology approach may be useful to predict novel GCBC-T cell interaction dynamics.

Ymir: A 3D structural affinity model for multi-epitope vaccine simulations

Vaccine development is challenged by the hierarchy of immunodominance between target antigen epitopes and the emergence of antigenic variants by pathogen mutation. The strength and breadth of antibody responses relies on selection and mutation in the germinal center and on the structural similarity between antigens. Computational methods for assessing the breadth of germinal center responses to multivalent antigens are critical to speed up vaccine development. Yet, such methods have poorly reflected the 3D antigen structure and antibody breadth. Here, we present Ymir, a new 3D-lattice-based framework that calculates in silico antibody-antigen affinities. Key physiological properties naturally emerge from Ymir such as affinity jumps, cross-reactivity, and differential epitope accessibility. We validated Ymir by replicating known features of germinal center dynamics. We show that combining antigens with mutated but structurally related epitopes enhances vaccine breadth. Ymir opens a new avenue for understanding vaccine potency based on the structural relationship between vaccine antigens.

Binding affinity landscapes constrain the evolution of broadly neutralizing anti-influenza antibodies

Over the past two decades, several broadly neutralizing antibodies (bnAbs) that confer protection against diverse influenza strains have been isolated. Structural and biochemical characterization of these bnAbs has provided molecular insight into how they bind distinct antigens. However, our understanding of the evolutionary pathways leading to bnAbs, and thus how best to elicit them, remains limited. Here, we measure equilibrium dissociation constants of combinatorially complete mutational libraries for two naturally isolated influenza bnAbs (CR9114, 16 heavy-chain mutations; CR6261, 11 heavy-chain mutations), reconstructing all possible evolutionary intermediates back to the unmutated germline sequences. We find that these two libraries exhibit strikingly different patterns of breadth: while many variants of CR6261 display moderate affinity to diverse antigens, those of CR9114 display appreciable affinity only in specific, nested combinations. By examining the extensive pairwise and higher order epistasis between mutations, we find key sites with strong synergistic interactions that are highly similar across antigens for CR6261 and different for CR9114. Together, these features of the binding affinity landscapes strongly favor sequential acquisition of affinity to diverse antigens for CR9114, while the acquisition of breadth to more similar antigens for CR6261 is less constrained. These results, if generalizable to other bnAbs, may explain the molecular basis for the widespread observation that sequential exposure favors greater breadth, and such mechanistic insight will be essential for predicting and eliciting broadly protective immune responses.

A molecular theory of germinal center B cell selection and division

The selection of B cells (BCs) in germinal centers (GCs) is pivotal to the generation of high-affinity antibodies and memory BCs, but it lacks global understanding. Based on the idea of a single Tfh-cell signal that controls BC selection and division, experiments appear contradictory. Here, we use the current knowledge on the molecular pathways of GC BCs to develop a theory of GC BC selection and division based on the dynamics of molecular factors. This theory explains the seemingly contradictory experiments by the separation of signals for BC fate decision from signals controlling the number of BC divisions. Three model variants are proposed and experiments are predicted that allow one to distinguish those. Understanding information processing in molecular BC states is critical for targeted immune interventions, and the proposed theory implies that selection and division can be controlled independently in GC reactions.

Coevolutionary transitions emerging from flexible molecular recognition and eco-evolutionary feedback

Highly mutable viruses evolve to evade host immunity that exerts selective pressure and adapts to viral dynamics. Here, we provide a framework for identifying key determinants of the mode and fate of viral-immune coevolution by linking molecular recognition and eco-evolutionary dynamics. We find that conservation level and initial diversity of antigen jointly determine the timing and efficacy of narrow and broad antibody responses, which in turn control the transition between viral persistence, clearance, and rebound. In particular, clearance of structurally complex antigens relies on antibody evolution in a larger antigenic space than where selection directly acts; viral rebound manifests binding-mediated feedback between ecology and rapid evolution. Finally, immune compartmentalization can slow viral escape but also delay clearance. This work suggests that flexible molecular binding allows a plastic phenotype that exploits potentiating neutral variations outside direct contact, opening new and shorter paths toward highly adaptable states.

•

A scale-crossing framework identifies key determinants of viral-immune coevolution

•

Fast specific response influences slow broad response by shaping antigen dynamics

•

Antibody footprint shift enables breadth acquisition and viral clearance

•

Model explains divergent kinetics and outcomes of HCV infection in humans

Administration of Multivalent Influenza Virus Recombinant Hemagglutinin Vaccine in Combination-Adjuvant Elicits Broad Reactivity Beyond the Vaccine Components

Combining variant antigens into a multivalent vaccine is a traditional approach used to provide broad coverage against antigenically variable pathogens, such as polio, human papilloma and influenza viruses. However, strategies for increasing the breadth of antibody coverage beyond the vaccine are not well understood, but may provide more anticipatory protection. Influenza virus hemagglutinin (HA) is a prototypic variant antigen. Vaccines that induce HA-specific neutralizing antibodies lose efficacy as amino acid substitutions accumulate in neutralizing epitopes during influenza virus evolution. Here we studied the effect of a potent combination adjuvant (CpG/MPLA/squalene-in-water emulsion) on the breadth and maturation of the antibody response to a representative variant of HA subtypes H1, H5 and H7. Using HA protein microarrays and antigen-specific B cell labelling, we show when administered individually, each HA elicits a cross-reactive antibody profile for multiple variants within the same subtype and other closely-related subtypes (homosubtypic and heterosubtypic cross-reactivity, respectively). Despite a capacity for each subtype to induce heterosubtypic cross-reactivity, broader coverage was elicited by simply combining the subtypes into a multivalent vaccine. Importantly, multiplexing did not compromise antibody avidity or affinity maturation to the individual HA constituents. The use of adjuvants to increase the breadth of antibody coverage beyond the vaccine antigens may help future-proof vaccines against newly-emerging variants.

In Silico Analysis of the Longevity and Timeline of Individual Germinal Center Reactions in a Primary Immune Response

Germinal centers (GCs) are transient structures in the secondary lymphoid organs, where B cells undergo affinity maturation to produce high affinity memory and plasma cells. The lifetime of GC responses is a critical factor limiting the extent of affinity maturation and efficiency of antibody responses. While the average lifetime of overall GC reactions in a lymphoid organ is determined experimentally, the lifetime of individual GCs has not been monitored due to technical difficulties in longitudinal analysis. In silico analysis of the contraction phase of GC responses towards primary immunization with sheep red blood cells suggested that if individual GCs had similar lifetimes, the data would be consistent only when new GCs were formed until a very late phase after immunization. Alternatively, there could be a large variation in the lifetime of individual GCs suggesting that both long and short-lived GCs might exist in the same lymphoid organ. Simulations predicted that such differences in the lifetime of GCs could arise due to variations in antigen availability and founder cell composition. These findings identify the potential factors limiting GC lifetime and contribute to an understanding of overall GC responses from the perspective of individual GCs in a primary immune response.

Biomaterials and nanomaterials for sustained release vaccine delivery

Vaccines are considered one of the most significant medical advancements in human history, as they have prevented hundreds of millions of deaths since their discovery; however, modern travel permits disease spread at unprecedented rates, and vaccine shortcomings like thermal sensitivity and required booster shots have been made evident by the COVID-19 pandemic. Approaches to overcoming these issues appear promising via the integration of vaccine technology with biomaterials, which offer sustained-release properties and preserve proteins, prevent conformational changes, and enable storage at room temperature. Sustained release and thermal stabilization of therapeutic biomacromolecules is an emerging area that integrates material science, chemistry, immunology, nanotechnology, and pathology to investigate different biocompatible materials. Biomaterials, including natural sugar polymers, synthetic polyesters produced from biologically derived monomers, hydrogel blends, protein-polymer blends, and metal-organic frameworks, have emerged as early players in the field. This overview will focus on significant advances of sustained release biomaterial in the context of vaccines against infectious disease and the progress made towards thermally stable "single-shot" formulations. This article is categorized under: Therapeutic Approaches and Drug Discovery > Nanomedicine for Infectious Disease.

Antigenic waves of virus-immune coevolution

The evolution of many microbes and pathogens, including circulating viruses such as seasonal influenza, is driven by immune pressure from the host population. In turn, the immune systems of infected populations get updated, chasing viruses even farther away. Quantitatively understanding how these dynamics result in observed patterns of rapid pathogen and immune adaptation is instrumental to epidemiological and evolutionary forecasting. Here we present a mathematical theory of coevolution between immune systems and viruses in a finite-dimensional antigenic space, which describes the cross-reactivity of viral strains and immune systems primed by previous infections. We show the emergence of an antigenic wave that is pushed forward and canalized by cross-reactivity. We obtain analytical results for shape, speed, and angular diffusion of the wave. In particular, we show that viral-immune coevolution generates an emergent timescale, the persistence time of the wave's direction in antigenic space, which can be much longer than the coalescence time of the viral population. We compare these dynamics to the observed antigenic turnover of influenza strains, and we discuss how the dimensionality of antigenic space impacts the predictability of the evolutionary dynamics. Our results provide a concrete and tractable framework to describe pathogen-host coevolution.

Roadmap on biology in time varying environments

Biological organisms experience constantly changing environments, from sudden changes in physiology brought about by feeding, to the regular rising and setting of the Sun, to ecological changes over evolutionary timescales. Living organisms have evolved to thrive in this changing world but the general principles by which organisms shape and are shaped by time varying environments remain elusive. Our understanding is particularly poor in the intermediate regime with no separation of timescales, where the environment changes on the same timescale as the physiological or evolutionary response. Experiments to systematically characterize the response to dynamic environments are challenging since such environments are inherently high dimensional. This roadmap deals with the unique role played by time varying environments in biological phenomena across scales, from physiology to evolution, seeking to emphasize the commonalities and the challenges faced in this emerging area of research.

Survival probability and size of lineages in antibody affinity maturation

Affinity maturation (AM) is the process through which the immune system is able to develop potent antibodies against new pathogens it encounters, and is at the base of the efficacy of vaccines. At its core AM is analogous to a Darwinian evolutionary process, where B cells mutate and are selected on the base of their affinity for an antigen (Ag), and Ag availability tunes the selective pressure. In cases when this selective pressure is high, the number of B cells might quickly decrease and the population might risk extinction in what is known as a population bottleneck. Here we study the probability for a B-cell lineage to survive this bottleneck scenario as a function of the progenitor affinity for the Ag. Using recursive relations and probability generating functions we derive expressions for the average extinction time and progeny size for lineages that go extinct. We then extend our results to the full population, both in the absence and presence of competition for T-cell help, and quantify the population survival probability as a function of Ag concentration and initial population size. Our study suggests the population bottleneck phenomenology might represent a limit case in the space of biologically plausible maturation scenarios, whose characterization could help guide the process of vaccine development.

Mechanisms underlying vaccination protocols that may optimally elicit broadly neutralizing antibodies against highly mutable pathogens

Effective prophylactic vaccines usually induce the immune system to generate potent antibodies that can bind to an antigen and thus prevent it from infecting host cells. B cells produce antibodies by a Darwinian evolutionary process called affinity maturation (AM). During AM, the B cell population evolves in response to the antigen to produce antibodies that bind specifically and strongly to the antigen. Highly mutable pathogens pose a major challenge to the development of effective vaccines because antibodies that are effective against one strain of the virus may not protect against a mutant strain. Antibodies that can protect against diverse strains of a mutable pathogen have high "breadth" and are called broadly neutralizing antibodies (bnAbs). In spite of extensive studies, an effective vaccination strategy that can generate bnAbs in humans does not exist for any highly mutable pathogen. Here we study a minimal model to explore the mechanisms underlying how the selection forces imposed by antigens can be optimally chosen to guide AM to maximize the evolution of bnAbs. For logistical reasons, only a finite number of antigens can be administered in a finite number of vaccinations; that is, guiding the nonequilibrium dynamics of AM to produce bnAbs must be accomplished nonadiabatically. The time-varying Kullback-Leibler divergence (KLD) between the existing B cell population distribution and the fitness landscape imposed by antigens is a quantitative metric of the thermodynamic force acting on B cells. If this force is too small, adaptation is minimal. If the force is too large, contrary to expectations, adaptation is not faster; rather, the B cell population is extinguished for reasons that we describe. We define the conditions necessary for the force to be set optimally such that the flux of B cells from low to high breadth states is maximized. Even in this case we show why the dynamics of AM prevent perfect adaptation. If two shots of vaccination are allowed, the optimal protocol is characterized by a relatively low optimal KLD during the first shot that appropriately increases the diversity of the B cell population so that the surviving B cells have a high chance of evolving into bnAbs upon subsequently increasing the KLD during the second shot. Phylogenetic tree analysis further reveals the evolutionary pathways that lead to bnAbs. The connections between the mechanisms revealed by our analyses and recent simulation studies of bnAb evolution, the problem of generalist versus specialist evolution, and learning theory are discussed.