An explainable language model for antibody specificity prediction using curated influenza hemagglutinin antibodies.

Despite decades of antibody research, it remains challenging to predict the specificity of an antibody solely based on its sequence. Two major obstacles are the lack of appropriate models and the inaccessibility of datasets for model training. In this study, we curated >5,000 influenza hemagglutinin (HA) antibodies by mining research publications and patents, which revealed many distinct sequence features between antibodies to HA head and stem domains. We then leveraged this dataset to develop a lightweight memory B cell language model (mBLM) for sequence-based antibody specificity prediction. Model explainability analysis showed that mBLM could identify key sequence features of HA stem antibodies. Additionally, by applying mBLM to HA antibodies with unknown epitopes, we discovered and experimentally validated many HA stem antibodies. Overall, this study not only advances our molecular understanding of the antibody response to the influenza virus but also provides a valuable resource for applying deep learning to antibody research.

Seasonal influenza A virus lineages exhibit divergent abilities to antagonize interferon induction and signaling.

The circulation of seasonal influenza A viruses (IAVs) in humans relies on effective evasion and subversion of the host immune response. While the evolution of seasonal H1N1 and H3N2 viruses to avoid humoral immunity is well characterized, relatively little is known about the evolution of innate immune antagonism phenotypes in these viruses. Numerous studies have established that only a small subset of infected cells is responsible for initiating the type I and type III interferon (IFN) response during IAV infection, emphasizing the importance of single cell studies to accurately characterize the IFN response during infection. We developed a flow cytometry-based method to examine transcriptional changes in IFN and interferon stimulated gene (ISG) expression at the single cell level. We observed that NS segments derived from seasonal H3N2 viruses are more efficient at antagonizing IFN signaling but less effective at suppressing IFN induction, compared to the pdm2009 H1N1 lineage. We compared a collection of NS segments spanning the natural history of the current seasonal IAV lineages and demonstrate long periods of stability in IFN antagonism potential, punctuated by occasional phenotypic shifts. Altogether, our data reveal significant differences in how seasonal and pandemic H1N1 and H3N2 viruses antagonize the human IFN response at the single cell level.

Epistasis mediates the evolution of the receptor binding mode in recent human H3N2 hemagglutinin.

The receptor-binding site of influenza A virus hemagglutinin partially overlaps with major antigenic sites and constantly evolves. In this study, we observe that mutations G186D and D190N in the hemagglutinin receptor-binding site have coevolved in two recent human H3N2 clades. X-ray crystallography results show that these mutations coordinately drive the evolution of the hemagglutinin receptor binding mode. Epistasis between G186D and D190N is further demonstrated by glycan binding and thermostability analyses. Immunization and neutralization experiments using mouse and human samples indicate that the evolution of receptor binding mode is accompanied by a change in antigenicity. Besides, combinatorial mutagenesis reveals that G186D and D190N, along with other natural mutations in recent H3N2 strains, alter the compatibility with a common egg-adaptive mutation in seasonal influenza vaccines. Overall, our findings elucidate the role of epistasis in shaping the recent evolution of human H3N2 hemagglutinin and substantiate the high evolvability of its receptor-binding mode.

Natural variation in neuraminidase activity influences the evolutionary potential of the seasonal H1N1 lineage hemagglutinin.

The antigenic evolution of the influenza A virus hemagglutinin (HA) gene poses a major challenge for the development of vaccines capable of eliciting long-term protection. Prior efforts to understand the mechanisms that govern viral antigenic evolution mainly focus on HA in isolation, ignoring the fact that HA must act in concert with the viral neuraminidase (NA) during replication and spread. Numerous studies have demonstrated that the degree to which the receptor-binding avidity of HA and receptor-cleaving activity of NA are balanced with each other influences overall viral fitness. We recently showed that changes in NA activity can significantly alter the mutational fitness landscape of HA in the context of a lab-adapted virus strain. Here, we test whether natural variation in relative NA activity can influence the evolutionary potential of HA in the context of the seasonal H1N1 lineage (pdmH1N1) that has circulated in humans since the 2009 pandemic. We observed substantial variation in the relative activities of NA proteins encoded by a panel of H1N1 vaccine strains isolated between 2009 and 2019. We comprehensively assessed the effect of NA background on the HA mutational fitness landscape in the circulating pdmH1N1 lineage using deep mutational scanning and observed pronounced epistasis between NA and residues in or near the receptor-binding site of HA. To determine whether NA variation could influence the antigenic evolution of HA, we performed neutralizing antibody selection experiments using a panel of monoclonal antibodies targeting different HA epitopes. We found that the specific antibody escape profiles of HA were highly contingent upon NA background. Overall, our results indicate that natural variation in NA activity plays a significant role in governing the evolutionary potential of HA in the currently circulating pdmH1N1 lineage.

Deep mutational scanning of influenza A virus NEP reveals pleiotropic mutations in its N-terminal domain.

The influenza A virus nuclear export protein (NEP) is a multifunctional protein that is essential for the viral life cycle and has very high sequence conservation. However, since the open reading frame of NEP largely overlaps with that of another influenza viral protein, non-structural protein 1, it is difficult to infer the functional constraints of NEP based on sequence conservation analysis. Besides, the N-terminal of NEP is structurally disordered, which further complicates the understanding of its function. Here, we systematically measured the replication fitness effects of >1,800 mutations of NEP. Our results show that the N-terminal domain has high mutational tolerance. Additional experiments demonstrate that N-terminal domain mutations pleiotropically affect viral transcription and replication dynamics, host cellular responses, and mammalian adaptation of avian influenza virus. Overall, our study not only advances the functional understanding of NEP, but also provides insights into its evolutionary constraints.

Functional and antigenic characterization of SARS-CoV-2 spike fusion peptide by deep mutational scanning.

The fusion peptide of SARS-CoV-2 spike protein is functionally important for membrane fusion during virus entry and is part of a broadly neutralizing epitope. However, sequence determinants at the fusion peptide and its adjacent regions for pathogenicity and antigenicity remain elusive. In this study, we perform a series of deep mutational scanning (DMS) experiments on an S2 region spanning the fusion peptide of authentic SARS-CoV-2 in different cell lines and in the presence of broadly neutralizing antibodies. We identify mutations at residue 813 of the spike protein that reduced TMPRSS2-mediated entry with decreased virulence. In addition, we show that an F823Y mutation, present in bat betacoronavirus HKU9 spike protein, confers resistance to broadly neutralizing antibodies. Our findings provide mechanistic insights into SARS-CoV-2 pathogenicity and also highlight a potential challenge in developing broadly protective S2-based coronavirus vaccines.

Evidence of antigenic drift in the fusion machinery core of SARS-CoV-2 spike.

Antigenic drift of SARS-CoV-2 is typically defined by mutations in the N-terminal domain and receptor binding domain of spike protein. In contrast, whether antigenic drift occurs in the S2 domain remains largely elusive. Here, we perform a deep mutational scanning experiment to identify S2 mutations that affect binding of SARS-CoV-2 spike to three S2 apex public antibodies. Our results indicate that spatially diverse mutations, including D950N and Q954H, which are observed in Delta and Omicron variants, respectively, weaken the binding of spike to these antibodies. Although S2 apex antibodies are known to be nonneutralizing, we show that they confer protection in vivo through Fc-mediated effector functions. Overall, this study indicates that the S2 domain of SARS-CoV-2 spike can undergo antigenic drift, which represents a potential challenge for the development of more universal coronavirus vaccines.

Natural variation in neuraminidase activity influences the evolutionary potential of the seasonal H1N1 lineage hemagglutinin.

The antigenic evolution of the influenza A virus hemagglutinin (HA) gene poses a major challenge for the development of vaccines capable of eliciting long-term protection. Prior efforts to understand the mechanisms that govern viral antigenic evolution mainly focus on HA in isolation, ignoring the fact that HA must act in concert with the viral neuraminidase (NA) during replication and spread. Numerous studies have demonstrated that the degree to which the receptor binding avidity of HA and receptor cleaving activity of NA are balanced with each other influences overall viral fitness. We recently showed that changes in NA activity can significantly alter the mutational fitness landscape of HA in the context of a lab-adapted virus strain. Here, we test whether natural variation in relative NA activity can influence the evolutionary potential of HA in the context of the seasonal H1N1 lineage (pdmH1N1) that has circulated in humans since the 2009 pandemic. We observed substantial variation in the relative activities of NA proteins encoded by a panel of H1N1 vaccine strains isolated between 2009 and 2019. We comprehensively assessed the effect of NA background on the HA mutational fitness landscape in the circulating pdmH1N1 lineage using deep mutational scanning and observed pronounced epistasis between NA and residues in or near the receptor binding site of HA. To determine whether NA variation could influence the antigenic evolution of HA, we performed neutralizing antibody selection experiments using a panel of monoclonal antibodies targeting different HA epitopes. We found that the specific antibody escape profiles of HA were highly contingent upon NA background. Overall, our results indicate that natural variation in NA activity plays a significant role in governing the evolutionary potential of HA in the currently circulating pdmH1N1 lineage.

Functional and antigenic characterization of SARS-CoV-2 spike fusion peptide by deep mutational scanning.

The fusion peptide of SARS-CoV-2 spike protein is functionally important for membrane fusion during virus entry and is part of a broadly neutralizing epitope. However, sequence determinants at the fusion peptide and its adjacent regions for pathogenicity and antigenicity remain elusive. In this study, we performed a series of deep mutational scanning (DMS) experiments on an S2 region spanning the fusion peptide of authentic SARS-CoV-2 in different cell lines and in the presence of broadly neutralizing antibodies. We identified mutations at residue 813 of the spike protein that reduced TMPRSS2-mediated entry with decreased virulence. In addition, we showed that an F823Y mutation, present in bat betacoronavirus HKU9 spike protein, confers resistance to broadly neutralizing antibodies. Our findings provide mechanistic insights into SARS-CoV-2 pathogenicity and also highlight a potential challenge in developing broadly protective S2-based coronavirus vaccines.

Stringent and complex sequence constraints of an IGHV1-69 broadly neutralizing antibody to influenza HA stem.

IGHV1-69 is frequently utilized by broadly neutralizing influenza antibodies to the hemagglutinin (HA) stem. These IGHV1-69 HA stem antibodies have diverse complementarity-determining region (CDR) H3 sequences. Besides, their light chains have minimal to no contact with the epitope. Consequently, sequence determinants that confer IGHV1-69 antibodies with HA stem specificity remain largely elusive. Using high-throughput experiments, this study reveals the importance of light-chain sequence for the IGHV1-69 HA stem antibody CR9114, which is the broadest influenza antibody known to date. Moreover, we demonstrate that the CDR H3 sequences from many other IGHV1-69 antibodies, including those to the HA stem, are incompatible with CR9114. Along with mutagenesis and structural analysis, our results indicate that light-chain and CDR H3 sequences coordinately determine the HA stem specificity of IGHV1-69 antibodies. Overall, this work provides molecular insights into broadly neutralizing antibody responses to influenza virus, which have important implications for universal influenza vaccine development.

Leveraging vaccination-induced protective antibodies to define conserved epitopes on influenza N2 neuraminidase.

There is growing appreciation for neuraminidase (NA) as an influenza vaccine target; however, its antigenicity remains poorly characterized. In this study, we isolated three broadly reactive N2 antibodies from the plasmablasts of a single vaccinee, including one that cross-reacts with NAs from seasonal H3N2 strains spanning five decades. Although these three antibodies have diverse germline usages, they recognize similar epitopes that are distant from the NA active site and instead involve the highly conserved underside of NA head domain. We also showed that all three antibodies confer prophylactic and therapeutic protection in vivo, due to both Fc effector functions and NA inhibition through steric hindrance. Additionally, the contribution of Fc effector functions to protection in vivo inversely correlates with viral growth inhibition activity in vitro. Overall, our findings advance the understanding of NA antibody response and provide important insights into the development of a broadly protective influenza vaccine.

Widespread impact of immunoglobulin V-gene allelic polymorphisms on antibody reactivity.

The ability of the human immune system to generate antibodies to any given antigen can be strongly influenced by immunoglobulin V-gene allelic polymorphisms. However, previous studies have provided only limited examples. Therefore, the prevalence of this phenomenon has been unclear. By analyzing >1,000 publicly available antibody-antigen structures, we show that many V-gene allelic polymorphisms in antibody paratopes are determinants for antibody binding activity. Biolayer interferometry experiments further demonstrate that paratope allelic polymorphisms on both heavy and light chains often abolish antibody binding. We also illustrate the importance of minor V-gene allelic polymorphisms with low frequency in several broadly neutralizing antibodies to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and influenza virus. Overall, this study not only highlights the pervasive impact of V-gene allelic polymorphisms on antibody binding but also provides mechanistic insights into the variability of antibody repertoires across individuals, which in turn have important implications for vaccine development and antibody discovery.

An explainable language model for antibody specificity prediction using curated influenza hemagglutinin antibodies.

Despite decades of antibody research, it remains challenging to predict the specificity of an antibody solely based on its sequence. Two major obstacles are the lack of appropriate models and inaccessibility of datasets for model training. In this study, we curated a dataset of >5,000 influenza hemagglutinin (HA) antibodies by mining research publications and patents, which revealed many distinct sequence features between antibodies to HA head and stem domains. We then leveraged this dataset to develop a lightweight memory B cell language model (mBLM) for sequence-based antibody specificity prediction. Model explainability analysis showed that mBLM captured key sequence motifs of HA stem antibodies. Additionally, by applying mBLM to HA antibodies with unknown epitopes, we discovered and experimentally validated many HA stem antibodies. Overall, this study not only advances our molecular understanding of antibody response to influenza virus, but also provides an invaluable resource for applying deep learning to antibody research.

Stringent and complex sequence constraints of an IGHV1-69 broadly neutralizing antibody to influenza HA stem.

IGHV1-69 is frequently utilized by broadly neutralizing influenza antibodies to the hemagglutinin (HA) stem. These IGHV1-69 HA stem antibodies have diverse complementarity-determining region (CDR) H3 sequences. Besides, their light chains have minimal to no contact with the epitope. Consequently, sequence determinants that confer IGHV1-69 antibodies with HA stem specificity remain largely elusive. Using high-throughput experiments, this study revealed the importance of light chain sequence for the IGHV1-69 HA stem antibody CR9114, which is the broadest influenza antibody known to date. Moreover, we demonstrated that the CDR H3 sequences from many other IGHV1-69 antibodies, including those to HA stem, were incompatible with CR9114. Along with mutagenesis and structural analysis, our results indicate that light chain and CDR H3 sequences coordinately determine the HA stem specificity of IGHV1-69 antibodies. Overall, this work provides molecular insights into broadly neutralizing antibody responses to influenza virus, which have important implications for universal influenza vaccine development.

Widespread impact of immunoglobulin V gene allelic polymorphisms on antibody reactivity.

The ability of human immune system to generate antibodies to any given antigen can be strongly influenced by immunoglobulin V gene (IGV) allelic polymorphisms. However, previous studies have provided only a limited number of examples. Therefore, the prevalence of this phenomenon has been unclear. By analyzing >1,000 publicly available antibody-antigen structures, we show that many IGV allelic polymorphisms in antibody paratopes are determinants for antibody binding activity. Biolayer interferometry experiment further demonstrates that paratope allelic mutations on both heavy and light chain often abolish antibody binding. We also illustrate the importance of minor IGV allelic variants with low frequency in several broadly neutralizing antibodies to SARS-CoV-2 and influenza virus. Overall, this study not only highlights the pervasive impact of IGV allelic polymorphisms on antibody binding, but also provides mechanistic insights into the variability of antibody repertoires across individuals, which in turn have important implications for vaccine development and antibody discovery.

Caprin-1 binding to the critical stress granule protein G3BP1 is influenced by pH.

G3BP is the central node within stress-induced protein-RNA interaction networks known as stress granules (SGs). The SG-associated proteins Caprin-1 and USP10 bind mutually exclusively to the NTF2 domain of G3BP1, promoting and inhibiting SG formation, respectively. Herein, we present the crystal structure of G3BP1-NTF2 in complex with a Caprin-1-derived short linear motif (SLiM). Caprin-1 interacts with His-31 and His-62 within a third NTF2-binding site outside those covered by USP10, as confirmed using biochemical and biophysical-binding assays. Nano-differential scanning fluorimetry revealed reduced thermal stability of G3BP1-NTF2 at acidic pH. This destabilization was counterbalanced significantly better by bound USP10 than Caprin-1. The G3BP1/USP10 complex immunoprecipated from human U2OS cells was more resistant to acidic buffer washes than G3BP1/Caprin-1. Acidification of cellular condensates by approximately 0.5 units relative to the cytosol was detected by ratiometric fluorescence analysis of pHluorin2 fused to G3BP1. Cells expressing a Caprin-1/FGDF chimera with higher G3BP1-binding affinity had reduced Caprin-1 levels and slightly reduced condensate sizes. This unexpected finding may suggest that binding of the USP10-derived SLiM to NTF2 reduces the propensity of G3BP1 to enter condensates.

High-throughput identification of prefusion-stabilizing mutations in SARS-CoV-2 spike.

Designing prefusion-stabilized SARS-CoV-2 spike is critical for the effectiveness of COVID-19 vaccines. All COVID-19 vaccines in the US encode spike with K986P/V987P mutations to stabilize its prefusion conformation. However, contemporary methods on engineering prefusion-stabilized spike immunogens involve tedious experimental work and heavily rely on structural information. Here, we establish a systematic and unbiased method of identifying mutations that concomitantly improve expression and stabilize the prefusion conformation of the SARS-CoV-2 spike. Our method integrates a fluorescence-based fusion assay, mammalian cell display technology, and deep mutational scanning. As a proof-of-concept, we apply this method to a region in the S2 domain that includes the first heptad repeat and central helix. Our results reveal that besides K986P and V987P, several mutations simultaneously improve expression and significantly lower the fusogenicity of the spike. As prefusion stabilization is a common challenge for viral immunogen design, this work will help accelerate vaccine development against different viruses.

Mutational fitness landscape of human influenza H3N2 neuraminidase.

Influenza neuraminidase (NA) has received increasing attention as an effective vaccine target. However, its mutational tolerance is not well characterized. Here, the fitness effects of >6,000 mutations in human H3N2 NA are probed using deep mutational scanning. Our result shows that while its antigenic regions have high mutational tolerance, there are solvent-exposed regions with low mutational tolerance. We also find that protein stability is a major determinant of NA mutational fitness. The deep mutational scanning result correlates well with mutational fitness inferred from natural sequences using a protein language model, substantiating the relevance of our findings to the natural evolution of circulating strains. Additional analysis further suggests that human H3N2 NA is far from running out of mutations despite already evolving for >50 years. Overall, this study advances our understanding of the evolutionary potential of NA and the underlying biophysical constraints, which in turn provide insights into NA-based vaccine design.

Molecular analysis of a public cross-neutralizing antibody response to SARS-CoV-2.

As severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concerns (VOCs) continue to emerge, cross-neutralizing antibody responses become key toward next-generation design of a more universal COVID-19 vaccine. By analyzing published data from the literature, we report here that the combination of germline genes IGHV2-5/IGLV2-14 represents a public antibody response to the receptor-binding domain (RBD) that potently cross-neutralizes a broad range of VOCs, including Omicron and its sub-lineages. Detailed molecular analysis shows that the complementarity-determining region H3 sequences of IGHV2-5/IGLV2-14-encoded RBD antibodies have a preferred length of 11 amino acids and a conserved HxIxxI motif. In addition, these antibodies have a strong allelic preference due to an allelic polymorphism at amino acid residue 54 of IGHV2-5, which is located at the paratope. These findings have important implications for understanding cross-neutralizing antibody responses to SARS-CoV-2 and its heterogenicity at the population level as well as the development of a universal COVID-19 vaccine.