key: cord-0866121-s69mcjv3 authors: Mast, Fred D; Fridy, Peter C; Ketaren, Natalia E; Wang, Junjie; Jacobs, Erica Y; Olivier, Jean Paul; Sanyal, Tanmoy; Molloy, Kelly R; Schmidt, Fabian; Rutkowska, Magdalena; Weisblum, Yiska; Rich, Lucille M; Vanderwall, Elizabeth R; Dambrauskas, Nicholas; Vigdorovich, Vladimir; Keegan, Sarah; Jiler, Jacob B; Stein, Milana E; Olinares, Paul Dominic B; Herlands, Louis; Hatziioannou, Theodora; Sather, D Noah; Debley, Jason S; Fenyö, David; Sali, Andrej; Bieniasz, Paul D; Aitchison, John D; Chait, Brian T; Rout, Michael P title: Highly synergistic combinations of nanobodies that target SARS-CoV-2 and are resistant to escape date: 2021-12-07 journal: eLife (Cambridge) DOI: 10.7554/elife.73027 sha: 400674b56113c100af27317face7be07e19c2974 doc_id: 866121 cord_uid: s69mcjv3 The emergence of SARS-CoV-2 variants threatens current vaccines and therapeutic antibodies and urgently demands powerful new therapeutics that can resist viral escape. We therefore generated a large nanobody repertoire to saturate the distinct and highly conserved available epitope space of SARS-CoV-2 spike, including the S1 receptor binding domain, N-terminal domain, and the S2 subunit, to identify new nanobody binding sites that may reflect novel mechanisms of viral neutralization. Structural mapping and functional assays show that indeed these highly stable monovalent nanobodies potently inhibit SARS-CoV-2 infection, display numerous neutralization mechanisms, are effective against emerging variants of concern, and are resistant to mutational escape. Rational combinations of these nanobodies that bind to distinct sites within and between spike subunits exhibit extraordinary synergy and suggest multiple tailored therapeutic and prophylactic strategies. SARS-CoV-2, the viral causative agent of COVID-19, is estimated to have infected some 10% of the world's population, killing a confirmed ~5 million but likely considerably more. Despite the great promise of vaccines, the pandemic is ongoing; inequities in vaccine distribution, waning immunity, the biological and behavioral diversity of the human population, the emergence of viral variants that compromise monoclonal therapies and vaccine efficacy, all challenge current and future containment (Diamond et al., 2021; Lavine et al., 2021; Fraser et al., 2004; Wang et al., 2021a; Wang et al., 2021b) . Thus, the best we can hope for now is an uneasy truce, in which multipronged containment strategies will be required for many years to keep SARS-CoV-2, future variants, and novel coronaviruses at bay (Phillips, 2021; McKenna, 2021; Steenhuysen and Kelland, 2021; Weisblum et al., 2020) . Spike (S), the major surface envelope glycoprotein of the SARS-CoV-2 virion, is key for infection as it attaches the virion to its cognate host surface receptor, angiotensin-converting enzyme 2 (ACE2) protein, and triggers fusion between the host and viral membranes, leading to viral entry into the cytoplasm Wrapp et al., 2020b; Walls et al., 2020) . The spike protein monomer is ~140 kDa, or ~180-200 kDa including its extensive glycosylation, and exists as a homotrimer on the viral surface. Spike is highly dynamic and is composed of two domains: S1, which contains the host receptor binding domain (RBD); and S2, which undergoes large conformational changes that enable fusion of the viral membrane with that of its host (Li et al., 2003; Li, 2016; Letko et al., 2020; Watanabe et al., 2020; Hsieh et al., 2020) . Based on its requirement for entry, the major target of immunotherapeutics has been the RBD (Hartenian et al., 2020; Wu et al., 2020; Baum et al., 2020; Finkelstein et al., 2021; Korber et al., 2020; Trigueiro-Louro et al., 2020; Barnes et al., 2020) . Major immunotherapeutic strategies to date have focused on immune sera and human monoclonal antibodies; however, these therapies now face the emergence of variants, particularly RBD point mutants, which have evolved to bypass the most potent neutralizing human antibodies (Wang et al., 2021b; Liu et al., 2021a; Weisblum et al., 2020; Garcia-Beltran et al., 2021; Starr et al., 2021) . A specific alternative class of single-chain monoclonal antibodies, commonly called nanobodies, are attractive alternatives to traditional monoclonal antibodies (Muyldermans, 2013) . Nanobodies are the smallest single-domain antigen binding proteins identified to date, possessing several potential advantages over conventional monoclonal antibodies. Nanobodies are derived from the variable domain (V H H) of variant heavy chain-only IgGs (HCAb) found in camelids (e.g., llamas, alpacas, and camels). They can bind in modes different from typical antibodies, covering more chemical space and binding with very high affinities (comparable to the very best antibodies) (Jovčevska and Muyldermans, 2020; Muyldermans, 2013) . Their small size (~15 kDa) allows them to bind tightly to otherwise inaccessible epitopes that may be obscured by the glycoprotein coat, as well as minimizing issues of steric hindrance of multiple antibodies binding to adjacent epitopes as observed with larger immunoglobulin G molecules (Corti et al., 2021) . Nanobodies are also highly soluble, very stable, lack glycans, and are readily cloned and produced in bacteria or yeast (Muyldermans, 2013) . They have low immunogenicity (Revets et al., 2005; Jovčevska and Muyldermans, 2020; Bannas et al., 2017) and can be readily 'humanized' (including Fc addition), modified to alter clearance rates, derivatized, combined for synergistic activity, and multimerized to improve characteristics (Chanier and Chames, 2019; Vincke et al., 2009; Duggan, 2018) . In the case of respiratory viruses like SARS-CoV-2, nanobodies' flexibility in drug delivery is a critical advantage. Beyond typical administration methods, a major advantage of nanobodies is their potential for direct delivery by nebulization deep into the lungs (Wölfel et al., 2020; Nambulli et al., 2021) . This route can provide a high local concentration in the airways and lungs to ensure rapid onset of therapeutic effects, while limiting the potential for unwanted systemic effects (Erreni et al., 2020) as exemplified by clinical trials (Van Heeke et al., 2017; Zare et al., 2021) . Moreover, with respect to deployment, nanobodies are relatively inexpensive and easy to reproducibly manufacture, with long shelf-lives and greater inherent stability compared to other biologicals, including monoclonals. Taken together nanobodies have great potential for the development of superior and differentiated therapeutics that would not only serve critically ill hospitalized patients, but also are especially well suited to the developing countries, most of which lack a reliable supply chain, or to stockpiling. To date, there are 453 nanobodies available against SARS-CoV-2 spike and those that are available primarily recognize regions of RBD with many subject to escape variation (Niu et al., 2021; Raybould et al., 2021; Schoof et al., 2020; Xiang et al., 2020; Koenig et al., 2021; Huo et al., 2020a; Pymm et al., 2021; Hanke et al., 2020; Custódio et al., 2020; Esparza et al., 2020; Wrapp, 2020a; Dong et al., 2020; Ye et al., 2021) . To address the urgent need for strongly neutralizing and escape resistant nanobodies, we generated a large repertoire of nanobodies that exploit the available epitope and vulnerability landscape of SARS-CoV-2 spike protein. The resulting repertoire provides a plethora of synergistically potent and escape resistant therapeutics. Maximizing the size and diversity of anti-SARS-CoV-2 spike nanobody repertoire We sought to isolate a large repertoire of highly diverse nanobodies against SARS-CoV-2 spike protein. Thus, we built on our existing nanobody generation pipeline (Fridy et al., 2014) , further optimizing each step, explicitly designing it to yield hundreds of high-quality, highly diverse nanobody candidates ( Figure 1A) . In this way, we took advantage both of the straightforward procedure of llama immunization and the powerful natural affinity maturation processes in vivo (Thompson et al., 2016) . To identify V H H domains that bind spike, we affinity-purified V H H domains from the immunized animals' sera against spike S1, S2, or RBD domains using independent domains in this purification step to maximize epitope accessibility. In parallel, lymphocyte RNA was taken from bone marrow aspirates and used to amplify V H H domain sequences by PCR, which were sequenced to generate an in silico library representative of all V H H sequences expressed in the individual animal. The affinity-purified V H H fragments were proteolyzed and the resulting peptides analyzed by LC-MS/MS. These data were searched against the V H H sequence library to identify and rank candidate nanobody sequences using our Llama-Magic software package (Fridy et al., 2014) with a series of key improvements (see Materials and methods). To maximize sequence diversity and thus the paratope space being explored, we clustered CDR sequences, revealing that many of the candidates form clusters likely to have similar antigen binding behavior. Here, partitioning of the clusters was performed by requiring that CDR3s in distinct clusters differ by a distance of more than three Damerau-Levenshtein edit operations (Bard, 2007) -that is, each operation being defined by insertion, deletion, or substitution of an amino acid residue, or transposition of two adjacent amino acid residues ( Figure 1B) . This partitioning was found to be effective, in that virtually no overlap was observed between those directed against S1 versus S2 (4 out of 183 clusters show overlap). The lengths of these CDR3 candidates also varied considerably, ranging from 3 to 22 amino acids in length. The use of two animals further expanded the paratope diversity in that only 4 out of 22 possible clusters from the second animal were observed to be shared with the first animal. In addition, we detected relatively little overlap between our CDR3 clusters and those observed by other groups; for example, only 1 out of 109 S1-specific clusters (Damerau-Levenshtein ≤3) were shared by Xiang et al., 2020 and the present work, indicating that our repertoires sampled extended regions of the available paratope space (see also below). Of the several hundred positives, 180 high-confidence candidates were selected for expression and screening. Of these, 66 were from S1 affinity purification, 63 from S2, and 51 from RBD, numbered S1-n, S2-n, and S1-RBD-n, respectively. These were then expressed with periplasmic secretion in bacteria, and crude periplasmic fractions were bound in large excess to the corresponding immobilized spike antigen to assay recombinant expression, specific binding, and degree of binding ( Figure 2-figure supplements 1 and 2) . 138 candidates were validated by this screen: 52 against S1-RBD-11 S1-RBD-4 S1-RBD-23 S2-10 S1-4 LaM2 Integrative structural modeling to map Nb binding sites [S1-RBD-15] (µM) 100 S1-1 S1-3 S1-6 S1-7 S1-10 S1-11 S1-12 S1-14 S1-15 S1-17 S1-19 S1-21 S1-22 S1-24 S1-9 S1-18 S1-25 S1-26 S1-27 S1-29 S1-31 S1-32 S1-33 S1-8 S1-28 S1-36 S1-38 S1-40 S1-30 S1-41 S1-43 S1-44 S1-45 S1-48 S1-49 S1-2 S1-50 S1-51 S2-2 S2-53 S1-RBD-1 S1-RBD-2 S1-RBD-3 S1-RBD-7 S1-RBD-9 S1-RBD-10 S1-RBD-12 S1-RBD-13 S1-RBD-14 S1-39 S1-RBD-15 S1-RBD-16 S1-RBD-19 S1-RBD-21 S1-4 S1-RBD-22 S1-RBD-24 S1-RBD-28 S1-RBD-29 S1-35 S1-RBD-5 S1-RBD-30 S1-RBD-6 S1-RBD-31 S1-13 S1-RBD-25 S1-RBD-32 S1-RBD-33 S1-RBD-8 S1-RBD-11 S1-RBD-35 S1-RBD-4 S1-RBD-46 S1-23 S1-RBD-51 S1-52 S1-53 S1-54 S1-55 S1-58 S1-59 S1-60 S1-61 S1-62 S2-56 S2-57 S2-59 S2-60 S2-61 S2-62 S1-66 S1-64 S1-65 S1-5 S1-16 S1-20 S1-37 S1-34 S1-42 S1-46 S1-47 S2-52 S1-RBD-17 S2-47 S1-RBD-37 S1-RBD-38 S1-RBD-40 S1-RBD-44 S1-RBD-45 S1-RBD-47 S1-RBD-18 S1-RBD-20 S1-RBD-23 S1-RBD-26 S1-RBD-27 S1-RBD-36 S1-RBD-34 S1-RBD-41 S1-RBD-39 S1-RBD-43 S1-RBD-48 S1-RBD-49 S1-RBD-50 S1-RBD-80 S1-56 S1-57 S1-63 Schematic of our strategy for generating, identifying, and characterizing large, diverse repertoires of nanobodies that bind the spike protein of SARS-CoV-2. The highest quality nanobodies were assayed for their ability to neutralize SARS-CoV-2 pseudovirus, SARS-CoV-2 virus, and viral entry into primary human airway epithelial cells. We also measured the activities of homodimers/homotrimers and mixtures. (B) A network visualization of 374 high-confidence CDR3 sequences identified from the mass spectrometry workflow. Nodes (CDR3 sequences) were connected by edges defined by a Damerau-Levenshtein distance of no more than 3, forming 183 isolated components. A thicker edge indicates a smaller distance value, that is, a closer relation. (C) Dendrogram showing sequence relationships between the 116 selected nanobodies, demonstrating that the repertoire generally retains significant diversity in both anti-S1 (green) and anti-S2 (blue) nanobodies, albeit with a few closely related members. Scale, 0.2 substitutions per residue. S1, 42 against S2, and 44 against RBD ( Figure 1B) . To eliminate candidates with weaker expression and binding affinity, only nanobodies in lysates with binding intensity >20% of the observed maximum across all those screened were chosen for follow-up study. This filtering identified the top 116 nanobodies that were purified for further characterization (Tables 1 and 2) . Note that these selections were designed to provide a strict cutoff in the interests of maximizing the quality of the repertoire selected for thorough characterization, but eliminated many additional nanobodies that nevertheless specifically bind to S1 and S2. While a few of these 116 nanobodies were chosen to share similar paratopes, overall, the group retained a high sequence and paratope diversity ( Figure 1C ). High-affinity nanobodies across the entire spike ectodomain that are refractory to common spike escape mutants Surface plasmon resonance (SPR) was used to detail the kinetic properties and affinities of the selected nanobodies ( Table 1 and 2). All bound with high affinity, with >60% binding with K D s < 1 nM, and two with single-digit picomolar affinities ( Figure 2) . While most S1-binding nanobodies bind RBD (71 nanobodies), 19 targeted non-RBD regions of S1 and 26 bind S2 ( Figure 2) . The lower number of non-RBD S1 and S2 nanobodies likely reflects the highly antigenic nature of the RBD and the occlusion of non-RBD S1 regions and S2 due to the glycan shield of SARS-CoV-2 spike (Grant et al., 2020; Watanabe et al., 2020) . At the same time, we observed no obvious bias in nanobody affinities for these different domains. While both high on rates and low off rates contributed to these high affinities, kinetic analyses underscore the generally fast association rates (many with k on ≥ 10 +6 ) of these nanobodies (likely due to their small size and proportionally large paratope surface area), with many surpassing the k on of high-performing monoclonal antibodies (k on ~ 10 +5 ) Figure 2 ), a property that would benefit translation of these nanobodies into rapid therapeutics and diagnostics (Carter, 2006) . For those nanobodies with apparently homologous paratopes ( Figure 1C ), we found no correlation in their kinetic properties (Tables 1 and 2), demonstrating that even small paratope changes can strongly alter behaviors (Fridy et al., 2014) . A worrying development is the continuing emergence of viral variants, including mutations in RBD that minimize or nullify binding of many currently available monoclonal antibodies and nanobodies, which solely target RBD Wang et al., 2021b; Diamond et al., 2021; Jangra et al., 2021; Garcia-Beltran et al., 2021; Liu et al., 2021b; Sun et al., 2021) . Indeed, in one study, the efficacy of 14 out of the 17 most potent monoclonal antibodies tested was compromised by such common RBD mutants (Wang et al., 2021b) . Here, based on the large size of our repertoire and its extensive binding across the available epitope space of spike, nanobodies or combinations thereof show great potential to be particularly resistant to these variants . RBD mutants represent a significant class of escape variants (Garcia-Beltran et al., 2021; Greaney et al., 2021) , leading us to employ two strategies to ensure the generation of numerous nanobodies whose binding (and virus-neutralizing activities) are resistant to emerging variants. First, we isolated a large diversity of high-quality anti-RBD nanobodies to maximize the probability of identifying ones that are refractory to escape. Second, to reveal additional nanobody-neutralizing potential, we deliberately targeted non-RBD regions of spike (see below) (Elshabrawy et al., 2012; Greaney et al., 2021) . To test the first strategy, we sampled RBD-binding nanobodies covering non-overlapping epitopes on RBD ( Figure 3 ) and examined their binding to SARS-CoV-2 variants B.1.1.7/20I/501Y.V1/alpha (United Kingdom) and B.1.351/20H/501Y.V2/beta (South Africa) (Wang et al., 2021a; Ho et al., 2021; Figure 2, Table 3) . Of the seven nanobodies tested, six of these (S1-1, S1-6, S1-RBD-9, S1-RBD-11, S1-RBD-15, and S1-RBD-35) retained their very strong binding to both variants, with only a modest reduction in affinity for S1-RBD-11 binding to variant B.1.351/20H/501Y.V2/beta (20-161 pM). For the seventh nanobody, S1-23, binding to variant B.1.1.7/20I/501Y.V1/alpha was only reduced from a K D of 17 pM to a still-respectable 230 pM, although its binding to variant B.1.351/20H/501Y.V2/beta was abolished ( Figure 2 ). As expected (VanCott et al., 1994; Magnus, 2013; Steckbeck et al., 2005) , it is the off rates that are most affected by these variants. Nevertheless, based on epitope mapping (below) and our identification of nanobodies that recognize epitopes not altered in the emerging variant strains, we expect that a high percentage of our nanobodies will remain resistant to these escape mutants; this would now include the B.1.617.2/21A/delta variant , making our collection a powerful resource for potential prophylactics and therapeutics. Nanobodies against S1 were determined to bind RBD or non-RBD epitopes by their affinity for recombinant full-length S1 and/or S1 RBD protein. Binding kinetics against these two recombinant proteins were determined by surface plasmon resonance (SPR), with on rates, off rates, and K D s determined by Langmuir fits to binding sensorgrams unless otherwise noted. Nanobody melting temperatures (T m ) were determined by differential scanning fluorimetry (DSF). Nanobodies were assayed for neutralization activity against a SARS-CoV-2 spike pseudotyped HIV-1 virus (PSV), with IC50s calculated from neutralization curves. Standard error of the mean (s.e.m.) is reported when available. The nanobody repertoire has favorable stability properties A key consideration for possible biological therapeutics and diagnostics for SARS-CoV-2 is their stability under potentially denaturing conditions (McConnell et al., 2014) . To address this, we performed differential scanning fluorimetry (DSF) experiments to determine the thermal stability (T m ) of each of our nanobodies. These studies revealed a thermal stability range between 50 and 80°C, similar to published results of other properly folded nanobodies and indicative of their generally high stability (Muyldermans, 2013) . In contrast to many conventional antibodies, nanobodies are also reported The data points highlighted in blue correspond to nanobodies that neutralize. The majority of nanobodies have high affinity for their antigen with K D s below 1 nm. 10 nanobodies are not included in this plot as they were unable to be analyzed successfully using surface plasmon resonance (SPR). (B) SPR sensorgrams for each of the three targets on SARS-CoV-2 spike protein of our nanobody repertoire, showing three representatives for each binding region. (C) The association rate of each nanobody (k on ) versus the corresponding dissociation rate (k off ). The majority of our nanobodies have fast association rates (~10 +5 -10 +7 M -1 s -1 ), with many surpassing the k on of Figure 2 continued on next page to remain fully active upon reconstitution after lyophilization, particularly in buffers lacking cryoprotectants (Schoof et al., 2020; Xiang et al., 2020) . A representative sample from our repertoire was thus freeze-dried without cryoprotectants, reconstituted, then analyzed via SPR and DSF to determine whether their properties were compromised due to lyophilization. The results revealed no significant effect on stability, kinetics, and affinity ( Figure 2E and F). Taken together, these data suggest that our nanobodies, like those published in other contexts (Xiang et al., 2020; Schoof et al., 2020) , are able to withstand various temperatures and storage conditions without affecting their stability and binding. These are essential requirements for downstream applications (e.g., use in a nebulizer) and ease of storage -important considerations if these are to be used for mass distribution, including in resourcepoor settings (Peeling and McNerney, 2014) . We applied a multifaceted approach to physically distinguish nanobodies that target common regions on the surface of the RBD. Using an eight-channel biolayer interferometer, we tested for pairwise competitive binding of nanobodies that bind the RBD, as well as for those that bind outside of the RBD (i.e., within the S1 non-RBD and S2 domains) ( Figure 3 ). Label-free binding of antibodies to antigens measured in a 'dip-and-read' mode provides a real-time analysis of affinity and the kinetics of the competitive binding of nanobody pairs and can distinguish between those that bind to similar or overlapping epitopes versus distinct, non-overlapping epitopes (Estep et al., 2013) . 56 anti-RBD nanobodies were screened in pairwise combinations. The response values were used to assist the discovery of nanobody groups that most likely bind non-overlapping epitopes by ensuring that the least response of pairwise nanobodies within the group was maximized. Eleven representative anti-RBD nanobodies were used as a foundation, selecting two or more representative nanobodies from each group to bin the remaining RBD nanobodies in our collection. Overlapping pairs from the foundation group and the remaining RBD binders were used to measure if a nanobody pair behaved similarly against other nanobodies measured in the dataset ( Figure 3A ), to comprehensively map nanobody competition and epitope bins ( Figure 3D ). Pearson's correlation coefficients were derived based on their binding characteristics, and the data were used to hierarchically cluster and group all RBD binders into bins. This approach revealed three large, mostly non-overlapping bins. However, each bin contained smaller, better-correlated clusters of nanobodies, reflected by the dendrogram, indicating the presence of numerous distinct sub-epitope bins present within each larger bin, that is, discrete epitopes that partially overlap with other discrete epitopes in the same bin. We calculated the gap statistic (Tibshirani et al., 2001) , to estimate the optimal cluster number, discerning at least eight epitope bins ( Figure 3A ). Nanobodies binding to regions outside of the RBD of S1 were binned in a similar fashion Figure 3B , C, E, and F. Using SPR, we binned 16 non-RBD S1-binding nanobodies and 19 S2-binding nanobodies in pairwise competition assays ( Figure 3B and C). Pearson's correlation coefficients were used to hierarchically cluster these nanobodies, revealing as many as four S1-non-RBD bins and five S2 bins. The binning data from pairwise combinations suggest numerous epitope bins, and thus it is reasonable to hypothesize that more than two nanobodies can bind a single domain at the same time. To test this hypothesis, we used mass photometry (MP) (Soltermann et al., 2020; Wu and Piszczek, 2021; Young et al., 2018) , which can accurately measure multiple binding events to a single antigen. This allowed us to determine which nanobodies share epitope space on spike S1 monomer through detection of additive mass accumulation of a nanobody (or nanobodies) on spike S1 depending on high-performing monoclonal antibodies (~10 +4 -10 +5 M -1 s -1 ). (D) Each nanobody plotted against their T m as measured by differential scanning fluorimetry (DSF), revealing all but two nanobodies fall within a T m range between 50 and 80°C, where the bulk of our nanobodies have a T m ≥ 60°C. No data could be collected for two nanobodies, and 10 nanobodies exhibited two dominant peaks in the thermal shift assay and were not included in this plot (a full summary of this data can be seen in Tables 1-3). The K D (E) and T m (F) of six nanobodies were assessed pre-and post-freeze-drying, revealing no significant change in affinity or T m after freeze-drying. (G) SPR sensorgrams comparing the kinetic and affinity analysis of seven nanobodies against wildtype spike S1 (Wuhan strain), spike 20I/S1 501Y.V1 (alpha variant), and 20H/spike S1 501Y.V2 (beta variant). The online version of this article includes the following figure supplement(s) for figure 2: 0 .00 0 .50 -0.50 S1-6 S1-39 S1-35 S1-4 S1-RBD-22 S1-1 S1-RBD-12 S1-RBD-14 S1-RBD-29 S1-RBD-4 S1-RBD-36 S1-RBD-9 S1-RBD-16 S1-RBD-10 S1-RBD-24 S1-RBD-30 S1-RBD-25 S1-RBD-32 S1-RBD-34 S1-RBD-19 S1-RBD-26 S1-RBD-44 S1-RBD-48 S1-RBD-51 S1-RBD-47 S1-RBD-49 S1-RBD-45 S1-RBD-46 S1-RBD-41 S1-RBD-37 S1-RBD-43 S1-RBD-21 S1-RBD-27 S1-RBD-20 S1-RBD-40 S1-RBD-15 S1-RBD-28 S1-RBD-3 S1-RBD-5 S1-RBD-6 S1-31 S1-27 S1-36 S1-23 S1-62 S1-RBD-18 S1-28 S1-20 S1-RBD-35 S1-27 S1-48 S1-29 S1-46 S1-RBD-23 S1-RBD-38 S1-RBD-11 S1-6 S1-39 S1-35 S1-4 S1-RBD-22 S1-1 S1-RBD-12 S1-RBD-14 S1-RBD-29 S1-RBD-4 S1-RBD-36 S1-RBD-9 S1-RBD-16 S1-RBD-10 S1-RBD-24 S1-RBD-30 S1-RBD-25 S1-RBD-32 S1-RBD-34 S1-RBD-19 S1-RBD-26 S1-RBD-44 S1-RBD-48 S1-RBD-51 S1-RBD-47 S1-RBD-49 S1-RBD-45 S1-RBD-46 S1-RBD-41 S1-RBD-37 S1-RBD-43 S1-RBD-21 S1-RBD-27 S1-RBD-20 S1-RBD-40 S1-RBD-15 S1-RBD-28 S1-RBD-3 S1-RBD-5 S1-RBD-6 S1-31 S1-27 S1-36 S1-23 S1-62 S1-RBD-18 S1-28 S1-20 S1-RBD-35 S1-27 S1-48 S1-29 S1-46 S1-RBD-23 S1-RBD-38 S1-RBD-11 S1 RBD epitope bins S1-25 S1-9 S1-10 S1-11 S1-58 S1-24 S1-3 S1-30 S1-41 S1-60 S1-2 S1-49 S1-17 S1-7 S1-32 S1-50 S1-25 S1-9 S1-10 S1-11 S1-58 S1-24 S1-3 S1-30 S1-41 S1-60 S1-2 S1-49 S1-17 S1-7 S1-32 S1-50 1.00 0.90 0.80 0.70 0.60 B S1 non-RBD epitope bins Edge opacity is proportional to the response value of the nanobody pair. A lighter color indicates a higher response value from the pair. Average shortest path length L = 1.64 Average clustering coefficient C = 0.831 Small-world coefficient σ = 1.031 S1 RBD clustering analysis S1-RBD-10 S1-RBD-11 S1-RBD-12 S1-RBD-14 S1-RBD-15 S1-RBD-16 S1-RBD-18 S1-RBD-19 S1-RBD-20 S1-RBD-21 S1-RBD-22 S1-RBD-23 S1-RBD-24 S1-RBD-25 S1-RBD-26 S1-RBD-27 S1-RBD-28 S1-RBD-29 S1-RBD-3 S1-RBD-30 S1-RBD-32 S1-RBD-34 S1-RBD-35 S1-RBD-36 S1-RBD-37 S1-RBD-38 S1-RBD-4 S1-RBD-40 S1-RBD-41 S1-RBD-43 S1-RBD-44 S1-RBD-45 S1-RBD-46 S1-RBD-47 S1-RBD-48 S1-RBD-49 S1-RBD-5 S1-RBD-51 S1-RBD-6 S1-RBD-9 S1-1 S1-20 S1-23 S1-27 S1-28 S1-29 S1-31 S1-35 S1-36 S1-37 S1-39 S1-4 S1-46 S1-48 S1-6 S1-62 E S1 non-RBD clustering analysis +S1-1 + S1-23 +S1-1 + S1-RBD-9 +S1-23 + S1-RBD-9 +S1-1 + S1-23 +S1-1 + S1-46 +S1-23 + S1-46 +S1-1 + S1-23 + S1-RBD-9 +S1-1 + S1-23 + S1-46 500 S1-10 S1-11 S1-17 S1-2 S1-24 S1-25 S1-3 S1-30 S1-32 S1-41 S1-49 S1-50 S1-58 S1-60 S1-7 S1-9 whether or not nanobodies share epitope space on spike S1. Several representative nanobodies that indicate that a nanobody pair responds divergently when measured against nanobodies in the representative panel and indicate binding to distinct or non-overlapping regions on the RBD. (B) As in (A), but for 16 S1 non-RBD-binding nanobodies. (C) As in (A), but for 19 S2-binding nanobodies. (D) A network visualization of anti-S1-RBD nanobodies. Each node is a nanobody and each edge is a response value measured by biolayer interferometry from pairwise cross-competition assays. Orange nodes represent 11 nanobodies used as a representative panel for clustering analysis in (A). Blue nodes represent the other nanobodies in the dataset. The average shortest distance between any nanobody pair in the dataset is 1.64. An average clustering coefficient of 0.831 suggests that the measurements are well distributed across the dataset. The small world coefficient of 1.031 indicates that the network is more connected than to be expected from random, but the average path length is what you would expect from a random network, together indicating that the relationship between nanobody pairs not actually measured can be inferred from the similar/neighboring nanobodies. (E, F) As in (D) but for S1 non-RBD and S2 nanobodies, respectively. These are complete networks with every nanobody measured against the others in the dataset. (G) Mass photometry (MP) analysis of spike S1 monomer incubated with different anti-spike S1 nanobodies. Two examples of an increase in mass as spike S1 monomers (black line) are incubated with 1-3 nanobodies. The accumulation in mass upon addition of each different nanobody on spike S1 monomer is due to each nanobody binding to non-overlapping space on spike S1, an observation consistent with Octet binning data. As a control, using MP, each individual nanobody was shown to bind spike S1 monomers on its own (data not shown). The online version of this article includes the following figure supplement(s) for figure 3: Table 3 . Nanobody binding activity against spike S1 variants; related to Figure 2 . Binding kinetics against wild-type spike S1 or two variants of concern were determined by surface plasmon resonance (SPR), with on rates, off rates, and K D s determined by Langmuir fits to binding sensorgrams. sample across the epitope space of our nanobody repertoire that bind the RBD were chosen for MP studies based on the epitope binning data ( Figure 3G ). These data confirmed the separation of our major epitope bins, and furthermore demonstrated that we can bind at least three different nanobodies simultaneously to the RBD, contrasting with the much larger conventional immunoglobulins, which may be too large to simultaneously bind either monomer or trimer S protein (Corti et al., 2021; Stewart et al., 1997; Xu et al., 2021) . This is a critical consideration for the design of complementary nanobody cocktails and multimers with synergistic-neutralizing activities (see below). We used a SARS-CoV-2 pseudovirus neutralization assay to screen and characterize our nanobody repertoire for antiviral activities ( Figure 4) . The lentiviral-based, single-round infection assay robustly measures the neutralization potential of a candidate nanobody and is a validated surrogate for replication competent SARS-CoV-2 (Riepler et al., 2020; Schmidt et al., 2020) . Because measured IC50s are dependent on assay conditions and so cannot be readily compared across laboratories (Cheng and Prusoff, 1973) , we included four other published nanobodies in this assay for comparison (Xiang et al., 2020; Wrapp, 2020a ; Figure 4G ). Overall, 36% of our monomeric nanobody repertoire neutralized with IC50s ≤ 100 nM, while 23% showed neutralization with IC50s < 50 nM and 17 potent neutralizers at 20 nM or lower ( Figure 4A) . Similarly, the four published nanobodies span the range of neutralization observed within our repertoire from potent (<20 nM) to relatively weak (between 1 and 10 µM). As a further comparison and validation of our IC50 values, we evaluated a subset of our nanobodies in a complementary neutralization assay , which revealed a strong correlation between these two assays with a Pearson's correlation coefficient of 0.98 and p-value < 0.0001 ( Figure 4 -figure supplement 1). Our most potent neutralizing nanobodies mapped to the RBD; neutralizing activity mapped to each of the major epitope bins of the RBD and were of similar efficacy to the most potent of the comparison nanobodies; importantly, nanobodies binding outside of the RBD also possess neutralizing activity (for example S1-64 and S1-65). Notably, nanobodies mapping outside of the RBD on S1 (anti-S1, non-RBD) and to S2 also neutralized the pseudovirus in our assay, albeit with somewhat higher IC50s ( Figure 4B and C). This is the first evidence of nanobody neutralization activity mapping outside of the RBD. As nanobodies are monomeric, the mechanism of this neutralization does not involve viral aggregation and likely reflects disruption of the virus binding or spike-driven fusion of viral and cellular membranes. Nanobodies, especially directed against relatively invariant regions of coronavirus spike proteins, may have broadly binding/neutralizing activities and are therefore important targets for optimization. Such optimization includes their use in cocktails and as oligomers. Oligomerization strongly enhances the affinity and neutralization activity of nanobodies A distinct advantage of nanobodies is the facility by which oligomers can be produced to improve their affinities and avidities (Fridy et al., 2014; Schoof et al., 2020; Xiang et al., 2020; Koenig et al., 2021) . Oligomerization of most of our nanobodies significantly improved their IC50s and measured affinities ( Table 4) . For example, dimerization and trimerization of S1-RBD-35 improved neutralization activity from IC50s of ~12 nM to ~150 pM and ~70 pM, respectively. Similar results were found with S1-23, improving neutralization from ~6 nM to ~220 pM and ~90 pM, respectively ( Figure 4D ). Dimerization of the anti-S1 non-RBD nanobody S1-49 improved IC50s from ~350 nM to ~9 nM, and trimerization improved its activity an additional ~10-fold. Multimerization of some nanobodies directed against regions outside of the RBD on both S1 and S2 led to nanomolar range IC50s ( Figure 4E and F). This includes S2-7, for which dimerization converted a nanobody that we considered to be a nonneutralizer to one having a respectable neutralizing activity (IC50 ~ 250 nM), with further potency achieved by trimerization and tetramerization, down to an IC50 of ~30 nM ( Figure 4F ). While these results show that multimerization can dramatically improve their activities, importantly this was not always the case ( Table 4 ), indicating that enhancement by multimerization is not a given, but must be determined empirically. % inhibition of pseudo-SARS-CoV-2 % inhibition of pseudo-SARS-CoV-2 % inhibition of pseudo-SARS-CoV-2 10 -7 10 -6 10 -5 10 -4 10 -3 10 -2 10 -1 10 0 10 1 10 2 0 50 100 (µM) [nanobody] 10 -7 10 -6 10 -5 10 -4 10 -3 10 -2 10 -1 10 0 10 1 10 2 0 50 100 10 -7 10 -6 10 -5 10 -4 10 -3 10 -2 10 -1 10 0 10 1 10 2 0 50 100 (µM) [nanobody] (µM) [nanobody] anti-S1 RBD with IC 50 <20 nM anti-S1 non-RBD a nti-S2 % inhibition of pseudo-SARS-CoV-2 % inhibition of pseudo-SARS-CoV-2 anti-S1 RBD oligomers 0 50 100 10 -7 10 -6 10 -5 10 -4 10 -3 10 -2 10 -1 10 0 10 1 10 2 % inhibition of pseudo-SARS-CoV-2 anti-S1 non-RBD oligomers 0 50 100 anti-S2 oligomers 10 -7 10 -6 10 -5 10 -4 10 -3 10 -2 10 -1 10 0 10 1 10 2 • S1-RBD-15 • S1-62 • S1-48 • S1-1 • S1-23 • S1-RBD-23 • S1-37 • S1-RBD-29 • S1-RBD-35 • S1-35 • S1-RBD-20 • S1-RBD-11 • S1-RBD-21 • S1-RBD-4 • S1-RBD-41 • S1-27 • S1-65 • S1-64 • S1-58 • S1-3 • S1-49 • S1-30 • S1-24 • S1-51 • S1-41 • S1-50 IC 50 10 pM 100 pM 1 nM 10 nM 100 nM 1 µM 10 µM 100 µM anti-S1 non-RBD anti-S2 anti-S1 RBD benchmarks anti-S1 RBD (Xiang et al., 2020 ) 10 -7 10 -6 10 -5 10 -4 10 -3 10 -2 10 -1 10 0 10 1 10 2 0 50 100 S1-49 1× 2× 3× monomers oligomers 0 50 100 0 50 100 S1-62 S1-1 S1-RBD-15 S1-48 S1-RBD-35 S1-RBD-21 S1-RBD-9 S1-RBD-16 S1-23 % inhibition of pseudo-SARS-CoV-2 10 -7 10 -6 10 -5 10 -4 10 -3 10 -2 10 -1 10 0 10 1 10 2 (µM) [nanobody] 0 50 100 % inhibition of pseudo-SARS-CoV-2 10 -7 10 -6 10 -5 10 -4 10 -3 10 -2 10 -1 10 0 10 1 10 2 (µM) [nanobody] anti-S1 RBD on 20H/501Y.V2 (beta) I S1-1 S1-RBD-15 S1-48 S1-RBD-9 S1-RBD-35 S1-RBD-21 S1-RBD-16 S1-23 S1-62 % inhibition of pseudo-SARS-CoV-2 10 -7 10 -6 10 -5 10 -4 10 -3 10 -2 10 -1 10 0 10 1 10 2 (µM) [nanobody] anti-S1 RBD on 20J/501Y.V3 (gamma) S1-RBD-15 S1-1 S1-RBD-16 S1-RBD-35 S1-RBD-9 S1-48 S1-RBD-21 S1-23 S1-62 K anti-S1 RBD on SARS-CoV-1 anti-S1 RBD on 20I/501Y.V1 (alpha) anti-S1 RBD on 20J/501Y.V3 (gamma) anti-S1 RBD on 20H/501Y.V2 (beta) 10 -7 10 -6 10 -5 10 -4 10 -3 10 -2 10 -1 10 0 10 1 10 2 (µM) • S1-1 • S1-39 • S1-RBD-6 • S1-6 • S1-RBD-9 • S1-4 • S1-35 • S1-37 • S1-23 • S1-48 • S1-RBD-21 • S1-RBD-47 Nanobodies neutralize SARS-CoV-2 variants and SARS-CoV-1 Certain mutations in spike, appearing in 'variants of concern' (VOC) associated with rapidly increasing case numbers in certain locales, have been demonstrated to reduce the neutralization potency of some monoclonal antibodies and polyclonal plasma, increase the frequency of serious illness, and are spreading rapidly (Wang et al., 2021b; Wibmer et al., 2021) . We therefore tested a subset of our neutralizing nanobodies against pseudovirus carrying the spike protein of the alpha variant (B.1.1.7 /20I/501Y.V1) ( Figure 4H ); the beta variant (B.1.351/20H/501Y.V2) (Tegally et al., 2021; Stamatatos et al., 2021; Figure 4I) ; and the gamma variant (P.1/20J/501Y.V3) ( Figure 4J , Table 5 ). These VOCs have mutations resulting in amino acid substitutions in spike, which impact the neutralization efficacy for some of the tested nanobodies. While all nanobodies neutralize the alpha variant, S1-23 showed an almost 14-fold drop in potency ( Figure 4A and I). S1-23 and S1-62 failed to neutralize the beta and gamma variants, while S1-1 and S1-RBD-15 were as efficacious against all three VOCs as they were against wild-type spike ( Figure 4I ). Both S1-RBD-21 and -35 also remained effective neutralizers of spike VOC pseudotypes, albeit with reduced IC50s compared to the wild-type spike ( Figure 4I ). Remarkably, S1-RBD-9 showed increased neutralization activity against all three VOC, improving ~2-fold against alpha, ~6-fold against beta, and ~10-fold against gamma ( Figure 4H -J). These results are in accord with our SPR studies, which showed that S1-1, S1-RBD-9, -15, and -35 retained very strong binding to the alpha and beta VOC, whereas binding of S1-23 was reduced against alpha and completely abolished against beta ( Figure 2G ). The B.1.617.2/21A/delta VOC has L452R and T478K as unique RBD amino acid substitutions , which, based on our epitope binning and escape mutants (below), we would predict to impact neutralization of S1-RBD-11 and S1-RBD-35 (T478K) or S1-RBD-23 and S1-36 (L452F). However, the great majority of nanobodies in our repertoire would be predicted to show similar high neutralization efficacy against all these VOCs as compared to the wild-type virus. Overall, these data indicate that comprehensive mining of our repertoire and multimerization can lead to nanobody-based therapies that remain fully effective against common and even potentially yet-to-emerge variants of SARS-CoV-2 and with broad-spectrum coronavirus inhibition activities. Both SARS-CoV-1 and SARS-CoV-2 share the same host receptor, ACE2, and the RBDs of the viruses share ~74% identity. As a result, some antibodies and nanobodies have been shown to be cross-neutralizing (Liu et al., 2021b; Wrapp, 2020a) . We therefore tested the ability of our nanobodies to neutralize SARS-CoV-1 in the pseudovirus assay. Of the nanobodies tested in this assay, several (7 of 27 tested) of our anti-RBD monomer nanobodies also displayed excellent neutralizing activities against SARS-CoV-1 spike pseudotyped virus ( Figure 4K ). S1-1, S1-39 and S1-RBD-6 had similar IC50s against both pseudotypes, while S1-35 and S1-6 showed reduced activity against SARS-CoV-1 pseudotypes. Notably, S1-23, S1-37, and S1-48 showed no activity against SARS-CoV-1 spike pseudotypes. These three nanobodies are highly correlated with one another in our epitope binning analysis, indicating their binding to proximal epitopes on the RBD ( Figure 3A ). Beyond nanobodies that bind to the RBD, 3 of 11 nanobodies that bind to non-RBD regions of S1 and S2 also neutralized SARS-CoV-1 spike pseudotypes (Table 5) . non-RBD, and S2 nanobodies significantly increases neutralization potency. (G) Summary scatter plot of all nanobody IC50s across the major domains of SARS-CoV-2 spike and where tested, across SARS-CoV-2 variant 20H/501Y.V2 and SARS-CoV-1. Representative published nanobodies were also tested in our neutralization assays and show similar potency towards SARS-CoV-2 pseudovirus. From Xiang et al., 2020: (1) Nb-21 (IC50 2.4 nM); (2) Nb-34 (IC50 5.6 nM); and (3) Nb-93 (IC50 123 nM). From Wrapp, 2020a: (4) VHH-72 (IC50 2.5 μM). (H-K) Representative SARS-CoV-2 RBD targeting nanobodies cross-neutralize the 20I/501Y.V1/alpha variant with H69-, V70-, Y144-amino acid deletions and N501Y, A570D, D614G, P681H, T716I, S982A, and D1118H amino acid substitutions in spike (H); the 20H/501Y.V2/beta variant with L18F, D80A, K417N, E484K, and N501Y amino acid substitutions in spike; (I) 20J/501Y.V3/gamma variant with L18F, T20N, P26S, D138Y, R190S, K417T, E484K, N501Y, D614G, H655Y, T1027I, and V1176F amino acid substitutions in spike (J); and SARS-CoV-1 spike pseudotyped lentivirus (K). In all cases, n ≥ 2 biological replicates of each nanobody monomer/oligomer with a representative biological replicate with n = 4 technical replicates per dilution are displayed. See also Table 1, Table 2, Table 4 , and Table 5 . The online version of this article includes the following figure supplement(s) for figure 4: Nanobodies effectively neutralize SARS-CoV-2 infection in human primary airway epithelium Nanobody and antibody neutralizations have been reported to yield similar results when performed with pseudovirus versus authentic virus (Schoof et al., 2020; Xiang et al., 2020; Schmidt et al., Table 4 . Characterization of oligomerized spike nanobodies; related to Figure 4. Nanobody oligomers (1-4 nanobody repeats) were assayed for neutralization activity against a SARS-CoV-2 spike pseudotyped HIV-1 virus (PSV), with IC50s calculated from neutralization curves. Standard error of the mean (s.e.m.) is reported where replicates were available. Epitopes were determined by relative affinity for recombinant S1 or S1 RBD protein. ID Epitope SARS-CoV-2 PSV IC50 (s.e.m.) (nM) SARS-CoV-1 PSV IC50 (s.e.m.) (nM) S1-1 RBD 6.7 (1.0) 8.6 (7.2) S1-1 dimer RBD 4.9 (0.1) -S1-1 trimer RBD 5.7 (0.1) -S1-23 RBD 5.7 (2.2) -S1-23 dimer RBD 0.22 (0.05) NA S1-23 trimer RBD 0.089 (0.019) NA S1-RBD-35 RBD 12.3 (2.4) NA S1-RBD-35 dimer RBD 0.15 (0.11) -S1-RBD-35 trimer RBD 0.068 (0.043) NA S1-3 S1 non-RBD 1,030 (666) 4161 S1-3 dimer S1 non-RBD 429 513 S1-3 trimer S1 non-RBD 411 -S1-30 S1 non-RBD 717 (388) -S1-30 dimer S1 non-RBD 18.3 662 S1-30 trimer S1 non-RBD 77.5 (3.6) -S1-7 S1 non-RBD NA -S1-7 dimer S1 non-RBD NA -S1-7 trimer S1 non-RBD NA -S1-17 S1 non-RBD 1271 (888) -S1-17 dimer S1 non-RBD 2144 -S1-49 S1 non-RBD 356 (32.8) NA S1-49 dimer S1 non-RBD 9.1 (1.2) NA S1-49 trimer S1 non-RBD 0.87 (0.08) NA Table 5 continued on next page 2020). However, discrepancies have also been reported, particularly for antibodies targeting regions outside the RBD (Chi et al., 2020; Huo et al., 2020b) . We therefore selected a panel of exemplar (monomeric) nanobodies, which target the RBD and domains outside of the RBD, to test for neutralization with authentic SARS-CoV-2. All nanobodies tested that neutralized pseudovirus also showed potent neutralization by plaque and focus reduction assays and correlated well with our pseudovirus assays ( Figure 5A , Table 6 ). To mimic human infection, we exploited human air-liquid interface (ALI) cultures of primary airway epithelium as an ex vivo model system of viral infection (Barrow et al., 2021) . This system mimics the lung environment as it contains pseudostratified, ciliated, and mucous-secreting cells that express ACE2 and has several advantages over animal models including representing the relevant physiological site of initial SARS-CoV-2 infection in humans (and associated innate responses), while enabling experimental control over infection, nanobody delivery, and quantification of viral RNA at the site of infection. We thus tested a subset of our nanobodies for their ability to block SARS-CoV-2 infection and spread in this model ( Figure 5B) . We treated the air-exposed apical surface of the culture with serial dilutions of S1-1 and S1-23 and then challenged them with SARS-CoV-2 at an MOI of 0.5. ALI cultures were then treated with nanobodies at 24 hr intervals for an additional 3 days before harvesting the cells, extracting RNA, and measuring SARS-CoV-2 levels by qPCR ( Figure 5B ). S1-1 potently neutralized SARS-CoV-2 at each concentration tested while S1-23 inhibited SARS-CoV-2 in a dose-dependent manner ( Figure 5C ). The efficacy of the S1-23 nanobody was strongly enhanced when provided to cells as a trimer, potently inhibiting viral replication ( Figure 5C ). As an additional comparator and as a control, we determined the inhibition of replication upon addition of recombinant competitor, ACE2. Nanobodies inhibited at lower doses than recombinant ACE2, reflective of our measured low K D of nanobody interactions with spike (<1 nM) compared to a reported K D of 14.7 nM or greater for ACE2 with spike (Huang and Chai, 2020; Shang et al., 2020; Chan et al., 2020; Rogers et al., 2020; Liu et al., 2020b; Cao et al., 2020) . These data highlight the potential for nanobodies to function as single-agent therapies against COVID-19, with efficacies comparable to monoclonal immunoglobulins. With the emerging variants of concern, our goal is to develop nanobody multimers and cocktails that are maximally refractory to escape by such variants. To do so, we used a previously employed method that drives the selection of antibody-resistant populations of rVSV/SARS-CoV-2 chimeric virus harboring variants of spike and measured the ability of the chimeric virus to escape nanobody-mediated neutralization . This approach simultaneously maps the escape potential of spike and the epitopes responsible for neutralization by nanobody binding (Figure 6-figure supplement 1) , with the goal of discovering spike variants that resist the neutralizing activity of individual nanobodies. Based on this information, we could then predict pairs of nanobodies whose escape mutants do not map to the same region of spike, the combination of which would thus likely prevent escape. Specifically, we prepared large and diversified populations (10 6 infectious units) of a recombinant UV inactivated S1-23 trimer S1-23 S1-1 SARS-CoV-2 copy number / ng RNA C S1-1 S1-23 S1-RBD-11 S1-RBD-4 S1-RBD-23 S2-10 S1-4 LaM2 with the indicated concentrations of anti-SARS-CoV-2 spike nanobodies. The infected cultures were maintained for 5 days with a daily 1 hr incubation of nanobodies before being harvested for RNA isolation and determination of the SARS-CoV-2 copy number by qPCR. SARS-CoV-2 copy number was normalized to total RNA measured by spectrophotometry. Mock-treated samples exposed to infectious and UVinactivated SARS-CoV-2 virions served as positive and negative controls. Recombinant soluble angiotensin converting enzyme 2 (rACE2) was used as a positive treatment control. The indicated nanobodies were used at 1, 10, and 100× their IC50 values determined in pseudovirus neutralization assays (Table 1 and Table 4 ). Table 6 . Nanobody neutralization activity against SARS-CoV-2; related to Figure 5 . Nanobodies were assayed for neutralization activity against authentic SARS-CoV-2, with IC50s calculated from neutralization curves. ID Epitope SARS-CoV-2 IC50 (nM) S1-1 RBD 1.1 S1-4 RBD 1310 S1-23 RBD 0.7 S1-RBD-4 RBD 5.4 S1-RBD-11 RBD 3.0 S1-RBD-23 RBD 6.1 S2-10 Non-RBD 91.2 Non-target ctrl NA NA, no activity. diverse spike epitopes ( Table 7) . For some of the non-RBD epitope nanobodies, we used dimeric or trimeric forms of the nanobodies to further enhance their activity, but in each case the selected viral isolates exhibited resistance to monomeric, dimeric, or trimeric forms. While some of the mutations that arose in the selection experiments were likely passenger mutations (Table 7) , a number of mutations clustered on the spike surface close to each other on RBD ( Figure 6-figure supplement 1 ; Weisblum et al., 2020; Muecksch et al., 2021; Wang et al., 2021b) . Some of the most potently neutralizing nanobodies selected resistant mutations at the same positions (e.g., E484K) as those selected by potent neutralizing antibodies that have been cloned from SARS-CoV-2 convalescents and vaccine recipients, confirming that the ACE2 binding site is a point of particular vulnerability for potent neutralization. Additionally, however, other nanobodies selected mutations that have not previously been encountered in human antibody selection experiments ( Table 7) . Nanobody cocktails are expected to be resistant to escape as they recognize multiple epitopes (Baum et al., 2020; De Gasparo et al., 2021; Weisblum et al., 2020) . As proof of principle, we generated sets of two-nanobody cocktails by combining specific nanobodies that selected spatially distinct resistance mutations on the RBD ( Figure 3A) . When rVSV/SARS-CoV-2/GFP was passaged in the presence of the single nanobodies, resistant mutants were rapidly selected, as before. Indeed, the yield of infectious virus obtained after two passages in the presence of the single nanobody was nearly indistinguishable from that when rVSV/SARS-CoV-2/GFP was passaged in the absence of nanobodies. In contrast, when nanobodies were combined in cocktails containing two nanobodies, at the same total concentration as was used for the individual nanobodies, in eight out of nine cases, no infectious rVSV/SARS-CoV-2/GFP was recovered after two passages ( Figure 6-figure supplement 1E) . In the ninth case in which S1-48 and RBD-15 were combined and virus was still recovered, sequence analysis revealed that this virus contained two amino acid substitutions, F490V and Y508H, in the RBD. These substitutions were similar or identical to the individual substitutions found in the selection experiments with the single S1-48 and S1-RBD-15 nanobodies, which gave escape variants carrying the substitutions F490S and Y508H, respectively ( Table 7) . These results show that simply combining two nanobodies imposed the requirement for a minimum of two amino acid substitutions to confer resistance to the nanobody cocktail, greatly elevating the genetic barrier for escape. Such mixtures or derived multimers may represent powerful escape-resistant therapeutics, and even more escape resistance should be possible using three or more carefully chosen nanobodies in cocktails or multimers. Integrative structural modeling reveals that the nanobody repertoire explores the available spike epitopes We have taken an integrative modeling approach to generate structural maps of representative nanobody-spike complexes from our repertoire, allowing us to infer likely mechanisms by which our different nanobodies and combinations inhibit the virus. We used the Integrative Modeling Platform (IMP) (Webb et al., 2018) to generate structures using multiple atomic resolution structures available for both spike and the invariant framework of nanobodies as our starting point. Spatial restraints for these calculations were based on our escape mutant data ( Table 7) because for each nanobody its escape mutants cluster around a highly restricted area of its binding epitope on spike (Garrett et al., 2021) ; additional residue-specific distance restraints were generated by cross-linking with MS readout (XL-MS) using the amine-specific bifunctional cross-linkers DSS and BS3 (Shi et al., 2015; Table 8 ). We also incorporated our epitope binning and MP findings (Figure 3, Figure 3 -figure supplement 1) to provide excluded volume validation data Webb et al., 2018) . We benchmarked this modeling approach using a published nanobody with escape mutant data and a solved cryo-EM structure ( Figure 6-figure supplement 2 ; Sun et al., 2021) . These models provide sufficient resolution to map the size and position of the epitopes bound by each nanobody; however, future higher-resolution studies using cryo-EM or crystallization are warranted for the highest priority nanobodies (Schoof et al., 2020; Pymm et al., 2021; Xiang et al., 2020; Wrapp, 2020a) . In sum, we solved integrative structures for 21 different neutralizing nanobodies that, based on our epitope binning data, appeared to collectively explore much of the spike surface, with 18 recognizing RBD, 1 recognizing the NTD of S1, and 2 recognizing S2 ( Figure 6 ). It should be noted that these represent only a small fraction of our total repertoire and so total coverage is greater than what is represented by these maps. Based on overlapping footprints, these 21 nanobodies are classified into 10 groups. Figure 6 summarizes the position of binding and the relative neutralization activity Table 7 . Nanobody neutralization of rVSV/SARS-CoV-2 and selected resistant mutants; related to Figure 6 . Neutralization assays were carried out using rVSV/SARS-CoV-2 and 293T/ACE2cl.22 target cells treated with the denoted nanobodies. Pseudovirus with either wild-type or variant spike (with escape mutants selected using the corresponding nanobody) was used. Escape mutants and IC50s are listed. Amino acid substitutions contributing to loss of neutralization activity are indicated in bold. W64R, Y170H, Q787H >1000 S1-30 trimer S1 non-RBD wt 150 T315I 2400 S1-49 S1 non-RBD wt 146 ± 53.8 S172G >1000 S1-49 dimer S1 non-RBD wt 3.38 ± 2.44 S172G >1000 S1-49 trimer S1 non-RBD wt 0.47 ± 0.00 S172G >1000 S2-10 S2 wt 6649 ± 2,545 W64R, S982R >100,000 (Casalino et al., 2020) , and are instead concentrated at the largely glycanfree RBD. Indeed, among our entire repertoire, epitope binning shows that neutralization activity, corresponding escape mutants, and the mapped epitopes are heavily concentrated on the RBD (Figures 3, 4 and 6) ; ~80% of our anti-RBD nanobodies are neutralizing, with many escape mutants mapping adjacent to the receptorbinding motif (RBM), the region of RBD that interacts directly with ACE2 and is most lightly glycosylated (Shajahan et al., 2020; Watanabe et al., 2020) , whereas ~20% of our anti-S2 nanobodies and ~60% of our non-RBD anti-S1 nanobodies are neutralizing. We note that, based on Table 8 . Cross-linked residues used in integrative modeling; related to Figure 6 . The indicated nanobodies were bound to RBD, NTD, or the spike ectodomain and cross-linked with disuccinimidyl suberate (DSS). Crosslinked complexes were excised from SDS-PAGE gels, reduced, alkylated, and digested with either trypsin or chymotrypsin. Peptides were extracted and analyzed by mass spectrometry. Cross-linked residues (listed) were identified using pLink, and spectra were manually validated to eliminate false positives. the fact that glycosylation obscures a considerable fraction of the spike surface (Watanabe et al., 2020; Zhang et al., 2020) , our repertoire explores much of the remaining available epitope space. The neutralization bias that we observe also likely reflects the most obvious mechanism of viral inhibition, namely, blocking the binding of spike's RBD domain to ACE2 on host membranes to preclude viral fusion, but the non-RBD-based neutralization also underscores that other important mechanisms for viral inhibition exist. Our RBD-binding nanobodies fall into at least seven groups ( Figure 6A ). Many of our mapped nanobodies bind epitopes that partially overlap with previously defined classes of IgG-binding epitopes, but are more compact due to the smaller nanobody paratopes (Corti et al., 2021; Barnes et al., 2020; Xu et al., 2021) ; however, many others define previously unreported binding sites. Overlapping with the RBM are groups 3 (S1-48, S1-RBD-40, S1-23, S1-37, S1-36) and 4 (S1-RBD-29, S1-RBD-21, S1-62); group 3 partially overlaps with previously defined site Ia/class I and group 4 partially overlaps with site Ib/class 2; however, S1-48 uniquely extends beyond the site Ia/class 1 sitting in the saddle of the RBM. Moving further out from the RBM, group 5 (S1-RBD-16 and S1-RBD-23) binds adjacent to the RBM (on the right of the RBD in Figure 6 ) and partially overlaps with site IV/class 3. Group 1 (S1-RBD-15, S1-1, S1-RBD-22, S1-RBD-24, S1-RBD-9) overlaps with site IIa/class 4 and partially overlaps with the opposite side of the RBM as groups 3 and 4. This site appears to be a common nanobody epitope and is shared by VHH-U, VHH-V, and WNb 10 (Koenig et al., 2021; Pymm et al., 2021) . However, a number of our nanobodies map to epitopes that appear distinct from those previously described: group 6 is represented by a single nanobody (S1-RBD-35) and also binds the RBM at a site partially overlapping with site Ib/class 2 (Corti et al., 2021; Tortorici et al., 2020; Dejnirattisai et al., 2021) , but is distinguished by its binding to the tip of the RBM at the left of this representation ( Figure 6 ). S1-46 represents group 7 and defines a unique binding site, adjacent to site IIc/class 4, but higher on the RBD and closer to the RBM. Adjacent and to the right of group 1 is a unique binding site we define as group 2, represented by S1-6. Outside of the RBD, group 8 (represented by S1-49), like many IgGs, binds to the NTD of spike (McCallum et al., 2021) . However, groups 9 (S2-10) and 10 (S2-40) are unique in binding to S2, with S2-10 binding to a region proximal to spike's heptad repeat 1 and S2-40 binding to a cleft between subunits on the stalk of S2 (Wrapp et al., 2020b; Walls et al., 2020) . Of the 10 groups into which our nanobodies were classified, 7 (groups 1, 2, 5, and 7 targeting RBD; group 8 targeting the S1-NTD; groups 9 and 10 targeting S2) do not overlap with the mutations that distinguish alpha, beta, gamma, and delta SARS-CoV-2 VOCs (Figure 6-figure supplement 3A) . Thus, to further validate our predictions for delta virus we tested two nanobodies against the virus in a plaque reduction neutralization assay. As predicted, the group 1 nanobody S1-1 effectively neutralized delta, whereas the group 6 nanobody S1-RBD-35 did not ( Figure 6-figure supplement 3B) . This confirms the potency of our repertoire against key VOC and underscoring the importance of generating large, diverse nanobody repertoires as presented here. A subset of the 21 modeled nanobodies bind sites that interfere with ACE2 binding, preventing the virus from initial binding to its host cell (Wrapp et al., 2020b; Walls et al., 2020; Starr et al., 2020) . Even here, more than one mechanism of inhibition can exist. The RBD is tethered to spike through a hinge, allowing it to fluctuate between either a 'down' conformation, hiding the vulnerable RBM seen in four different orientations. All 21 nanobodies are categorized into 10 groups based on their footprint on spike where groups 1-7 are anti-RBD nanobodies; group 8 contains an anti-S1-NTD nanobody and groups 9 and 10 contain anti-S2 nanobodies. (B) Heatmap of neutralizing epitopes on the structure of SARS-CoV-2 full spike (PDB ID: 6VYB). Epitopes are colored from pale yellow (epitopes with weak neutralization against SARS-CoV-2) to dark red (strong neutralization against SARS-CoV-2). The online version of this article includes the following figure supplement(s) for figure 6: Figure supplement 3. Nanobody groups resistant or predicted to be resistant against SARS-CoV-2 variants of concern (VOC). from the host immune system, or an 'up' conformation, exposing the RBM for potential ACE2 binding and so spike activation/disassembly (Corti et al., 2021; Cai et al., 2020) . Group 1, with overlapping and adjacent epitopes to the RBM, should only bind to the RBD when it is in the 'up' conformation ( Figure 6 ). In contrast, groups 3, 4, and 6 also have overlapping and adjacent epitopes to the RBM, but appear agnostic to RBD conformation and may bind to both 'up' and 'down' positions of the RBD. For nanobodies whose epitopes actually overlap significantly with the RBM (e.g., S1-48, S1-62, S1-RBD-15, or S1-RBD-35) (Figure 6 ), neutralization may occur by directly competing with ACE2 and preventing its binding. For nanobodies within these groups with epitopes more adjacent to the RBM (e.g., S1-36, S1-1, or S1-RBD-9), neutralization may occur through sterically occluding ACE2 binding. However, interestingly, several nanobodies sharing similar epitope bins as S1-RBD-9 (group 1), such as S1-RBD-34, S1-RBD-19, S1-RBD-25, S1-RBD-32, and S1-RBD-36, do not neutralize spike ( Table 1) , suggesting that neutralization may occur via an additional mechanism. Alternatively, different nanobodies can engage a shared epitope with different binding orientations that may or may not hinder ACE2 binding. Additionally, the binding of nanobodies in groups 1, 3, 4, and 6 may mimic ACE2 binding, thus trapping the RBD in its 'up' position to either catalyze the spike trimer rearrangements that prematurely convert spike into a post-fusion state, suppressing viral fusion, or destabilizing the trimer to cause its premature disassembly (Huo et al., 2020b; Walls et al., 2019; Liu et al., 2020a; Turoňová et al., 2020; Benton et al., 2020; Cai et al., 2020; Koenig et al., 2021) . Nanobodies in groups 2, 5, and 7 bind sites distal to the RBM, and therefore are unlikely to neutralize spike through direct ACE2 competition ( Figure 6 ). However, we speculate that they neutralize spike via similar mechanisms. Group 5 (S1-RBD-16, S1-RBD-23) binds to the exposed face of the RBD, and its binding is likely to be agnostic to RBD conformation, and group 2 (S1-6) binds at a position adjacent to group 1, which is enough of a shift that it also may bind both 'up' and 'down' RBD conformations. Group 7 (S1-46) has a peculiar epitope that is only exposed in the 'up' conformation, so nanobody binding likely sterically blocks additional RBDs from accessing the 'up' position (although S1-46 should also be able to bind more than one simultaneously 'up' RBD). In each of these cases, nanobody binding is expected to stabilize the fluctuating RBD in its 'up,' ACE2-engaging, position, potentially destabilizing the trimer similar to mimics of ACE2. S1-49, which is a member of group 8 (Figure 6 ), binds to the NTD for which neutralization activities remain unclear (McCallum et al., 2021) . Human monoclonal antibodies specific to the NTD have been shown not to inhibit ACE2 binding and are instead proposed to inhibit viral infection by blocking membrane fusion, interaction with a different receptor, or proteolytic activation of spike (McCallum et al., 2021) . It remains to be determined if these mechanisms of neutralization hold for our nanobodies that bind non-RBD domains of S1, or if S1-49 suppresses ACE2 binding. The human monoclonals that neutralize the virus by binding outside of the RBD, and their yet to be discovered orthogonal mechanisms of neutralization, emphasize the potential and need for further characterization of our large repertoire of nanobodies. The S2 domain is also a prime, but largely unexplored, therapeutic target (Elshabrawy et al., 2012; Shah et al., 2021) . It is also not where the great majority of mutants in the current VOCs map, making it a particularly exciting target for potentially universal and VOC-resistant therapeutics. Here, we present the first neutralizing nanobodies that bind to S2 (groups 9 and 10) (Figures 2, 4 and 6 , Figure 6-figure supplement 1) . Some monoclonal antibodies that target S2 have been identified and shown to have neutralizing activity, but to our knowledge none have been structurally mapped (Andreano et al., 2021; Wang et al., 2021c; Poh et al., 2020; Li et al., 2020; Song et al., 2021) . Because S2's function is primarily membrane fusion rather than receptor binding, the nanobodies' neutralization mechanisms must differ from those discussed above. For example, the S2-10 escape mutant S982R (Figure 6-figure supplement 1) indicates binding at S982 of spike, positioned at the end of the highly conserved heptad repeat 1, within a region of the S2 that undergoes large dynamic changes as the protein adopts a post fusion conformation; this suggests that S2-10 may restrict this conformational change, thereby inhibiting viral fusion (Cai et al., 2020; Pierri, 2020; Turoňová et al., 2020; Walls et al., 2020) . Notably, the region proximal to S982 appears accessible through an ~30 Å portal, even in the prefusion form with the RBDs in the 'up' position. This is a size not inconsistent with the binding face of a diminutive nanobody but likely inaccessible to conventional antibodies, as has been suggested by others . S2-40 uniquely sits at the interface between spike subunits, which raises fascinating possibilities for its neutralizing activity, perhaps involving the alteration of spike's quaternary structure or dynamics. Nanobodies, as monomeric proteins, can provide a unique opportunity to define possible mechanisms of activity that may otherwise be difficult to distinguish. For example, the dimeric nature of conventional antibodies can introduce ambiguities regarding the mechanisms of neutralization because they can operate either as individual or pairwise binders. In the latter case, they may operate, for example, by aggregation (Thomas et al., 1986) , increased avidity, enhanced steric hindrance via the larger binding entity, or by simultaneously binding and locking two separate moieties within a viral particle. In some cases, for example, S1-7 and S1-25 (which are non-neutralizing as monomers), dimerization did not convert them into neutralizers. In other cases, dimerization and trimerization can engender several folds to orders of magnitude increased neutralization potency (e.g., S1-RBD-35 and S1-23, respectively) ( Figure 4, Table 4) . We even have a curious case where a nanobody such as S2-7 that is essentially non-neutralizing as a monomer becomes strongly neutralizing upon dimerization (Figure 4) . In this latter case, aggregation is a possible contributory mechanism, both between virions -which would lower effective virion concentration -or within a virion, with adjacent spike trimers being cross-linked to each other, inhibiting their function. Although a tremendous range in neutralization improvements by oligomerization is observed both by us and others, there is likely a limit to how much improvement can be induced by oligomerization as, for example, the trimers of S1-23 and S1-RBD-35 do not show a similar fold improvement as to what was observed for the monomer to dimer transition ( Figure 4D , Table 4 ; Schoof et al., 2020; Xiang et al., 2020; Koenig et al., 2021; Xu et al., 2021; Ma et al., 2021) . Drugs are often combined to improve single-agent therapies and dramatically enhance the therapeutic potential of either drug alone while reducing the drug concentrations to be administered. Synergy occurs when the combination of drugs has a greater effect than the sum of the individual effects of each drug. For example, tixagevimab and cilgavimab, two human monoclonal antibodies that target non-overlapping regions of the RBD, function synergistically and show promise as prophylactic and therapeutic agents against COVID-19 (Dong et al., 2020; Zost et al., 2020) . A major advantage of a large repertoire of nanobodies that bind to different epitopes on spike is their strong potential for cooperative activity among nanobody pairs (or higher-order combinations), leading to synergistic viral-neutralizing effects. The small size of nanobodies also provides a great advantage over much larger immunoglobulins in this context as the binding of a nanobody has a lower chance of sterically occluding the binding of a second nanobody to a distinct epitope and because they are monovalent. Moreover, as discussed above, nanobodies binding to the RBD may stabilize the otherwise 'up'-'down' fluctuating RBD in its 'up,' ACE2-engaging, position ( Figure 6 ; Xiang et al., 2020; Schoof et al., 2020; Bracken et al., 2021) . This can have three effects, all of which can potentially promote nanobody synergy: first, it will increase the effective on rate for spike trimer to that measured for monomer, and therefore, make it easier for complementary nanobodies to bind and inhibit; second, by stabilizing this 'up' for any one of the three RBDs in each spike trimer, it destabilizes the 'down' position for the remaining two RBDs, again making nanobody binding from the second class more likely; and third, the 'up' position exposes additional nanobody epitopes that would otherwise be buried (Figure 6 ; Xiang et al., 2020; Sun et al., 2021) . Using an automated platform, we titrated pairwise combinations of nanobodies in a 2D dilution format and measured their IC50s in the pseudovirus assay. IC50s were modeled using a multifaceted synergy framework (Wooten et al., 2021) , including a parameterized version of the equivalent dose model (Zimmer et al., 2016) , the Bivariate Response to Additive Interacting Doses (BRAID) model (Twarog et al., 2016) , and the multidimensional synergy of combinations (MuSyC) model, which models a two-dimensional (2D) Hill equation and extends it to a 2D surface plot (Meyer et al., 2019) . Synergy is evidenced by the parameters of the respective models (Table 9 ). To select nanobody pairs to test for synergy, we took advantage of epitope mapping, structural data, and biophysical characterization. We tested pairwise combinations of nanobodies that bind to similar epitopes, to different epitopes on RBD, and to regions outside and within the RBD (Figure 7, Figure 7-figure supplement 1, Table 9 ). Based on our structural mapping (Figure 6 ), we were able to infer some of the major molecular mechanisms by which these synergistic effects may occur. Table 9 . Nanobody synergy of neutralization activity; related to Parameters from modeling the synergy observed for the indicated nanobody pairs. Multidimensional synergy of combinations (MuSyC), equivalent dose, and Bivariate Response to Additive Interacting Doses (BRAID) models were used to determine if statistically significant synergy was evident from the neutralization response in a 2D grid of nanobody concentrations. 2D neutralization response synergy 10 -3 10 -2 10 -1 10 0 10 2 10 1 10 -3 10 -2 10 -1 10 0 10 2 10 1 10 - 10 -5 10 -4 10 -3 10 -2 10 0 10 -1 10 -5 10 -4 10 -3 10 -2 10 0 10 -1 synergy of neutralization within the S1-RBD E synergy of neutralization within the S1-RBD F S1-23 IC 50 25 nM 50 20 nM S1-23 + S1-27 IC 50 9 nM (IC 50 11 10 -5 10 0 10 -1 10 -2 10 -3 10 -4 10 -5 10 -4 10 -3 10 -2 10 -1 10 0 n o i t a z i l a r t u e n n i e g n a h c % S1-23 promotes "up" S1-RBD-15 stabilizes "up" synergy of neutralization between the S1-RBD and S1-NTD synergy of neutralization between the S1-RBD and S1-NTD 90°S 1-23 promotes "up" S1-46 stabilizes "up" & blocks "down" RBDs 90°S 1-1 stabilizes "up" S1-49 promotes "up" 90°S 1-23 promotes "up" S2-10-dimer requires "up" to bind 90°S 1-23 promotes "up" S1-1 stabilizes "up" 90°S 1-RBD-23 promotes "up" S1-RBD-15 stabilizes "up" 90°S 1-RBD-15 stabilizes "up" S1-46 stabilizes "up" & blocks "down" RBDs 90°S 1-RBD-15 stabilizes "up" S1-49 promotes "up" 90°1 20°S 1-RBD-15 stabilizes "up" S2-10-dimer requires "up" to bind 10 0 0 -1 M) D S1-23 promotes S1-46 stabilizes " & blocks "d " own" RB " 90°1 0 0 0 -1 M) D S1-RBD-1 stabilizes S1-46 stabilizes & blocks "d " own" RB " 90°1 0 0 D S1-23 promotes "u S1-RBD-15 stabilizes "u 90°1 0 0 D S1-RBD-23 promotes " S1-RBD-15 stabilizes "u 90°S 1-23 promotes "u S1-1 stabilizes "u 10 0 90°B D and S S1-RBD-15 stabilizes "u S1-49 promotes "u BD and S1-NTD 10 2 BD and S1 S1-1 stabilizes "u S1-49 promotes " BD and S1-NTD 90°1 0 0 BD a S1-23 promotes "u S2-10-dime requires "up o bind to 90°) RBD a S1-RBD-1 stabilizes S2-10-dim requires " ) 10 0 10 -1 up 120°r equires to bind Figure 7 . Synergistic neutralization of spike with nanobody cocktails. (A) An example of additive effects between two anti-SARS-CoV-2 spike nanobodies. S1-23 and S1-27 were prepared in a two-dimensional serial dilution matrix and then incubated with SARS-CoV-2 pseudovirus for 1 hr before adding the mixture to cells. After 56 hr, the expression of luciferase in each well was measured by addition of Steady-Glo reagent and read out on a spectrophotometer. The left panel shows a heatmap of pseudovirus neutralization by a two-dimensional serial dilution of combinations of S1-23 and S1-27. Lines and red numbers demarcate the % inhibition, that is, inhibitory concentration where X% of the virus is neutralized, e.g., IC50. Dark Combinations of S1-27 and S1-23 showed simple additive effects ( Figure 7A ). These nanobodies belong to the same epitope bin ( Figure 3A) ; their additive effect is as expected for two nanobodies accessing the same site on S1-RBD, but, for example, in equal concentrations, effectively doubling the concentration of a single nanobody. The potential for synergy resides instead in nanobodies that bind to different epitopes and can bind to spike monomers simultaneously. We therefore tested combinations that bind to different epitopes, first focusing on the RBD. Indeed, powerfully synergistic effects were observed between numerous nanobody pairs. For example, the combination of S1-23 and S1-1, which bind to opposite sides of the RBD, dramatically increased the potency of both nanobodies by ~32-and ~ 21-fold, respectively ( Figure 7B, Table 9 ). S1-1 is expected to bind to RBD in its 'up' position, while S1-23 can bind to both 'up' and 'down' RBDs. We interpret this reciprocal synergy observed as S1-23 promoting the 'up' position and S1-1 stabilizing the 'up' position. S1-RBD-15, which binds to a similar epitope as S1-1, shows corresponding synergy with S1-23, suggesting that synergy may be predictable based on epitope mapping ( Figure 7C ). In this case, however, S1-RBD-15 had a greater influence on S1-23, promoting its potency by ~300-fold. S1-RBD-15 showed a comparable synergy profile against S1-RBD-23, which binds to the opposite side of the RBD, and adjacent to the site occupied by S1-23 ( Figures 3A, 6A and 7D) . The synergy profiles observed by these pairs of nanobodies highlight how stabilizing the RBD in the 'up' position can have a dramatic effect on their ability to neutralize spike activity. However, it should not be taken for granted that simply binding to distinct epitopes on RBD simultaneously will always be sufficient to generate a strongly synergistic response. For example, S1-46 failed to show synergy with either S1-23 or S1-RBD-15 ( Figure 7E and F). Indeed, synergy modeling indicates that S1-46 actually mildly antagonizes both S1-23 and S1-RBD-15 (Table 9 ). S1-46 binds an epitope on RBD only in the 'up' position ( Figure 6 ), but in this case binding hinders the movement of adjacent RBDs, and therefore reduces the activity of nanobodies depending on the dynamics of adjacent RBDs. S1-49, which binds to the NTD of spike (Figure 6) , substantially improved the neutralization potency of either S1-1 or S1-RBD-15 ( Figure 7G and H) . The synergy observed with S1-49 with S1-1 elicited a >1000-fold increase in potency. Interpretation of the mechanism underlying this remarkable synergistic effect will require a greater understanding of the mechanism of NTD activity on spike (McCallum et al., 2021) . We also tested for synergy between nanobodies targeting the S1-RBD and S2, which revealed remarkable results. In this case, we used a dimer of S2-10 to increase its potency to be closer to that of the nanobodies to which it was paired. Interestingly, among all the nanobody pairs that we tested, the synergy was greatest with the S2-10-dimer, which showed >4000-fold increase in potency when combined with either S1-23 or S1-RBD-15 ( Figure 7I and J) . S2-10 recognizes a site occluded by the S1-RBD when it is in its 'down'' position, but is revealed and becomes accessible when in its 'up' position. Thus, we interpret the mechanism of synergy as one of cooperativity where these RBD-binders promote the 'up' state and provide increased access for S2-10. The strong synergy observed may also reflect the distinct mechanisms by which either the RBD-binders or S2-10 operate individually. blue regions are concentrations that potently neutralize the pseudovirus, as per the heatmap legend. The right panel shows neutralization curves (with 90% confidence interval bands) and the calculated IC50 of each nanobody alone, or in a 1:1 combination was determined along with a calculated IC50 based on the theoretical additive mixture model of the pair (curve with dotted gray line). The inset shows a difference (synergy) map calculated as the difference between the parameterized 2D neutralization response and that expected in a null model of only additive effects. Here, no difference is observed. (B) S1-1 synergizes with S1-23 in neutralizing SARS-CoV-2 pseudovirus. The left panel shows the heatmap of pseudovirus neutralization observed by a two-dimensional serial dilution of combinations of S1-1 and S1-23. The middle panel shows a heatmap mapping the synergy of neutralization observed for this pair. The lines bounding the darker purple areas demarcate regions in the heatmap where the observed neutralization is greater than additive by the indicated percentages (yellow numbers), as per the heatmap legend. The right panel shows two representations of spike with the accessible S1-1 (salmon) and S1-23 (steel blue) epitopes (PDB ID: 6VYB). (C-J) Examples of synergy between nanobodies binding the S1-RBD, or between the S1-RBD and S1-NTD or S2 domains of spike. The layout is as found in (B), but comparing S1-RBD-15 with S1-23 (C), S1-RBD-15 with S1-RBD-23 (D), S1-23 with S1-46 (E), S1-RBD-15 with S1-46 (F), S1-49 with S1-1 (G), S1-49 with S1-RBD-15 (H), S1-23 with S2-10-dimer (I), and S1-RBD-15 with S2-10-dimer (J). The online version of this article includes the following figure supplement(s) for figure 7: Our repeated observation of strong synergistic effects between nanobodies is especially noteworthy, reflecting the unique properties of nanobodies such as their small size that are not shared by, for example, human monoclonal antibodies. For example, we structurally aligned human IgGs with the S1-23 and S1-1 paratopes and found the same binding characteristics would not be predicted to act synergistically because they could not bind the same monomer simultaneously as they would clash with other RBDs on the spike trimer. In the case of S1-23 and S1-RBD-15 epitopes, structural alignment of IgGs with nanobody paratopes suggested a strong inter-IgG steric clash in addition to a clash with other RBDs in the down position. While native nanobodies' small size and lack of an Fc domain affect their pharmacokinetic behavior in comparison to standard IgGs when used therapeutically, numerous studies have demonstrated that various modifications, such as albumin-binding domains or a synthetic Fc, are available to tune half-life and other behaviors upon either intravenous or intranasal delivery (Tijink et al., 2008; Nambulli et al., 2021; Pymm et al., 2021; Shen et al., 2021) . The data presented here demonstrate the power of raising large and diverse repertoires of nanobodies against the entire ectodomain of SARS-CoV-2 spike to maximize the likelihood of generating potent reagents for prophylactics and therapeutics. Moreover, several neutralizing nanobodies in our current repertoire are, or are predicted to be, effective against current circulating variants (Figure 6figure supplement 3) . However, our escape experiments support the idea that the current circulating variants are not yet necessarily exploring the full potential of the virus to escape our current and emerging therapeutic arsenals, and that even if antibodies or nanobodies are resistant to the current variants, they will not necessarily be resistant to variants as they continually emerge. To counteract this eventuality, we show that judicious choice of nanobody combinations that can synergize and have orthogonal and complementary neutralization mechanisms have the potential to result in potent and broadly neutralizing reagents that are resistant to viral escape. Collectively, this unique and readily modifiable repertoire has the potential to complement vaccines, drugs and single epitope reagents, and guard against single-molecule failure in human trials even in the face of emerging variants. Most urgently, it paves the way to develop therapeutics for hospitalized patients with acute disease, and address the unmet needs of patients in the developing world, many of which will not see a COVID-19 vaccine before 2023 (Padma, 2021) . See Appendix 1-key resources table. To maximize the purity of the serum HCAb sample, we explored different binding conditions to select for the tightest V H H binders -a key step not generally available to display panning methods (Fridy et al., 2014) . We also used an additional HCAb purification step to deplete VH IgG by incubation with immobilized Protein M, a mycoplasma protein specific for IgG light chain (Grover et al., 2014) . To further enrich the V H H sample for MS analysis and remove Fc, we performed a digest with IdeS, a protease that cleaves the V H H domain from the HCAb with higher specificity than conventionally used papain (von Pawel-Rammingen et al., 2002) . Greater peptide coverage for LC-MS was attained by using complementary digestion with trypsin and chymotrypsin , augmented by partial SDS-PAGE gel-based separation of different V H Hs to reduce the V H H complexity and to give more complete peptide coverage and candidate selection. We redesigned PCR primers to maximize coverage of V H H sequences for our cDNA libraries. Also, to increase the reliability of the library, singletons were not considered as candidates and priority was given to sequences with high counts. Finally, we refined our Llama-Magic software package (Fridy et al., 2014) to include improved scoring functions, weighting the length, uniqueness, and quality of the MS data especially for complementaritydetermining regions. This optimized protocol allowed us to identify 374 unique CDR3 sequences (from 847 unique V H H candidates). Details are provided below. Recombinant Fc-tagged SARS-CoV-2 spike S1 and S2 proteins purified from HEK293 cells were used for llama immunization (The Native Antigen Company; REC31806 and REC31807). For affinity isolation, binding screens, SPR analysis, and MP, recombinant spike S1-His, untagged RBD, or S2-His proteins expressed in HEK293 (S1 and RBD), or insect cells (S2) were used (Sino Biological; 40591-V08H, 40592-VNAH, and 40590-V08B). Native mass spectrometry (Olinares et al., 2016; Olinares et al., 2021) was used to confirm the quality of these proteins and determine their glycosylation state, with S1 and S2 observed to be heavily glycosylated (at least 10 kDa of attached glycans). RBD was observed to be monomeric, S1-His likely monomeric, and S2-His, a mix of monomer and trimer. We used a pre-screening protocol to select llamas with naturally strong immune responses, as determined by activity against standard animal vaccines (Thompson et al., 2016) . Two llamas, Marley (9-year-old male) and Rocky (5-year-old male), were immunized with recombinant SARS-CoV-2 spike S1 and SARS-CoV-2 spike S2 expressed in HEK293 cells as Fc fusion proteins. Llamas were injected subcutaneously with 0.25 mg of each antigen with CFA, then boosted with the same amounts with IFA at intervals of 14, 7, 21, and 10 days. Serum bleeds and bone marrow aspirates were obtained 9 days after the final injection. From the production serum bleeds, HCAb fractions of IgG were obtained by purification with immobilized Protein A and Protein G as previously described (Fridy et al., 2014) . Residual light-chain-containing IgG was removed from this fraction by incubating with 25 µl of 10 mg/ ml Protein M-Sepharose per mg of HCAb (Grover et al., 2014) . After a 30 min incubation, the HCAb flow-through was collected. 12 mg of this HCAb fraction was then incubated with Sepharoseconjugated recombinant SARS-CoV-2 spike S1-His, RBD, or S2-His protein. This resin was washed with (1) 20 mM sodium phosphate, pH 7.4 + 500 mM NaCl; (2) 2 M MgCl 2 in 20 mM Tris, pH 7.5; (3) PBS + 0.5 % Triton X-100; and (4) PBS. The resin was then resuspended in a 200 μl solution of 2 U/ µl IdeS enzyme (Genovis) in PBS and digested for 3.5 hr at 37°C on an orbital shaker. The resin was then washed with (1) PBS, (2) PBS plus 0.1% Tween-20, and (3) PBS. Bound protein was eluted by incubating 10 min at 72°C in 1× NuPAGE LDS sample buffer (Thermo Fisher). The samples were reduced with DTT and alkylated with iodoacetamide, then run on a 4-12% Bis-Tris gel. Bands at ~15 kDa and ~20 kDa corresponding to digested V H H region were then cut out and prepared for MS. Bone marrow aspirates were obtained from immunized llamas concurrent with production serum bleeds. Bone marrow plasma cells were isolated on a Ficoll gradient using Ficoll-Paque (Cytiva). RNA was isolated from approximately 3-4 × 10 7 cells using TRIzol reagent (Thermo Fisher), according to the manufacturer's instructions. cDNA was synthesized using SuperScript IV reverse transcriptase (Thermo Fisher). A PCR was then performed with V H H IgG specific primers and Deep Vent polymerase (New England Biolabs). Forward primers 6N_CALL001 5′-NNNN NNGT CCTG GCTG CTCT TCTA CAAGG-3′ and 6N_CALL001B 5′-NNNN NNGT CCTG GCTG CTCT TTTA CAAGG-3′ target the leader sequence (Conrath et al., 2001) while reverse primers 6N_VHH_SH_rev 5′-NNNN NNCT GGGG TCTT CGCT GTGGTGC-3′ and 6N_VHH_LH_rev 5′-NNNN NNGT GGTT GTGG TTTT GGTG TCTTGGG-3′ target short and long hinge sequences at the 3′ side of V H H. Primers included six random bases (N) to aid cluster identification. The approximately 350-450 bp product of this reaction was gel purified, then ligated to Illumina adaptors before library preparation using Illumina kits, before MiSeq sequencing using two 300 bp paired end reads. Trypsin (Roche) or chymotrypsin (Promega) solution was added to previously reduced, alkylated, diced, destained, and dehydrated gel pieces at ~1:4-3:1 enzyme to substrate mass ratios. Gel pieces were allowed to rehydrate with enzyme solution for 10 min on ice. 45 μl of digestion buffer (trypsin: 50 mM ammonium bicarbonate, 10% acetonitrile; chymotrypsin: 100 mM Tris pH 7.8, 10 mM CaCl 2 ) were then added, and samples were incubated for 6 hr at 37°C (trypsin) or 25°C (chymotrypsin). Supernatant was then removed from gel pieces and transferred to a new tube. 150 μl of a 1.67% FA, 67% ACN, 0.05% TFA solution were added to gel pieces, and shaken at 4°C for ~6 hr. Supernatant was removed from gel pieces, transferred to the tube with previous supernatant, and evaporated in a speedvac until dry. Samples were resuspended in 5% formic acid, 0.1% TFA, and cleaned on StageTips (Rappsilber, 2012) . Samples were analyzed with a nano-LC 1200 (Thermo Fisher) using an EASYspray PepMap RSLC C18 3 µm, 100 Å, 75 µm × 15 cm column coupled to an Orbitrap Fusion Lumos Tribrid mass spectrometer (Thermo Fisher). An Active Background Ion Reduction Device (ABIRD, ESI Source Solutions) was used to reduce background. The Lumos was operated in data-dependent mode, and top intensity ions were fragmented by high-energy collisional dissociation (normalized collision energy 28). Ions with charge states 2-5 were selected for fragmentation. Orbitrap resolution was 120,000. The quadrupole isolation window was 1.4, and the MS/MS used a maximum injection time of 250 ms with one microscan. The initial identification of nanobody sequences was performed as described (Fridy et al., 2014) using the program Llama-Magic (https:// github. com/ FenyoLab/ llama-magic) with a few added features (including being able to deal with chymotryptic proteolysis and to rank V H Hs by corresponding read counts in high-throughput sequencing data), where 23 MS datasets (concatenated from all MS acquisition data according to antigens, animal individuals, gel band positions and proteases) were independently searched. The results were filtered with criteria including read counts, uniqueness score, and quality and coverage of MS/MS fragments to generate a collection of high-confidence nanobody sequences. A CDR3 network graph was created by connecting nodes (unique high-confidence CDR3 sequences) by edges where a CDR3 pair has a Damerau-Levenshtein distance of no more than three by using NetworkX 2.5 (https:// networkx. org) and pyxDamerauLevenshtein (https:// github. com/ gfairchild/ pyxD amer auLe vens htein ; Fairchild, 2013) . The diversity of nanobodies for screening was maximized by selecting CDR3 sequences from isolated components of the network graph, together with varying CDR3 lengths and animal individual origin. Nanobody sequences were codon-optimized for expression in Escherichia coli and synthesized as gene fragments (IDT), incorporating BamHI and XhoI restriction sites at 5′ and 3′ ends, respectively. Nanobody sequences were then subcloned into pET21-pelB using BamHI and XhoI restriction sites as previously described (Fridy et al., 2014) . pelB-fused nanobodies were expressed and purified using Arctic Express (DE3) cells (Agilent) as previously described using TALON metal affinity resin (Takara) (Fridy et al., 2014) . Nanobodies to be oligomerized were ordered from IDT as minigenes incorporating at the 5′ end a SalI site followed by codon optimized sequence for the linker GGGGSGGGGSGGGGSGGGGS upstream of the start codon of the nanobody cDNA, and at the 3′ end of the nanobody the coding sequence a XhoI site was added. The minigene was cut with SalI and XhoI, the linker-nanobody insert was gel purified and ligated with the XhoI linearized recipient nanobody expression vector (pET21-pelB + nanobody). Restriction digests and sequencing was performed to identify two (dimer) and three (trimer) oligomers. To validate nanobody candidates, pelB-fused nanobodies were expressed in 50 ml cultures of Arctic Express (DE3) cells, and the periplasmic fractions were isolated by osmotic shock as previously described (Fridy et al., 2014) . Spike S1-His, RBD, or S2-His proteins (Sino Biological 40591-V08H, 40592-VNAH, and 40590-V08B) were conjugated to cyanogen bromide-activated Sepharose 4 Fast Flow resin (Cytiva) according to the manufacturer's instructions using 100 µg protein per mg of resin. Periplasm was incubated with 15 µl of the corresponding antigen-conjugated Sepharose for 30 min while rotating at room temperature (RT). The resin was then transferred to a spin column and washed twice with buffer TBT-100 (20 mM HEPES pH 7.4, 100 mM NaCl, 110 mM KOAc, 2 mM MgCl 2 , 0.1% Tween 20). Bound protein was eluted with 1.2× NuPAGE LDS sample buffer (Thermo Fisher) for 10 min at 72°C, then reduced with 50 mM DTT (10 min at 72°C). Input and elution samples were separated by SDS-PAGE, and Coomassie-stained bands were quantified using ImageJ software. Surface plasmon resonance K D s were determined via SPR experiments. Measurements were either taken on a Proteon XPR36 Protein Interaction Array System (Bio-Rad) or a Biacore 8k (Cytiva). Recombinant spike S1, RBD, and spike S2 were immobilized at 5 µg/ml, 5 µg/ml, and 12.5 µg/ml, respectively, using the ProteOn Amine Coupling Kit (EDC/NHS coupling chemistry, Bio-Rad) according to the respective manufacturer's guidelines either on a ProteOn GLC sensor chip or a Series S CM5 sensor chip. All purified nanobodies in a final buffer containing 20 mM HEPES pH 7.4, 150 mM NaCl, 0.02% Tween, were prepared in 5-8 concentrations. For experiments performed on the Proteon XPR36, protein was then injected at 50 μl/min for 120 s, followed by a dissociation time of 600 s. Residual bound proteins were removed by regenerating the chip surface using glycine pH 3 + 1 M MgCl 2 . Data were processed and analyzed using the ProteOn Manager software. For experiments performed on the Biacore 8k, protein was injected at 60 μl/min for 120 s, followed by a dissociation time of either 1200 s or 2400s. Residual bound proteins were removed by regenerating the chip surface using glycine pH 2.5 + 1 M MgCl 2 . Data were processed and analyzed using the Biacore Insight Evaluation software. Nanobody melting temperatures (T m ) were measured by DSF using a CFX96 Real-Time PCR Detection System (Bio-Rad, Hercules, CA). A 96-well thin-wall hard-shell PCR plate (Bio-Rad) was set up with each well containing 10-40 µM of protein in 20 mM HEPES, 150 mM NaCl buffer (pH 7.4), 5× SyproOrange Protein Gel Stain (Sigma-Aldrich). Fluorescence variation was measured from 25 to 95°C at a ramp rate of 0.5°C/5 s. Excitation was between 515 and 535 nm, and emission was monitored between 560 and 580 nm. T m was the transition midpoint value, calculated using the manufacturer's software (Niesen et al., 2007) . Nanobodies in 20 mM HEPES, 150 mM NaCl, pH 7.4 at concentrations between 0.33 mg/ml and 0.63 mg/ml were snap-frozen in liquid nitrogen and dried in a speed-vac to replicate lyophilization conditions. Nanobodies were then reconstituted in dd H 2 O and characterized using SPR and DSF. Epitope mapping studies were carried out using the Octet system (ForteBio, USA, version 7) that measures biolayer interferometry (BLI). All steps were performed at 30°C with shaking at 1300 rpm in a black 96-well plate containing 300 μl kinetics buffer (PBS; 0.2% BSA; 0.02% sodium azide) in each well. AMC-coated biosensors were loaded with mFc tagged RBD (Sino Biological) at 40 μg/ ml to reach >1 nm wavelength shift following binding and washing. The sensors were then reacted for ~300 s with reference nanobodies and then transferred to kinetics buffer-containing wells for another 180 s. A new baseline was set, sensors were then reacted for 180 s with analyte nanobodies (association phase), and then transferred to buffer-containing wells for another 180 s (dissociation phase). Binding and dissociation were measured as changes over time in light interference after subtraction of parallel measurements from unloaded biosensors. Sensorgrams of analyte association/dissociation responses were analyzed using the Octet data analysis software 7.1 (Fortebio, USA, 2015) . Analyte binding to mFc RBD was also measured in parallel to get response levels in the absence of the reference nanobodies. Octet response values were used to compute a Pearson's correlation coefficient for pairwise combinations of nanobodies using Pandas (McKinney, 2010) in Python 3.7.6 (https://www. python. org/). These coefficients were then used to hierarchically cluster the nanobodies and were visualized as a heatmap (Pedregosa, 2011) . The undirected unweighted network graph of Octet response values was constructed by treating each nanobody as a node, adding an edge to each measured pair of different nanobodies, and setting the maximum response value of a nanobody pair as an attribute to the edge, by using NetworkX 2.5 (https:// networkx. org). The least responses of pairwise nanobodies within all fully measured nanobody subsets were computed by iterating through all network cliques of size 2-14 by using NetworkX's 'find_cliques' function, and taking the minimum value of edge attributes within each clique. Network coefficients (average shortest path length, average clustering coefficient, and small-world coefficient sigma) were computed using NetworkX's 'average_shortest_path_length,' 'average_clustering,' and 'sigma' functions. Network visualization was created by using D3. js (https:// d3js. org). Select nanobodies were binned using MP. Experiments were performed on a Refeyn OneMP instrument (Refeyn Ltd). The instrument was calibrated with a mix of BSA (Sigma-Aldrich), thyroglobulin (Sigma-Aldrich), and beta-amylase (Sigma-Aldrich). Coverslips (Thorlabs) and gaskets (Grace Bio-Labs) were prepared by washing with 100% IPA followed by dd H 2 O, repeated three times, followed by drying with HEPA filtered air. 12 μl of buffer was added to each well to focus the instrument after which 8 μl of protein solution was added and pipetted up and down to briefly mix after which movies/frame acquisition was promptly started. The final concentration in each experiment of recombinant spike S1 monomer (Sino Biological) and each nanobody was 30 nM and between 25 and 40 nM, respectively. Movies were acquired for 60 s (6000 frames) using AcquireMP (version 2.3.0; Refeyn Ltd) using standard settings. All movies were processed, analyzed, and masses estimated by fitting a Gaussian distribution to the data using DiscoverMP (version 2.3.0; Refeyn Ltd). Epitope mapping of anti-S2 and non-RBD anti-S1 nanobodies SPR was utilized to perform epitope binning experiments using a Biacore 8k (Cytiva) supplemented with the Biacore Insight Epitope Binning Extension. All nanobodies' concentrations were ≥20× the concentration of their K D for binning experiments, with the majority surpassing their K D by 50×. For non-RBD anti-S1 nanobodies, experiments were performed utilizing either the tandem method or dual-tandem method for epitope binning, whereas for anti-S2 nanobodies, only the tandem method was utilized. Series S CM5 sensor chips immobilized with spike S1 and spike S2 were used (see 'Surface plasmon resonance' section above for full details). For the tandem method, nanobody '1' was injected at 10 μl/min for 240 s, followed by a brief wash, after which nanobody '2' was injected at 10 μl/min for 240 s and dissociated for 30 s. Residual bound proteins were removed by washing the chip surface four times with 10 mM glycine pH 2 + 1 M MgCl 2 at 60 μl/min for 60 s. For the dual-tandem method, nanobody '1' was injected at 10 μl/min for 120 s, followed by nanobody '2,' which was injected at 10 μl/min for 150 s and dissociated for 30 s. Residual bound proteins were removed by washing the chip surface three times with 10 mM glycine pH 2 + 1 M MgCl 2 at 60 μl/min for 60 s. For the anti-S2 nanobody binning experiments, residual bound proteins were removed by washing the chip surface first with 0.1 M HCl at 60 μl/min for 60 s, followed by a second wash with 3 M MgCl 2 at 60 μl/min for 60 s. Data were processed and analyzed using the Biacore Insight Evaluation software utilizing the Epitope Binning Extension. Vero E6 cells (ATCC) were cultured at 30°C in the presence of 5% CO 2 in medium composed of high-glucose Dulbecco's modified Eagle's medium (DMEM, Gibco) supplemented with 5% (v/v) heatinactivated fetal bovine serum (FBS) (VWR). TMPRSS expressing Vero E6 cells (gift from Rhea Coler) were cultured in DMEM supplemented with 5% (v/v) FBS and 1 mg/ml geneticin. 293T/17 and 293T-hACE2 (Cawford et al., 2020) cells (Life Technologies; Cat# R70007; RRID:CVCL_6911) were cultured in DMEM (Gibco) supplemented with 10% 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-1, SARS-CoV-2, and SARS-CoV-2 variant pseudotyped lentiviral reporter particles Pseudovirus stocks were prepared using a modified protocol published by Cawford et al., 2020; Qing et al., 2020 . Briefly, pseudovirus stocks were prepared by cotransfecting 4.75 µg pHAGE-CMV-Luc2-IRES-ZsGreen-W (BEI Cat # NR-52516) (Cawford et al., 2020) , 3.75 µg psPAX and 1.5 µg spike containing plasmid using lipofectamine 3000 (Thermo Fisher). 4 × 10 6 cells were plated 16-24 hr prior to transfection. 60 hr post transfection, pseudovirus containing media was collected, cleared by centrifugation at 1000 × g, and filtered through a 0.45 µm syringe filter to clear debris. 1 ml aliquots were frozen at -80°C. Pseudovirus was titered by threefold serial dilution on 293T-hACE2 cells (Cawford et al., 2020) , treated with 2 µg/ml polybrene (Sigma). 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. All periplasmic purified nanobodies were treated with Triton X-114 to remove any residual endotoxins so as to not have endotoxins contribute to the effective neutralization (Aida and Pabst, 1990) , and residual detergent was removed using Pierce Detergent Removal Resin according to the manufacturer's instructions (Thermo Fisher). 293-hACE2 cells were plated at 2500-4000 cells per well on 384 solid white TC-treated plates. Threefold serially diluted nanobodies (10 dilutions in total) were incubated with 40,000-60,000 RLU equivalents of pseudotyped SARS-CoV-2-Luc for 1 hr at 37°C. Mock treatment and a sham treatment with LaM2 nanobodies (Fridy et al., 2014) that do not bind to spike were included as negative controls while untreated wells were used to monitor background levels. 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 2 . Infected cells were processed between 52 and 60 hr as described above. Data were processed using Prism 7 (GraphPad), using fourparameter nonlinear regression (least-squares regression method without weighting). All nanobodies were tested at least two times and with more than one pseudovirus preparation. Experiments were performed as per our pseudovirus neutralization assay. A robotic liquid handler was used to prepare 2D matrices of serial dilutions of two nanobodies and then mix these with SARS-CoV-2 pseudovirus 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 synergy software (Wooten et al., 2021) . Integrative structural modeling proceeded through the standard four-stage protocol (Russel et al., 2012; Kim et al., 2018; Rout and Sali, 2019; Saltzberg et al., 2021) , which was scripted using the Python Modeling Interface package, a library for modeling macromolecular complexes based on the Integrative Modeling Platform software (Russel et al., 2012) , version develop-31a0ad09b4 (https:// integrativemodeling. org). Separate models were computed for rigid-receptor-rigid ligand-type binary docking of (1) 18 nanobodies (S1-RBD- [9, 15, 16, 21, 22, 23, 24, 29, 35, 40] and S1- [1, 6, 23, 36, 37, 46, 48, 62] ) on a monomeric S1-RBD domain, (2) the S1-49 nanobody on a monomeric S1-NTD domain, and (3) S2-10 and S2-40 nanobodies on the trimeric S2 domain of the spike protein. Monomeric S1-RBD, spanning amino acids T333-G526, was represented using the crystal structure of the co-complex of ACE2 bound RBD (PDB ID: 6M0J; Lan et al., 2020) . Monomeric S1-NTD, spanning amino acids V16-S305, was represented using the crystal structure of the S2M28 Fab bound NTD (PDB ID: 7LY3; McCallum et al., 2021) . Trimeric S2 was represented using the residues S689-P1162 (for each monomer) from a 2.9 Å cryo-EM structure with PDB ID: 6XR8 (Cai et al., 2020 , McCallum et al., 2021 . Comparative models of all 21 nanobodies were built from the crystal structure of the human Vsig4 targeting nanobody Nb119 (PDB ID: 5IML; Wen et al., 2017) as template using MODELLER (Sali and Blundell, 1993) , and their CDR3 regions were refined using MODELLER's loop modeling algorithm (Fiser et al., 2000) . To maximize the efficiency of structural sampling while avoiding too much information loss, 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 alternate 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 select a random subsample of 30,000 models, which were clustered by the similarity of their interfaces to the receptor. The fraction of common contacts (fcc) between receptor and nanobody was used to characterize interface similarity between alternate nanobody poses (Fiser et al., 2000; Rodrigues et al., 2012) . Binding poses belonging to only 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. However, no significant changes were observed in the satisfaction of restraints. 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; Rodrigues et al., 2012) . Files containing input data, scripts and output results are available at https:// github. com/ integrativemodeling/ nbspike (Sanyal, 2021 ; copy archived at swh:1:rev:2607a97503e1d764107 9641142734f4075d334e2). SARS-related coronavirus 2, isolate USA-WA1/2020, NR-52281, was deposited by the Centers for Disease Control and Prevention and obtained through BEI Resources, NIAID, NIH. SARS-CoV-2, isolate USA (B.1.617.2), was a kind gift from Rhea Coler. Viral stocks were propagated in Vero E6 cells. 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. 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, which were prominent after 48 hr of infection. After 72 hr, infectious supernatants were collected, cleared of cellular debris by centrifugation, and stored at -80°C until use. Viral titers were determined by plaque assay using a liquid overlay and fixation-staining method, as described Case et al., 2020) . Briefly, serially diluted virus stocks were used to infect confluent monolayers of Vero E6 cells (~1.2 × 10 6 cells per well) cultured in six-well plates. After a 90 min incubation, the virus was removed, and the cell monolayer overlaid with a medium composed of 3% (w/v) carboxymethylcellulose and 4% (v/v) FBS in phenyl-free Opti-MEM. 96 hr post infection, the viscous carboxymethylcellulose medium was removed and the cells were washed once with Dulbecco's phosphate buffered saline (DPBS; Gibco) before being fixed with 4% (w/v) paraformaldehyde in DPBS. After a 30 min incubation, the fixative was removed, and the cells were rinsed with DBPS before being stained with 1% (w/v) crystal violet in 20% (v/v) ethanol. Contrast was enhanced by successive washes with DPBS, and clear plaques representing individual viral infections were visualized as spots lacking the stain. Plaques were enumerated by first identifying the dilution factor of the well containing 10-100 plaques. After counting the plaques, the average number of plaque forming units (pfus) from three experiments was used to determine the viral titer by dividing the average by the dilution factor and volume of virus delivered per well. Focus forming reduction assay with authentic SARS-CoV-2 Nanobody neutralization of infectious SARS-CoV-2 was performed using a focus forming reduction assay. Briefly, eight threefold serial dilutions of nanobodies were incubated with ~7.5 × 10 4 focus forming units of SARS-CoV-2 for 1 hr at RT. The mixture was then added to a confluent monolayer of Vero E6 cells or 293-ACE2 (Cawford et al., 2020) plated at ~1.5 × 10 5 cells per well and seeded in 48-well plates. 24 hr post infection, the cells were washed once with DPBS, trypsinized with 0.05% trypsin (Gibco), and fixed for 30 min with 4% paraformaldehyde in DPBS. After fixation, the cells were permeabilized with 1% (w/v) Triton X-100 (Sigma Aldrich) for 30 min. After permeabilization, the cells were incubated with a blocking buffer (1% [w/v] bovine serum albumin [Calbiochem] and 0.5% [w/v] Triton X-100 in DBPS) for 60 min, and then stained with primary anti-spike CR3022 (Absolute Antibody) monoclonal antibodies (1:1000), and secondary anti-human IgG antibodies (1:2000) conjugated to Alexa Fluor 488 (Invitrogen). Cells staining positive for spike were measured by flow cytometry on a Becton Dickinson BD LSR II Special Order System Flow Cytometer With HTS Sampler. The percentage of spike-positive cells from triplicate wells for each dilution was used to determine the half maximal inhibitory concentrations (IC50) using a parametric 1D Hill fitting algorithm with synergy (Wooten et al., 2021) . A mock treatment, sham treatment with LaM2 nanobodies (Fridy et al., 2014) , and untreated cells were used as controls. Nanobody neutralization of infectious SARS-CoV-2 was performed using a plaque reduction neutralization test (PRNT) assay. Briefly, 10 threefold serial dilutions of nanobodies were incubated with ~100-300 pfus of SARS-CoV-2 for 1 hr at RT. The mixture was then added to a confluent monolayer of TMPRSS2+ Vero E6 cells (~6 × 10 5 cells per well) in 12-well plates. After a 90 min incubation, the virus was removed, and the cell monolayer overlaid with a medium composed of 3% (w/v) carboxymethylcellulose and 4% (v/v) FBS in phenyl-free Opti-MEM. 96 hr post infection, the viscous carboxymethylcellulose medium was removed and the cells were washed once with DPBS (Gibco) before being fixed with 4% (w/v) paraformaldehyde in DPBS. After a 30 min incubation, the fixative was removed, and the cells were rinsed with DBPS before being stained with 1% (w/v) crystal violet in 20% (v/v) ethanol. Contrast was enhanced by washing with DPBS, and clear plaques representing individual viral infections were visualized as spots lacking the stain. The number of plaques at each dilution was used to determine the IC50s of each nanobody. Assays with primary airway epithelial cell cultures were performed as described (Barrow et al., 2021) . Briefly, bronchial AECs were obtained under study #12490 approved by the Seattle Children's Institutional Review Board, with investigations carried out following the rules of the Declaration of Helsinki of 1975. AECs were differentiated for 21 days at an ALI on 12-well collagen-coated Corning plates with permeable transwells in PneumaCult ALI media (Stemcell, Vancouver, BC, Canada). Differentiated AECs were treated with nanobodies diluted in PBS, or PBS alone for 60 min, the liquid was removed, and the AECs were then infected with SARS-CoV-2 at a multiplicity of infection (MOI) of 0.5. At 24 hr intervals, the cells were treated with nanobodies or PBS for 60 min. After 96 hr of infection, SARS-CoV-2 viral replication was measured in AEC cultures by quantitative PCR, with triplicate assays of harvested RNA from each SARS-CoV-2-infected AEC donor cell line (Genesig Coronavirus Strain 2019-nCoV Advanced PCR Kit, Primerdesign, Southampton, UK). The concentration of RNA harvested from AECs was used to normalize the qPCR data and was measured on a spectrophotometer (NanoDrop). Nanobodies were fivefold serially diluted and then incubated with rVSV/SARS-CoV-2/GFP wt 2E1 or plaque-purified selected variants for 1 hr at 37°C. The nanobody/recombinant virus mixture was then added to 293T/ACE2.cl22 cells. After 16 hr, cells were harvested, and GFP-positive cells quantified by flow cytometry. The percentage of GFP-positive cells was normalized to that derived from cells infected with rVSV/SARS-CoV-2 in the absence of nanobodies. The half-maximal inhibitory concentrations (IC50) for the nanobodies were determined using four-parameter nonlinear regression (leastsquares regression method without weighting) (GraphPad Prism). To identify putative nanobody resistance mutations, RNA was isolated from aliquots of supernatant containing selected viral populations or individual plaque purified variants using NucleoSpin 96 Virus Core Kit (Macherey-Nagel). The purified RNA was subjected to reverse transcription using random hexamer primers and SuperScript VILO cDNA Synthesis Kit (Thermo Fisher Scientific). The cDNA was amplified using KOD Xtreme Hot Start DNA 396 Polymerase (MilliporeSigma) flanking the spike encoding sequences. The PCR products were gel-purified and sequenced using Sanger sequencing. For selection of spike variants that were resistant to nanobodies, rVSV/SARS-CoV-2/GFP wt 2E1 was passaged to generate diversity and populations containing 10 6 infectious particles were used. The rVSV/SARS-CoV-2/GFP wt 2E1 populations were incubated with dilutions of nanobodies (10× to 100× the IC50 excess) for 1 hr at 37°C. Then, the virus-nanobody mixtures were incubated with 5 × 10 5 293T/ACE2.22 cells in six-well plates. Two days later, the cells were imaged and supernatant were harvested from cultures that showed evidence of viral replication (GFP-positive foci) or large numbers of GFP-positive cells. A 100 μl of the cleared supernatant was incubated with the same dilution of nanobodies and then used to infect 5 × 10 5 293T/ACE2.22 cells in six-well plates, as before. rVSV/ SARS-CoV-2/GFP wt 2E1 were passaged in the present combination of nanobodies two times before complete escape was evaluated. To isolate individual mutant viruses, selected rVSV/SARS-CoV-2/GFP wt 2E1 populations were serially diluted in medium without nanobodies and individual viral variants isolated by visualizing single GFPpositive plaques at limiting dilutions in 96-well plates containing 1 × 10 4 293T/ACE2.22 cells. These plaque-purified viruses were expanded and further characterized using sequencing and nanobody neutralization assays. Nanobodies and antigens were incubated together at a 2× molar excess of nanobody at RT for 10 min in 20 mM HEPES pH 7.4 and 150 mM NaCl. Cross-linker was then added to a final concentration of 5 mM bissulfosuccinimidyl suberate (BS3) or 1 mM disuccinimidyl suberate (DSS), and samples were cross-linked for 30 min (RBD, NTD) or 18 min (ectodomain trimer) at RT. Reactions were quenched, reduced, and alkylated, and run on an SDS-PAGE gel. The band corresponding to the cross-linked nanobody-antigen complex was then excised from the gel and subjected to in-gel digestion at 37°C with trypsin (Roche, 1 μg, 4 hr) or chymotrypsin (Roche, 0.5 μg, 1.5 hr). Peptides were extracted and analyzed with a nano-LC 1200 (Thermo Fisher) with an EASYspray PepMap RSLC C18 3 µm, 100 Å, 75 µm × 15 cm column coupled to an Orbitrap Fusion Lumos Tribrid mass spectrometer (Thermo Fisher). An ABIRD (ESI Source Solutions) was used to reduce background. The Lumos was operated in a data-dependent mode, and ions were fragmented by high-energy collisional dissociation (normalized collision energy 28). Separate LC runs were used to analyze the +3 and the +4 through +7 charge states. Higher charge species were prioritized for selection for fragmentation when analyzing the 4-7 species. Orbitrap resolution was 30,000 for MS and 15,000 for MS/MS analyses. The quadrupole isolation window was 1.4, and the MS/MS used a maximum injection time of 800 ms with four microscans. Data were then searched by pLink 2.3 (Chen et al., 2019) to identify cross-linked peptides. The mass accuracy in pLink was set to 10 ppm for MS and 20 ppm for MS/MS. Cysteine carbamidomethylation was included as a fixed modification and methionine oxidation as a variable modification. For trypsin, up to three missed cleavages were permitted. For chymotrypsin, the enzyme setting was 'nonspecific.' Spectra were manually checked to ensure correct identifications of cross-linked peptides. Human subjects: Bronchial airway epithelial cells were obtained under study #12490 approved by the Seattle Children's Institutional Review Board, with investigations carried out following the rules of the Declaration of Helsinki of 1975. Removal of endotoxin from protein solutions by phase separation using Triton X-114 Extremely potent human monoclonal antibodies from COVID-19 convalescent patients Nanobodies and nanobody-based human heavy chain antibodies as antitumor therapeutics Spelling-Error Tolerant, Order-Independent Pass-Phrases via the Damerau-Levenshtein String-Edit Distance Metric SARS-CoV-2 neutralizing antibody structures inform therapeutic strategies Inactivation of material from SARS-CoV-2-infected primary airway epithelial cell cultures Antibody cocktail to SARS-CoV-2 spike protein prevents rapid mutational escape seen with individual antibodies Receptor binding and priming of the spike protein of SARS-CoV-2 for membrane fusion Bi-paratopic and multivalent VH domains block ACE2 binding and neutralize SARS-CoV-2 Distinct conformational states of SARS-CoV-2 spike protein Increased transmissibility and global spread of SARS-CoV-2 variants of concern as at Potent Neutralizing Antibodies against SARS-CoV-2 identified by high-throughput single-cell sequencing of convalescent patients B cells Potent antibody therapeutics by design Beyond shielding: the roles of glycans in the SARS-CoV-2 spike protein Growth, detection, quantification, and inactivation of SARS-CoV-2 Protocol and reagents for pseudotyping lentiviral particles with SARS-CoV-2 spike protein for neutralization assays Engineering human ACE2 to optimize binding to the spike protein of SARS coronavirus 2 Nanobody Engineering: Toward Next Generation Immunotherapies and Immunoimaging of Cancer A high-speed search engine pLink 2 with systematic evaluation for proteomescale identification of cross-linked peptides Relationship between the inhibition constant (K1) and the concentration of inhibitor which causes 50 per cent inhibition (I50) of an enzymatic reaction A neutralizing human antibody binds to the N-terminal domain of the Spike protein of SARS-CoV-2 Betalactamase inhibitors derived from single-domain antibody fragments elicited in the camelidae Tackling COVID-19 with neutralizing monoclonal antibodies 2020. Selection, biophysical and structural analysis of synthetic nanobodies that effectively neutralize SARS-CoV-2 Bispecific IgG neutralizes SARS-CoV-2 variants and prevents escape in mice The antigenic anatomy of SARS-CoV-2 receptor binding domain SARS-CoV-2 variants show resistance to neutralization by many monoclonal and serum-derived polyclonal antibodies Development of multi-specific humanized llama antibodies blocking SARS-CoV-2/ACE2 interaction with high affinity and avidity Caplacizumab: First Global Approval Human monoclonal antibodies against highly conserved HR1 and HR2 domains of the SARS-CoV spike protein are more broadly neutralizing Nanobodies as versatile tool for multiscale imaging modalities High affinity nanobodies block SARS-CoV-2 spike receptor binding domain interaction with human angiotensin converting enzyme High throughput solution-based measurement of antibody-antigen affinity and epitope binning Structural analysis of neutralizing epitopes of the SARS-CoV-2 spike to guide therapy and vaccine design strategies Modeling of loops in protein structures Factors that make an infectious disease outbreak controllable A robust pipeline for rapid production of versatile nanobody repertoires Multiple SARS-CoV-2 variants escape neutralization by vaccine-induced humoral immunity High-resolution profiling of pathways of escape for SARS-CoV-2 spike-binding antibodies Analysis of the SARS-CoV-2 spike protein glycan shield reveals implications for immune recognition Comprehensive mapping of mutations in the SARS-CoV-2 receptor-binding domain that affect recognition by polyclonal human plasma antibodies A structurally distinct human mycoplasma protein that generically blocks antigen-antibody union An alpaca nanobody neutralizes SARS-CoV-2 by blocking receptor interaction The molecular virology of coronaviruses Increased Resistance of SARS-CoV-2 Variants B.1.351 and B.1.1.7 to Antibody Neutralization Structure-based design of prefusion-stabilized SARS-CoV-2 spikes Door to the cell for COVID-19 opened, leading way to therapies Matplotlib: A 2D Graphics Environment Neutralizing nanobodies bind SARS-CoV-2 spike RBD and block interaction with ACE2 Neutralization of SARS-CoV-2 by Destruction of the Prefusion Spike SARS-CoV-2 spike E484K mutation reduces antibody neutralisation The Therapeutic Potential of Nanobodies Integrative structure and functional anatomy of a nuclear pore complex Structure-guided multivalent nanobodies block SARS-CoV-2 infection and suppress mutational escape Tracking Changes in SARS-CoV-2 Spike: Evidence that D614G Increases Infectivity of the COVID-19 Virus Structure of the SARS-CoV-2 spike receptor-binding domain bound to the ACE2 receptor Immunological characteristics govern the transition of COVID-19 to endemicity Functional assessment of cell entry and receptor usage for SARS-CoV-2 and other lineage B betacoronaviruses Angiotensin-converting enzyme 2 is a functional receptor for the SARS coronavirus Structure, Function, and Evolution of Coronavirus Spike Proteins Linear epitopes of SARS-CoV-2 spike protein elicit neutralizing antibodies in COVID-19 patients The Architecture of Inactivated SARS-CoV-2 with Postfusion Spikes Revealed by Cryo-EM and Cryo-ET Potent neutralizing antibodies against multiple epitopes on SARS-CoV-2 spike Identification of SARS-CoV-2 spike mutations that attenuate monoclonal and serum antibody neutralization 501Y.V2 and 501Y.V3 variants of SARS-CoV-2 lose binding to bamlanivimab in vitro Potent Neutralization of SARS-CoV-2 by Hetero-Bivalent Alpaca Nanobodies Targeting the Spike Receptor-Binding Domain Virus neutralisation: new insights from kinetic neutralisation curves N-terminal domain antigenic mapping reveals a site of vulnerability for SARS-CoV-2 A general approach to antibody thermostabilization Vaccines Need Not Completely Stop COVID Transmission to Curb the Pandemic Data Structures for Statistical Computing in Python. Python in Science Conference Two Detailed Plaque Assay Protocols for the Quantification of Infectious SARS-CoV-2 Quantifying Drug Combination Synergy along Potency and Efficacy Axes Affinity maturation of SARS-CoV-2 neutralizing antibodies confers potency, breadth, and resilience to viral escape mutations Effects of Asthma and Human Rhinovirus A16 on the Expression of SARS-CoV-2 Entry Factors in Human Airway Epithelium Nanobodies: natural single-domain antibodies Inhalable Nanobody (PiN-21) prevents and treats SARS-CoV-2 infections in Syrian hamsters at ultra-low doses The use of differential scanning fluorimetry to detect ligand interactions that promote protein stability A Structural Landscape of Neutralizing Antibodies Against SARS-CoV-2 Receptor Binding Domain A Robust Workflow for Native Mass Spectrometric Analysis of Affinity-Isolated Endogenous Protein Assemblies Native mass spectrometry-based screening for optimal sample preparation in single-particle Cryo-EM COVID vaccines to reach poorest countries in 2023 -despite recent pledges Scikit-learn: Machine Learning in Python Emerging technologies in point-of-care molecular diagnostics for resourcelimited settings UCSF ChimeraX: Structure visualization for researchers, educators, and developers The coronavirus is here to stay -here's what that means SARS-CoV-2 spike protein: flexibility as a new target for fighting infection Two linear epitopes on the SARS-CoV-2 spike protein that elicit neutralising antibodies in COVID-19 patients Nanobody cocktails potently neutralize SARS-CoV-2 D614G N501Y variant and protect mice Evaluating MERS-CoV Entry Pathways Cross-linking/mass spectrometry as a new field and the proteomics information mountain of tomorrow CoV-AbDab: the coronavirus antibody database Nanobodies as novel agents for cancer therapy Comparison of Four SARS-CoV-2 Neutralization Assays Clustering biomolecular complexes by residue contacts similarity Isolation of potent SARS-CoV-2 neutralizing antibodies and protection from disease in a small animal model Principles for Integrative Structural Biology Studies Putting the pieces together: integrative modeling platform software for structure determination of macromolecular assemblies Comparative protein modelling by satisfaction of spatial restraints Using integrative modeling platform to compute, validate, and archive a model of a protein complex structure Integrative models of neutralizing nanobody epitopes on the SARS-CoV-2 spike protein. 2607a97 Measuring SARS-CoV-2 neutralizing antibody activity using pseudotyped and chimeric viruses An ultrapotent synthetic nanobody neutralizes SARS-CoV-2 by stabilizing inactive Spike The Case for S2: The Potential Benefits of the S2 Subunit of the SARS-CoV-2 Spike Protein as an Immunogen in Fighting the COVID-19 Pandemic Deducing the N-and O-glycosylation profile of the spike protein of novel coronavirus SARS-CoV-2 Structural basis of receptor recognition by SARS-CoV-2 A resource of high-quality and versatile nanobodies for drug delivery A strategy for dissecting the architectures of native macromolecular assemblies Quantifying proteinprotein interactions by molecular counting with mass photometry Cross-reactive serum and memory B-cell responses to spike protein in SARS-CoV-2 and endemic coronavirus infection 2021. mRNA vaccination boosts cross-variant neutralizing antibodies elicited by SARS-CoV-2 infection Deep Mutational Scanning of SARS-CoV-2 Receptor Binding Domain Reveals Constraints on Folding and ACE2 Binding Prospective mapping of viral mutations that escape antibodies used to treat COVID-19 Kinetic rates of antibody binding correlate with neutralization sensitivity of variant simian immunodeficiency virus strains When will it end?': How a changing virus is reshaping scientists' views on COVID-19 Cryo-EM visualization of an exposed RGD epitope on adenovirus that escapes antibody neutralization Potent neutralizing nanobodies resist convergent circulating variants of SARS-CoV-2 by targeting diverse and conserved epitopes Detection of a SARS-CoV-2 variant of concern in South Africa Relationship between poliovirus neutralization and aggregation Optimizing selection of large animals for antibody production by screening immune response to standard vaccines Potent binding of 2019 novel coronavirus spike protein by a SARS coronavirus-specific human monoclonal antibody Estimating the number of clusters in a data set via the gap statistic Improved tumor targeting of anti-epidermal growth factor receptor Nanobodies through albumin binding: taking advantage of modular Nanobody technology Ultrapotent human antibodies protect against SARS-CoV-2 challenge via multiple mechanisms Unlocking COVID therapeutic targets: A structure-based rationale against SARS-CoV-2, SARS-CoV and MERS-CoV Spike In situ structural analysis of SARS-CoV-2 spike reveals flexibility mediated by three hinges BRAID: A unifying paradigm for the analysis of combined drug action Nanobodies as inhaled biotherapeutics for lung diseases Python 3 Reference Manual. CreateSpace Dissociation rate of antibody-gp120 binding interactions is predictive of V3-mediated neutralization of HIV-1 General strategy to humanize a camelid single-domain antibody and identification of a universal humanized nanobody scaffold SciPy 1.0: fundamental algorithms for scientific computing in Python IdeS, a novel streptococcal cysteine proteinase with unique specificity for immunoglobulin G Unexpected receptor functional mimicry elucidates activation of coronavirus fusion Structure, Function, and Antigenicity of the SARS-CoV-2 Spike Glycoprotein Antibody resistance of SARS-CoV-2 variants B.1.351 and B.1.1.7 mRNA vaccine-elicited antibodies to SARS-CoV-2 and circulating variants A conserved immunogenic and vulnerable site on the coronavirus spike protein delineated by cross-reactive monoclonal antibodies 2021. seaborn: statistical data visualization Site-specific glycan analysis of the SARS-CoV-2 spike Integrative structure modeling with the Integrative Modeling Platform Protein structure modeling with modeller Escape from neutralizing antibodies by SARS-CoV-2 spike protein variants Structural evaluation of a nanobody targeting complement receptor Vsig4 and its cross reactivity SARS-CoV-2 501Y.V2 escapes neutralization by South African COVID-19 donor plasma Virological assessment of hospitalized patients with COVID-2019 synergy: a Python library for calculating, analyzing and visualizing drug combination synergy Structural basis for potent neutralization of betacoronaviruses by single-domain camelid antibodies Cryo-EM structure of the 2019-nCoV spike in the prefusion conformation A noncompeting pair of human neutralizing antibodies block COVID-19 virus binding to its receptor ACE2 Rapid Determination of Antibody-Antigen Affinity by Mass Photometry Versatile and multivalent nanobodies efficiently neutralize SARS-CoV-2 Integrative proteomics identifies thousands of distinct, multi-epitope, and high-affinity nanobodies Nanobodies from camelid mice and llamas neutralize SARS-CoV-2 variants The development of Nanosota-1 as anti-SARS-CoV-2 nanobody drug candidates Quantitative mass imaging of single biological macromolecules Nanobodies, the potent agents to detect and treat the Coronavirus infections: A systematic review Site-specific N-glycosylation Characterization of Recombinant SARS-CoV-2 Spike Proteins A pneumonia outbreak associated with a new coronavirus of probable bat origin Prediction of multidimensional drug dose responses based on measurements of drug pairs Potently neutralizing and protective human antibodies against SARS-CoV-2 We are very grateful to Sonya Paske and the staff at Capralogics, Inc for raising antisera; the High-Throughput and Spectroscopy Resource Center (HTSRC), Rockefeller University, and especially Lavoisier Ramos-Espiritu for technical and data analysis support; Brenda Watt from Refeyn Ltd for technical and data analysis support; Jason Wendler and Maxwell Neal for help with statistical analyses; Thomas Gallagher and Enya Qing, Loyola University Chicago, for the kind gift of the SARS-CoV-1 and SARS-CoV-2 spike expression plasmids; Andrew McGuire and Leonidas Stamatatos (Fred Hutchinson Cancer Research Center) for their kind gift of the SARS-CoV-2 20H/501Y.V2 spike variant expression plasmid; Rhea Coler (Seattle Children's Research Institute) for the kind gift of the SARS-CoV-2 delta variant and the TMPRSS2+ Vero E6 cells; and the rest of the Aitchison, Chait and Rout laboratories, as well as Fred Cross for technical and intellectual support. The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication. Biological sample (Lama glama) From two male llamas immunized with SARS-CoV-2 spike S1, RBD, and S2.Biological sample (L. From two male llamas immunized with SARS-CoV-2 spike S1, RBD, and S2.