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

  1. Angela M Phillips
  2. Katherine R Lawrence
  3. Alief Moulana
  4. Thomas Dupic
  5. Jeffrey Chang
  6. Milo S Johnson
  7. Ivana Cvijovic
  8. Thierry Mora
  9. Aleksandra M Walczak
  10. Michael M Desai  Is a corresponding author
  1. Harvard University, United States
  2. Massachusetts Institute of Technology, United States
  3. Stanford University, United States
  4. École Normale Supérieure, France

Abstract

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

Data availability

Data and code used for this study are available at https://github.com/klawrence26/bnab-landscapes. CR9114 data are also available in an interactive data browser at https://yodabrowser.netlify.app/yoda_browser/.FASTQ files from high-throughput sequencing have been deposited in the NCBI BioProject database with accession number PRJNA741613, and will be publicly released upon acceptance.

The following data sets were generated

Article and author information

Author details

  1. Angela M Phillips

    Harvard University, Cambridge, United States
    Competing interests
    No competing interests declared.
  2. Katherine R Lawrence

    Physics, Massachusetts Institute of Technology, Cambridge, United States
    Competing interests
    No competing interests declared.
  3. Alief Moulana

    Harvard University, Cambridge, United States
    Competing interests
    No competing interests declared.
  4. Thomas Dupic

    Harvard University, Cambridge, United States
    Competing interests
    No competing interests declared.
  5. Jeffrey Chang

    Harvard University, Cambridge, United States
    Competing interests
    No competing interests declared.
  6. Milo S Johnson

    Harvard University, Cambridge, United States
    Competing interests
    No competing interests declared.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-0169-2494
  7. Ivana Cvijovic

    Stanford University, Stanford, United States
    Competing interests
    No competing interests declared.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-6272-2979
  8. Thierry Mora

    École Normale Supérieure, Paris, France
    Competing interests
    No competing interests declared.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-5456-9361
  9. Aleksandra M Walczak

    École Normale Supérieure, Paris, France
    Competing interests
    Aleksandra M Walczak, Senior editor, eLife.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-2686-5702
  10. Michael M Desai

    Harvard University, Cambridge, United States
    For correspondence
    mdesai@oeb.harvard.edu
    Competing interests
    No competing interests declared.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-9581-1150

Funding

Howard Hughes Medical Institute (Hanna H. Gray Postdoctoral Fellowship)

  • Angela M Phillips

International Human Frontier Science Program Organization

  • Thomas Dupic

The NSF-Simons Center for Mathematical and Statistical Analysis of Biology at Harvard (1764269)

  • Katherine R Lawrence

Hertz Foundation (Graduate Fellowship Award)

  • Katherine R Lawrence

National Science Foundation (Graduate Research Fellowship Program)

  • Katherine R Lawrence

National Science Foundation (Graduate Research Fellowship Program)

  • Milo S Johnson

National Science Foundation (Graduate Research Fellowship Program)

  • Jeffrey Chang

European Research Council (COG 724208)

  • Thierry Mora
  • Aleksandra M Walczak

National Institutes of Health (GM104239)

  • Michael M Desai

National Science Foundation (PHY-1914916)

  • Michael M Desai

Stanford University (Stanford Science Fellowship)

  • Ivana Cvijovic

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Copyright

© 2021, Phillips et al.

This article is distributed under the terms of the Creative Commons Attribution License permitting unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 4,880
    views
  • 548
    downloads
  • 52
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Angela M Phillips
  2. Katherine R Lawrence
  3. Alief Moulana
  4. Thomas Dupic
  5. Jeffrey Chang
  6. Milo S Johnson
  7. Ivana Cvijovic
  8. Thierry Mora
  9. Aleksandra M Walczak
  10. Michael M Desai
(2021)
Binding affinity landscapes constrain the evolution of broadly neutralizing anti-influenza antibodies
eLife 10:e71393.
https://doi.org/10.7554/eLife.71393

Share this article

https://doi.org/10.7554/eLife.71393

Further reading

    1. Evolutionary Biology
    2. Microbiology and Infectious Disease
    Zach Hensel
    Short Report

    Accurate estimation of the effects of mutations on SARS-CoV-2 viral fitness can inform public-health responses such as vaccine development and predicting the impact of a new variant; it can also illuminate biological mechanisms including those underlying the emergence of variants of concern. Recently, Lan et al. reported a model of SARS-CoV-2 secondary structure and its underlying dimethyl sulfate reactivity data (Lan et al., 2022). I investigated whether base reactivities and secondary structure models derived from them can explain some variability in the frequency of observing different nucleotide substitutions across millions of patient sequences in the SARS-CoV-2 phylogenetic tree. Nucleotide basepairing was compared to the estimated ‘mutational fitness’ of substitutions, a measurement of the difference between a substitution’s observed and expected frequency that is correlated with other estimates of viral fitness (Bloom and Neher, 2023). This comparison revealed that secondary structure is often predictive of substitution frequency, with significant decreases in substitution frequencies at basepaired positions. Focusing on the mutational fitness of C→U, the most common type of substitution, I describe C→U substitutions at basepaired positions that characterize major SARS-CoV-2 variants; such mutations may have a greater impact on fitness than appreciated when considering substitution frequency alone.

    1. Evolutionary Biology
    Yiheng Zhang, Xing Wang ... Xiaoguang Yang
    Research Article

    Although fossil evidence suggests the existence of an early muscular system in the ancient cnidarian jellyfish from the early Cambrian Kuanchuanpu biota (ca. 535 Ma), south China, the mechanisms underlying the feeding and respiration of the early jellyfish are conjectural. Recently, the polyp inside the periderm of olivooids was demonstrated to be a calyx-like structure, most likely bearing short tentacles and bundles of coronal muscles at the edge of the calyx, thus presumably contributing to feeding and respiration. Here, we simulate the contraction and expansion of the microscopic periderm-bearing olivooid Quadrapyrgites via the fluid-structure interaction computational fluid dynamics (CFD) method to investigate their feeding and respiratory activities. The simulations show that the rate of water inhalation by the polyp subumbrella is positively correlated with the rate of contraction and expansion of the coronal muscles, consistent with the previous feeding and respiration hypothesis. The dynamic simulations also show that the frequent inhalation/exhalation of water through the periderm polyp expansion/contraction conducted by the muscular system of Quadrapyrgites most likely represents the ancestral feeding and respiration patterns of Cambrian sedentary medusozoans that predated the rhythmic jet-propelled swimming of the modern jellyfish. Most importantly for these Cambrian microscopic sedentary medusozoans, the increase of body size and stronger capacity of muscle contraction may have been indispensable in the stepwise evolution of active feeding and subsequent swimming in a higher flow (or higher Reynolds number) environment.