A complex epistatic network limits the mutational reversibility in the influenza hemagglutinin receptor-binding site

Predicting pathogen evolution and immune evasion in the age of artificial intelligence

The genomic diversification of viral pathogens during viral epidemics and pandemics represents a major adaptive route for infectious agents to circumvent therapeutic and public health initiatives. Historically, strategies to address viral evolution have relied on responding to emerging variants after their detection, leading to delays in effective public health responses. Because of this, a long-standing yet challenging objective has been to forecast viral evolution by predicting potentially harmful viral mutations prior to their emergence. The promises of artificial intelligence (AI) coupled with the exponential growth of viral data collection infrastructures spurred by the COVID-19 pandemic, have resulted in a research ecosystem highly conducive to this objective. Due to the COVID-19 pandemic accelerating the development of pandemic mitigation and preparedness strategies, many of the methods discussed here were designed in the context of SARS-CoV-2 evolution. However, most of these pipelines were intentionally designed to be adaptable across RNA viruses, with several strategies already applied to multiple viral species. In this review, we explore recent breakthroughs that have facilitated the forecasting of viral evolution in the context of an ongoing pandemic, with particular emphasis on deep learning architectures, including the promising potential of language models (LM). The approaches discussed here employ strategies that leverage genomic, epidemiologic, immunologic and biological information.

T4 Phage Displaying Dual Antigen Clusters Against H3N2 Influenza Virus Infection

BACKGROUND

The current H3N2 influenza subunit vaccine exhibits weak immunogenicity, which limits its effectiveness in preventing and controlling influenza virus infections.

METHODS

In this study, we aimed to develop a T4 phage-based nanovaccine designed to enhance the immunogenicity of two antigens by displaying the HA1 and M2e antigens of the H3N2 influenza virus on each phage nanoparticle. Specifically, we fused the Soc protein with the HA1 antigen and the Hoc protein with the M2e antigen, assembling them onto a T4 phage that lacks Soc and Hoc proteins (Soc-Hoc-T4), thereby constructing a nanovaccine that concurrently presents both HA1 and M2e antigens.

RESULTS

The analysis of the optical density of the target protein bands indicated that each particle could display approximately 179 HA1 and 68 M2e antigen molecules. Additionally, animal experiments demonstrated that this nanoparticle vaccine displaying dual antigen clusters induced a stronger specific immune response, higher antibody titers, a more balanced Th1/Th2 immune response, and enhanced CD4+ and CD8+ T cell effects compared to immunization with HA1 and M2e antigen molecules alone. Importantly, mice immunized with the T4 phage displaying dual antigen clusters achieved full protection (100% protection) against the H3N2 influenza virus, highlighting its robust protective efficacy.

CONCLUSIONS

In summary, our findings indicate that particles based on a T4 phage displaying antigen clusters exhibit ideal immunogenicity and protective effects, providing a promising strategy for the development of subunit vaccines against various viruses beyond influenza.

2018-2019 human seasonal H3N2 influenza A virus spillovers into swine with demonstrated virus transmission in pigs were not sustained in the pig population

Human seasonal H3 clade 3C3a influenza A viruses (IAV) were detected four times in U.S. pigs from commercial swine farms in Michigan, Illinois, and Virginia in 2019. To evaluate the relative risk of this spillover to the pig population, whole genome sequencing and phylogenetic characterization were conducted, and the results revealed that all eight viral gene segments were closely related to 2018-2019 H3N2 human seasonal IAV. Next, a series of in vitro viral kinetics, receptor binding, and antigenic characterization studies were performed using a representative A/swine/Virginia/A02478738/2018(H3N2) (SW/VA/19) isolate. Viral replication kinetic studies of SW/VA/19 demonstrated less efficient replication curves than all 10 swine H3N2 viruses tested but higher than three human H3N2 strains. Serial passaging experiments of SW/VA/19 in swine cells did not increase virus replication, but changes at HA amino acid positions 9 and 159 occurred. In swine transmission studies, wild-type SW/VA/19 was shed in nasal secretions and transmitted to all indirect contact pigs, whereas the human seasonal strain A/Switzerland/9715293/2013(H3N2) from the same 3C3a clade failed to transmit. SW/VA/19 induced minimal macroscopic and microscopic lung lesions. Collectively, these findings demonstrate that these human seasonal H3N2 3C3a-like viruses did not require reassortment with endemic swine IAV gene segments for virus shedding and transmission in pigs. Limited detections in the U.S. pig population in the subsequent period of time suggest a yet-unknown restriction factor likely limiting the spread of these viruses in the U.S. pig population.IMPORTANCEInterspecies human-to-swine IAV transmission occurs globally and contributes to increased IAV diversity in pig populations. We present data that a swine isolate from a 2018-2019 human-to-swine transmission event was shed for multiple days in challenged and contact pigs. By characterizing this introduction through bioinformatic, molecular, and animal experimental approaches, these findings better inform animal health practices and vaccine decision-making. Since wholly human seasonal H3N2 viruses in the United States were not previously identified as being transmissible in pigs (i.e., reverse zoonosis), these findings reveal that the interspecies barriers for transmission to pigs may not require significant changes to all human seasonal H3N2, although additional changes may be required for sustained transmission in swine populations.

Hemagglutinin glycosylation pattern-specific effects: implications for the fitness of H9.4.2.5-branched H9N2 avian influenza viruses

Since 2007, h9.4.2.5 has emerged as the most predominant branch of H9N2 avian influenza viruses (AIVs) that affects the majority of the global poultry population. The spread of this viral branch in vaccinated chicken flocks has not been considerably curbed despite numerous efforts. The evolutionary fitness of h9.4.2.5-branched AIVs must consequently be taken into consideration. The glycosylation modifications of hemagglutinin (HA) play a pivotal role in regulating the balance between receptor affinity and immune evasion for influenza viruses. Sequence alignment showed that five major HA glycosylation patterns have evolved over time in h9.4.2.5-branched AIVs. Here, we compared the adaptive phenotypes of five virus mutants with different HA glycosylation patterns. According to the results, the mutant with 6 N-linked glycans displayed the best acid and thermal stability and a better capacity for multiplication, although having a relatively lower receptor affinity than 7 glycans. The antigenic profile between the five mutants revealed a distinct antigenic distance, indicating that variations in glycosylation level have an impact on antigenic drift. These findings suggest that changes in the number of glycans on HA can not only modulate the receptor affinity and antigenicity of H9N2 AIVs, but also affect their stability and multiplication. These adaptive phenotypes may underlie the biological basis for the dominant strain switchover of h9.4.2.5-branched AIVs. Overall, our study provides a systematic insight into how changes in HA glycosylation patterns regulate the evolutionary fitness and epidemiological dominance drift of h9.4.2.5-branched H9N2 AIVs, which will be of great benefit for the glycosylation-dependent vaccine design.

Epistasis in the receptor binding domain of contemporary H3N2 viruses that reverted to bind sialylated diLacNAc repeats

Since the introduction of H3N2 influenza A viruses in the human population, these viruses have continuously evolved to escape human immunity, with mutations occurring in and around the receptor binding site. This process, called antigenic drift, recently resulted in viruses that recognize elongated glycans that are not abundantly displayed in the human respiratory tract. Such receptor specificities hampered our ability to pick and propagate vaccine strains. Using ELISA, glycan array, tissue staining, flow cytometry, and hemagglutinin assays, this study revealed that the most recent H3N2 viruses have expanded receptor specificity by regaining effective recognition to shorter glycans. In recent H3 strains, Y159 and T160 are responsible for restricted binding to elongated glycans; in contemporary strains, however, Y159N and T160I dominate with a consequent loss of strength in receptor binding. Yet, effective receptor interaction is rescued by a remote mutation in the 190-helix, Y195F. The results demonstrate epistasis of critical residues in three of the four structural elements composing the HA receptor-binding site (the 130-loop, 150-loop, and 190-helix), which synergistically contribute to shape receptor binding specificity. Interestingly, a positive correlation exists between binding to an asymmetrical N-glycan containing an α2,6 sialylated tri-LacNAc arm and binding to human and ferret respiratory tract tissues. Together, these results elucidate the epistatic nature of receptor binding specificity during influenza A virus H3N2 evolution.

Deep mutational scanning of H5 hemagglutinin to inform influenza virus surveillance

H5 influenza is considered a potential pandemic threat. Recently, H5 viruses belonging to clade 2.3.4.4b have caused large outbreaks in avian and multiple nonhuman mammalian species. Previous studies have identified molecular phenotypes of the viral hemagglutinin (HA) protein that contribute to pandemic potential in humans, including cell entry, receptor preference, HA stability, and reduced neutralization by polyclonal sera. However, prior experimental work has only measured how these phenotypes are affected by a handful of the >10,000 different possible amino-acid mutations to HA. Here, we use pseudovirus deep mutational scanning to measure how all mutations to a 2.3.4.4b H5 HA affect each phenotype. We identify mutations that allow HA to better bind α2-6-linked sialic acids and show that some viruses already carry mutations that stabilize HA. We also measure how all HA mutations affect neutralization by sera from mice and ferrets vaccinated against or infected with 2.3.4.4b H5 viruses. These antigenic maps enable rapid assessment of when new viral strains have acquired mutations that may create mismatches with candidate vaccine virus, and we show that a mutation present in some recent H5 HAs causes a large antigenic change. Overall, the systematic nature of deep mutational scanning combined with the safety of pseudoviruses enables comprehensive measurements of the phenotypic effects of mutations that can inform real-time interpretation of viral variation observed during surveillance of H5 influenza.

Avian influenza A H5N1 hemagglutinin protein models have distinct structural patterns re-occurring across the 1959-2023 strains

Influenza A H5N1 hemagglutinin (HA) plays a crucial role in viral pathogenesis and changes in the HA receptor binding domain (RBD) have been attributed to alterations in viral pathogenesis. Mutations often occur within the HA which in-turn results in HA structural changes that consequently contribute to protein evolution. However, the possible occurrence of mutations that results to reversion of the HA protein (going back to an ancestral protein conformation) which in-turn creates distinct HA structural patterns across the 1959-2023 H5N1 viral evolution has never been investigated. Here, we generated and verified the quality of the HA models, identified similar HA structural patterns, and elucidated the possible variations in HA RBD structural dynamics. Our results show that there are 7 distinct structural patterns occurring among the 1959-2023 H5N1 HA models which suggests that reversion of the HA protein putatively occurs during viral evolution. Similarly, we found that the HA RBD structural dynamics vary among the 7 distinct structural patterns possibly affecting viral pathogenesis.

Exploiting the Affimer platform against influenza A virus

Influenza A virus (IAV) is well known for its pandemic potential. While current surveillance and vaccination strategies are highly effective, therapeutic approaches are often short-lived due to the high mutation rates of IAV. Recently, monoclonal antibodies (mAbs) have emerged as a promising therapeutic approach, both against current strains and future IAV pandemics. In addition to mAbs, several antibody-like alternatives exist, which aim to improve upon mAbs. Among these, Affimers stand out for their short development time, high expression levels in Escherichia coli, and animal-free production. In this study, we utilized the Affimer platform to isolate and produce specific and potent inhibitors of IAV. Using a monomeric version of the IAV trimeric hemagglutinin (HA) fusion protein, we isolated 12 Affimers that inhibit IAV infection in vitro. Two of these Affimers were characterized in detail and exhibited nanomolar-binding affinities to the target H3 HA protein, specifically binding to the HA1 head domain. Cryo-electron microscopy (cryo-EM), employing a novel spray approach to prepare cryo-grids, allowed us to image HA-Affimer complexes. Combined with functional assays, we determined that these Affimers inhibit IAV by blocking the interaction of HA with the host-cell receptor, sialic acid. Furthermore, these Affimers inhibited IAV strains closely related to the one used for their isolation. Overall, our results support the use of Affimers as a viable alternative to existing targeted therapies for IAV and highlight their potential as diagnostic reagents.

IMPORTANCE

Influenza A virus is one of the few viruses that can cause devastating pandemics. Due to the high mutation rates of this virus, annual vaccination is required, and antivirals are short-lived. Monoclonal antibodies present a promising approach to tackle influenza virus infections but are associated with some limitations. To improve on this strategy, we explored the Affimer platform, which are antibody-like proteins made in bacteria. By performing phage-display against a monomeric version of influenza virus fusion protein, an established viral target, we were able to isolate Affimers that inhibit influenza virus infection in vitro. We characterized the mechanism of inhibition of the Affimers by using assays targeting different stages of the viral replication cycle. We additionally characterized HA-Affimer complex structure, using a novel approach to prepare samples for cryo-electron microscopy. Overall, these results show that Affimers are a promising tool against influenza virus infection.

Deep mutational scanning of H5 hemagglutinin to inform influenza virus surveillance

H5 influenza is a potential pandemic threat. Previous studies have identified molecular phenotypes of the viral hemagglutinin (HA) protein that contribute to pandemic risk, including cell entry, receptor preference, HA stability, and reduced neutralization by polyclonal sera. Here we use pseudovirus deep mutational scanning to measure how all mutations to a clade 2.3.4.4b H5 HA affect each phenotype. We identify mutations that allow HA to better bind a2-6-linked sialic acids, and show that some viruses already carry mutations that stabilize HA. We also identify recent viral strains with reduced neutralization to sera elicited by candidate vaccine virus. Overall, the systematic nature of deep mutational scanning combined with the safety of pseudoviruses enables comprehensive characterization of mutations to inform surveillance of H5 influenza.

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.

Probing altered receptor specificities of antigenically drifting human H3N2 viruses by chemoenzymatic synthesis, NMR, and modeling

Prototypic receptors for human influenza viruses are N-glycans carrying α2,6-linked sialosides. Due to immune pressure, A/H3N2 influenza viruses have emerged with altered receptor specificities that bind α2,6-linked sialosides presented on extended N-acetyl-lactosamine (LacNAc) chains. Here, binding modes of such drifted hemagglutinin's (HAs) are examined by chemoenzymatic synthesis of N-glycans having 13C-labeled monosaccharides at strategic positions. The labeled glycans are employed in 2D STD-1H by 13C-HSQC NMR experiments to pinpoint which monosaccharides of the extended LacNAc chain engage with evolutionarily distinct HAs. The NMR data in combination with computation and mutagenesis demonstrate that mutations distal to the receptor binding domain of recent HAs create an extended binding site that accommodates with the extended LacNAc chain. A fluorine containing sialoside is used as NMR probe to derive relative binding affinities and confirms the contribution of the extended LacNAc chain for binding.

Asymmetrical Biantennary Glycans Prepared by a Stop-and-Go Strategy Reveal Receptor Binding Evolution of Human Influenza A Viruses

Glycan binding properties of respiratory viruses have been difficult to probe due to a lack of biologically relevant glycans for binding studies. Here, a stop-and-go chemoenzymatic methodology is presented that gave access to a panel of 32 asymmetrical biantennary N-glycans having various numbers of N-acetyl lactosamine (LacNAc) repeating units capped by α2,3- or α2,6-sialosides resembling structures found in airway tissues. It exploits that the branching enzymes MGAT1 and MGAT2 can utilize unnatural UDP-2-deoxy-2-trifluoro-N-acetamido-glucose (UDP-GlcNTFA) as donor. The TFA moiety of the resulting glycans can be hydrolyzed to give GlcNH2 at one of the antennae, which temporarily blocks extension by glycosyl transferases. The N-glycans were printed as a microarray that was probed for receptor binding specificities of the evolutionary distinct human A(H3N2) and A(H1N1)pdm09 viruses. It was found that not only the sialoside type but also the length of the LacNAc chain and presentation at the α1,3-antenna of N-glycans are critical for binding. Early A(H3N2) viruses bound to 2,6-sialosides at a single LacNAc moiety at the α1,3-antenna whereas later viruses required the sialoside to be presented at a tri-LacNAc moiety. Surprisingly, most of the A(H3N2) viruses that appeared after 2021 regained binding capacity to sialosides presented at a di-LacNAc moiety. As a result, these viruses again agglutinate erythrocytes, commonly employed for antigenic characterization of influenza viruses. Human A(H1N1)pdm09 viruses have similar receptor binding properties as recent A(H3N2) viruses. The data indicate that an asymmetric N-glycan having 2,6-sialoside at a di-LacNAc moiety is a commonly employed receptor by human influenza A viruses.

Co-evolution of immunity and seasonal influenza viruses

Seasonal influenza viruses cause recurring global epidemics by continually evolving to escape host immunity. The viral constraints and host immune responses that limit and drive the evolution of these viruses are increasingly well understood. However, it remains unclear how most of these advances improve the capacity to reduce the impact of seasonal influenza viruses on human health. In this Review, we synthesize recent progress made in understanding the interplay between the evolution of immunity induced by previous infections or vaccination and the evolution of seasonal influenza viruses driven by the heterogeneous accumulation of antibody-mediated immunity in humans. We discuss the functional constraints that limit the evolution of the viruses, the within-host evolutionary processes that drive the emergence of new virus variants, as well as current and prospective options for influenza virus control, including the viral and immunological barriers that must be overcome to improve the effectiveness of vaccines and antiviral drugs.

Asymmetrical Bi-antennary Glycans Prepared by a Stop-and-Go Strategy Reveal Receptor Binding Evolution of Human Influenza A Viruses

Glycan binding properties of respiratory viruses have been difficult to probe due to a lack of biological relevant glycans for binding studies. Here, a stop-and-go chemoenzymatic methodology is presented that gave access to a panel of 32 asymmetrical bi-antennary N-glycans having various numbers of N-acetyl lactosamine (LacNAc) repeating units capped by α2,3- or α2,6-sialosides resembling structures found in airway tissues. It exploits that the branching enzymes MGAT1 and MGAT2 can utilize unnatural UDP-2-deoxy-2-trifluoro-N-acetamido-glucose (UDP-GlcNTFA) as donor. The TFA moiety of the resulting glycans can be hydrolyzed to give GlcNH2 at one of the antennae that temporarily blocks extension by glycosyl transferases. The N-glycans were printed as a microarray that was probed for receptor binding specificities of evolutionary distinct human A(H3N2) and A(H1N1)pdm09 viruses. It was found that not only the sialoside type but also the length of the LacNAc chain and presentation at the α1,3-antenna of N-glycans is critical for binding. Early A(H3N2) viruses bound to 2,6-sialosides at a single LacNAc moiety at the α1,3-antenna whereas later viruses required the sialoside to be presented at a tri-LacNAc moiety. Surprisingly, most of the A(H3N2) viruses that appeared after 2021 regained binding capacity to sialosides presented at a di-LacNAc moiety. As a result, these viruses agglutinate erythrocytes again, commonly employed for antigenic characterization of influenza viruses. Human A(H1N1)pdm09 viruses have similar receptor binding properties as recent A(H3N2) viruses. The data indicates that an asymmetric N-glycan having 2,6-sialoside at a di-LacNAc moiety is a commonly employed receptor by human influenza A viruses.

Paired associated SARS-CoV-2 spike variable positions: a network analysis approach to emerging variants

Amino acids in variable positions of proteins may be correlated, with potential structural and functional implications. Here, we apply exact tests of independence in R × C contingency tables to examine noise-free associations between variable positions of the SARS-CoV-2 spike protein, using as a paradigm sequences from Greece deposited in GISAID (N = 6,683/1,078 full length) for the period 29 February 2020 to 26 April 2021 that essentially covers the first three pandemic waves. We examine the fate and complexity of these associations by network analysis, using associated positions (exact P ≤ 0.001 and Average Product Correction ≥ 2) as links and the corresponding positions as nodes. We found a temporal linear increase of positional differences and a gradual expansion of the number of position associations over time, represented by a temporally evolving intricate web, resulting in a non-random complex network of 69 nodes and 252 links. Overconnected nodes corresponded to the most adapted variant positions in the population, suggesting a direct relation between network degree and position functional importance. Modular analysis revealed 25 k-cliques comprising 3 to 11 nodes. At different k-clique resolutions, one to four communities were formed, capturing epistatic associations of circulating variants (Alpha, Beta, B.1.1.318), but also Delta, which dominated the evolutionary landscape later in the pandemic. Cliques of aminoacidic positional associations tended to occur in single sequences, enabling the recognition of epistatic positions in real-world virus populations. Our findings provide a novel way of understanding epistatic relationships in viral proteins with potential applications in the design of virus control procedures. IMPORTANCE Paired positional associations of adapted amino acids in virus proteins may provide new insights for understanding virus evolution and variant formation. We investigated potential intramolecular relationships between variable SARS-CoV-2 spike positions by exact tests of independence in R × C contingency tables, having applied Average Product Correction (APC) to eliminate background noise. Associated positions (exact P ≤ 0.001 and APC ≥ 2) formed a non-random, epistatic network of 25 cliques and 1-4 communities at different clique resolutions, revealing evolutionary ties between variable positions of circulating variants and a predictive potential of previously unknown network positions. Cliques of different sizes represented theoretical combinations of changing residues in sequence space, allowing the identification of significant aminoacidic combinations in single sequences of real-world populations. Our analytic approach that links network structural aspects to mutational aminoacidic combinations in the spike sequence population offers a novel way to understand virus epidemiology and evolution.

Re-Emergence of Circulation of Seasonal Influenza during COVID-19 Pandemic in Russia and Receptor Specificity of New and Dominant Clade 3C.2a1b.2a.2 A(H3N2) Viruses in 2021-2022

The circulation of seasonal influenza in 2020-2021 around the world was drastically reduced after the start of the COVID-19 pandemic and the implementation of mitigation strategies. The influenza virus circulation reemerged in 2021-2022 with the global spread of the new genetic clade 3C.2a1b.2a.2 of A(H3N2) viruses. The purpose of this study was to characterize influenza viruses in the 2021-2022 season in Russia and to analyze the receptor specificity properties of the 3C.2a1b.2a.2 A(H3N2) viruses. Clinical influenza samples were collected at the local Sanitary-and-Epidemiological Centers of Rospotrebnadzor. Whole genome sequencing was performed using NGS. The receptor specificity of hemagglutinin was evaluated using molecular modeling and bio-layer interferometry. Clinical samples from 854 cases of influenza A and B were studied; A(H3N2) viruses were in the majority of the samples. All genetically studied A(H3N2) viruses belonged to the new genetic clade 3C.2a1b.2a.2. Molecular modeling analysis suggested a higher affinity of hemagglutinin of 3C.2a1b.2a.2. A(H3N2) viruses to the α2,6 human receptor. In vitro analysis using a trisaccharide 6'-Sialyl-N-acetyllactosamine receptor analog did not resolve the differences in the receptor specificity of 3C.2a1b.2a.2 clade viruses from viruses belonging to the 3C.2a1b.2a.1 clade. Further investigation of the A(H3N2) viruses is required for the evaluation of their possible adaptive advantages. Constant monitoring and characterization of influenza are critical for epidemiological analysis.

Prevalence and mechanisms of evolutionary contingency in human influenza H3N2 neuraminidase

Neuraminidase (NA) of human influenza H3N2 virus has evolved rapidly and been accumulating mutations for more than half-century. However, biophysical constraints that govern the evolutionary trajectories of NA remain largely elusive. Here, we show that among 70 natural mutations that are present in the NA of a recent human H3N2 strain, >10% are deleterious for an ancestral strain. By mapping the permissive mutations using combinatorial mutagenesis and next-generation sequencing, an extensive epistatic network is revealed. Biophysical and structural analyses further demonstrate that certain epistatic interactions can be explained by non-additive stability effect, which in turn modulates membrane trafficking and enzymatic activity of NA. Additionally, our results suggest that other biophysical mechanisms also contribute to epistasis in NA evolution. Overall, these findings not only provide mechanistic insights into the evolution of human influenza NA and elucidate its sequence-structure-function relationship, but also have important implications for the development of next-generation influenza vaccines.

The evolutionary potential of influenza A virus hemagglutinin is highly constrained by epistatic interactions with neuraminidase

Antigenic evolution of the influenza A virus (IAV) hemagglutinin (HA) gene limits efforts to effectively control the spread of the virus in the population. Efforts to understand the mechanisms governing HA antigenic evolution typically examine the HA gene in isolation. This can ignore the importance of balancing HA receptor binding activities with the receptor-destroying activities of the viral neuraminidase (NA) to maintain viral fitness. We hypothesize that the need to maintain functional balance with NA significantly constrains the evolutionary potential of the HA. We use deep mutational scanning and show that variation in NA activity significantly reshapes the HA fitness landscape by modulating the overall mutational robustness of HA. Consistent with this, we observe that different NA backgrounds support the emergence of distinct repertoires of HA escape variants under neutralizing antibody pressure. Our results reveal a critical role for intersegment epistasis in influencing the evolutionary potential of the HA gene.

Egg-adaptive mutations of human influenza H3N2 virus are contingent on natural evolution

Egg-adaptive mutations in influenza hemagglutinin (HA) often emerge during the production of egg-based seasonal influenza vaccines, which contribute to the largest share in the global influenza vaccine market. While some egg-adaptive mutations have minimal impact on the HA antigenicity (e.g. G186V), others can alter it (e.g. L194P). Here, we show that the preference of egg-adaptive mutation in human H3N2 HA is strain-dependent. In particular, Thr160 and Asn190, which are found in many recent H3N2 strains, restrict the emergence of L194P but not G186V. Our results further suggest that natural amino acid variants at other HA residues also play a role in determining the preference of egg-adaptive mutation. Consistently, recent human H3N2 strains from different clades acquire different mutations during egg passaging. Overall, these results demonstrate that natural mutations in human H3N2 HA can influence the preference of egg-adaptation mutation, which has important implications in seed strain selection for egg-based influenza vaccine.

Predicting Egg Passage Adaptations to Design Better Vaccines for the H3N2 Influenza Virus

Seasonal H3N2 influenza evolves rapidly, leading to an extremely poor vaccine efficacy. Substitutions employed during vaccine production using embryonated eggs (i.e., egg passage adaptation) contribute to the poor vaccine efficacy (VE), but the evolutionary mechanism remains elusive. Using an unprecedented number of hemagglutinin sequences (n = 89,853), we found that the fitness landscape of passage adaptation is dominated by pervasive epistasis between two leading residues (186 and 194) and multiple other positions. Convergent evolutionary paths driven by strong epistasis explain most of the variation in VE, which has resulted in extremely poor vaccines for the past decade. Leveraging the unique fitness landscape, we developed a novel machine learning model that can predict egg passage substitutions for any candidate vaccine strain before the passage experiment, providing a unique opportunity for the selection of optimal vaccine viruses. Our study presents one of the most comprehensive characterizations of the fitness landscape of a virus and demonstrates that evolutionary trajectories can be harnessed for improved influenza vaccines.

Membrane-Tethered Mucin 1 Is Stimulated by Interferon and Virus Infection in Multiple Cell Types and Inhibits Influenza A Virus Infection in Human Airway Epithelium

Influenza A virus (IAV) causes significant morbidity and mortality in the human population. Tethered mucin 1 (MUC1) is highly expressed in airway epithelium, the primary site of IAV replication, and also by other cell types that influence IAV infection, including macrophages. MUC1 has the potential to influence infection dynamics through physical interactions and/or signaling activity, yet MUC1 modulation and its impact during viral pathogenesis remain unclear. Thus, we investigated MUC1-IAV interactions in an in vitro model of human airway epithelium (HAE). Our data indicate that a recombinant IAV hemagglutinin (H3) and H3N2 virus can bind endogenous HAE MUC1. Notably, infection of HAE with H1N1 or H3N2 IAV strains does not trigger MUC1 shedding but instead stimulates an increase in cell-associated MUC1 protein. We observed a similar increase after type I or III interferon (IFN) stimulation; however, inhibition of IFN signaling during H1N1 infection only partially abrogated this increase, indicating that multiple soluble factors contribute to MUC1 upregulation during the antiviral response. In addition to HAE, primary human monocyte-derived macrophages also upregulated MUC1 protein in response to IFN treatment and conditioned media from IAV-infected HAE. Then, to determine the impact of MUC1 on IAV pathogenesis, we developed HAE genetically depleted of MUC1 and found that MUC1 knockout cultures exhibited enhanced viral growth compared to control cultures for several IAV strains. Together, our data support a model whereby MUC1 inhibits productive uptake of IAV in HAE. Infection then stimulates MUC1 expression on multiple cell types through IFN-dependent and -independent mechanisms that further impact infection dynamics.

The 2009 Pandemic H1N1 Hemagglutinin Stalk Remained Antigenically Stable after Circulating in Humans for a Decade

An H1N1 influenza virus caused a pandemic in 2009, and descendants of this virus continue to circulate seasonally in humans. Upon infection with the 2009 H1N1 pandemic strain (pH1N1), many humans produced antibodies against epitopes in the hemagglutinin (HA) stalk. HA stalk-focused antibody responses were common among pH1N1-infected individuals because HA stalk epitopes were conserved between the pH1N1 strain and previously circulating H1N1 strains. Here, we completed a series of experiments to determine if the pH1N1 HA stalk has acquired substitutions since 2009 that prevent the binding of human antibodies. We identified several amino acid substitutions that accrued in the pH1N1 HA stalk from 2009 to 2019. We completed enzyme-linked immunosorbent assays, absorption-based binding assays, and surface plasmon resonance experiments to determine if these substitutions affect antibody binding. Using sera collected from 230 humans (aged 21 to 80 years), we found that pH1N1 HA stalk substitutions that have emerged since 2009 do not affect antibody binding. Our data suggest that the HA stalk domain of pH1N1 viruses remained antigenically stable after circulating in humans for a decade. IMPORTANCE In 2009, a new pandemic H1N1 (pH1N1) virus began circulating in humans. Many individuals mounted hemagglutinin (HA) stalk-focused antibody responses upon infection with the 2009 pH1N1 strain, since the HA stalk of this virus was relatively conserved with other seasonal H1N1 strains. Here, we completed a series of studies to determine if the 2009 pH1N1 strain has undergone antigenic drift in the HA stalk domain over the past decade. We found that serum antibodies from 230 humans could not antigenically distinguish the 2009 and 2019 HA stalk. These data suggest that the HA stalk of pH1N1 has remained antigenically stable, despite the presence of high levels of HA stalk antibodies within the human population.

Conserved topology of virus glycoepitopes presents novel targets for repurposing HIV antibody 2G12

Complex glycans decorate viral surface proteins and play a critical role in virus-host interactions. Viral surface glycans shield vulnerable protein epitopes from host immunity yet can also present distinct "glycoepitopes" that can be targeted by host antibodies such as the potent anti-HIV antibody 2G12 that binds high-mannose glycans on gp120. Two recent publications demonstrate 2G12 binding to high mannose glycans on SARS-CoV-2 and select Influenza A (Flu) H3N2 viruses. Previously, our lab observed 2G12 binding and functional inhibition of a range of Flu viruses that include H3N2 and H1N1 lineages. In this manuscript, we present these data alongside structural analyses to offer an expanded picture of 2G12-Flu interactions. Further, based on the remarkable breadth of 2G12 N-glycan recognition and the structural factors promoting glycoprotein oligomannosylation, we hypothesize that 2G12 glycoepitopes can be defined from protein structure alone according to N-glycan site topology. We develop a model describing 2G12 glycoepitopes based on N-glycan site topology, and apply the model to identify viruses within the Protein Data Bank presenting putative 2G12 glycoepitopes for 2G12 repurposing toward analytical, diagnostic, and therapeutic applications.

Inferring the intrinsic mutational fitness landscape of influenzalike evolving antigens from temporally ordered sequence data

There still are no effective long-term protective vaccines against viruses that continuously evolve under immune pressure such as seasonal influenza, which has caused, and can cause, devastating epidemics in the human population. To find such a broadly protective immunization strategy, it is useful to know how easily the virus can escape via mutation from specific antibody responses. This information is encoded in the fitness landscape of the viral proteins (i.e., knowledge of the viral fitness as a function of sequence). Here we present a computational method to infer the intrinsic mutational fitness landscape of influenzalike evolving antigens from yearly sequence data. We test inference performance with computer-generated sequence data that are based on stochastic simulations mimicking basic features of immune-driven viral evolution. Although the numerically simulated model does create a phylogeny based on the allowed mutations, the inference scheme does not use this information. This provides a contrast to other methods that rely on reconstruction of phylogenetic trees. Our method just needs a sufficient number of samples over multiple years. With our method, we are able to infer single as well as pairwise mutational fitness effects from the simulated sequence time series for short antigenic proteins. Our fitness inference approach may have potential future use for the design of immunization protocols by identifying intrinsically vulnerable immune target combinations on antigens that evolve under immune-driven selection. In the future, this approach may be applied to influenza and other novel viruses such as SARS-CoV-2, which evolves and, like influenza, might continue to escape the natural and vaccine-mediated immune pressures.

Key Substitutions in the Spike Protein of SARS-CoV-2 Variants Can Predict Resistance to Monoclonal Antibodies, but Other Substitutions Can Modify the Effects

Mutations in the spike protein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants can compromise the effectiveness of therapeutic antibodies. Most clinical-stage therapeutic antibodies target the spike receptor binding domain (RBD), but variants often have multiple mutations in several spike regions. To help predict antibody potency against emerging variants, we evaluated 25 clinical-stage therapeutic antibodies for neutralization activity against 60 pseudoviruses bearing spikes with single or multiple substitutions in several spike domains, including the full set of substitutions in B.1.1.7 (alpha), B.1.351 (beta), P.1 (gamma), B.1.429 (epsilon), B.1.526 (iota), A.23.1, and R.1 variants. We found that 14 of 15 single antibodies were vulnerable to at least one RBD substitution, but most combination and polyclonal therapeutic antibodies remained potent. Key substitutions in variants with multiple spike substitutions predicted resistance, but the degree of resistance could be modified in unpredictable ways by other spike substitutions that may reside outside the RBD. These findings highlight the importance of assessing antibody potency in the context of all substitutions in a variant and show that epistatic interactions in spike can modify virus susceptibility to therapeutic antibodies. IMPORTANCE Therapeutic antibodies are effective in preventing severe disease from SARS-CoV-2 infection (COVID-19), but their effectiveness may be reduced by virus variants with mutations affecting the spike protein. To help predict resistance to therapeutic antibodies in emerging variants, we profiled resistance patterns of 25 antibody products in late stages of clinical development against a large panel of variants that include single and multiple substitutions found in the spike protein. We found that the presence of a key substitution in variants with multiple spike substitutions can predict resistance against a variant but that other substitutions can affect the degree of resistance in unpredictable ways. These findings highlight complex interactions among substitutions in the spike protein affecting virus neutralization and, potentially, virus entry into cells.

Antigenic evolution of human influenza H3N2 neuraminidase is constrained by charge balancing

As one of the main influenza antigens, neuraminidase (NA) in H3N2 virus has evolved extensively for more than 50 years due to continuous immune pressure. While NA has recently emerged as an effective vaccine target, biophysical constraints on the antigenic evolution of NA remain largely elusive. Here, we apply combinatorial mutagenesis and next-generation sequencing to characterize the local fitness landscape in an antigenic region of NA in six different human H3N2 strains that were isolated around 10 years apart. The local fitness landscape correlates well among strains and the pairwise epistasis is highly conserved. Our analysis further demonstrates that local net charge governs the pairwise epistasis in this antigenic region. In addition, we show that residue coevolution in this antigenic region is correlated with the pairwise epistasis between charge states. Overall, this study demonstrates the importance of quantifying epistasis and the underlying biophysical constraint for building a model of influenza evolution.

Glycan remodeled erythrocytes facilitate antigenic characterization of recent A/H3N2 influenza viruses

During circulation in humans and natural selection to escape antibody recognition for decades, A/H3N2 influenza viruses emerged with altered receptor specificities. These viruses lost the ability to agglutinate erythrocytes critical for antigenic characterization and give low yields and acquire adaptive mutations when cultured in eggs and cells, contributing to recent vaccine challenges. Examination of receptor specificities of A/H3N2 viruses reveals that recent viruses compensated for decreased binding of the prototypic human receptor by recognizing α2,6-sialosides on extended LacNAc moieties. Erythrocyte glycomics shows an absence of extended glycans providing a rationale for lack of agglutination by recent A/H3N2 viruses. A glycan remodeling approach installing functional receptors on erythrocytes, allows antigenic characterization of recent A/H3N2 viruses confirming the cocirculation of antigenically different viruses in humans. Computational analysis of HAs in complex with sialosides having extended LacNAc moieties reveals that mutations distal to the RBD reoriented the Y159 side chain resulting in an extended receptor binding site.

An evolution-based high-fidelity method of epistasis measurement: Theory and application to influenza

Linkage effects in a multi-locus population strongly influence its evolution. The models based on the traveling wave approach enable us to predict the average speed of evolution and the statistics of phylogeny. However, predicting statistically the evolution of specific sites and pairs of sites in the multi-locus context remains a mathematical challenge. In particular, the effects of epistasis, the interaction of gene regions contributing to phenotype, is difficult to predict theoretically and detect experimentally in sequence data. A large number of false-positive interactions arises from stochastic linkage effects and indirect interactions, which mask true epistatic interactions. Here we develop a proof-of-principle method to filter out false-positive interactions. We start by demonstrating that the averaging of haplotype frequencies over multiple independent populations is necessary but not sufficient for epistatic detection, because it still leaves high numbers of false-positive interactions. To compensate for the residual stochastic noise, we develop a three-way haplotype method isolating true interactions. The fidelity of the method is confirmed analytically and on simulated genetic sequences evolved with a known epistatic network. The method is then applied to a large sequence database of neurominidase protein of influenza A H1N1 obtained from various geographic locations to infer the epistatic network responsible for the difference between the pre-pandemic virus and the pandemic strain of 2009. These results present a simple and reliable technique to measure epistatic interactions of any sign from sequence data.

Interactions between genomic sites create a fitness landscape. The knowledge of topology and strength of interactions is vital for predicting the escape of viruses from drugs and immune response and their passing through fitness valleys. Many efforts have been invested into measuring these interactions from DNA sequence sets. Unfortunately, reproducibility of the results remains low due partly to a very small fraction of interaction pairs and partly to stochastic linkage noise masking true interactions. Here we propose a method to separate stochastic linkage and indirect interactions from epistatic interactions and apply it to influenza virus sequence data.

Glycans in Virus-Host Interactions: A Structural Perspective

Many interactions between microbes and their hosts are driven or influenced by glycans, whose heterogeneous and difficult to characterize structures have led to an underappreciation of their role in these interactions compared to protein-based interactions. Glycans decorate microbe glycoproteins to enhance attachment and fusion to host cells, provide stability, and evade the host immune system. Yet, the host immune system may also target these glycans as glycoepitopes. In this review, we provide a structural perspective on the role of glycans in host-microbe interactions, focusing primarily on viral glycoproteins and their interactions with host adaptive immunity. In particular, we discuss a class of topological glycoepitopes and their interactions with topological mAbs, using the anti-HIV mAb 2G12 as the archetypical example. We further offer our view that structure-based glycan targeting strategies are ready for application to viruses beyond HIV, and present our perspective on future development in this area.

Contingency and chance erase necessity in the experimental evolution of ancestral proteins

The roles of chance, contingency, and necessity in evolution are unresolved because they have never been assessed in a single system or on timescales relevant to historical evolution. We combined ancestral protein reconstruction and a new continuous evolution technology to mutate and select proteins in the B-cell lymphoma-2 (BCL-2) family to acquire protein–protein interaction specificities that occurred during animal evolution. By replicating evolutionary trajectories from multiple ancestral proteins, we found that contingency generated over long historical timescales steadily erased necessity and overwhelmed chance as the primary cause of acquired sequence variation; trajectories launched from phylogenetically distant proteins yielded virtually no common mutations, even under strong and identical selection pressures. Chance arose because many sets of mutations could alter specificity at any timepoint; contingency arose because historical substitutions changed these sets. Our results suggest that patterns of variation in BCL-2 sequences – and likely other proteins, too – are idiosyncratic products of a particular and unpredictable course of historical events.

One of the most fundamental and unresolved questions in evolutionary biology is whether the outcomes of evolution are predictable. Is the diversity of life we see today the expected result of organisms adapting to their environment throughout history (also known as natural selection) or the product of random chance? Or did chance events early in history shape the paths that evolution could take next, determining the biological forms that emerged under natural selection much later?

These questions are hard to study because evolution happened only once, long ago. To overcome this barrier, Xie, Pu, Metzger et al. developed an experimental approach that can evolve reconstructed ancestral proteins that existed deep in the past. Using this method, it is possible to replay evolution multiple times, from various historical starting points, under conditions similar to those that existed long ago. The end products of the evolutionary trajectories can then be compared to determine how predictable evolution actually is.

Xie, Pu, Metzger et al. studied proteins belonging to the BCL-2 family, which originated some 800 million years ago. These proteins have diversified greatly over time in both their genetic sequences and their ability to bind to specific partner proteins called co-regulators. Xie, Pu, Metzger et al. synthesized BCL-2 proteins that existed at various times in the past. Each ancestral protein was then allowed to evolve repeatedly under natural selection to acquire the same co-regulator binding functions that evolved during history.

At the end of each evolutionary trajectory, the genetic sequence of the resulting BCL-2 proteins was recorded. This revealed that the outcomes of evolution were almost completely unpredictable: trajectories initiated from the same ancestral protein produced proteins with very different sequences, and proteins launched from different ancestral starting points were even more dissimilar.

Further experiments identified the mutations in each trajectory that caused changes in coregulator binding. When these mutations were introduced into other ancestral proteins, they did not yield the same change in function. This suggests that early chance events influenced each protein’s evolution in an unpredictable way by opening and closing the paths available to it in the future.

This research expands our understanding of evolution on a molecular level whilst providing a new experimental approach for studying evolutionary drivers in more detail. The results suggest that BCL-2 proteins, in all their various forms, are unique products of a particular, unpredictable course of history set in motion by ancient chance events.

Applications of Deep Mutational Scanning in Virology

Several recently developed high-throughput techniques have changed the field of molecular virology. For example, proteomics studies reveal complete interactomes of a viral protein, genome-wide CRISPR knockout and activation screens probe the importance of every single human gene in aiding or fighting a virus, and ChIP-seq experiments reveal genome-wide epigenetic changes in response to infection. Deep mutational scanning is a relatively novel form of protein science which allows the in-depth functional analysis of every nucleotide within a viral gene or genome, revealing regions of importance, flexibility, and mutational potential. In this review, we discuss the application of this technique to RNA viruses including members of the Flaviviridae family, Influenza A Virus and Severe Acute Respiratory Syndrome Coronavirus 2. We also briefly discuss the reverse genetics systems which allow for analysis of viral replication cycles, next-generation sequencing technologies and the bioinformatics tools that facilitate this research.

A cross-neutralizing antibody between HIV-1 and influenza virus

Incessant antigenic evolution enables the persistence and spread of influenza virus in the human population. As the principal target of the immune response, the hemagglutinin (HA) surface antigen on influenza viruses continuously acquires and replaces N-linked glycosylation sites to shield immunogenic protein epitopes using host-derived glycans. Anti-glycan antibodies, such as 2G12, target the HIV-1 envelope protein (Env), which is even more extensively glycosylated and contains under-processed oligomannose-type clusters on its dense glycan shield. Here, we illustrate that 2G12 can also neutralize human seasonal influenza A H3N2 viruses that have evolved to present similar oligomannose-type clusters on their HAs from around 20 years after the 1968 pandemic. Using structural biology and mass spectrometric approaches, we find that two N-glycosylation sites close to the receptor binding site (RBS) on influenza hemagglutinin represent the oligomannose cluster recognized by 2G12. One of these glycan sites is highly conserved in all human H3N2 strains and the other emerged during virus evolution. These two N-glycosylation sites have also become crucial for fitness of recent H3N2 strains. These findings shed light on the evolution of the glycan shield on influenza virus and suggest 2G12-like antibodies can potentially act as broad neutralizers to target human enveloped viruses.

Site-Specific Evolutionary Rate Shifts in HIV-1 and SIV

Site-specific evolutionary rate shifts are defined as protein sites, where the rate of substitution has changed dramatically across the phylogeny. With respect to a given clade, sites may either undergo a rate acceleration or a rate deceleration, reflecting a site that was conserved and became variable, or vice-versa, respectively. Sites displaying such a dramatic evolutionary change may point to a loss or gain of function at the protein site, reflecting adaptation, or they may indicate epistatic interactions among sites. Here, we analyzed full genomes of HIV and SIV-1 and identified 271 rate-shifting sites along the HIV-1/SIV phylogeny. The majority of rate shifts occurred at long branches, often corresponding to cross-species transmission branches. We noted that in most proteins, the number of rate accelerations and decelerations was equal, and we suggest that this reflects epistatic interactions among sites. However, several accessory proteins were enriched for either accelerations or decelerations, and we suggest that this may be a signature of adaptation to new hosts. Interestingly, the non-pandemic HIV-1 group O clade exhibited a substantially higher number of rate-shift events than the pandemic group M clade. We propose that this may be a reflection of the height of the species barrier between gorillas and humans versus chimpanzees and humans. Our results provide a genome-wide view of the constraints operating on proteins of HIV-1 and SIV.

Host Receptors of Influenza Viruses and Coronaviruses-Molecular Mechanisms of Recognition

Among the four genera of influenza viruses (IVs) and the four genera of coronaviruses (CoVs), zoonotic αIV and βCoV have occasionally caused airborne epidemic outbreaks in humans, who are immunologically naïve, and the outbreaks have resulted in high fatality rates as well as social and economic disruption and losses. The most devasting influenza A virus (IAV) in αIV, pandemic H1N1 in 1918, which caused at least 40 million deaths from about 500 million cases of infection, was the first recorded emergence of IAVs in humans. Usually, a novel human-adapted virus replaces the preexisting human-adapted virus. Interestingly, two IAV subtypes, A/H3N2/1968 and A/H1N1/2009 variants, and two lineages of influenza B viruses (IBV) in βIV, B/Yamagata and B/Victoria lineage-like viruses, remain seasonally detectable in humans. Both influenza C viruses (ICVs) in γIV and four human CoVs, HCoV-229E and HCoV-NL63 in αCoV and HCoV-OC43 and HCoV-HKU1 in βCoV, usually cause mild respiratory infections. Much attention has been given to CoVs since the global epidemic outbreaks of βSARS-CoV in 2002-2004 and βMERS-CoV from 2012 to present. βSARS-CoV-2, which is causing the ongoing COVID-19 pandemic that has resulted in 890,392 deaths from about 27 million cases of infection as of 8 September 2020, has provoked worldwide investigations of CoVs. With the aim of developing efficient strategies for controlling virus outbreaks and recurrences of seasonal virus variants, here we overview the structures, diversities, host ranges and host receptors of all IVs and CoVs and critically review current knowledge of receptor binding specificity of spike glycoproteins, which mediates infection, of IVs and of zoonotic, pandemic and seasonal CoVs.

Computationally Designed Cyclic Peptides Derived from an Antibody Loop Increase Breadth of Binding for Influenza Variants

The influenza hemagglutinin (HA) glycoprotein is the target of many broadly neutralizing antibodies. However, influenza viruses can rapidly escape antibody recognition by mutation of hypervariable regions of HA that overlap with the binding epitope. We hypothesized that by designing peptides to mimic antibody loops, we could enhance breadth of binding to HA antigenic variants by reducing contact with hypervariable residues on HA that mediate escape. We designed cyclic peptides that mimic the heavy-chain complementarity-determining region 3 (CDRH3) of anti-influenza broadly neutralizing antibody C05 and show that these peptides bound to HA molecules with <100 nM affinity, comparable with that of the full-length parental C05 IgG. In addition, these peptides exhibited increased breadth of recognition to influenza H4 and H7 subtypes by eliminating clashes between the hypervariable antigenic regions and the antibody CDRH1 loop. This approach can be used to generate antibody-derived peptides against a wide variety of targets.

Structural Biology of Influenza Hemagglutinin: An Amaranthine Adventure

Hemagglutinin (HA) glycoprotein is an important focus of influenza research due to its role in antigenic drift and shift, as well as its receptor binding and membrane fusion functions, which are indispensable for viral entry. Over the past four decades, X-ray crystallography has greatly facilitated our understanding of HA receptor binding, membrane fusion, and antigenicity. The recent advances in cryo-EM have further deepened our comprehension of HA biology. Since influenza HA constantly evolves in natural circulating strains, there are always new questions to be answered. The incessant accumulation of knowledge on the structural biology of HA over several decades has also facilitated the design and development of novel therapeutics and vaccines. This review describes the current status of the field of HA structural biology, how we got here, and what the next steps might be.

Influenza Hemagglutinin Structures and Antibody Recognition

Hemagglutinin (HA) is most abundant glycoprotein on the influenza virus surface. Influenza HA promotes viral entry by engaging the receptor and mediating virus-host membrane fusion. At the same time, HA is the major antigen of the influenza virus. HA antigenic shift can result in pandemics, whereas antigenic drift allows human circulating strains to escape herd immunity. Most antibody responses against HA are strain-specific. However, antibodies that have neutralizing activities against multiple strains or even subtypes have now been discovered and characterized. These broadly neutralizing antibodies (bnAbs) target conserved regions on HA, such as the receptor-binding site and the stem domain. Structural studies of such bnAbs have provided important insight into universal influenza vaccine and therapeutic design. This review discusses the HA functions as well as HA-antibody interactions from a structural perspective.

Immunodominance and Antigenic Variation of Influenza Virus Hemagglutinin: Implications for Design of Universal Vaccine Immunogens

Influenza viruses routinely acquire mutations in their hemagglutinin (HA) and neuraminidase (NA) glycoproteins that abrogate binding of pre-existing antibodies in a process known as antigenic drift. Most human antibodies against HA and NA are directed against epitopes that are hypervariable and not against epitopes that are conserved among different influenza virus strains. Universal influenza vaccines are currently being developed to elicit protective responses against functionally conserved sites on influenza proteins where viral escape mutations can result in large fitness costs [1]. Universal vaccine targets include the highly conserved HA stem domain [2-12], the less conserved HA receptor-binding site (RBS) [13-16], as well as conserved sites on NA [17-19]. One central challenge of universal vaccine efforts is to steer human antibody responses away from immunodominant, variable epitopes and towards subdominant, functionally conserved sites. Overcoming this challenge will require further understanding of the structural basis of broadly neutralizing HA and NA antibody binding epitopes and factors that influence immunodominance hierarchies of human antibody responses.

Major antigenic site B of human influenza H3N2 viruses has an evolving local fitness landscape

Antigenic drift of influenza virus hemagglutinin (HA) is enabled by facile evolvability. However, HA antigenic site B, which has become immunodominant in recent human H3N2 influenza viruses, is also evolutionarily constrained by its involvement in receptor binding. Here, we employ deep mutational scanning to probe the local fitness landscape of HA antigenic site B in six different human H3N2 strains spanning from 1968 to 2016. We observe that the fitness landscape of HA antigenic site B can be very different between strains. Sequence variants that exhibit high fitness in one strain can be deleterious in another, indicating that the evolutionary constraints of antigenic site B have changed over time. Structural analysis suggests that the local fitness landscape of antigenic site B can be reshaped by natural mutations via modulation of the receptor-binding mode. Overall, these findings elucidate how influenza virus continues to explore new antigenic space despite strong functional constraints.

Influenza A Virus Hemagglutinin-Neuraminidase-Receptor Balance: Preserving Virus Motility

Influenza A viruses (IAVs) occasionally cross the species barrier and adapt to novel host species. This requires readjustment of the functional balance of the sialic acid receptor-binding hemagglutinin (HA) and the receptor-destroying neuraminidase (NA) to the sialoglycan-receptor repertoire of the new host. Novel techniques have revealed mechanistic details of this HA-NA-receptor balance, emphasizing a previously underappreciated crucial role for NA in driving the motility of receptor-associated IAV particles. Motility enables virion penetration of the sialylated mucus layer as well as attachment to, and uptake into, underlying epithelial cells. As IAVs are essentially irreversibly bound in the absence of NA activity, the fine-tuning of the HA-NA-receptor balance rather than the binding avidity of IAV particles per se is an important factor in determining host species tropism.

Enzyme Evolution: An Epistatic Ratchet versus a Smooth Reversible Transition

Evolutionary trajectories are deemed largely irreversible. In a newly diverged protein, reversion of mutations that led to the functional switch typically results in loss of both the new and the ancestral functions. Nonetheless, evolutionary transitions where reversions are viable have also been described. The structural and mechanistic causes of reversion compatibility versus incompatibility therefore remain unclear. We examined two laboratory evolution trajectories of mammalian paraoxonase-1, a lactonase with promiscuous organophosphate hydrolase (OPH) activity. Both trajectories began with the same active-site mutant, His115Trp, which lost the native lactonase activity and acquired higher OPH activity. A neo-functionalization trajectory amplified the promiscuous OPH activity, whereas the re-functionalization trajectory restored the native activity, thus generating a new lactonase that lacks His115. The His115 revertants of these trajectories indicated opposite trends. Revertants of the neo-functionalization trajectory lost both the evolved OPH and the original lactonase activity. Revertants of the trajectory that restored the original lactonase function were, however, fully active. Crystal structures and molecular simulations show that in the newly diverged OPH, the reverted His115 and other catalytic residues are displaced, thus causing loss of both the original and the new activity. In contrast, in the re-functionalization trajectory, reversion compatibility of the original lactonase activity derives from mechanistic versatility whereby multiple residues can fulfill the same task. This versatility enables unique sequence-reversible compositions that are inaccessible when the active site was repurposed toward a new function.

Why Glycosylation Matters in Building a Better Flu Vaccine

Low vaccine efficacy against seasonal influenza A virus (IAV) stems from the ability of the virus to evade existing immunity while maintaining fitness. Although most potent neutralizing antibodies bind antigenic sites on the globular head domain of the IAV envelope glycoprotein hemagglutinin (HA), the error-prone IAV polymerase enables rapid evolution of key antigenic sites, resulting in immune escape. Significantly, the appearance of new N-glycosylation consensus sequences (sequons, NXT/NXS, rarely NXC) on the HA globular domain occurs among the more prevalent mutations as an IAV strain undergoes antigenic drift. The appearance of new glycosylation shields underlying amino acid residues from antibody contact, tunes receptor specificity, and balances receptor avidity with virion escape, all of which help maintain viral propagation through seasonal mutations. The World Health Organization selects seasonal vaccine strains based on information from surveillance, laboratory, and clinical observations. Although the genetic sequences are known, mature glycosylated structures of circulating strains are not defined. In this review, we summarize mass spectrometric methods for quantifying site-specific glycosylation in IAV strains and compare the evolution of IAV glycosylation to that of human immunodeficiency virus. We argue that the determination of site-specific glycosylation of IAV glycoproteins would enable development of vaccines that take advantage of glycosylation-dependent mechanisms whereby virus glycoproteins are processed by antigen presenting cells.

Influenza Virus-Specific Human Antibody Repertoire Studies

The diversity of Ag-specific adaptive receptors on the surface of B cells and in the population of secreted Abs is enormous, but increasingly, we are acquiring the technical capability to interrogate Ab repertoires in great detail. These Ab technologies have been especially pointed at understanding the complex issues of immunity to infection and disease caused by influenza virus, one of the most common and vexing medical problems in man. Influenza immunity is particularly interesting as a model system because the antigenic diversity of influenza strains and proteins is high and constantly evolving. Discovery of canonical features in the subset of the influenza repertoire response that is broadly reactive for diverse influenza strains has spurred the recent optimism for creating universal influenza vaccines. Using new technologies for sequencing Ab repertoires at great depth is helping us to understand the central features of influenza immunity.

Molecular Docking of Broad-Spectrum Antibodies on Hemagglutinins of Influenza A Virus

Influenza A has caused several deadly pandemics throughout human history. The virus is often resistant to developed treatments because of its genetic drift or shift property. Broad-spectrum antibodies show a promising potential to overcome the resistance of influenza viruses. In silico studies on broad-reactive antibodies and their interactions with hemagglutinins might shed light on the rational design of a universal vaccine. In this study, 11 broad-spectrum antibodies (or antigen-binding fragments) and 14 hemagglutinins of H3N2 and H5N1 strains were docked and analyzed to provide information about the construction of the scaffold for using universal antibodies against the influenza A virus. Antigen-binding fragments that have high number of appearances in the top 3 within each H3 and H5 subtypes were chosen for protein-protein interaction analysis. The results show that while the hydrogen bond is important for Ab/Fab binding to H3, the H5-Ab/Fab system may need cation-pi interaction for a strong interaction.

Mapping person-to-person variation in viral mutations that escape polyclonal serum targeting influenza hemagglutinin

A longstanding question is how influenza virus evolves to escape human immunity, which is polyclonal and can target many distinct epitopes. Here, we map how all amino-acid mutations to influenza's major surface protein affect viral neutralization by polyclonal human sera. The serum of some individuals is so focused that it selects single mutations that reduce viral neutralization by over an order of magnitude. However, different viral mutations escape the sera of different individuals. This individual-to-individual variation in viral escape mutations is not present among ferrets that have been infected just once with a defined viral strain. Our results show how different single mutations help influenza virus escape the immunity of different members of the human population, a phenomenon that could shape viral evolution and disease susceptibility.

Preventing an Antigenically Disruptive Mutation in Egg-Based H3N2 Seasonal Influenza Vaccines by Mutational Incompatibility

Egg-based seasonal influenza vaccines are the major preventive countermeasure against influenza virus. However, their effectiveness can be compromised when antigenic changes arise from egg-adaptive mutations on influenza hemagglutinin (HA). The L194P mutation is commonly observed in egg-based H3N2 vaccine seed strains and significantly alters HA antigenicity. An approach to prevent L194P would therefore be beneficial. We show that emergence of L194P during egg passaging can be impeded by preexistence of a G186V mutation, revealing strong incompatibility between these mutations. X-ray structures illustrate that individual G186V and L194P mutations have opposing effects on the HA receptor-binding site (RBS), and when both G186V and L194P are present, the RBS is severely disrupted. Importantly, wild-type HA antigenicity is maintained with G186V, but not L194P. Our results demonstrate that these epistatic interactions can be used to prevent the emergence of mutations that adversely alter antigenicity during egg adaptation.

•

Most H3N2 egg isolates carry hemagglutinin mutation G186V or L194P, but not both

•

Hemagglutinin double mutation G186V/L194P is highly deleterious to the virus

•

Hemagglutinin double mutation G186V/L194P disrupts the receptor-binding site

•

Wild-type hemagglutinin antigenicity is maintained in G186V, but not in L194P

Targeting Hemagglutinin: Approaches for Broad Protection against the Influenza A Virus

Influenza A viruses are dynamically epidemic and genetically diverse. Due to the antigenic drift and shift of the virus, seasonal vaccines are required to be reformulated annually to match with current circulating strains. However, the mismatch between vaccinal strains and circulating strains occurs frequently, resulting in the low efficacy of seasonal vaccines. Therefore, several “universal” vaccine candidates based on the structure and function of the hemagglutinin (HA) protein have been developed to meet the requirement of a broad protection against homo-/heterosubtypic challenges. Here, we review recent novel constructs and discuss several important findings regarding the broad protective efficacy of HA-based universal vaccines.

Epistatic interactions can moderate the antigenic effect of substitutions in haemagglutinin of influenza H3N2 virus

We previously showed that single amino acid substitutions at seven positions in haemagglutinin determined major antigenic change of influenza H3N2 virus. Here, the impact of two such substitutions was tested in 11 representative H3 haemagglutinins to investigate context-dependence effects. The antigenic effect of substitutions introduced at haemagglutinin position 145 was fully independent of the amino acid context of the representative haemagglutinins. Antigenic change caused by substitutions introduced at haemagglutinin position 155 was variable and context-dependent. Our results suggest that epistatic interactions with contextual amino acids in the haemagglutinin can moderate the magnitude of antigenic change.

Antibody Determinants of Influenza Immunity

Understanding antigenic variation in influenza virus strains and how the human immune system recognizes strains are central challenges for vaccinologists. Antibodies directed to the 2 major viral surface membrane proteins, hemagglutinin (HA) and neuraminidase (NA), mediate protection against reinfection following natural infection or vaccination, but HA and NA protein sequences in field strains are highly variable. The central questions are how to achieve protective antibody responses in a higher proportion of individuals and how to induce responses with more breadth and durability. Studies using isolation of human monoclonal antibodies followed by structural and functional characterization revealed conserved antigenic sites recognized by broadly cross-reactive antibodies. The antigenic landscape on HA and NA proteins is coming into focus to inform studies of the correlates and mechanisms of immunity. Understanding the antibody determinants of influenza immunity points the way toward development and testing of next-generation vaccines with potential to confer broadly protective immunity.

Antibody responses to viral infections: a structural perspective across three different enveloped viruses

Antibodies serve as critical barriers to viral infection. Humoral immunity to a virus is achieved through the dual role of antibodies in communicating the presence of invading pathogens in infected cells to effector cells and interfering with processes essential to the viral lifecycle, chiefly entry into the host cell. For individuals that successfully control infection, virus-elicited antibodies can provide lifelong surveillance and protection from future insults. One approach to understand the nature of a successful immune response has been to utilize structural biology to uncover the molecular details of the antibodies derived from vaccines or natural infection and how they interact with their cognate microbial antigens. The ability to isolate antigen specific B-cells and rapidly solve structures of functional, monoclonal antibodies in complex with viral glycoprotein surface antigens has greatly expanded our knowledge of the sites of vulnerability on viruses. In this review, we compare the adaptive humoral immune responses to HIV, influenza, and filoviruses, with a particular focus on neutralizing antibodies. The pathogenesis of each of these viruses is quite different, providing an opportunity for comparison of immune responses: HIV causes a persistent, chronic infection; influenza an acute infection with multiple exposures during a lifetime and annual vaccination; and filoviruses, a virulent, acute infection. Neutralizing antibodies that develop under these different constraints are therefore sentinels that can provide insight into the underlying humoral immune responses and important lessons to guide future development of vaccines and immunotherapeutics.