Implications of localized charge for human influenza A H1N1 hemagglutinin evolution: Insights from deep mutational scans

Seasonal influenza A viruses of humans evolve rapidly due to strong selection pressures from host immune responses, principally on the hemagglutinin (HA) viral surface protein. Based on mouse transmission experiments, a proposed mechanism for immune evasion consists of increased avidity to host cellular receptors, mediated by electrostatic charge interactions with negatively charged cell surfaces. In support of this, the HA charge of the globally circulating H3N2 has increased over time since its pandemic. However, the same trend was not seen in H1N1 HA sequences. This is counter-intuitive, since immune escape due to increased avidity (due itself to an increase in charge) was determined experimentally. Here, we explore whether patterns of local charge of H1N1 HA can explain this discrepancy and thus further associate electrostatic charge with immune escape and viral evolutionary dynamics. Measures of site-wise functional selection and expected charge computed from deep mutational scan data on an early H1N1 HA yield a striking division of residues into three groups, separated by charge. We then explored evolutionary dynamics of these groups from 1918 to 2008. In particular, one group increases in net charge over time and consists of sites that are evolving the fastest, that are closest to the receptor binding site (RBS), and that are exposed to solvent (i.e., on the surface). By contrast, another group decreases in net charge and consists of sites that are further away from the RBS and evolving slower, but also exposed to solvent. The last group consists of those sites in the HA core, with no change in net charge and that evolve very slowly. Thus, there is a group of residues that follows the same trend as seen for the entire H3N2 HA. It is possible that the H1N1 HA is under other biophysical constraints that result in compensatory decreases in charge elsewhere on the protein. Our results implicate localized charge in HA interactions with host cells, and highlight how deep mutational scan data can inform evolutionary hypotheses.

[1]  Steven J. Gamblin,et al.  Influenza Hemagglutinin and Neuraminidase Membrane Glycoproteins , 2010, The Journal of Biological Chemistry.

[2]  J. Rodrigues,et al.  pdb-tools: a swiss army knife for molecular structures , 2018, bioRxiv.

[3]  Nicole M. Bouvier,et al.  The biology of influenza viruses. , 2008, Vaccine.

[4]  L. Kasturi,et al.  The Amino Acid at the X Position of an Asn-X-Ser Sequon Is an Important Determinant of N-Linked Core-glycosylation Efficiency (*) , 1996, The Journal of Biological Chemistry.

[5]  S. Hensley,et al.  Immune history and influenza virus susceptibility. , 2017, Current opinion in virology.

[6]  Rommie E. Amaro,et al.  Human Influenza A Virus Hemagglutinin Glycan Evolution Follows a Temporal Pattern to a Glycan Limit , 2019, mBio.

[7]  J. Bloom.  An Experimentally Informed Evolutionary Model Improves Phylogenetic Fit to Divergent Lactamase Homologs , 2014, bioRxiv.

[8]  N. Ben-Tal,et al.  Comparison of site-specific rate-inference methods for protein sequences: empirical Bayesian methods are superior. , 2004, Molecular biology and evolution.

[9]  Michael B. Doud,et al.  Accurate Measurement of the Effects of All Amino-Acid Mutations on Influenza Hemagglutinin , 2016, Viruses.

[10]  Trevor Bedford,et al.  Deep mutational scanning of hemagglutinin helps predict evolutionary fates of human H3N2 influenza variants , 2018, Proceedings of the National Academy of Sciences.

[11]  Nimalan Arinaminpathy,et al.  Dynamics of Glycoprotein Charge in the Evolutionary History of Human Influenza , 2010, PloS one.

[12]  Sara Linse,et al.  Salting the charged surface: pH and salt dependence of protein G B1 stability. , 2006, Biophysical journal.

[13]  E. C. Pielou.  The measurement of diversity in different types of biological collections , 1966 .

[14]  Walter Krämer,et al.  Review of Modern applied statistics with S, 4th ed. by W.N. Venables and B.D. Ripley. Springer-Verlag 2002 , 2003 .

[15]  Austin G. Meyer,et al.  Geometric constraints dominate the antigenic evolution of 1 influenza H 3 N 2 hemagglutinin 2 , 2015 .

[16]  D. Fowler,et al.  Deep mutational scanning: assessing protein function on a massive scale. , 2011, Trends in biotechnology.

[17]  Yuji Ikegaya,et al.  Genetically Encoded Green Fluorescent Ca2+ Indicators with Improved Detectability for Neuronal Ca2+ Signals , 2012, PloS one.

[18]  Jonathan Dushoff,et al.  Ecology and evolution of the flu , 2002 .

[19]  Michael W. Deem,et al.  Selective Pressure to Increase Charge in Immunodominant Epitopes of the H3 Hemagglutinin Influenza Protein , 2010, Journal of Molecular Evolution.

[20]  Trevor Bedford,et al.  Global circulation patterns of seasonal influenza viruses vary with antigenic drift , 2015, Nature.

[21]  R. Webster,et al.  Evolution and ecology of influenza A viruses. , 1992, Current topics in microbiology and immunology.

[22]  O. Pybus,et al.  Unifying the Epidemiological and Evolutionary Dynamics of Pathogens , 2004, Science.

[23]  Xiaomei Li,et al.  Influenza Research Database: An integrated bioinformatics resource for influenza virus research , 2016, Nucleic Acids Res..

[24]  Jesse D. Bloom,et al.  Software for the analysis and visualization of deep mutational scanning data , 2015, bioRxiv.

[25]  William N. Venables,et al.  Modern Applied Statistics with S , 2010 .

[26]  Ivo M Foppa,et al.  Annual estimates of the burden of seasonal influenza in the United States: A tool for strengthening influenza surveillance and preparedness , 2018, Influenza and other respiratory viruses.

[27]  Rahul Raman,et al.  Hemagglutinin Receptor Binding Avidity Drives Influenza A Virus Antigenic Drift , 2009, Science.

[28]  C. M. Saad-Roy,et al.  Dynamic Perspectives on the Search for a Universal Influenza Vaccine. , 2019, The Journal of infectious diseases.

[29]  J. H. Ward.  Hierarchical Grouping to Optimize an Objective Function , 1963 .

[30]  Fionn Murtagh,et al.  Ward’s Hierarchical Agglomerative Clustering Method: Which Algorithms Implement Ward’s Criterion? , 2011, Journal of Classification.

[31]  Yoshiyuki Suzuki,et al.  Compensatory Evolution of Net-Charge in Influenza A Virus Hemagglutinin , 2012, PloS one.

[32]  P. Wilson,et al.  The influenza virus hemagglutinin head evolves faster than the stalk domain , 2018, Scientific Reports.

[33]  Nicole M. Bouvier,et al.  Antigenic sites in influenza H1 hemagglutinin display species-specific immunodominance , 2018, The Journal of clinical investigation.

[34]  David F. Burke,et al.  Koel Antigenic Change During Influenza Virus Evolution Substitutions Near the Receptor Binding Site Determine Major , 2013 .

[35]  Hadley Wickham,et al.  ggplot2 - Elegant Graphics for Data Analysis (2nd Edition) , 2017 .

[36]  J. Yewdell,et al.  The antigenic structure of the influenza virus A/PR/8/34 hemagglutinin (H1 subtype) , 1982, Cell.

[37]  C. Bender,et al.  Differences in the biological phenotype of low-yielding (L) and high-yielding (H) variants of swine influenza virus A/NJ/11/76 are associated with their different receptor-binding activity. , 1998, Virology.

[38]  D. J. Stevens,et al.  The Structure and Receptor Binding Properties of the 1918 Influenza Hemagglutinin , 2004, Science.

[39]  O. AltmanMeghan,et al.  Antibody Immunodominance: The Key to Understanding Influenza Virus Antigenic Drift , 2018 .

[40]  Matthew Angel,et al.  Defining B Cell Immunodominance to Viruses , 2017, Nature Immunology.