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.

Reviewing Editor

  1. Sarel Jacob Fleishman, Weizmann Institute of Science, Israel

Version history

  1. Preprint posted: May 25, 2021 (view preprint)
  2. Received: June 17, 2021
  3. Accepted: September 5, 2021
  4. Accepted Manuscript published: September 7, 2021 (version 1)
  5. Version of Record published: September 27, 2021 (version 2)
  6. Version of Record updated: November 3, 2021 (version 3)

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

  • 3,574
    views
  • 428
    downloads
  • 30
    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. Immunology and Inflammation
    Mark S Lee, Peter J Tuohy ... Michael S Kuhns
    Research Advance

    CD4+ T cell activation is driven by five-module receptor complexes. The T cell receptor (TCR) is the receptor module that binds composite surfaces of peptide antigens embedded within MHCII molecules (pMHCII). It associates with three signaling modules (CD3γε, CD3δε, and CD3ζζ) to form TCR-CD3 complexes. CD4 is the coreceptor module. It reciprocally associates with TCR-CD3-pMHCII assemblies on the outside of a CD4+ T cells and with the Src kinase, LCK, on the inside. Previously, we reported that the CD4 transmembrane GGXXG and cytoplasmic juxtamembrane (C/F)CV+C motifs found in eutherian (placental mammal) CD4 have constituent residues that evolved under purifying selection (Lee et al., 2022). Expressing mutants of these motifs together in T cell hybridomas increased CD4-LCK association but reduced CD3ζ, ZAP70, and PLCγ1 phosphorylation levels, as well as IL-2 production, in response to agonist pMHCII. Because these mutants preferentially localized CD4-LCK pairs to non-raft membrane fractions, one explanation for our results was that they impaired proximal signaling by sequestering LCK away from TCR-CD3. An alternative hypothesis is that the mutations directly impacted signaling because the motifs normally play an LCK-independent role in signaling. The goal of this study was to discriminate between these possibilities. Using T cell hybridomas, our results indicate that: intracellular CD4-LCK interactions are not necessary for pMHCII-specific signal initiation; the GGXXG and (C/F)CV+C motifs are key determinants of CD4-mediated pMHCII-specific signal amplification; the GGXXG and (C/F)CV+C motifs exert their functions independently of direct CD4-LCK association. These data provide a mechanistic explanation for why residues within these motifs are under purifying selection in jawed vertebrates. The results are also important to consider for biomimetic engineering of synthetic receptors.

    1. Evolutionary Biology
    Robert Horvath, Nikolaos Minadakis ... Anne C Roulin
    Research Article

    Understanding how plants adapt to changing environments and the potential contribution of transposable elements (TEs) to this process is a key question in evolutionary genomics. While TEs have recently been put forward as active players in the context of adaptation, few studies have thoroughly investigated their precise role in plant evolution. Here, we used the wild Mediterranean grass Brachypodium distachyon as a model species to identify and quantify the forces acting on TEs during the adaptation of this species to various conditions, across its entire geographic range. Using sequencing data from more than 320 natural B. distachyon accessions and a suite of population genomics approaches, we reveal that putatively adaptive TE polymorphisms are rare in wild B. distachyon populations. After accounting for changes in past TE activity, we show that only a small proportion of TE polymorphisms evolved neutrally (<10%), while the vast majority of them are under moderate purifying selection regardless of their distance to genes. TE polymorphisms should not be ignored when conducting evolutionary studies, as they can be linked to adaptation. However, our study clearly shows that while they have a large potential to cause phenotypic variation in B. distachyon, they are not favored during evolution and adaptation over other types of mutations (such as point mutations) in this species.