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, Fridy 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.
This paper by Aitchison and colleagues describes nanobody neutralizing and binding activity against various SARS-CoV-2 variants of concern. The findings are important in that the described nanobodies may have broad therapeutic relevance against current and future variants of concern and may be able to avoid significant resistance. The claims are incomplete: while the study is well-executed and uses a nice balance of biochemical and cellular assays, the efficacy of the proposed nanobody library against VOCs is not completely supported as IC50 values appear to increase against newer variants and are higher than previously used therapeutic bNAbs, animal data showing in vivo efficacy is lacking, and protection against future possible variants is not proven.
SARS-CoV-2 has infected >40% of the world’s population (Collaborators 2022) resulting in a devastating loss of life. As the SARS-CoV-2 pandemic enters its endemic phase (Meng, Irwin 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 (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, Anderson et al. 2020, Krammer 2020, Letko, Marzi et al. 2020, Polack, Thomas 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, Park et al. 2020, Jackson, Farzan et al. 2022). Glycosylation is most extensive on the NTD and the S2 domain, whereas the RBD is largely glycan free (Watanabe, Allen et al. 2020, Zhao, Praissman et al. 2020). Consequently, it is unsurprising that the most antigenic domain on spike is the RBD, where the vast majority of neutralizing antibodies have been shown to bind. A comprehensive mapping of the epitopes from 1,640 neutralizing monoclonal antibodies (mAbs), all targeting the RBD, revealed 12 epitope groups (Cao, Yisimayi et al. 2022). This data, combined with previous studies mapping the epitopes of antibodies targeting spike, reveal a total of 19 mAb epitope groups, including 7 on the NTD (Wang, Muecksch et al. 2022). Very few anti-S2 antibodies have been shown to be effective therapeutic options (Wec, Wrapp et al. 2020), likely due to the shielding effect of S2 glycans (Grant, Montgomery 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, Loes et al. 2021, Starr, Greaney et al. 2022, Dadonaite, Crawford 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, Saville et al. 2022). The omicron BA.4/BA.5 variants, which have identical spike proteins (Tegally, Moir 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, Yisimayi et al. 2022, Hachmann, Miller et al. 2022, Wang, Guo et al. 2022). Of the monoclonal antibodies (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, Yamayoshi et al. 2022, Focosi, Quiroga et al. 2023, Imai, Ito et al. 2023). As a result, no mAb therapy is currently approved by the FDA for treatment of SARS-CoV-2 infection (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, Li et al. 2014, Mast, Fridy 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, Fridy et al. 2021), retains efficacy against circulating variants of concern (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 remain efficacious against delta, omicron BA.1 and omicron BA.4/5, targeting all major regions of spike
From the original nanobody repertoire that we generated against SARS-CoV-2 wild-type spike protein (Mast, Fridy et al. 2021), representative nanobodies from all 10 structurally mapped epitope groups that we previously identified (Mast, Fridy et al. 2021, Cross, Fridy et al. 2023), were selected for SARS-CoV-2 pseudovirus (PSV) neutralization assays against the SARS-CoV-2 delta and omicron BA.1 strains (Fig. 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.
These observations are consistent with data from mAbs approved by the FDA, where nanobody epitopes from groups I, II and IV (Fig. 1B) overlap with the epitopes of five mAbs (grouped into three classes on RBD: class I, class II and class III) effective against delta: three class I (etesivimab, casirivimab and amubarvimab) (Greaney, Starr et al. 2021, Li, Xue et al. 2021, Planas, Veyer et al. 2021, Takashita, Yamayoshi et al. 2022, Cox, Peacock et al. 2023) and two class III (imdemivab and bebtelovimab) (Cox, Peacock et al. 2023). Additionally, the majority of our group I and II nanobodies neutralized delta better than wild-type (Fig.1), as observed with the mAb etesivimab (Wang, Li et al. 2022, Cox, Peacock 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 decreased compared to wild-type. This decrease in neutralization potency largely correlates with the accumulation of omicron BA.1 specific mutations throughout the RBD, which likely alters the epitope-binding site of these nanobodies, weakening their interaction with BA.1 spike (Fig. 1B). Concomitantly, groups I, I/II, I/IV and the anti-S2 nanobody 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.
Of the nanobodies that neutralized both delta and omicron BA.1, representatives from each of the nanobody epitope groups were selected for SPR analysis, where S1 binders with mapped epitopes that neutralized one or both variants well, were prioritized. The binding data as determined by SPR experiments performed against either spike S1 domain or RBD from delta and omicron BA.1, largely correlated with the neutralization data, where nanobodies that retained binding to the spike variant also neutralized the virus. 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, Fridy 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 (Fig. S1) (Focosi, Quiroga et al. 2023). The remaining six nanobody groups tested, including the major nanobody groups III and IV, showed no detectable binding to omicron BA.4/BA.5 RBD, save for RBD-40 (group III) and RBD-46 (group IV), which exhibited a ∼40-fold and ∼20-fold decrease in affinity respectively, compared to wild-type. The likely reason behind this loss of binding is the concentration of omicron BA.4/BA.5 specific mutations that overlap with the epitope regions of the affected nanobodies, resulting in epitopes altered enough to abolish or significantly reduce binding (Fig. 2B). Interestingly, one nanobody, S1-46 (group VII, Fig 1B and 2B) retained wild-type binding affinity to RBDs from delta, omicron BA.1, and omicron BA.4/BA.5 (Fig. 2, Table 1). S1-46 binds a region on spike that is conserved across all variants to date and is not targeted by any of the mAbs that previously received EUA by the FDA (Cox, Peacock 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 (PRNT) as previously described (Mast, Fridy et al. 2021). The results showed that all 5 nanobodies tested neutralized omicron BA.5 live virus (Fig. 3), supporting our SPR observations, to reveal our nanobody repertoire generated against wildtype spike has retained efficacy against omicron BA.4/BA.5.
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, has 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, Yamayoshi et al., Cox, Peacock et al., Focosi, Quiroga et al.). These structural differences have also somewhat 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 (Fig. 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, Liu et al. 2022); 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 (Fig. 1B and 2B). 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 (Fig. 1B and 2B), 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, Marson et al. 2021, Starr, Greaney et al. 2021, Starr, Greaney et al. 2021). 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, Toyoda et al. 2021, Yan, Hou 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, Zhou 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 (Fig. 1B and 2B). 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 (Fig. 1A and 2A). 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 I), S1-4 (Group I/II, Class III), S1-RBD-5 (Group I/IV, Class III), 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 sub-lineages, 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 group II nanobodies against omicron variants BA.4/5 may be because of the T376A, D405N and R408S substitutions, which lie within and near the epitope regions of groups I and II respectively (Fig. 2B). Importantly, at least 5 nanobodies from groups I, I/II, II and V retained neutralization activity against omicron BA.5 (Fig. 3), 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 against both delta and omicron BA.1, where the significantly decreased binding affinity correlated to decreased neutralization potency, which is likely due to the L452R mutation against delta, and the large number of amino acid substitutions compared to wild-type for omicron BA.1 (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 (Fig. 2A), yet none neutralized the variant in the pseudovirus assay (Fig. 1A). These nanobodies are in groups III, II and IV respectively (Fig. 2B), 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 vivo. 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 (Fig. 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 bound both delta and omicron BA.1 with similar affinities (>100-fold decrease compared to wild-type), yet only neutralized omicron BA.1 (Fig. 1A). As mentioned above, the epitope of S1-37 overlaps with the delta L452R mutation (Fig. 1B), which has impacted the effectiveness of numerous neutralizing antibodies (Bian, Gao 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 (Fig. 1B and 2B), showed binding to delta with equal affinity to wild-type (Fig. 2A), yet does not neutralize delta (Fig. 1A). The delta mutation T478K is present within the epitope region of group IV (Fig. 1B and 2B), 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, Krieger 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 one-third of our original repertoire of 116 nanobodies (Mast, Fridy et al.) generated against wild-type SARS-CoV-2 remained effective binders/neutralizers of the variants tested, where nanobodies from 11 (including the 10 mapped epitopes) of the 18 nanobody groups are efficacious against the delta, omicron BA.1, and omicron BA.4/BA.5 lineages. The ability of nanobodies within a group to retain efficacy against different variants, coupled with the structural modeling data, allowed us to further refine our original six structurally modeled RBD epitope groups, for a total of 12 RBD epitope groups (Fig. 4A) compared to the three classes used to group mAbs (all RBD binders) previously approved for EUA by the FDA. When comparing our 12 RBD nanobody epitope groups to the three mAb classes, groups I, III, IV and VII overlap with class I binders; group II, III, IV and V overlap with class II binders; and group II, III and V overlap with class III binders (Fig. S1).
Interestingly, the epitope footprint of many of our nanobody classes extends beyond those of the three mAb classes. 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, Mellado-Sanchez et al. 2022, Cao, Yisimayi et al. 2022); instead, these epitopes extend away from the ACE2 binding site (Fig 1.B and 2B), 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 (Fig. 4B) where: of the 26 nanobodies that showed binding/neutralization to delta, 6 were specific only for delta; of the 21 nanobodies that showed binding/neutralization to omicron BA.1, two were specific for omicron BA.1; none of the 10 nanobodies that retained binding to omicron BA.4 were specific only for omicron BA.4. Strikingly, we have in our cohort 11 nanobodies able to bind delta, omicron BA.1, and omicron BA.4/BA.5 (Fig. 4B), which we further predict will be effective binders against current circulating strains of the virus including omicron XBB and BQ1.1 (Fig. 4C).
Nanobody synergy involving a non-neutralizing nanobody
Previously, we showed that cocktails of nanobodies were more resistant to mutational escape (Mast, Fridy 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 predicted from the neutralization by either nanobody alone (Mast, Fridy et al. 2021). Our present observations of nanobodies that retained binding to variants of the spike RBD (Table 1) despite losing neutralization efficacy (Fig. 1A; Table 1) 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, Fridy et al. 2021) (Fig. 1A, 2A, and 5A 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 (Fig. 1A, 2A, and 5A 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-3 μM (Fig. 5A 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 (Fig. 1A, 2A, & 5A panel iii).
We also tested the broadly neutralizing nanobody S1-RBD-22 in combination with S1-36 (Fig. 5B). Like S1-23, the epitope of S1-36 is opposite that of S1-RBD-22, permitting simultaneous binding to a single RBD (Mast, Fridy 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 (Fig. 1A, 2A, 5B). 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 PSB (Fig. 5B panels ii, and iii).
The synergies observed between S1-1 and S1-23, and between S1-RBD-22 and S1-36, are likely not due to pleiotropic effects, as providing the nonspecific LaM2 nanobody in combination with S1-1 had no impact on the neutralization efficacy of S1-1 against delta PSV (Fig. 5C). 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 (Fig. 5D) (Mast, Fridy et al. 2021).
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 (Fig. 6) (Mast, Fridy 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 (Fig 5 and 6), 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 (Fig. 1B and 2B).
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, De Vlieger et al. 2020, Mast, Fridy 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.
In this study, we identified nanobodies that specifically recognize only certain SARS-CoV-2 variants (Fig. 4B), 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
Series S CM5 sensor chips (Cytiva) were immobilized with recombinant delta Spike S1, delta RBD and omicron RBD at 12.5 ug/mL, 10 ug/mL and 15 ug/mL respectively using EDC/NHS coupling chemistry according to the manufacturer’s instructions. KD measurements were performed on a Biacore 8k (Cytiva) as previously described (Mast, Fridy et al. 2021), using 5-8 concentrations of each nanobody. Data was processed and analyzed using the Biacore Insight Evaluation software.
Integrative structure modeling of nanobody epitopes on the different variant spike proteins, proceeded through the standard four-stage protocol (Kim, Fernandez-Martinez et al. 2018, Webb, Viswanath et al. 2018, Rout and Sali 2019, Sali 2021, Saltzberg, Viswanath 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, Ge et al. 2020)) for the original virus, the 4.30 Å cryo-EM structure of up-RBD pre-fusion spike (7SBO.A; (Zhang, Xiao et al. 2021)) for the delta variant, the ACE2 bound 2.45 Å cryo-EM spike structure (7T9K.B; (Mannar, Saville et al. 2022)) for the omicron BA.1 variant and a structure predicted using AlphaFold-2 for the omicron BA.4/BA.5 variant. 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, De Marco 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, Zhang 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, Evans et al. 2021, Mirdita, Schutze et al. 2022). The protocol included automatic refinement of the CDR loops through an all-atom energy minimization of the AlphaFold2-predicted structure using the AMBER molecular mechanics force field (Hornak, Abel et al. 2006). We verified that the predicted nanobody structures are within 3-4 Å backbone RMSD from the comparative models of these nanobody sequences published previously (Mast, Fridy 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, Trellet 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 don’t 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 (Kd) and neutralization potential (IC50) 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 and normalized to the range from 0-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, Goddard et al. 2021). Files containing input data, scripts and output results are available at https://github.com/integrativemodeling/nbspike/tree/main/integrative_modeling_VOC. 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.
TMPRSS2-expressing Vero E6 cells, 293T/17 cells and 293T-hACE2 cells were cultured as described previously (Mast, Fridy et al. 2021). Briefly, TMPRSS2+ Vero E6 cells were cultured at 37°C in the presence of 5% CO2 in medium composed of in high-glucose Dulbecco’s modified Eagle’s medium (DMEM, Gibco) supplemented with 10% (v/v) FBS and 1 mg/ml geneticin. 293T/17 were cultured at 37°C in the presence of 5% CO2 in a medium composed of DMEM supplemented with 10% (v/v) FBS and penicillin/streptomycin. 293t-hACE2 cells were cultured at 37°C in the presence of 5% CO2 in medium composed of DMEM supplemented with 10% (v/v) FBS, penicillin/streptomycin, 10 mM HEPES, and with 0.1 mM 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
Pseudovirus stocks were prepared and titered as described previously (Mast, Fridy et al. 2021). Variant spike containing plasmids were combined with pHAGE-CMV-Luc2-IRES-ZsGreen-W (BEI Cat # NR-52516) (Crawford, Eguia 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, Fridy et al. 2021).
SARS-CoV-2 pseudovirus neutralization assay
Nanobodies were tested for their neutralization properties as described previously (Mast, Fridy et al. 2021). Briefly, threefold serial dilutions of nanobodies were incubated with pseudotyped SARS-CoV-2 for 1 h 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% CO2. Infected cells were processed between 52 and 60 hr as described above. 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 nonlinear regression (least-squares regression method without weighting). All nanobodies were tested at least two times and with more than one pseudovirus preparation.
Synergy experiments were performed as described previously, (Mast, Fridy 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 one h. After incubation with the virus, the mixture was overlaid on a monolayer of 293-hACE2 cells and left to incubate for 56 h. Luminescence was quantified as described above. Data were processed using the Bivariate Response to Additive Interacting Doses (BRAID) model (Twarog, Stewart et al. 2016) as implemented in the synergy software package for python (Wooten and Albert 2021).
SARS-CoV-2 stocks and titers
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% (v/v) 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 was 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, Manguiat et al. 2020, Mast, Fridy et al. 2021).
We are very grateful to The Fisher Drug Discovery Resource Center (DDRC), Rockefeller University, and the rest of the Aitchison, Chait and Rout laboratories for intellectual support.
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)
National Institutes of Health grant R01GM083960 (AS)
Conceptualization: JDA, BTC, MPR, NEK, PCF, FDM, JPO, TS
Methodology: JDA, BTC, MPR, NEK, FDM, JPO, TS, PCF
Investigation: NEK, FDM, JPO, TS, PCF
Visualization: JDA, BTC, MPR, NEK, PCF, FDM, JPO, TS
Funding acquisition: JDA, BTC, MPR, AS
Project administration: JDA, BTC, MPR
Supervision: JDA, BTC, MPR, AS
Writing – original draft: NEK, FDM, PCF, JPO, TS
Writing – review & editing: JDA, BTC, MPR, NEK, PCF, FDM, JPO, TS, AS
JDA, BTC, MPR, PCF, NEK, FDM, and JPO are inventors on a provisional patent describing the anti-spike nanobodies described in this manuscript.
Supplementary Table 1 – Neutralization data from synergy experiments
- Evolution of Anti-SARS-CoV-2 Therapeutic AntibodiesInt J Mol Sci 23
- Impact of the Delta variant on vaccine efficacy and response strategiesExpert Rev Vaccines 20:1201–1209
- Distinct conformational states of SARS-CoV-2 spike proteinScience 369:1586–1592
- BA.2.12.1, BA.4 and BA.5 escape antibodies elicited by Omicron infectionNature 608:593–602
- Impact of new variants on SARS-CoV-2 infectivity and neutralization: A molecular assessment of the alterations in the spike-host protein interactionsiScience 25
- Estimating global, regional, and national daily and cumulative infections with SARS-CoV-2 through Nov 14, 2021: a statistical analysisLancet 399:2351–2380
- SARS-CoV-2 variant evasion of monoclonal antibodies based on in vitro studiesNat Rev Microbiol 21:112–124
- Protocol and Reagents for Pseudotyping Lentiviral Particles with SARS-CoV-2 Spike Protein for Neutralization AssaysViruses 12
- Expanding and improving nanobody repertoires using a yeast display method: Targeting SARS-CoV-2J Biol Chem 299
- A pseudovirus system enables deep mutational scanning of the full SARS-CoV-2 spikeCell 186:1263–1278
- Convergent Evolution in SARS-CoV-2 Spike Creates a Variant Soup from Which New COVID-19 Waves EmergeInt J Mol Sci 24
- A robust pipeline for rapid production of versatile nanobody repertoiresNature methods 11:1253–1260
- Analysis of the SARS-CoV-2 spike protein glycan shield reveals implications for immune recognitionSci Rep 10
- Comprehensive mapping of mutations in the SARS-CoV-2 receptor-binding domain that affect recognition by polyclonal human plasma antibodiesCell Host Microbe 29:463–476
- Mapping mutations to the SARS-CoV-2 RBD that escape binding by different classes of antibodiesNat Commun 12
- Neutralization Escape by SARS-CoV-2 Omicron Subvariants BA.2.12.1, BA.4, and BA.5N Engl J Med 387:86–88
- Comparison of multiple Amber force fields and development of improved protein backbone parametersProteins 65:712–725
- Efficacy of Antiviral Agents against Omicron Subvariants BQ.1.1 and XBBN Engl J Med 388:89–91
- Mechanisms of SARS-CoV-2 entry into cellsNat Rev Mol Cell Biol 23:3–20
- An mRNA Vaccine against SARS-CoV-2 – Preliminary ReportN Engl J Med 383:1920–1931
- Highly accurate protein structure prediction with AlphaFoldNature 596:583–589
- Integrative structure and functional anatomy of a nuclear pore complexNature 555:475–482
- SARS-CoV-2 vaccines in developmentNature 586:516–527
- Structure of the SARS-CoV-2 spike receptor-binding domain bound to the ACE2 receptorNature 581:215–220
- Molecular rationale for SARS-CoV-2 spike circulating mutations able to escape bamlanivimab and etesevimab monoclonal antibodiesSci Rep 11
- Functional assessment of cell entry and receptor usage for SARS-CoV-2 and other lineage B betacoronavirusesNat Microbiol 5:562–569
- Cross-neutralizing antibodies bind a SARS-CoV-2 cryptic site and resist circulating variantsNat Commun 12
- SARS-CoV-2 Omicron variant: Antibody evasion and cryo-EM structure of spike protein-ACE2 complexScience 375:760–764
- Highly synergistic combinations of nanobodies that target SARS-CoV-2 and are resistant to escapeElife 10
- N-terminal domain antigenic mapping reveals a site of vulnerability for SARS-CoV-2Cell 184:2332–2347
- Two Detailed Plaque Assay Protocols for the Quantification of Infectious SARS-CoV-2Curr Protoc Microbiol 57
- Ecology of SARS-CoV-2 in the post-pandemic eraLancet Microbe 4
- ColabFold: making protein folding accessible to allNat Methods 19:679–682
- SARS-CoV-2 spike L452R variant evades cellular immunity and increases infectivityCell Host Microbe 29:1124–1136
- Nanobodies: natural single-domain antibodiesAnnu Rev Biochem 82:775–797
- UCSF ChimeraX: Structure visualization for researchers, educators, and developersProtein Sci 30:70–82
- Does natural and hybrid immunity obviate the need for frequent vaccine boosters against SARS-CoV-2 in the endemic phase?Eur J Clin Invest 53
- Reduced sensitivity of SARS-CoV-2 variant Delta to antibody neutralizationNature 596:276–280
- Safety and Efficacy of the BNT162b2 mRNA Covid-19 VaccineN Engl J Med 383:2603–2615
- Effects of temporal lobe surgery on behaviorAdv Neurol 55:279–292
- Clustering biomolecular complexes by residue contacts similarityProteins 80:1810–1817
- Principles for Integrative Structural Biology StudiesCell 177:1384–1403
- From integrative structural biology to cell biologyJournal of Biological Chemistry 296
- Using Integrative Modeling Platform to compute, validate, and archive a model of a protein complex structureProtein Sci 30:250–261
- Prospective mapping of viral mutations that escape antibodies used to treat COVID-19Science 371:850–854
- Complete map of SARS-CoV-2 RBD mutations that escape the monoclonal antibody LY-CoV555 and its cocktail with LY-CoV016Cell Rep Med 2
- Shifting mutational constraints in the SARS-CoV-2 receptor-binding domain during viral evolutionScience 377:420–424
- Efficacy of Antibodies and Antiviral Drugs against Omicron BA.2.12.1, BA.4, and BA.5 SubvariantsN Engl J Med 387:468–470
- Emergence of SARS-CoV-2 Omicron lineages BA.4 and BA.5 in South AfricaNat Med 28:1785–1790
- BRAID: A Unifying Paradigm for the Analysis of Combined Drug ActionSci Rep 6
- Structure, Function, and Antigenicity of the SARS-CoV-2 Spike GlycoproteinCell 181:281–292
- Etesevimab in combination with JS026 neutralizing SARS-CoV-2 and its variantsEmerg Microbes Infect 11:548–551
- Antibody evasion by SARS-CoV-2 Omicron subvariants BA.2.12.1, BA.4 and BA.5Nature 608:603–608
- Structural basis for SARS-CoV-2 Delta variant recognition of ACE2 receptor and broadly neutralizing antibodiesNat Commun 13
- Conserved Neutralizing Epitopes on the N-Terminal Domain of Variant SARS-CoV-2 Spike Proteins
- Site-specific glycan analysis of the SARS-CoV-2 spikeScience 369:330–333
- Integrative structure modeling with the Integrative Modeling PlatformProtein Sci 27:245–258
- Broad neutralization of SARS-related viruses by human monoclonal antibodiesScience 369:731–736
- synergy: a Python library for calculating, analyzing and visualizing drug combination synergyBioinformatics 37:1473–1474
- Structural basis for potent neutralization of betacoronaviruses by single-domain camelid antibodies (preprint)
- Shared IGHV1-69-encoded neutralizing antibodies contribute to the emergence of L452R substitution in SARS-CoV-2 variantsEmerg Microbes Infect 11:2749–2761
- Membrane fusion and immune evasion by the spike protein of SARS-CoV-2 Delta variantScience 374:1353–1360
- Virus-Receptor Interactions of Glycosylated SARS-CoV-2 Spike and Human ACE2 ReceptorCell Host Microbe 28:586–601
- Omicron SARS-CoV-2 mutations stabilize spike up-RBD conformation and lead to a non-RBM-binding monoclonal antibody escapeNat Commun 13