Research Advance

Nanobody repertoire generated against the spike protein of ancestral SARS-CoV-2 remains efficacious against the rapidly evolving virus

Laboratory of Cellular and Structural Biology, The Rockefeller University, United States
Center for Global Infectious Disease Research, Seattle Children's Research Institute, United States
Department of Bioengineering and Therapeutic Sciences, Department of Pharmaceutical Chemistry, California Institute for Quantitative Biosciences, Byers Hall, University of California, San Francisco, United States
Laboratory of Mass Spectrometry and Gaseous Ion Chemistry, The Rockefeller University, United States
Department of Pediatrics, University of Washington, United States
Department of Biochemistry, University of Washington, United States

May 7, 2024

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

Open access
Copyright information

Version of Record

The authors declare this version of their article to be the Version of Record.

About eLife's process

Download
Cite
Share
CommentOpen annotations (there are currently 0 annotations on this page).

Version of Record published: May 7, 2024 (This version)
Reviewed preprint version 2: April 15, 2024 (Go to version)
Reviewed preprint version 1: November 9, 2023 (Go to version)
Preprint posted: July 14, 2023 (Go to version)
Sent for peer review: June 21, 2023

1. Builds upon
Highly synergistic combinations of nanobodies that target SARS-CoV-2 and are resistant to escape

Fred D Mast, Peter C Fridy ... Michael P Rout

Research Article Dec 7, 2021
Further reading

Abstract
eLife assessment
Introduction
Results and discussion
Materials and methods
Data availability
References
Article and author information
Metrics

Abstract

To date, all major modes of monoclonal antibody therapy targeting SARS-CoV-2 have lost significant efficacy against the latest circulating variants. As SARS-CoV-2 omicron sublineages account for over 90% of COVID-19 infections, evasion of immune responses generated by vaccination or exposure to previous variants poses a significant challenge. A compelling new therapeutic strategy against SARS-CoV-2 is that of single-domain antibodies, termed nanobodies, which address certain limitations of monoclonal antibodies. Here, we demonstrate that our high-affinity nanobody repertoire, generated against wild-type SARS-CoV-2 spike protein (Mast et al., 2021), remains effective against variants of concern, including omicron BA.4/BA.5; a subset is predicted to counter resistance in emerging XBB and BQ.1.1 sublineages. Furthermore, we reveal the synergistic potential of nanobody cocktails in neutralizing emerging variants. Our study highlights the power of nanobody technology as a versatile therapeutic and diagnostic tool to combat rapidly evolving infectious diseases such as SARS-CoV-2.

eLife assessment

This study presents important insights on the impact of SARS-CoV-2 variants on the binding and neutralization of a small library of nanobodies. The authors should be applauded for their comprehensive in vitro and in silico analyses of nanobody targeting of SARS-CoV-2 variants. The evidence supporting the claims of the authors is now convincing. This work will be of great interest to researchers in the fields of antibody/nanobody engineering and SARS-CoV-2 therapeutics.

https://doi.org/10.7554/eLife.89423.3.sa0

About eLife assessments

Introduction

SARS-CoV-2 has infected >40% of the world’s population (COVID-19 Cumulative Infection Collaborators, 2022)resulting in a devastating loss of life. As the SARS-CoV-2 pandemic enters its endemic phase (Meng et al., 2023; Pilz and Ioannidis, 2023), multiple new variants continue to circulate. Since its initial spread, the rapid adaptation of the virus to selective pressures continues to produce variants of concern (VoC), of which the omicron variants presently account for over 90% of current SARS-CoV-2 infections (https://www.cdc.gov). SARS-CoV-2 displays three structural proteins that are potential targets for therapeutic intervention, but the primary focal point of vaccine development and many therapeutic strategies is the spike surface glycoprotein, which the virus uses to gain cell entry by attaching to the host cell angiotensin-converting enzyme 2 (ACE2) receptor (Jackson et al., 2020; Krammer, 2020; Letko et al., 2020; Polack et al., 2020). The spike protein trimer consists of three domains: the receptor-binding domain (RBD) on S1 that binds ACE2, the S1 N-terminal domain (NTD) that has a poorly defined function, and the S2 domain that is involved in virus–host cell membrane fusion (Walls et al., 2020; Jackson et al., 2022). Glycosylation is most extensive on the NTD and the S2 domain, whereas the RBD is largely glycan free (Watanabe et al., 2020; Zhao et al., 2020). Consequently, it is unsurprising that the most antigenic domain on spike is the RBD, where the majority of neutralizing antibodies have been shown to bind. A comprehensive mapping of the epitopes from 1640 neutralizing monoclonal antibodies (mAbs), all targeting the RBD, revealed 12 epitope groups (Cao et al., 2022). These data, combined with previous studies mapping the epitopes of antibodies targeting spike, reveal a total of 19 mAb epitope groups, including seven on the NTD (Wang et al., 2022d). Very few anti-S2 antibodies have been shown to be effective therapeutic options (Wec et al., 2020), likely due to the shielding effect of S2 glycans (Grant et al., 2020).

The mechanism by which omicron variants of SARS-CoV-2 (e.g., BA.1, BA.4, BA.5, and XBB) escape the neutralizing abilities of antibodies generated against spike proteins from preceding variants, whether by vaccination or infection, is largely attributed to the extensive number of mutations accumulated in spike (Greaney et al., 2021a; Starr et al., 2022; Dadonaite et al., 2023). Compared to wild-type SARS-CoV-2 spike, omicron BA.1 spike has 37 amino acid residue differences, with almost half located in the RBD domain (Mannar et al., 2022). The omicron BA.4/BA.5 variants, which have identical spike proteins (Tegally et al., 2022), have additional mutations (including the L452R substitution first seen in the delta variant), that render many previously broadly neutralizing antibodies ineffective (Cao et al., 2022; Hachmann et al., 2022; Wang et al., 2022b). Of the mAbs that previously received emergency use authorization (EUA) by the FDA for the treatment of SARS-CoV-2 infection, even cilgavimab and bebtelovimab that were, respectively, moderately and highly efficacious against omicron BA.5, are no longer effective against the current circulating variants XBB, BQ.1.1, and related sublineages (Takashita et al., 2022; Focosi et al., 2023; Imai et al., 2023). As a result, no mAb therapy is currently approved by the FDA for treatment of SARS-CoV-2 infection (https://www.fda.gov).

Nanobodies, single-domain antibodies derived from a unique heavy chain-only class of llama antibodies, present numerous therapeutic benefits compared to mAbs. Their smaller size and increased stability make them more resistant to denaturation, simpler to produce, and easier to modify in order to adjust properties such as immunogenicity and half-life (Muyldermans, 2013). One potential advantage for neutralizing the spike protein is their compact size and distinctive binding attributes, which allow them to access and bind to epitopes that mAbs cannot reach. Consequently, while the antigenic evolution of the spike protein in response to antibodies has largely rendered mAbs ineffective in a therapeutic context, it remains uncertain how this applies to nanobodies. Moreover, the diminutive size of nanobodies enables them to bind concurrently to a single antigen through non-overlapping epitopes, making them well suited for creating nanobody mixtures with the potential for highly synergistic effects (Fridy et al., 2014; Mast et al., 2021).

Here, we demonstrate that a subset of our previously published repertoire of nanobodies, generated against spike from the ancestral SARS-CoV-2 virus (Mast et al., 2021), retains binding and in vitro neutralization efficacy against circulating VoC, including omicron BA.4/BA.5. We show the power of nanobodies when working synergistically to create potent neutralizing mixtures against the different VoCs. We also predict that a subset of these nanobodies will remain efficacious against the circulating XBB and BQ.1.1 sublineages. Our study underscores the importance and versatility of large, diverse repertoires of nanobodies, in their potential to create long-term therapeutic options against rapidly evolving infectious agents such as the SARS-CoV-2 virus.

Results and discussion

Nanobodies generated against wild-type SARS-CoV-2 spike remains efficacious against delta, and omicron lineages BA.1, BA.4/BA.5, XBB, and BQ.1.1

From the original nanobody repertoire that we generated against SARS-CoV-2 wild-type spike protein (Mast et al., 2021), representative nanobodies from all 10 structurally mapped epitope groups that we previously identified (Mast et al., 2021; Cross et al., 2023), were selected for SARS-CoV-2 pseudovirus (PSV) neutralization assays against the SARS-CoV-2 delta and omicron BA.1 strains (Figure 1; Figure 2—source data 1). Of the 41 nanobodies tested, 35 remained efficacious against at least one variant, where 28 neutralized delta, 23 neutralized omicron BA.1, and 15 neutralized both. The RBD groups I, I/II, II, I/IV, and IV – where nanobodies whose epitopes could not be distinguished between two groups are demarcated I/II and I/IV – and the anti-S2 groups (groups IX and X) contain a high number of nanobodies that neutralized delta. We note that IC₅₀s are not directly comparable across different experimental setups because measured values are highly dependent on the experimental conditions. For this reason, we included other published nanobodies as benchmarks in our original publication (Mast et al., 2021) and have subsequently maintained standard experimental conditions. Additionally, for computational epitope modeling, we selected nanobody candidates using a series of experimentally obtained structural restraints, as described in Mast et al., 2021.

Figure 1

Download asset Open asset

Nanobody repertoires generated against wild-type SARS-CoV-2 remain efficacious.

Nanobodies targeting the S1-RBD, S1 non-RBD, and S2 regions of spike effectively neutralize lentivirus pseudotyped with delta and omicron BA.1 SARS-CoV-2 spikes (PSV) from infecting angiotensin-converting enzyme 2 (ACE2)-expressing HEK293T cells. (A) The half-maximal inhibitory concentration (IC₅₀) is reported for the indicated nanobodies against wild-type (Mast et al., 2021), delta, and omicron BA.1 PSV. These values are summarized in Figure 2—source data 1. Nanobodies are grouped by epitope and arranged within each epitope by neutralization efficacy against the wild-type PSV. n ≥ 4 (B) The structural differences in the receptor-binding domain (RBD) of the delta (PDB ID: 7SBO) and omicron BA.1 (PDB ID: 7T9K) variants are depicted. Nanobody epitopes are heat-mapped ranging from pale white (epitopes with weak neutralization against SARS-CoV-2) to dark red (indicating strong neutralization). Boxed in gray are mutations specific to each variant mapped in blue on the aforementioned structures. The nanobodies that contributed to epitope mapping are in bold in panel A. The color bar scale for each epitope is the neutralizing strength of each nanobody epitope, calculated as the normalized −log10 ratio of nanobody binding (IC₅₀) to variant versus wild-type SARS-CoV-2 Spike S1. For groups with multiple nanobodies, the average −log10 (IC₅₀) is first calculated for the nanobodies within that group, then normalized to a neutralization score within the 0–100 range using the min and max average −log10 (IC₅₀) for that group. A higher score indicates more potent neutralization of the variant relative to the wild-type. All structural representations were created on ChimeraX (Pettersen et al., 2021).

These observations align with data from mAbs approved by the FDA, where nanobody epitopes from groups I, II, and IV (Figure 1B) overlap with the epitopes of five mAbs effective against delta, categorized into two of three classes on RBD: class 1 (etesivimab, casirivimab, and amubarvimab; Greaney et al., 2021b; Li et al., 2021; Planas et al., 2021; Takashita et al., 2022; Cox et al., 2023) and class 3 (imdemivab and bebtelovimab; Cox et al., 2023; Figure 2). Notably, most of our group I and II nanobodies neutralized delta better than wild-type (Figure 1), paralleling observations made with the mAb etesivimab (Wang et al., 2022c, Cox et al., 2023). In comparison, groups I, I/II, I/IV, V, VII, VIII, and the anti-S2 nanobodies contained the majority of omicron BA.1 neutralizers, though here the neutralization potency of many nanobodies was generally decreased tenfold compared to wild-type. This decrease in neutralization potency largely correlates with the accumulation of omicron BA.1-specific mutations throughout the RBD, which potentially alter the nanobodies’ binding sites, weakening their interaction with the BA.1 spike (Figure 1B). Concomitantly, groups I, I/II, I/IV, and the anti-S2 groups contain nanobodies able to neutralize both delta and omicron BA.1. These results demonstrate the effectiveness of our original nanobody cohort against delta and omicron BA.1, targeting all major regions of spike.

Figure 2

Download asset Open asset

Nanobody epitope groups and mAb epitope classes mapped on receptor-binding domain (RBD).

(A) Nanobody epitope groups overlapping with the three mAb epitope classes (classes 1, 2, and 3). Nanobody groups are highlighted in gold, while mAb class footprints are outlined in black. mAb epitopes are taken from Cox et al. (B) A single RBD subunit with the angiotensin-converting enzyme 2 (ACE2) footprint/RBM mapped in green. All epitopes are represented on the structure of wild-type RBD (PDB ID: 6M0J). All structure representations were generated using ChimeraX (Pettersen et al.).

Figure 2—source data 1 Nanobody binding and neutralization characterization; related to Figures 2 and 3.: https://cdn.elifesciences.org/articles/89423/elife-89423-fig2-data1-v1.xlsx
Download elife-89423-fig2-data1-v1.xlsx

Of the nanobodies that neutralized both delta and omicron BA.1, representatives from each of the nanobody epitope groups were selected for surface plasmon resonance (SPR) analysis, where S1 binders with mapped epitopes that neutralized one or both variants well were prioritized. SPR-binding assessments to the spike S1 domain or RBD of delta revealed a pattern: nanobodies maintaining binding affinity generally also neutralized the virus with a statistically significant correlation between binding affinity and neutralization efficacy (Pearson’s correlation coefficient: 0.71, p-value: 0.01; Spearman’s rho: 0.63, p-value: 0.07). However, this correlation was not statistically significant for omicron BA.1 (Pearson’s correlation coefficient: 0.27, p-value: 0.31) (Figure 3A). Notably, while some nanobodies bound to the variants, they did not consistently neutralize them, suggesting additional factors influence neutralization beyond mere binding. We also tested this cohort for binding to omicron BA.4/BA.5 RBD using SPR, revealing that almost all the nanobodies in group I, I/II, II, and I/IV retained binding to omicron BA.4/BA.5. Based on our previously mapped nanobody epitopes on spike (Mast et al., 2021), these four groups appear to at least partially overlap with the spike epitope of the potent omicron BA.4/5 neutralizing mAb formerly FDA-approved, bebtelovimab (Figure 2; Focosi et al., 2023). The nanobodies that retained binding to omicron BA.4/5 were further tested for binding against the omicron XBB and BQ.1.1 lineages in addition to S1-1 and S1-23, revealing nanobodies that bound omicron BA.4/BA.5, aside from S1-39, also bound omicron variants XBB and BQ.1.1 (Figure 3A).

Figure 3

Download asset Open asset

Affinities of the nanobody repertoire against SARS-CoV-2 variants.

(A) Each nanobody is plotted against their affinity (K_D) for SARS-CoV-2 Spike S1 from wild-type, delta and omicron BA.1, BA.4, XBB, and BQ.1.1 strains. Kinetic values are summarized in . Nanobodies are characterized into their respective epitope groups as described previously (Mast et al., 2021). (B) Displayed are the structures of the receptor-binding domain (RBD) of spike delta (PDB ID: 7SBO), omicron BA.1 (PDB ID: 7T9K), omicron BA.4 RBD modeled using AlphaFold (Jumper et al., 2021), omicron XBB (PDB ID: 8IOU), and omicron BQ.1.1 (PDB ID: 8FXC). The structures feature heat-mapped epitopes of binding, ranging from pale white (weak binding to SARS-CoV-2) to dark red (strong binding to SARS-CoV-2). In the gray box, mutations specific to each variant are highlighted in blue. The nanobodies that contributed to epitope mapping are in bold in panel A. The color bar scale indicates each epitope’s binding affinity strength, represented as the normalized −log10 ratio of nanobody binding (K_D) of variant versus wild-type SARS-CoV-2 Spike S1. For groups with multiple nanobodies, the average −log10 (K_D) for the nanobodies within that group was calculated, then normalized to an affinity score ranging from 0 to 100 using the min and max average −log10 (K_D) for that group. Higher −log10 ratios indicate stronger binding of the nanobody to the variant versus wild-type. S1-RBD-16 bound omicron BA.1 and BA.4/5 in ELISA. S1-RBD-11 was not tested against omicron BA.4. S1-65 was not tested against BA.1. Only S1-1, S1-RBD-22, S1-RBD-9, S1-4, S1-35, S1-39, S1-RBD-6, S1-RBD-5, S1-23, S1-RBD-40, S1-RBD-46, and S1-46 were tested against omicron XBB and BQ.1.1. All structure representations were generated using ChimeraX (Pettersen et al., 2021).

The remaining six nanobody groups tested, including the major nanobody groups III and IV, showed no detectable binding to omicron BA.4/BA.5/XBB/BQ.1.1 RBD, save for RBD-40 (group III) and RBD-46 (group IV), which exhibited a ~40- and ~20-fold decrease in affinity, respectively, compared to wild-type. The diminished binding likely results from the concentration of omicron BA.4/BA.5/XBB/BQ.1.1-specific mutations that overlap with the epitope regions of the affected nanobodies. These mutations sufficiently alter the epitopes, to either abolish or significantly reduce binding (Figure 3B). Interestingly, one nanobody, S1-46 (group VII, Figures 1B and 3B) retained wild-type binding affinity to RBDs from delta and all four omicron variants tested (Figure 2). S1-46 binds a region on spike that is conserved across all variants to date, but which may be relatively inaccessible unless the RBDs are in the ‘up’ conformation. The epitope of S1-46 is not targeted by any of the mAbs that previously received EUA by the FDA (Cox et al., 2023). Guided by the SPR results against omicron BA.4/BA.5 RBD, nanobody neutralization of live BA.5 was performed using the plaque reduction neutralization test as previously described (Mast et al., 2021). All five nanobodies tested neutralized the live omicron BA.5 live virus (Figure 4), corroborating our SPR observations. This indicates that our nanobody repertoire generated against wild-type spike has retained efficacy against omicron BA.4/BA.5.

Figure 4

Download asset Open asset

Potent neutralization by broadly active nanobodies.

Nanobodies targeting the S1-RBD of spike and raised against the original wild-type sequence remain highly efficacious in neutralizing an evolved variant of SARS-CoV-2, omicron variant BA.5. The derived neutralization curves are plotted from the results of a plaque-forming reduction neutralization test with the indicated nanobodies. Serial dilutions of each nanobody were incubated with ~200 SARS-CoV-2 virions for 60 min and then overlaid on a monolayer of TMPRSS2-expressing Vero E6 cells. After 72 hr, cells were fixed and stained with crystal violet stain (1% wt/vol in 20% ethanol) allowing for the enumeration of viral plaques. The percent plaque inhibition for each nanobody dilution, summarized in Figure 4—source data 1, was used to fit the neutralization curves depicted in the figure. The colored shaded areas denote 90% confidence intervals for each fitted curve. n ≥ 3.

Figure 4—source data 1 Neutralization data from the plaque reduction neutralization test (PRNT) assay with omicron BA.5 virus.: https://cdn.elifesciences.org/articles/89423/elife-89423-fig4-data1-v1.csv
Download elife-89423-fig4-data1-v1.csv

Impact of spike structural differences across variants on nanobody binding and neutralization potency

The numerous structural differences observed in the spike protein of the delta and omicron sublineages compared to wild-type, have enabled these variants to escape the neutralizing effect of most mAbs, with only one clinically approved mAb retaining potency against omicron BA.4/BA.5 (Takashita et al., Cox et al., Focosi et al.). These structural differences have also affected our nanobody cohort, likely playing a role in the differential binding and neutralizing abilities observed against the three tested variants (delta, omicron BA.1, and omicron BA.4/5). Notably, all but one of our nanobodies in groups I, I/II, and II, displayed better neutralizing ability against delta compared to wild-type (Figure 1). It is possible that this difference is due to the delta spike trimer’s preference for the up orientation of its RBDs, resulting from the increased dynamics in its S1 domain compared to wild-type (Wang et al., 2022a); this more ‘open’ state likely allows greater accessibility to the epitopes of these nanobody groups, for net stronger binding. Additionally, the absence of delta-specific mutations within the epitope regions of groups I and II preserves their integrity for nanobody binding (Figures 1B and 3B). For the remaining RBD nanobody groups, we are likely observing the impact of the two delta mutations T478K and L452R on nanobody binding and neutralization. These two mutations lie within the mapped epitope regions of group III, IV, and V nanobodies (Figures 1B and 3B), which may alter the epitopes of many of these nanobodies enough to negatively impact both binding and neutralization. Importantly, the L452R mutation seems to be a key substitution that contributes to reduced or abolished neutralizing abilities of many mAbs (Laurini et al., 2021; Starr et al., 2021a, Starr et al., 2021b). Antibodies that rely on L452 to create hydrophobic interactions within their epitope will most likely have their binding greatly disrupted with a substituted arginine. Coupled with the stronger affinity between ACE2 and spike caused by L452R (Motozono et al., 2021; Yan et al., 2022), this substitution can greatly lessen the neutralization ability of nanobodies and antibodies targeting this region.

Unlike wild-type, the spike trimer of omicron BA.1 favors a one-RBD up confirmation (Zhao et al., 2022). Though this conformation may facilitate access of our nanobodies to their epitopes on RBD, unlike delta, omicron BA.1 contains many unique mutations distributed throughout the RBD domain that overlap or flank the epitopes of almost all our nanobody groups (Figures 1B and 3B). These mutations likely alter the epitopes of our nanobodies and may be the major contributing factor to the observed decrease in binding affinity and neutralizing potency of our nanobodies against omicron BA.1 (Figures 1A and 3A). However, the following cohort of nanobodies retained binding and neutralizing ability similar to wild-type against omicron BA.1: S1-RBD-22, S1-RBD-9 (both group I, class 1), S1-4 (group I/II, class 3), S1-RBD-5 (group I/IV, class 3), S1-46 (group VII), and members of group X. The additional mutations on the spike of omicron BA.4/BA.5, many of which overlap with our nanobody epitopes, are predicted to further impact the binding of our nanobody repertoire to this variant. For example, omicron BA.4/BA.5 re-introduces the key L452R substitution from delta, which when combined with the accumulated mutations on spike from omicron precursor sublineages, could be responsible for the observed loss of binding of numerous nanobodies in groups III, IV, and V. Furthermore, the differential binding and neutralizing abilities of group I and II nanobodies against omicron variants BA.4/BA.5/XBB/BQ.1.1 may be because of the T376A, D405N, and R408S substitutions, which lie within and near the epitope regions of groups I and II, respectively (Figure 3B). Importantly, at least five nanobodies from groups I, I/II, II, and V retained neutralization activity against omicron BA.5 (Figure 4), demonstrating the broad specificity of this set of nanobodies.

Nanobodies that lose neutralization ability can still bind spike

Our SPR experiments largely correlated with the neutralization data: nanobodies that showed binding to spike also neutralized the virus, and where the binding affinity decreased significantly, a loss of neutralization ability was observed. This was seen with the group III nanobodies S1-23 and S1-37 that both demonstrated significantly reduced binding affinity and neutralization potency against both delta and omicron BA.1. This reduction in efficacy is likely attributable to the L452R mutation against delta, and the extensive amino acid changes in omicron BA.1 compared to wild-type (see above). However, we observed instances where nanobodies retained binding to spike yet no longer neutralized the virus. Nanobodies S1-36, S1-39, and S1-RBD-29 showed binding to omicron BA.1 (Figure 3A), yet none neutralized the variant in the pseudovirus assay (Figure 1A). These nanobodies are in groups III, II, and IV, respectively (Figure 3B) – epitope regions that contain numerous omicron BA. 1 mutations. For all three nanobodies, the mutations peppered throughout their epitope space have likely altered the binding landscape, resulting in decreased affinity (~100-fold for S1-36 and S1-39) and ineffective neutralization in vitro. Additionally, S1-RBD-29 binds an epitope that overlaps with the ACE2-binding site, and likely neutralizes the wild-type strain by blocking the ACE2 interaction (Figure 1B). Alterations in the omicron BA.1 epitope may have changed the orientation of the nanobody as it binds, negating effective blocking of ACE2 binding and thus neutralization. In contrast, S1-37 binds both delta and omicron BA.1 with similar reduced affinities (>100-fold decrease compared to wild-type), yet only neutralized omicron BA.1 (Figure 1A). As mentioned above, the epitope of S1-37 overlaps with the delta L452R mutation (Figure 1B), which has impacted the effectiveness of numerous neutralizing antibodies (Bian et al., 2021). It is possible that this mutation alone drastically weakens S1-37 binding to delta and consequently virus neutralization. Lastly, S1-RBD-43, whose epitope as a group IV binder is predicted to overlap with the ACE2 binding footprint (Figures 1B and 3B), showed binding to delta with equal affinity to wild-type (Figure 3A), yet does not neutralize delta (Figure 1A). The delta mutation T478K is present within the epitope region of group IV (Figures 1B and 3B), and was shown to significantly increase the interaction of delta spike for ACE2 by creating a new salt bridge at the RBD/ACE2 interface (Cheng et al., 2022). It may be that, as suggested for S1-RBD-29, the binding orientation of S1-RBD-43 has been altered allowing the RBD to maintain an interaction with ACE2 despite the presence of S1-RBD-43, thus rendering the nanobody ineffective in neutralizing delta.

Identification of variant-specific epitopes and broadly neutralizing epitope groups

The results of our neutralization assays and affinity measurements revealed that ~1/3 of our original repertoire of 116 nanobodies (Mast et al., 2021) generated against wild-type SARS-CoV-2 remain effective binders/neutralizers of the variants tested. Specifically, nanobodies from 11 of the 18 nanobody groups (inclusive of the 10 mapped epitopes) demonstrated efficacy against one or more of the delta, omicron BA.1, omicron BA.4/BA.5, omicron XBB, and omicron BQ.1.1 lineages. The varied efficacy of nanobodies within each group, along with structural modeling, enabled us to expand and further refine our original six structurally modeled RBD epitope groups to a total of 12 (Figure 5A). These 12 groups were then contrasted with the three mAb classes containing mAbs previously approved for EUA by the FDA (Figure 2). This comparison revealed the following overlaps between our 12 RBD nanobody epitope groups and the three mAb classes: groups I, III, IV, and VII overlap with class 1 binders; groups II, III, IV, and V overlap with class 2 binders; and groups II, III, and V overlap with class 3 binders (Figure 2).

Figure 5 with 1 supplement see all

Download asset Open asset

Refining epitopes of the nanobody repertoire.

All epitopes are mapped on the structure of wild-type receptor-binding domain (RBD) (PDB ID: 6M0J). (A) The original six RBD nanobody epitope groups capable of binding and/or neutralizing one or more variants are highlighted in dark gray. Further refinement of four groups: I, III, IV, and V led to the identification of six additional epitope groups – resulting in a total of 12 epitope groups able to bind one or more variants of concern. (B) Summary of variant specific and broadly reactive nanobodies. (C) Nanobody groups predicted to bind/neutralize the circulating omicron variants EG.5 and HV.1. All structure representations were created on ChimeraX (Pettersen et al., 2021).

Interestingly, the epitope footprint of many of our nanobody classes extends beyond those of the three mAb classes (Figure 2). For example, the majority of the epitope regions of the group I, II, V, and VII nanobodies do not overlap with the mAb classes and are not a binding/neutralizing hotspot for mAbs (Almagro et al., 2022; Cao et al., 2022); instead, these epitopes extend away from the ACE2-binding site (Figures 1B—3B), as seen in particular with groups I, II, and V nanobodies. These regions may be inaccessible to mAbs, possibly due to steric limitations, a property nanobodies readily overcome due to their small size. Additionally, our data allowed us to identify variant-specific epitope groups (Figure 5B), where we define variant-specific nanobodies as nanobodies that bind a single additional variant alongside the original Wuhan strain, summarized as follows: of the 26 nanobodies that showed binding/neutralization to delta, 5 were specific only for delta; of the 21 nanobodies that showed binding/neutralization to omicron BA.1, 2 were specific for omicron BA.1. Strikingly, we have in our cohort eight nanobodies able to bind delta and the omicron lineages BA.1/BA.4/BA.5/XBB/BQ.1.1 (Figure 5B). We further predict these eight nanobodies will be effective binders against current circulating strains of the virus including omicron EG.5 and HV.1 as the epitope regions (or predicted epitopes) of these nanobodies do not vary significantly from omicron lineages XBB and BQ.1.1 (Figure 5C; Figure 5—figure supplement 1).

Nanobody synergy involving a non-neutralizing nanobody

We previously established that cocktails of nanobodies exhibit enhanced resistance to mutational escape (Mast et al., 2021). Excitingly, not only was the barrier to mutational escape extremely enhanced, but for certain combinations of nanobodies, their mechanisms of neutralization were synergistic, providing far more potent neutralization in combination than expected from the neutralization by either nanobody alone (Mast et al., 2021). Our present observations of nanobodies that retained binding to variants of the spike RBDdespite losing neutralization efficacy (Figure 1A) afforded us an opportunity to test whether the synergy observed for certain nanobody combinations was dependent on their ability to neutralize.

The synergistic S1-1 and S1-23 pair effectively neutralized the wild-type PSV with their epitopes on opposing surfaces of the RBD, permitting simultaneous binding and enhanced neutralization when delivered as a cocktail (Mast et al., 2021; Figures 1A, 3A,, 6A, panel i). While S1-1 remained efficacious against the delta variant of SARS-CoV-2, the L452R mutation in the delta RBD likely negatively impacted S1-23 (discussed above), weakening its binding affinity by ~1000-fold (12 nM) and negating its neutralization efficacy at concentrations <10 μM (Figures 1A, 3A,, 6A, panel ii). Surprisingly, when provided in combination with S1-1, which displays increased binding affinity and enhanced neutralization against the delta variant, S1-23 was able to further enhance the neutralization capabilities of S1-1, synergistically, at concentrations above 10⁻³ μM (Figure 6A, panel ii), by up to 42-fold. This synergistic interaction, however, did not apply to situations where extensive mutations are present in the RBD, such as in the omicron sublineages of SARS-CoV-2, which ablated the binding and neutralization efficacy of both S1-1 and S1-23 (Figures 1A, 3A,, 6A, panel iii).

Figure 6

Download asset Open asset

Persistence of synergistic neutralization with nanobody cocktails against SARS-CoV-2 variants of concern.

(A) S1-1 synergizes with S1-23 in neutralizing SARS-CoV-2 PSV. The upper left panel shows two representations of spike with the accessible S1-1 (dark goldenrod) and S1-23 (sky blue) epitopes (PDB ID: 6VYB). The measured affinities for S1-1 and S1-23for the receptor-binding domain (RBD) of wild-type, delta, and omicron BA.1 are displayed. Both S1-1 and S1-23 neutralize wild-type (i), whereas only S1-1 neutralizes delta at the concentrations shown (ii). In spite of a lack of neutralization at these concentrations, S1-23, synergizes with S1-1 and enhances its neutralization of delta SARS-CoV-2 PSV (ii). As neither S1-1 nor S1-23 are able to bind to the RBD of omicron BA.1, neither nanobody neutralizes omicron BA.1 SARS-CoV-2 PSV (**iii**). (B) S1-36 synergizes with S1-RBD-22 in neutralizing SARS-CoV-2 PSV. As in A, the upper left panel shows two representations of spike with the accessible S1-36 (cornflower blue) and S1-RBD-22 (sandy brown) epitopes. The measured affinities for S1-36 and S1-RBD-22 are displayed. Both S1-36 and S1-RBD-22 neutralize wild-type (i), whereas only S1-RBD-22 effectively neutralizes delta and omicron BA.1 SARS-CoV-2 PSV at the concentrations shown (ii and **iii**, respectively). However, S1-36 synergizes with S1-RBD-22 and enhances its neutralization of the three depicted SARS-CoV-2 pseudoviruses (i, ii, and **iii**). (C) An example of no interactions (synergistic or antagonistic) between S1-1 and LaM2 (Fridy et al., 2014), a non-specific nanobody that does not bind the RBD of delta. (D) An example of antagonism, where higher concentrations of S1-23 interferes with the ability of S1-RBD-16 to neutralize delta SARS-CoV-2 PSV. These nanobodies have adjacent epitopes on the RBD of spike and were previously shown to interfere with each other’s binding to their respective epitope (Mast et al., 2021). n = 4. Source data in .

Figure 6—source data 1 Neutralization data from synergy experiment.: https://cdn.elifesciences.org/articles/89423/elife-89423-fig6-data1-v1.xlsx
Download elife-89423-fig6-data1-v1.xlsx

We also tested the broadly neutralizing nanobody S1-RBD-22 in combination with S1-36 (Figure 6B). Like S1-23, the epitope of S1-36 is opposite that of S1-RBD-22, permitting simultaneous binding to a single RBD (Mast et al., 2021). However, while its neutralization efficacy dropped off when delivered to either delta or omicron BA.1 PSV alone, its ability to bind to its epitope was only marginally impacted (Figures 1A, 3A,, 6B). When provided in combination with S1-RBD-22, S1-36 synergistically enhanced by up to 80-fold the neutralization efficacy of S1-RBD-22 against both delta and omicron BA.1 PSV (Figure 6B, panels ii, iii). The synergies observed between S1-1 and S1-23, and between S1-RBD-22 and S1-36, appear to be specific rather than due to pleiotropic effects. This is evidenced by the lack of enhanced neutralization when S1-1 is combined with the non-specific LaM2 nanobody against delta PSV (Figure 6C). Furthermore, in the case of S1-23, binding to its non-neutralizing epitope on the RBD of delta PSV was able to induce dose-dependent antagonistic effects on the neutralizing efficacy of S1-RBD-16, which binds to a neighboring epitope that can be competitively blocked by S1-23 (Figure 6D; Mast et al., 2021).

Conclusions

Collectively, our binding and neutralization data allowed us to identify the regions of spike in multiple VoC that remained vulnerable to our original repertoire of nanobodies raised against wild-type SARS-CoV-2 spike (Figure 7; Mast et al., 2021). Unlike the epitope groups defining mAb-binding sites on spike, many of our nanobody epitope groups remained efficacious in the neutralization of different VoC (Figures 6 and 7), possibly due to each nanobody’s smaller epitope footprint allowing their access to regions of spike inaccessible to mAbs. However, the substantial changes on the surface of spike that has occurred as SARS-CoV-2 has evolved from one variant to the next has also negatively impacted many nanobodies by abolishing or weakening their binding and/or neutralization activity. This weakening is most evident for nanobodies directed against the NTD, and against the receptor-binding motif that engages ACE2 (Figures 1B and 2B).

Figure 7

Download asset Open asset

Nanobodies effective against circulation variants of concern.

The nanobody epitopes (red) that retain effectiveness against wild-type (PDB ID: 7KNB), delta (PDB ID: 7V7O), omicron BA.1, omicron BA.4 (the nanobody epitopes were mapped on to PDB ID: 7XO5 for both BA.1 and BA.4 due to the lack of a suitable BA.4 structure),omicron XBB.1 and BQ.1.1 (the nanobody epitopes were mapped on to PDB ID: 8IOU for both XBB.1 and BQ.1.1 due to the lack of a suitable BQ.1.1 structure). The spike trimer (*silver*) where glycans are represented as *light blue* sticks for each variant is displayed in three views: top = side view of spike trimer; middle = birds eye view looking down on the S1 domain and; bottom = birds eye view (same as middle view), with the S1 domain rendered transparent to enable visualization of the S2 domain. All structure representations were created on ChimeraX (Pettersen et al., 2021).

Intriguingly, we discovered nanobodies that maintained binding capabilities while losing their neutralizing properties. This finding paves the way for engineering these nanobodies through approaches like oligomerization, which has proven effective in boosting neutralization (Wrapp et al., 2020; Mast et al., 2021). Remarkably, we also showed that such non-neutralizing binders can nevertheless retain effectiveness as components of synergistic nanobody cocktails, suggesting their potential for widespread antiviral applications. This finding indicates that the mechanism of synergy can operate through epistatic interactions from binding alone, not solely through direct neutralization. Furthermore, these synergistic pairings hold promise as therapeutic options when formulated as multivalent hetero-oligomers and we propose that any of the nanobodies we have demonstrated to show pan-VoC activity would be prime candidates for further optimization.

In this study, we identified nanobodies that specifically recognize only certain SARS-CoV-2 variants (Figure 5B), allowing for the possibility of distinguishing between different VoC. By utilizing these nanobodies as molecular probes in diagnostic tests, a unique ‘molecular fingerprint’ could define each variant based on the combinations of nanobodies that bind to and recognize the specific virus particle. Consequently, this approach could enable the accurate and rapid diagnosis of SARS-CoV-2 infections as well as provide real-time identification of the specific variant causing the infection, thus enhancing surveillance and tailoring treatment strategies accordingly to optimize patient outcomes and control the spread of the virus.

Together, our findings highlight the strength and variety of the heavy chain-only category of llama antibody immune responses and, as a result, the extensive repertoires of high-affinity nanobodies. This result emphasizes the distinct benefits of nanobody technology, which enables extensive coverage of antigenic regions while simultaneously targeting multiple unique epitope sites. This approach paves the way for investigating a wide range of therapeutic options against rapidly evolving proteins like the SARS-CoV-2 spike protein, ultimately aiding in our preparedness and defense against future pandemics or major outbreaks.

Materials and methods

K_D measurements

Request a detailed protocol

Series S CM5 sensor chips (Cytiva) were immobilized with recombinant delta Spike S1, delta RBD, and omicron RBD at 12.5, 10, and 15 µg/ml, respectively, using EDC/NHS coupling chemistry according to the manufacturer’s instructions. K_D measurements were performed on a Biacore 8k (Cytiva) as previously described (Mast et al., 2021), using 5–8 concentrations of each nanobody. Data were processed and analyzed using the Biacore Insight Evaluation software.

Structural modeling

Request a detailed protocol

Integrative structure modeling of nanobody epitopes on the different variant spike proteins, proceeded through the standard four-stage protocol (Kim et al., 2018; Webb et al., 2018; Rout and Sali, 2019; Sali, 2021; Saltzberg et al., 2021). This protocol was implemented using the Python Modeling Interface package, a library for modeling macromolecular complexes based on the open-source Integrative Modeling Platform software, version 2.15.0 (https://integrativemodeling.org). Only (a subset of) nanobodies with pre-determined experimental escape mutations on the Wuhan spike structure were selected for modeling. Separate models were computed for rigid-receptor–rigid ligand-type binary docking of representative nanobodies from group 1 (S1-1, S1-RBD-[9, 15, 22, 24]), group 2 (S1-6), group 3 (S1-[23, 36]), group 4 (S1-RBD-[21, 29]), group 5 (S1-RBD-[16, 23]), group 7 (S1-46), group 8 (S1-49), and group 9 (S2-10), to the variant spike structures. S1-49 was docked to a monomeric S1-NTD domain, S2-10 was docked to the trimeric S2 (ecto-) domain, while the remaining nanobodies were docked to the S1-RBD domain. Monomeric S1-RBD spanning amino acid residues 333–526 was represented using the 2.45 Å crystal structure of the ACE2-bound RBD (6M0J.E; Lan et al., 2020) for the original virus, the 4.30 Å cryo-EM structure of up-RBD pre-fusion spike (7SBO.A; Zhang et al., 2021) for the delta variant, the ACE2 bound 2.45 Å cryo-EM spike structure (7T9K.B; Mannar et al., 2022) for the omicron BA.1 variant and a structure predicted using AlphaFold-2 for the omicron BA.4/BA.5 variant, the 3.18 Å cryo-EM structure of ACE2-bound up-RBD pre-fusion spike (8IOU:A; Tamura et al., 2023) for the XBB.1 variant, and the 3.20 Å cryo-EM structure of the BQ.1.1 variant RBD (8FXC:E; Tamura et al., 2023) in complex with ACE2 and the S309-neutralizing antibody Fab fragment. Missing residues in 8FXC:E were filled in using Modeller (Webb and Sali, 2016).

For the original virus, monomeric S1-NTD, spanning amino acid residues 16–305, was represented using the crystal structure of the S2M28 Fab bound NTD (7LY3.A; McCallum et al., 2021), while trimeric S2 was represented using the amino residues 689–1162 (for each monomer) from the 2.9 Å cryo-EM structure 6XR8 (Cai et al., 2020). NTD and S2 structures for delta and omicron variants were extracted from the corresponding whole spike structures.

Structural models for all 15 nanobodies and the omicron BA.4 RBD were built with the ColabFold implementation of AlphaFold2 (Jumper et al., 2021; Mirdita et al., 2022). The protocol included automatic refinement of the complementarity-determining region (CDR) loops through an all-atom energy minimization of the AlphaFold2-predicted structure using the AMBER molecular mechanics force field (Hornak et al., 2006). We verified that the predicted nanobody structures are within 3–4 Å backbone root-mean-square deviation (RMSD) from the comparative models of these nanobody sequences published previously (Mast et al., 2021). The CDR region boundaries in the nanobody structures were assigned using the FREAD algorithm as implemented within the SabPred web server (Rausch, 1991).

To make structural sampling sufficiently efficient, the system was represented at a resolution of one bead per residue, and the receptors and all nanobodies were treated as rigid bodies. For each nanobody, alternate binding modes were scored using spatial restraints enforcing receptor–ligand shape complementarity, cross-link satisfaction and proximity of CDR3 loops on the nanobodies to escape mutant residues on the corresponding receptor. With the receptor fixed in space, 1,200,000 alternative docked nanobody models were produced through 20 independent runs of replica exchange Gibbs sampling based on the Metropolis Monte Carlo algorithm, where each Monte Carlo step consisted of a series of random rotations and translations of rigid nanobodies. The initial set of models was filtered to obtain a random subsample of 30,000 models, which were clustered by the structural similarity of their interfaces to the receptor; this similarity was quantified by the fraction of common contacts (fcc) between receptor and nanobody was used to characterize interface similarity between alternate nanobody poses (Rodrigues et al., 2012). Binding poses from the most populated cluster were selected for further analysis. Five independent random subsamples of 30,000 models each were generated from the set of all models post-structural sampling, and the entire protocol of interface similarity-based clustering and top cluster selection was repeated each time. No significant differences among these five subsamples were observed in the satisfaction of restraints. Structural differences among the variants, as well as between AlphaFold2 models and previously published comparative models of nanobodies, lead to differences in binding modes of the same nanobody to different spike variants. Thus, for the sake of consistency we limit our comparison to the receptor epitopes, which are defined as all receptor atoms that are within 6 Å of the framework and CDR regions of the nanobodies (excluding the flexible N- and C-terminal regions). Although we do not include the nanobody paratopes in our analysis, we verified that all binding modes are primarily through CDR3, except for the CDR1 contribution to the binding of S1-1 and S1-RBD-15 to the RBD, for all variants. Relative differences in binding affinity (K_D) and neutralization potential (IC₅₀) between the original virus and other variants (delta and omicron) were projected onto the Wuhan epitopes to create the heatmaps in Figures 1B and 2B. Relative differences were reported as and/or normalized to the range from 0 to 100.

Integrative models of nanobody epitopes on the spike protein were computed on the Wynton HPC cluster at UCSF. Receptor epitopes were visualized in UCSF ChimeraX (Pettersen et al., 2021). Files containing input data, scripts, and output results are available at https://github.com/integrativemodeling/nbspike/tree/main/integrative_modeling_VOC (copy archived at Sanyal, 2024). Structure predictions using ColabFold utilized the ‘AlphaFold2_batch’ notebook, with the default settings. All modeled structures were subjected to molecular-mechanics-based relaxation, followed by using the model with the top pLDDT score was selected for the integrative modeling pipeline.

Cell lines

Request a detailed protocol

TMPRSS2-expressing Vero E6 cells, 293T/17 cells and 293T-hACE2 cells were cultured as described previously (Mast et al., 2021). Briefly, TMPRSS2 + Vero E6 cells were cultured at 37°C in the presence of 5% CO₂ in medium composed of in high-glucose Dulbecco’s modified Eagle’s medium (DMEM, Gibco) supplemented with 10% (vol/vol) fetal bovine serum (FBS) and 1 mg/ml geneticin. 293T/17 were cultured at 37°C in the presence of 5% CO₂ in a medium composed of DMEM supplemented with 10% (vol/vol) FBS and penicillin/streptomycin. 293T-hACE2 cells were cultured at 37°C in the presence of 5% CO₂ in medium composed of DMEM supplemented with 10% (vol/vol) FBS, penicillin/streptomycin, 10 mM HEPES (4-(2-hydroxyethyl)-1-piperazineethanesulfonic acid), and with 0.1 mM modified Eagle’s medium (MEM) non-essential amino acids (Thermo Fisher). All experiments were performed with cells passaged less than 15 times. The identities of cell lines were confirmed by chromosomal marker analysis and tested negative for mycoplasma using a MycoStrip (InvivoGen).

Production of SARS-CoV-2 variant pseudotyped lentiviral reporter particles

Request a detailed protocol

Pseudovirus stocks were prepared and ittered as described previously (Mast et al., 2021). Variant spike containing plasmids were combined with pHAGE-CMV-Luc2-IRES-ZsGreen-W (BEI Cat # NR-52516) (Crawford et al., 2020), and psPAX using lipofectamine 3000 and cotransfected into 293T/17 cells. Pseudovirus was titered by threefold serial dilution on 293T-hACE2 cells, as described previously (Mast et al., 2021).

SARS-CoV-2 pseudovirus neutralization assay

Request a detailed protocol

Nanobodies were tested for their neutralization properties as described previously (Mast et al., 2021). Briefly, threefold serial dilutions of nanobodies were incubated with pseudotyped SARS-CoV-2 for 1 hr at 37°C. The nanobody–pseudovirus mixtures were then added in quadruplicate to 293T-hACE2 cells along with 2 µg/ml polybrene (Sigma). Cells were incubated at 37°C with 5% CO₂. Infected cells were processed between 52 and 60 hr by adding equal volume of Steady-Glo (Promega), and firefly luciferase signal was measured using the Biotek Model N4 with integration at 0.5 ms. Data were processed using Prism 7 (GraphPad), using four-parameter non-linear regression (least-squares regression method without weighting). All nanobodies were tested at least two times and with more than one pseudovirus preparation.

Plaque reduction neutralization assay with SARS-CoV-2 BA.5

Request a detailed protocol

Briefly, 10 threefold serial dilutions of representative nanobodies of select epitope groups in Opti-MEM (Gibco) were incubated with approximately 100–200 pfus of SARS-CoV-2 BA.5 for 1 hr at RT. The nanobody/virus mixture was then added to a confluent monolayer of TMPRSS2 + Vero E6 cells in 12-well plates and incubated at RT for 90 min. One well was overlaid with virus only while another well was uninfected. Virus/nanobody mixture was removed, and the cell monolayer overlaid with a medium composed of 3% (wt/vol) carboxymethycellulose and 4% (vol/vol) FBS in Opti-MEM. 96 hr post infection, the overlay was removed and cell monolayer was washed with Dulbecco’s phosphate-buffered saline (DPBS, Gibco) before being fixed with 4% (wt/vol) paraformaldehyde in DPBS for 30 min. Fixative was removed and the cells were rinsed with DPBS before being stained with 1% (wt/vol) crystal violet in 20% (vol/vol) ethanol. Contrast was enhanced by washing with DPBS, and clear plaques representing individual viral infections were visualized as spots lacking crystal violet stain. Plaques were quantified and the ratio of plaques at each nanobody dilution to ‘virus only’ well was used to determine the IC₅₀s of each nanobody.

Nanobody synergy

Request a detailed protocol

Synergy experiments were performed as described previously (Mast et al., 2021). Briefly, a robotic liquid handler was used to prepare 2D matrices of threefold serial dilutions of two nanobodies and then mix these combinations with different variant pseudotyped SARS-CoV-2 for 1 hr. After incubation with the virus, the mixture was overlaid on a monolayer of 293-hACE2 cells and left to incubate for 56 hr. Luminescence was quantified as described above. Data were processed using the Bivariate Response to Additive Interacting Doses (BRAID) model (Twarog et al., 2016) as implemented in the synergy software package for python (Wooten and Albert, 2021).

SARS-CoV-2 stocks and titers

Request a detailed protocol

All experimental work involving live SARS-CoV-2 was performed at Seattle Children’s Research Institute (SCRI) in compliance with SCRI guidelines for BioSafety Level 3 (BSL-3) containment. SARS-CoV-2 isolate CGIDR_SARS2 omicron BA.5 was obtained from an infected individual. An initial inoculum was diluted in Opti-MEM (Gibco) at 1:1000, overlaid on a monolayer of Vero E6 and incubated for 90 min. Following the incubation, the supernatant was removed and replaced with 2% (vol/vol) FBS in Opti-MEM medium. The cultures were inspected for cytopathic effects, and infectious supernatants were collected, cleared of cellular debris by centrifugation, and stored at −80°C until use. Whole viral genome sequencing and variant analysis were performed by the University of Washington Department of Laboratory Medicine & Pathology. Viral titers were determined by plaque assay using a liquid overlay and fixation-staining method, as described previously (Mendoza et al., 2020; Mast et al., 2021).

Data availability

All data generated or analyzed during this study are included in the manuscript and supporting files.

References

(2022) Evolution of anti-SARS-CoV-2 therapeutic antibodies
International Journal of Molecular Sciences 23:17.
https://doi.org/10.3390/ijms23179763
- PubMed
- Google Scholar
1. Bian L
2. Gao Q
3. Gao F
4. Wang Q
5. He Q
6. Wu X
7. Mao Q
8. Xu M
9. Liang Z
(2021) Impact of the Delta variant on vaccine efficacy and response strategies
Expert Review of Vaccines 20:1201–1209.
https://doi.org/10.1080/14760584.2021.1976153
- PubMed
- Google Scholar
1. Cai Y
2. Zhang J
3. Xiao T
4. Peng H
5. Sterling SM
6. Walsh RM Jr
7. Rawson S
8. Rits-Volloch S
9. Chen B
(2020) Distinct conformational states of SARS-CoV-2 spike protein
Science 369:1586–1592.
https://doi.org/10.1126/science.abd4251
- PubMed
- Google Scholar
1. Cao Y
2. Yisimayi A
3. Jian F
4. Song W
5. Xiao T
6. Wang L
7. Du S
8. Wang J
9. Li Q
10. Chen X
11. Yu Y
12. Wang P
13. Zhang Z
14. Liu P
15. An R
16. Hao X
17. Wang Y
18. Wang J
19. Feng R
20. Sun H
21. Zhao L
22. Zhang W
23. Zhao D
24. Zheng J
25. Yu L
26. Li C
27. Zhang N
28. Wang R
29. Niu X
30. Yang S
31. Song X
32. Chai Y
33. Hu Y
34. Shi Y
35. Zheng L
36. Li Z
37. Gu Q
38. Shao F
39. Huang W
40. Jin R
41. Shen Z
42. Wang Y
43. Wang X
44. Xiao J
45. Xie XS
(2022) BA.2.12.1, BA.4 and BA.5 escape antibodies elicited by Omicron infection
Nature 608:593–602.
https://doi.org/10.1038/s41586-022-04980-y
- PubMed
- Google Scholar
1. Cheng MH
2. Krieger JM
3. Banerjee A
4. Xiang Y
5. Kaynak B
6. Shi Y
7. Arditi M
8. Bahar I
(2022) Impact of new variants on SARS-CoV-2 infectivity and neutralization: A molecular assessment of the alterations in the spike-host protein interactions
iScience 25:103939.
https://doi.org/10.1016/j.isci.2022.103939
- PubMed
- Google Scholar
1. COVID-19 Cumulative Infection Collaborators
(2022) Estimating global, regional, and national daily and cumulative infections with SARS-CoV-2 through Nov 14, 2021: a statistical analysis
Lancet (London, England) 399:2351–2380.
https://doi.org/10.1016/S0140-6736(22)00484-6
- PubMed
- Google Scholar
(2023) SARS-CoV-2 variant evasion of monoclonal antibodies based on in vitro studies
Nature Reviews. Microbiology 21:112–124.
https://doi.org/10.1038/s41579-022-00809-7
- PubMed
- Google Scholar
1. Crawford KHD
2. Eguia R
3. Dingens AS
4. Loes AN
5. Malone KD
6. Wolf CR
7. Chu HY
8. Tortorici MA
9. Veesler D
10. Murphy M
11. Pettie D
12. King NP
13. Balazs AB
14. Bloom JD
(2020) Protocol and reagents for Pseudotyping lentiviral particles with SARS-CoV-2 spike protein for neutralization assays
Viruses 12:513.
https://doi.org/10.3390/v12050513
- PubMed
- Google Scholar
1. Cross FR
2. Fridy PC
3. Ketaren NE
4. Mast FD
5. Li S
6. Olivier JP
7. Pecani K
8. Chait BT
9. Aitchison JD
10. Rout MP
(2023) Expanding and improving nanobody repertoires using a yeast display method: Targeting SARS-CoV-2
The Journal of Biological Chemistry 299:102954.
https://doi.org/10.1016/j.jbc.2023.102954
- PubMed
- Google Scholar
1. Dadonaite B
2. Crawford KHD
3. Radford CE
4. Farrell AG
5. Yu TC
6. Hannon WW
7. Zhou P
8. Andrabi R
9. Burton DR
10. Liu L
11. Ho DD
12. Chu HY
13. Neher RA
14. Bloom JD
(2023) A pseudovirus system enables deep mutational scanning of the full SARS-CoV-2 spike
Cell 186:1263–1278.
https://doi.org/10.1016/j.cell.2023.02.001
- PubMed
- Google Scholar
(2023) Convergent evolution in SARS-CoV-2 spike creates a variant soup from which new COVID-19 waves emerge
International Journal of Molecular Sciences 24:2264.
https://doi.org/10.3390/ijms24032264
- PubMed
- Google Scholar
1. Fridy PC
2. Li Y
3. Keegan S
4. Thompson MK
5. Nudelman I
6. Scheid JF
7. Oeffinger M
8. Nussenzweig MC
9. Fenyö D
10. Chait BT
11. Rout MP
(2014) A robust pipeline for rapid production of versatile nanobody repertoires
Nature Methods 11:1253–1260.
https://doi.org/10.1038/nmeth.3170
- PubMed
- Google Scholar
1. Grant OC
2. Montgomery D
3. Ito K
4. Woods RJ
(2020) Analysis of the SARS-CoV-2 spike protein glycan shield reveals implications for immune recognition
Scientific Reports 10:14991.
https://doi.org/10.1038/s41598-020-71748-7
- PubMed
- Google Scholar
1. Greaney AJ
2. Loes AN
3. Crawford KHD
4. Starr TN
5. Malone KD
6. Chu HY
7. Bloom JD
(2021a) Comprehensive mapping of mutations in the SARS-CoV-2 receptor-binding domain that affect recognition by polyclonal human plasma antibodies
Cell Host & Microbe 29:463–476.
https://doi.org/10.1016/j.chom.2021.02.003
- PubMed
- Google Scholar
1. Greaney AJ
2. Starr TN
3. Barnes CO
4. Weisblum Y
5. Schmidt F
6. Caskey M
7. Gaebler C
8. Cho A
9. Agudelo M
10. Finkin S
11. Wang Z
12. Poston D
13. Muecksch F
14. Hatziioannou T
15. Bieniasz PD
16. Robbiani DF
17. Nussenzweig MC
18. Bjorkman PJ
19. Bloom JD
(2021b) Mapping mutations to the SARS-CoV-2 RBD that escape binding by different classes of antibodies
Nature Communications 12:4196.
https://doi.org/10.1038/s41467-021-24435-8
- PubMed
- Google Scholar
1. Hachmann NP
2. Miller J
3. Collier A-RY
4. Ventura JD
5. Yu J
6. Rowe M
7. Bondzie EA
8. Powers O
9. Surve N
10. Hall K
11. Barouch DH
(2022) Neutralization escape by SARS-CoV-2 Omicron subvariants BA.2.12.1, BA.4, and BA.5
The New England Journal of Medicine 387:86–88.
https://doi.org/10.1056/NEJMc2206576
- PubMed
- Google Scholar
1. Hornak V
2. Abel R
3. Okur A
4. Strockbine B
5. Roitberg A
6. Simmerling C
(2006) Comparison of multiple Amber force fields and development of improved protein backbone parameters
Proteins 65:712–725.
https://doi.org/10.1002/prot.21123
- PubMed
- Google Scholar
1. Imai M
2. Ito M
3. Kiso M
4. Yamayoshi S
5. Uraki R
6. Fukushi S
7. Watanabe S
8. Suzuki T
9. Maeda K
10. Sakai-Tagawa Y
11. Iwatsuki-Horimoto K
12. Halfmann PJ
13. Kawaoka Y
(2023) Efficacy of Antiviral Agents against Omicron Subvariants BQ.1.1 and XBB
The New England Journal of Medicine 388:89–91.
https://doi.org/10.1056/NEJMc2214302
- PubMed
- Google Scholar
1. Jackson LA
2. Anderson EJ
3. Rouphael NG
4. Roberts PC
5. Makhene M
6. Coler RN
7. McCullough MP
8. Chappell JD
9. Denison MR
10. Stevens LJ
11. Pruijssers AJ
12. McDermott A
13. Flach B
14. Doria-Rose NA
15. Corbett KS
16. Morabito KM
17. O’Dell S
18. Schmidt SD
19. Swanson PA II
20. Padilla M
21. Mascola JR
22. Neuzil KM
23. Bennett H
24. Sun W
25. Peters E
26. Makowski M
27. Albert J
28. Cross K
29. Buchanan W
30. Pikaart-Tautges R
31. Ledgerwood JE
32. Graham BS
33. Beigel JH
(2020) An mRNA Vaccine against SARS-CoV-2 — Preliminary Report
New England Journal of Medicine 383:1920–1931.
https://doi.org/10.1056/NEJMoa2022483
- Google Scholar
1. Jackson CB
2. Farzan M
3. Chen B
4. Choe H
(2022) Mechanisms of SARS-CoV-2 entry into cells
Nature Reviews. Molecular Cell Biology 23:3–20.
https://doi.org/10.1038/s41580-021-00418-x
- PubMed
- Google Scholar
1. Jumper J
2. Evans R
3. Pritzel A
4. Green T
5. Figurnov M
6. Ronneberger O
7. Tunyasuvunakool K
8. Bates R
9. Žídek A
10. Potapenko A
11. Bridgland A
12. Meyer C
13. Kohl SAA
14. Ballard AJ
15. Cowie A
16. Romera-Paredes B
17. Nikolov S
18. Jain R
19. Adler J
20. Back T
21. Petersen S
22. Reiman D
23. Clancy E
24. Zielinski M
25. Steinegger M
26. Pacholska M
27. Berghammer T
28. Bodenstein S
29. Silver D
30. Vinyals O
31. Senior AW
32. Kavukcuoglu K
33. Kohli P
34. Hassabis D
(2021) Highly accurate protein structure prediction with AlphaFold
Nature 596:583–589.
https://doi.org/10.1038/s41586-021-03819-2
- PubMed
- Google Scholar
1. Kim SJ
2. Fernandez-Martinez J
3. Nudelman I
4. Shi Y
5. Zhang W
6. Raveh B
7. Herricks T
8. Slaughter BD
9. Hogan JA
10. Upla P
11. Chemmama IE
12. Pellarin R
13. Echeverria I
14. Shivaraju M
15. Chaudhury AS
16. Wang J
17. Williams R
18. Unruh JR
19. Greenberg CH
20. Jacobs EY
21. Yu Z
22. de la Cruz MJ
23. Mironska R
24. Stokes DL
25. Aitchison JD
26. Jarrold MF
27. Gerton JL
28. Ludtke SJ
29. Akey CW
30. Chait BT
31. Sali A
32. Rout MP
(2018) Integrative structure and functional anatomy of a nuclear pore complex
Nature 555:475–482.
https://doi.org/10.1038/nature26003
- PubMed
- Google Scholar
1. Krammer F
(2020) SARS-CoV-2 vaccines in development
Nature 586:516–527.
https://doi.org/10.1038/s41586-020-2798-3
- PubMed
- Google Scholar
1. Lan J
2. Ge J
3. Yu J
4. Shan S
5. Zhou H
6. Fan S
7. Zhang Q
8. Shi X
9. Wang Q
10. Zhang L
11. Wang X
(2020) Structure of the SARS-CoV-2 spike receptor-binding domain bound to the ACE2 receptor
Nature 581:215–220.
https://doi.org/10.1038/s41586-020-2180-5
- PubMed
- Google Scholar
1. Laurini E
2. Marson D
3. Aulic S
4. Fermeglia A
5. Pricl S
(2021) Molecular rationale for SARS-CoV-2 spike circulating mutations able to escape bamlanivimab and etesevimab monoclonal antibodies
Scientific Reports 11:20274.
https://doi.org/10.1038/s41598-021-99827-3
- PubMed
- Google Scholar
(2020) Functional assessment of cell entry and receptor usage for SARS-CoV-2 and other lineage B betacoronaviruses
Nature Microbiology 5:562–569.
https://doi.org/10.1038/s41564-020-0688-y
- PubMed
- Google Scholar
1. Li T
2. Xue W
3. Zheng Q
4. Song S
5. Yang C
6. Xiong H
7. Zhang S
8. Hong M
9. Zhang Y
10. Yu H
11. Zhang Y
12. Sun H
13. Huang Y
14. Deng T
15. Chi X
16. Li J
17. Wang S
18. Zhou L
19. Chen T
20. Wang Y
21. Cheng T
22. Zhang T
23. Yuan Q
24. Zhao Q
25. Zhang J
26. McLellan JS
27. Zhou ZH
28. Zhang Z
29. Li S
30. Gu Y
31. Xia N
(2021) Cross-neutralizing antibodies bind a SARS-CoV-2 cryptic site and resist circulating variants
Nature Communications 12:5652.
https://doi.org/10.1038/s41467-021-25997-3
- Google Scholar
1. Mannar D
2. Saville JW
3. Zhu X
4. Srivastava SS
5. Berezuk AM
6. Tuttle KS
7. Marquez AC
8. Sekirov I
9. Subramaniam S
(2022) SARS-CoV-2 Omicron variant: Antibody evasion and cryo-EM structure of spike protein-ACE2 complex
Science 375:760–764.
https://doi.org/10.1126/science.abn7760
- PubMed
- Google Scholar
1. Mast FD
2. Fridy PC
3. Ketaren NE
4. Wang J
5. Jacobs EY
6. Olivier JP
7. Sanyal T
8. Molloy KR
9. Schmidt F
10. Rutkowska M
11. Weisblum Y
12. Rich LM
13. Vanderwall ER
14. Dambrauskas N
15. Vigdorovich V
16. Keegan S
17. Jiler JB
18. Stein ME
19. Olinares PDB
20. Herlands L
21. Hatziioannou T
22. Sather DN
23. Debley JS
24. Fenyö D
25. Sali A
26. Bieniasz PD
27. Aitchison JD
28. Chait BT
29. Rout MP
(2021) Highly synergistic combinations of nanobodies that target SARS-CoV-2 and are resistant to escape
eLife 10:e73027.
https://doi.org/10.7554/eLife.73027
- PubMed
- Google Scholar
1. McCallum M
2. De Marco A
3. Lempp FA
4. Tortorici MA
5. Pinto D
6. Walls AC
7. Beltramello M
8. Chen A
9. Liu Z
10. Zatta F
11. Zepeda S
12. di Iulio J
13. Bowen JE
14. Montiel-Ruiz M
15. Zhou J
16. Rosen LE
17. Bianchi S
18. Guarino B
19. Fregni CS
20. Abdelnabi R
21. Foo S-YC
22. Rothlauf PW
23. Bloyet L-M
24. Benigni F
25. Cameroni E
26. Neyts J
27. Riva A
28. Snell G
29. Telenti A
30. Whelan SPJ
31. Virgin HW
32. Corti D
33. Pizzuto MS
34. Veesler D
(2021) N-terminal domain antigenic mapping reveals a site of vulnerability for SARS-CoV-2
Cell 184:2332–2347.
https://doi.org/10.1016/j.cell.2021.03.028
- PubMed
- Google Scholar
1. Mendoza EJ
2. Manguiat K
3. Wood H
4. Drebot M
(2020) Two Detailed Plaque Assay Protocols for the Quantification of Infectious SARS-CoV-2
Current Protocols in Microbiology 57:ecpmc105.
https://doi.org/10.1002/cpmc.105
- PubMed
- Google Scholar
1. Meng Y
2. Irwin DM
3. Shen Y
(2023) Ecology of SARS-CoV-2 in the post-pandemic era
The Lancet. Microbe 4:e208.
https://doi.org/10.1016/S2666-5247(22)00361-5
- PubMed
- Google Scholar
(2022) ColabFold: making protein folding accessible to all
Nature Methods 19:679–682.
https://doi.org/10.1038/s41592-022-01488-1
- PubMed
- Google Scholar
1. Motozono C
2. Toyoda M
3. Zahradnik J
4. Saito A
5. Nasser H
6. Tan TS
7. Ngare I
8. Kimura I
9. Uriu K
10. Kosugi Y
11. Yue Y
12. Shimizu R
13. Ito J
14. Torii S
15. Yonekawa A
16. Shimono N
17. Nagasaki Y
18. Minami R
19. Toya T
20. Sekiya N
21. Fukuhara T
22. Matsuura Y
23. Schreiber G
24. Genotype to Phenotype Japan (G2P-Japan) Consortium
25. Ikeda T
26. Nakagawa S
27. Ueno T
28. Sato K
(2021) SARS-CoV-2 spike L452R variant evades cellular immunity and increases infectivity
Cell Host & Microbe 29:1124–1136.
https://doi.org/10.1016/j.chom.2021.06.006
- PubMed
- Google Scholar
1. Muyldermans S
(2013) Nanobodies: natural single-domain antibodies
Annual Review of Biochemistry 82:775–797.
https://doi.org/10.1146/annurev-biochem-063011-092449
- PubMed
- Google Scholar
1. Pettersen EF
2. Goddard TD
3. Huang CC
4. Meng EC
5. Couch GS
6. Croll TI
7. Morris JH
8. Ferrin TE
(2021) UCSF ChimeraX: Structure visualization for researchers, educators, and developers
Protein Science 30:70–82.
https://doi.org/10.1002/pro.3943
- PubMed
- Google Scholar
1. Pilz S
2. Ioannidis JPA
(2023) Does natural and hybrid immunity obviate the need for frequent vaccine boosters against SARS-CoV-2 in the endemic phase?
European Journal of Clinical Investigation 53:e13906.
https://doi.org/10.1111/eci.13906
- PubMed
- Google Scholar
1. Planas D
2. Veyer D
3. Baidaliuk A
4. Staropoli I
5. Guivel-Benhassine F
6. Rajah MM
7. Planchais C
8. Porrot F
9. Robillard N
10. Puech J
11. Prot M
12. Gallais F
13. Gantner P
14. Velay A
15. Le Guen J
16. Kassis-Chikhani N
17. Edriss D
18. Belec L
19. Seve A
20. Courtellemont L
21. Péré H
22. Hocqueloux L
23. Fafi-Kremer S
24. Prazuck T
25. Mouquet H
26. Bruel T
27. Simon-Lorière E
28. Rey FA
29. Schwartz O
(2021) Reduced sensitivity of SARS-CoV-2 variant Delta to antibody neutralization
Nature 596:276–280.
https://doi.org/10.1038/s41586-021-03777-9
- PubMed
- Google Scholar
1. Polack FP
2. Thomas SJ
3. Kitchin N
4. Absalon J
5. Gurtman A
6. Lockhart S
7. Perez JL
8. Pérez Marc G
9. Moreira ED
10. Zerbini C
11. Bailey R
12. Swanson KA
13. Roychoudhury S
14. Koury K
15. Li P
16. Kalina WV
17. Cooper D
18. Frenck RW Jr
19. Hammitt LL
20. Türeci Ö
21. Nell H
22. Schaefer A
23. Ünal S
24. Tresnan DB
25. Mather S
26. Dormitzer PR
27. Şahin U
28. Jansen KU
29. Gruber WC
30. C4591001 Clinical Trial Group
(2020) Safety and Efficacy of the BNT162b2 mRNA Covid-19 Vaccine
The New England Journal of Medicine 383:2603–2615.
https://doi.org/10.1056/NEJMoa2034577
- PubMed
- Google Scholar
1. Rausch R
(1991)
Effects of temporal lobe surgery on behavior

Advances in Neurology 55:279–292.
- PubMed
- Google Scholar
(2012) Clustering biomolecular complexes by residue contacts similarity
Proteins 80:1810–1817.
https://doi.org/10.1002/prot.24078
- PubMed
- Google Scholar
1. Rout MP
2. Sali A
(2019) Principles for integrative structural biology studies
Cell 177:1384–1403.
https://doi.org/10.1016/j.cell.2019.05.016
- PubMed
- Google Scholar
1. Sali A
(2021) From integrative structural biology to cell biology
The Journal of Biological Chemistry 296:100743.
https://doi.org/10.1016/j.jbc.2021.100743
- PubMed
- Google Scholar
(2021) Using Integrative Modeling Platform to compute, validate, and archive a model of a protein complex structure
Protein Science 30:250–261.
https://doi.org/10.1002/pro.3995
- PubMed
- Google Scholar
Software
1. Sanyal T
(2024) Integrativemodeling/Nbspike, version swh:1:rev:cedbb9bfe83695187ccab0b416b1e85f4b09759f
Software Heritage.

https://archive.softwareheritage.org/swh:1:dir:951f79f25238d840c9c3527c1666ea3ba4f54ca0;origin=https://github.com/integrativemodeling/nbspike;visit=swh:1:snp:3236ecf8876b8002943690d6b885c39a14ff32ec;anchor=swh:1:rev:cedbb9bfe83695187ccab0b416b1e85f4b09759f
1. Starr TN
2. Greaney AJ
3. Addetia A
4. Hannon WW
5. Choudhary MC
6. Dingens AS
7. Li JZ
8. Bloom JD
(2021a) Prospective mapping of viral mutations that escape antibodies used to treat COVID-19
Science 371:850–854.
https://doi.org/10.1126/science.abf9302
- PubMed
- Google Scholar
(2021b) Complete map of SARS-CoV-2 RBD mutations that escape the monoclonal antibody LY-CoV555 and its cocktail with LY-CoV016
Cell Reports. Medicine 2:100255.
https://doi.org/10.1016/j.xcrm.2021.100255
- PubMed
- Google Scholar
1. Starr TN
2. Greaney AJ
3. Hannon WW
4. Loes AN
5. Hauser K
6. Dillen JR
7. Ferri E
8. Farrell AG
9. Dadonaite B
10. McCallum M
11. Matreyek KA
12. Corti D
13. Veesler D
14. Snell G
15. Bloom JD
(2022) Shifting mutational constraints in the SARS-CoV-2 receptor-binding domain during viral evolution
Science 377:420–424.
https://doi.org/10.1126/science.abo7896
- PubMed
- Google Scholar
1. Takashita E
2. Yamayoshi S
3. Simon V
4. van Bakel H
5. Sordillo EM
6. Pekosz A
7. Fukushi S
8. Suzuki T
9. Maeda K
10. Halfmann P
11. Sakai-Tagawa Y
12. Ito M
13. Watanabe S
14. Imai M
15. Hasegawa H
16. Kawaoka Y
(2022) Efficacy of antibodies and antiviral drugs against omicron BA.2.12.1, BA.4, and BA.5 subvariants
The New England Journal of Medicine 387:468–470.
https://doi.org/10.1056/NEJMc2207519
- PubMed
- Google Scholar
1. Tamura T
2. Ito J
3. Uriu K
4. Zahradnik J
5. Kida I
6. Anraku Y
7. Nasser H
8. Shofa M
9. Oda Y
10. Lytras S
11. Nao N
12. Itakura Y
13. Deguchi S
14. Suzuki R
15. Wang L
16. Begum MM
17. Kita S
18. Yajima H
19. Sasaki J
20. Sasaki-Tabata K
21. Shimizu R
22. Tsuda M
23. Kosugi Y
24. Fujita S
25. Pan L
26. Sauter D
27. Yoshimatsu K
28. Suzuki S
29. Asakura H
30. Nagashima M
31. Sadamasu K
32. Yoshimura K
33. Yamamoto Y
34. Nagamoto T
35. Schreiber G
36. Maenaka K
37. Ito H
38. Misawa N
39. Kimura I
40. Suganami M
41. Chiba M
42. Yoshimura R
43. Yasuda K
44. Iida K
45. Ohsumi N
46. Strange AP
47. Takahashi O
48. Ichihara K
49. Shibatani Y
50. Nishiuchi T
51. Kato M
52. Ferdous Z
53. Mouri H
54. Shishido K
55. Sawa H
56. Hashimoto R
57. Watanabe Y
58. Sakamoto A
59. Yasuhara N
60. Suzuki T
61. Kimura K
62. Nakajima Y
63. Nakagawa S
64. Wu J
65. Shirakawa K
66. Takaori-Kondo A
67. Nagata K
68. Kazuma Y
69. Nomura R
70. Horisawa Y
71. Tashiro Y
72. Kawai Y
73. Irie T
74. Kawabata R
75. Motozono C
76. Toyoda M
77. Ueno T
78. Hashiguchi T
79. Ikeda T
80. Fukuhara T
81. Saito A
82. Tanaka S
83. Matsuno K
84. Takayama K
85. Sato K
86. The Genotype to Phenotype Japan (G2P-Japan) Consortium
(2023) Virological characteristics of the SARS-CoV-2 XBB variant derived from recombination of two Omicron subvariants
Nature Communications 14:2800.
https://doi.org/10.1038/s41467-023-38435-3
- Google Scholar
1. Tegally H
2. Moir M
3. Everatt J
4. Giovanetti M
5. Scheepers C
6. Wilkinson E
7. Subramoney K
8. Makatini Z
9. Moyo S
10. Amoako DG
11. Baxter C
12. Althaus CL
13. Anyaneji UJ
14. Kekana D
15. Viana R
16. Giandhari J
17. Lessells RJ
18. Maponga T
19. Maruapula D
20. Choga W
21. Matshaba M
22. Mbulawa MB
23. Msomi N
24. Naidoo Y
25. Pillay S
26. Sanko TJ
27. San JE
28. Scott L
29. Singh L
30. Magini NA
31. Smith-Lawrence P
32. Stevens W
33. Dor G
34. Tshiabuila D
35. Wolter N
36. Preiser W
37. Treurnicht FK
38. Venter M
39. Chiloane G
40. McIntyre C
41. O’Toole A
42. Ruis C
43. Peacock TP
44. Roemer C
45. Kosakovsky Pond SL
46. Williamson C
47. Pybus OG
48. Bhiman JN
49. Glass A
50. Martin DP
51. Jackson B
52. Rambaut A
53. Laguda-Akingba O
54. Gaseitsiwe S
55. von Gottberg A
56. de Oliveira T
57. NGS-SA consortium
(2022) Emergence of SARS-CoV-2 Omicron lineages BA.4 and BA.5 in South Africa
Nature Medicine 28:1785–1790.
https://doi.org/10.1038/s41591-022-01911-2
- PubMed
- Google Scholar
(2016) BRAID: A unifying paradigm for the analysis of combined drug action
Scientific Reports 6:25523.
https://doi.org/10.1038/srep25523
- PubMed
- Google Scholar
1. Walls AC
2. Park YJ
3. Tortorici MA
4. Wall A
5. McGuire AT
6. Veesler D
(2020) Structure, function, and antigenicity of the SARS-CoV-2 spike glycoprotein
Cell 181:281–292.
https://doi.org/10.1016/j.cell.2020.02.058
- PubMed
- Google Scholar
1. Wang Y
2. Liu C
3. Zhang C
4. Wang Y
5. Hong Q
6. Xu S
7. Li Z
8. Yang Y
9. Huang Z
10. Cong Y
(2022a) Structural basis for SARS-CoV-2 Delta variant recognition of ACE2 receptor and broadly neutralizing antibodies
Nature Communications 13:871.
https://doi.org/10.1038/s41467-022-28528-w
- Google Scholar
1. Wang Q
2. Guo Y
3. Iketani S
4. Nair MS
5. Li Z
6. Mohri H
7. Wang M
8. Yu J
9. Bowen AD
10. Chang JY
11. Shah JG
12. Nguyen N
13. Chen Z
14. Meyers K
15. Yin MT
16. Sobieszczyk ME
17. Sheng Z
18. Huang Y
19. Liu L
20. Ho DD
(2022b) Antibody evasion by SARS-CoV-2 Omicron subvariants BA.2.12.1, BA.4 and BA.5
Nature 608:603–608.
https://doi.org/10.1038/s41586-022-05053-w
- Google Scholar
1. Wang F
2. Li L
3. Dou Y
4. Shi R
5. Duan X
6. Liu H
7. Zhang J
8. Liu D
9. Wu J
10. He Y
11. Lan J
12. Lu B
13. Feng H
14. Yan J
(2022c) Etesevimab in combination with JS026 neutralizing SARS-CoV-2 and its variants
Emerging Microbes & Infections 11:548–551.
https://doi.org/10.1080/22221751.2022.2032374
- Google Scholar
1. Wang Z
2. Muecksch F
3. Cho A
4. Gaebler C
5. Hoffmann HH
6. Ramos V
7. Zong S
8. Cipolla M
9. Johnson B
10. Schmidt F
11. DaSilva J
12. Bednarski E
13. Ben Tanfous T
14. Raspe R
15. Yao K
16. Lee Y
17. Chen T
18. Turroja M
19. Milard KG
20. Dizon J
21. Kaczynska A
22. Gazumyan A
23. Oliveira TY
24. Rice CM
25. Caskey M
26. Bieniasz P
27. Hatziioannou T
28. Barnes CO
29. Nussenzweig M
(2022d) Conserved neutralizing epitopes on the N-Terminal domain of variant SARS-CoV-2 spike proteins
SSRN Electronic Journal 01:4023678.
https://doi.org/10.2139/ssrn.4023678
- Google Scholar
1. Watanabe Y
2. Allen JD
3. Wrapp D
4. McLellan JS
5. Crispin M
(2020) Site-specific glycan analysis of the SARS-CoV-2 spike
Science 369:330–333.
https://doi.org/10.1126/science.abb9983
- PubMed
- Google Scholar
1. Webb B
2. Sali A
(2016) Comparative protein structure modeling using MODELLER
Current Protocols in Bioinformatics 54:5.
https://doi.org/10.1002/cpbi.3
- PubMed
- Google Scholar
1. Webb B
2. Viswanath S
3. Bonomi M
4. Pellarin R
5. Greenberg CH
6. Saltzberg D
7. Sali A
(2018) Integrative structure modeling with the Integrative Modeling Platform
Protein Science 27:245–258.
https://doi.org/10.1002/pro.3311
- PubMed
- Google Scholar
1. Wec AZ
2. Wrapp D
3. Herbert AS
4. Maurer DP
5. Haslwanter D
6. Sakharkar M
7. Jangra RK
8. Dieterle ME
9. Lilov A
10. Huang D
11. Tse LV
12. Johnson NV
13. Hsieh C-L
14. Wang N
15. Nett JH
16. Champney E
17. Burnina I
18. Brown M
19. Lin S
20. Sinclair M
21. Johnson C
22. Pudi S
23. Bortz R III
24. Wirchnianski AS
25. Laudermilch E
26. Florez C
27. Fels JM
28. O’Brien CM
29. Graham BS
30. Nemazee D
31. Burton DR
32. Baric RS
33. Voss JE
34. Chandran K
35. Dye JM
36. McLellan JS
37. Walker LM
(2020) Broad neutralization of SARS-related viruses by human monoclonal antibodies
Science 369:731–736.
https://doi.org/10.1126/science.abc7424
- Google Scholar
1. Wooten DJ
2. Albert R
(2021) synergy: a Python library for calculating, analyzing and visualizing drug combination synergy
Bioinformatics 37:1473–1474.
https://doi.org/10.1093/bioinformatics/btaa826
- PubMed
- Google Scholar
Preprint
(2020) Structural Basis for Potent Neutralization of Betacoronaviruses by Single-Domain Camelid Antibodies
bioRxiv.
https://doi.org/10.1101/2020.03.26.010165
- Google Scholar
1. Yan Q
2. Hou R
3. Huang X
4. Zhang Y
5. He P
6. Zhang Y
7. Liu B
8. Wang Q
9. Rao H
10. Chen X
11. Zhao X
12. Niu X
13. Zhao J
14. Xiong X
15. Chen L
(2022) Shared IGHV1-69-encoded neutralizing antibodies contribute to the emergence of L452R substitution in SARS-CoV-2 variants
Emerging Microbes & Infections 11:2749–2761.
https://doi.org/10.1080/22221751.2022.2140611
- PubMed
- Google Scholar
1. Zhang J
2. Xiao T
3. Cai Y
4. Lavine CL
5. Peng H
6. Zhu H
7. Anand K
8. Tong P
9. Gautam A
10. Mayer ML
11. Walsh RM Jr
12. Rits-Volloch S
13. Wesemann DR
14. Yang W
15. Seaman MS
16. Lu J
17. Chen B
(2021) Membrane fusion and immune evasion by the spike protein of SARS-CoV-2 Delta variant
Science 374:1353–1360.
https://doi.org/10.1126/science.abl9463
- PubMed
- Google Scholar
1. Zhao P
2. Praissman JL
3. Grant OC
4. Cai Y
5. Xiao T
6. Rosenbalm KE
7. Aoki K
8. Kellman BP
9. Bridger R
10. Barouch DH
11. Brindley MA
12. Lewis NE
13. Tiemeyer M
14. Chen B
15. Woods RJ
16. Wells L
(2020) Virus-receptor interactions of glycosylated SARS-CoV-2 spike and human ACE2 receptor
Cell Host & Microbe 28:586–601.
https://doi.org/10.1016/j.chom.2020.08.004
- PubMed
- Google Scholar
1. Zhao Z
2. Zhou J
3. Tian M
4. Huang M
5. Liu S
6. Xie Y
7. Han P
8. Bai C
9. Han P
10. Zheng A
11. Fu L
12. Gao Y
13. Peng Q
14. Li Y
15. Chai Y
16. Zhang Z
17. Zhao X
18. Song H
19. Qi J
20. Wang Q
21. Wang P
22. Gao GF
(2022) Omicron SARS-CoV-2 mutations stabilize spike up-RBD conformation and lead to a non-RBM-binding monoclonal antibody escape
Nature Communications 13:4958.
https://doi.org/10.1038/s41467-022-32665-7
- PubMed
- Google Scholar

Article and author information

Author details

Natalia E Ketaren

Laboratory of Cellular and Structural Biology, The Rockefeller University, New York, United States

Contribution
Conceptualization, Formal analysis, Validation, Investigation, Visualization, Methodology, Writing – original draft, Writing – review and editing

Contributed equally with
Fred D Mast, Peter C Fridy, Jean Paul Olivier and Tanmoy Sanyal

Competing interests
Inventor on a provisional patent (US20230331824A1) describing the anti-spike nanobodies described in this manuscript

"This ORCID iD identifies the author of this article:" 0000-0002-7869-6162
Fred D Mast

Center for Global Infectious Disease Research, Seattle Children's Research Institute, Seattle, United States

Contribution
Conceptualization, Formal analysis, Validation, Investigation, Visualization, Methodology, Writing – original draft, Writing – review and editing

Contributed equally with
Natalia E Ketaren, Peter C Fridy, Jean Paul Olivier and Tanmoy Sanyal

Competing interests
Inventor on a provisional patent (US20230331824A1) describing the anti-spike nanobodies described in this manuscript

"This ORCID iD identifies the author of this article:" 0000-0002-2177-6647
Peter C Fridy

Laboratory of Cellular and Structural Biology, The Rockefeller University, New York, United States

Contribution
Conceptualization, Formal analysis, Validation, Investigation, Visualization, Methodology, Writing – original draft, Writing – review and editing

Contributed equally with
Natalia E Ketaren, Fred D Mast, Jean Paul Olivier and Tanmoy Sanyal

Competing interests
Inventor on a provisional patent (US20230331824A1) describing the anti-spike nanobodies described in this manuscript

"This ORCID iD identifies the author of this article:" 0000-0002-8208-9154
Jean Paul Olivier

Center for Global Infectious Disease Research, Seattle Children's Research Institute, Seattle, United States

Contribution
Conceptualization, Formal analysis, Validation, Investigation, Visualization, Methodology, Writing – original draft, Writing – review and editing

Contributed equally with
Natalia E Ketaren, Fred D Mast, Peter C Fridy and Tanmoy Sanyal

Competing interests
Inventor on a provisional patent (US20230331824A1) describing the anti-spike nanobodies described in this manuscript

"This ORCID iD identifies the author of this article:" 0000-0002-2197-271X
Tanmoy Sanyal

Department of Bioengineering and Therapeutic Sciences, Department of Pharmaceutical Chemistry, California Institute for Quantitative Biosciences, Byers Hall, University of California, San Francisco, San Francisco, United States

Contribution
Conceptualization, Formal analysis, Validation, Investigation, Visualization, Methodology, Writing – original draft, Writing – review and editing

Contributed equally with
Natalia E Ketaren, Fred D Mast, Peter C Fridy and Jean Paul Olivier

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-6009-9431
Andrej Sali

Department of Bioengineering and Therapeutic Sciences, Department of Pharmaceutical Chemistry, California Institute for Quantitative Biosciences, Byers Hall, University of California, San Francisco, San Francisco, United States

Contribution
Supervision, Funding acquisition, Writing – review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0003-0435-6197
Brian T Chait

Laboratory of Mass Spectrometry and Gaseous Ion Chemistry, The Rockefeller University, New York, United States

Contribution
Conceptualization, Supervision, Visualization, Methodology, Project administration, Writing – review and editing

For correspondence
chait@rockefeller.edu

Competing interests
Inventor on a provisional patent (US20230331824A1) describing the anti-spike nanobodies described in this manuscript
Michael P Rout

Laboratory of Cellular and Structural Biology, The Rockefeller University, New York, United States

Contribution
Conceptualization, Supervision, Visualization, Methodology, Project administration, Writing – review and editing

For correspondence
rout@rockefeller.edu

Competing interests
Inventor on a provisional patent (US20230331824A1) describing the anti-spike nanobodies described in this manuscript

"This ORCID iD identifies the author of this article:" 0000-0003-2010-706X
John D Aitchison
1. Center for Global Infectious Disease Research, Seattle Children's Research Institute, Seattle, United States
2. Department of Pediatrics, University of Washington, Seattle, United States
3. Department of Biochemistry, University of Washington, Seattle, United States
Contribution
Conceptualization, Supervision, Visualization, Methodology, Project administration, Writing – review and editing

For correspondence
John.Aitchison@seattlechildrens.org

Competing interests
Inventor on a provisional patent (US20230331824A1) describing the anti-spike nanobodies described in this manuscript

"This ORCID iD identifies the author of this article:" 0000-0002-9153-6497

Funding

G. Harold and Leila Y. Mathers Foundation

Brian T Chait
Michael P Rout
John D Aitchison

Robertson Therapeutic Development Fund

Brian T Chait
Michael P Rout
John D Aitchison

Jain Foundation

Brian T Chait
Michael P Rout
John D Aitchison

National Institutes of Health (P41GM109824)

Andrej Sali
Brian T Chait
Michael P Rout
John D Aitchison

National Institutes of Health (R01GM083960)

Andrej Sali

American Lung Association (923055)

Brian T Chait
Michael P Rout
John D Aitchison

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

Acknowledgements

We are very grateful to The Fisher Drug Discovery Resource Center (DDRC), Rockefeller University (RRID:SCR_020985), and the rest of the Aitchison, Chait and Rout laboratories for intellectual support. Funding. American Lung Association (JDA, BTC, MPR), G Harold and Leila Y Mathers Charitable Foundation (JDA, BTC, MPR), Robertson Therapeutic Development Fund (JDA, BTC, MPR) Jain Foundation (JDA, BTC, MPR), National Institutes of Health grant P41GM109824 (JDA, BTC, AS, MPR), and National Institutes of Health grant R01GM083960 (AS).

Version history

Sent for peer review: June 21, 2023
Preprint posted: July 14, 2023 (view preprint)
Preprint posted: November 9, 2023 (view preprint)
Preprint posted: April 15, 2024 (view preprint)
Version of Record published: May 7, 2024 (version 1)

Cite all versions

You can cite all versions using the DOI https://doi.org/10.7554/eLife.89423. This DOI represents all versions, and will always resolve to the latest one.

Copyright

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