Delineating the functional activity of antibodies with cross-reactivity to SARS-CoV-2, SARS-CoV-1 and related sarbecoviruses.

The recurring spillover of pathogenic coronaviruses and demonstrated capacity of sarbecoviruses, such SARS-CoV-2, to rapidly evolve in humans underscores the need to better understand immune responses to this virus family. For this purpose, we characterized the functional breadth and potency of antibodies targeting the receptor binding domain (RBD) of the spike glycoprotein that exhibited cross-reactivity against SARS-CoV-2 variants, SARS-CoV-1 and sarbecoviruses from diverse clades and animal origins with spillover potential. One neutralizing antibody, C68.61, showed remarkable neutralization breadth against both SARS-CoV-2 variants and viruses from different sarbecovirus clades. C68.61, which targets a conserved RBD class 5 epitope, did not select for escape variants of SARS-CoV-2 or SARS-CoV-1 in culture nor have predicted escape variants among circulating SARS-CoV-2 strains, suggesting this epitope is functionally constrained. We identified 11 additional SARS-CoV-2/SARS-CoV-1 cross-reactive antibodies that target the more sequence conserved class 4 and class 5 epitopes within RBD that show activity against a subset of diverse sarbecoviruses with one antibody binding every single sarbecovirus RBD tested. A subset of these antibodies exhibited Fc-mediated effector functions as potent as antibodies that impact infection outcome in animal models. Thus, our study identified antibodies targeting conserved regions across SARS-CoV-2 variants and sarbecoviruses that may serve as therapeutics for pandemic preparedness as well as blueprints for the design of immunogens capable of eliciting cross-neutralizing responses.

Broadly neutralizing antibodies against emerging delta-coronaviruses.

Porcine deltacoronavirus (PDCoV) spillovers were recently detected in children with acute undifferentiated febrile illness, underscoring recurrent zoonoses of divergent coronaviruses. To date, no vaccines or specific therapeutics are approved for use in humans against PDCoV. To prepare for possible future PDCoV epidemics, we isolated human spike (S)-directed monoclonal antibodies from transgenic mice and found that two of them, designated PD33 and PD41, broadly neutralized a panel of PDCoV variants. Cryo-electron microscopy structures of PD33 and PD41 in complex with the PDCoV receptor-binding domain and S ectodomain trimer provide a blueprint of the epitopes recognized by these mAbs, rationalizing their broad inhibitory activity. We show that both mAbs inhibit PDCoV by competitively interfering with host APN binding to the PDCoV receptor-binding loops, explaining the mechanism of viral neutralization. PD33 and PD41 are candidates for clinical advancement, which could be stockpiled to prepare for possible future PDCoV outbreaks.

Deep mutational scans of XBB.1.5 and BQ.1.1 reveal ongoing epistatic drift during SARS-CoV-2 evolution.

Substitutions that fix between SARS-CoV-2 variants can transform the mutational landscape of future evolution via epistasis. For example, large epistatic shifts in mutational effects caused by N501Y underlied the original emergence of Omicron, but whether such epistatic saltations continue to define ongoing SARS-CoV-2 evolution remains unclear. We conducted deep mutational scans to measure the impacts of all single amino acid mutations and single-codon deletions in the spike receptor-binding domain (RBD) on ACE2-binding affinity and protein expression in the recent Omicron BQ.1.1 and XBB.1.5 variants, and we compared mutational patterns to earlier viral strains that we have previously profiled. As with previous deep mutational scans, we find many mutations that are tolerated or even enhance binding to ACE2 receptor. The tolerance of sites to single-codon deletion largely conforms with tolerance to amino acid mutation. Though deletions in the RBD have not yet been seen in dominant lineages, we observe tolerated deletions including at positions that exhibit indel variation across broader sarbecovirus evolution and in emerging SARS-CoV-2 variants of interest, most notably the well-tolerated Δ483 deletion in BA.2.86. The substitutions that distinguish recent viral variants have not induced as dramatic of epistatic perturbations as N501Y, but we identify ongoing epistatic drift in SARS-CoV-2 variants, including interaction between R493Q reversions and mutations at positions 453, 455, and 456, including F456L that defines the XBB.1.5-derived EG.5 lineage. Our results highlight ongoing drift in the effects of mutations due to epistasis, which may continue to direct SARS-CoV-2 evolution into new regions of sequence space.

Broad receptor tropism and immunogenicity of a clade 3 sarbecovirus.

Although Rhinolophus bats harbor diverse clade 3 sarbecoviruses, the structural determinants of receptor tropism along with the antigenicity of their spike (S) glycoproteins remain uncharacterized. Here, we show that the African Rhinolophus bat clade 3 sarbecovirus PRD-0038 S has a broad angiotensin-converting enzyme 2 (ACE2) usage and that receptor-binding domain (RBD) mutations further expand receptor promiscuity and enable human ACE2 utilization. We determine a cryo-EM structure of the PRD-0038 RBD bound to Rhinolophus alcyone ACE2, explaining receptor tropism and highlighting differences with SARS-CoV-1 and SARS-CoV-2. Characterization of PRD-0038 S using cryo-EM and monoclonal antibody reactivity reveals its distinct antigenicity relative to SARS-CoV-2 and identifies PRD-0038 cross-neutralizing antibodies for pandemic preparedness. PRD-0038 S vaccination elicits greater titers of antibodies cross-reacting with vaccine-mismatched clade 2 and clade 1a sarbecoviruses compared with SARS-CoV-2 S due to broader antigenic targeting, motivating the inclusion of clade 3 antigens in next-generation vaccines for enhanced resilience to viral evolution.

Deep mutational scans of XBB.1.5 and BQ.1.1 reveal ongoing epistatic drift during SARS-CoV-2 evolution.

Substitutions that fix between SARS-CoV-2 variants can transform the mutational landscape of future evolution via epistasis. For example, large epistatic shifts in mutational effects caused by N501Y underlied the original emergence of Omicron variants, but whether such large epistatic saltations continue to define ongoing SARS-CoV-2 evolution remains unclear. We conducted deep mutational scans to measure the impacts of all single amino acid mutations and single-codon deletions in the spike receptor-binding domain (RBD) on ACE2-binding affinity and protein expression in the recent Omicron BQ.1.1 and XBB.1.5 variants, and we compared mutational patterns to earlier viral strains that we have previously profiled. As with previous RBD deep mutational scans, we find many mutations that are tolerated or even enhance binding to ACE2 receptor. The tolerance of sites to single-codon deletion largely conforms with tolerance to amino acid mutation. Though deletions in the RBD have not yet been seen in dominant lineages, we observe many tolerated deletions including at positions that exhibit indel variation across broader sarbecovirus evolution and in emerging SARS-CoV-2 variants of interest, most notably the well-tolerated Δ483 deletion in BA.2.86. The substitutions that distinguish recent viral variants have not induced as dramatic of epistatic perturbations as N501Y, but we identify ongoing epistatic drift in SARS-CoV-2 variants, including interaction between R493Q reversions and mutations at positions 453, 455, and 456, including mutations like F456L that define the newly emerging EG.5 lineage. Our results highlight ongoing drift in the effects of mutations due to epistasis, which may continue to direct SARS-CoV-2 evolution into new regions of sequence space.

Broad receptor tropism and immunogenicity of a clade 3 sarbecovirus.

Although Rhinolophus bats harbor diverse clade 3 sarbecoviruses, the structural determinants of receptor tropism along with the antigenicity of their spike (S) glycoproteins remain uncharacterized. Here, we show that the African Rinolophus bat clade 3 sarbecovirus PRD-0038 S has a broad ACE2 usage and that RBD mutations further expand receptor promiscuity and enable human ACE2 utilization. We determined a cryoEM structure of the PRD-0038 RBD bound to R. alcyone ACE2, explaining receptor tropism and highlighting differences with SARS-CoV-1 and SARS-CoV-2. Characterization of PRD-0038 S using cryoEM and monoclonal antibody reactivity revealed its distinct antigenicity relative to SARS-CoV-2 and identified PRD-0038 cross-neutralizing antibodies for pandemic preparedness. PRD-0038 S vaccination elicited greater titers of antibodies cross-reacting with vaccine-mismatched clade 2 and clade 1a sarbecoviruses compared to SARS-CoV-2 S due to broader antigenic targeting, motivating the inclusion of clade 3 antigens in next-generation vaccines for enhanced resilience to viral evolution.

Potent pan huACE2-dependent sarbecovirus neutralizing monoclonal antibodies isolated from a BNT162b2-vaccinated SARS survivor.

The emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concern such as Omicron hampered efforts in controlling the ongoing coronavirus disease 2019 pandemic due to their ability to escape neutralizing antibodies induced by vaccination or prior infection, highlighting the need to develop broad-spectrum vaccines and therapeutics. Most human monoclonal antibodies (mAbs) reported to date have not demonstrated true pan-sarbecovirus neutralizing breadth especially against animal sarbecoviruses. Here, we report the isolation and characterization of highly potent mAbs targeting the receptor binding domain (RBD) of huACE2-dependent sarbecovirus from a SARS-CoV survivor vaccinated with BNT162b2. Among the six mAbs identified, one (E7) showed better huACE2-dependent sarbecovirus neutralizing potency and breadth than any other mAbs reported to date. Mutagenesis and cryo-electron microscopy studies indicate that these mAbs have a unique RBD contact footprint and that E7 binds to a quaternary structure-dependent epitope.

Evolution of antibody immunity following Omicron BA.1 breakthrough infection.

Understanding the longitudinal dynamics of antibody immunity following heterologous SAR-CoV-2 breakthrough infection will inform the development of next-generation vaccines. Here, we track SARS-CoV-2 receptor binding domain (RBD)-specific antibody responses up to six months following Omicron BA.1 breakthrough infection in six mRNA-vaccinated individuals. Cross-reactive serum neutralizing antibody and memory B cell (MBC) responses decline by two- to four-fold through the study period. Breakthrough infection elicits minimal de novo Omicron BA.1-specific B cell responses but drives affinity maturation of pre-existing cross-reactive MBCs toward BA.1, which translates into enhanced breadth of activity across other variants. Public clones dominate the neutralizing antibody response at both early and late time points following breakthough infection, and their escape mutation profiles predict newly emergent Omicron sublineages, suggesting that convergent antibody responses continue to shape SARS-CoV-2 evolution. While the study is limited by our relatively small cohort size, these results suggest that heterologous SARS-CoV-2 variant exposure drives the evolution of B cell memory, supporting the continued development of next-generation variant-based vaccines.

Identification of broad, potent antibodies to functionally constrained regions of SARS-CoV-2 spike following a breakthrough infection.

The antiviral benefit of antibodies can be compromised by viral escape especially for rapidly evolving viruses. Therefore, durable, effective antibodies must be both broad and potent to counter newly emerging, diverse strains. Discovery of such antibodies is critically important for SARS-CoV-2 as the global emergence of new variants of concern (VOC) has compromised the efficacy of therapeutic antibodies and vaccines. We describe a collection of broad and potent neutralizing monoclonal antibodies (mAbs) isolated from an individual who experienced a breakthrough infection with the Delta VOC. Four mAbs potently neutralize the Wuhan-Hu-1 vaccine strain, the Delta VOC, and also retain potency against the Omicron VOCs through BA.4/BA.5 in both pseudovirus-based and authentic virus assays. Three mAbs also retain potency to recently circulating VOCs XBB.1.5 and BQ.1.1 and one also potently neutralizes SARS-CoV-1. The potency of these mAbs was greater against Omicron VOCs than all but one of the mAbs that had been approved for therapeutic applications. The mAbs target distinct epitopes on the spike glycoprotein, three in the receptor-binding domain (RBD) and one in an invariant region downstream of the RBD in subdomain 1 (SD1). The escape pathways we defined at single amino acid resolution with deep mutational scanning show they target conserved, functionally constrained regions of the glycoprotein, suggesting escape could incur a fitness cost. Overall, these mAbs are unique in their breadth across VOCs, their epitope specificity, and include a highly potent mAb targeting a rare epitope outside of the RBD in SD1.

The landscape of antibody binding affinity in SARS-CoV-2 Omicron BA.1 evolution.

The Omicron BA.1 variant of SARS-CoV-2 escapes convalescent sera and monoclonal antibodies that are effective against earlier strains of the virus. This immune evasion is largely a consequence of mutations in the BA.1 receptor binding domain (RBD), the major antigenic target of SARS-CoV-2. Previous studies have identified several key RBD mutations leading to escape from most antibodies. However, little is known about how these escape mutations interact with each other and with other mutations in the RBD. Here, we systematically map these interactions by measuring the binding affinity of all possible combinations of these 15 RBD mutations (215=32,768 genotypes) to 4 monoclonal antibodies (LY-CoV016, LY-CoV555, REGN10987, and S309) with distinct epitopes. We find that BA.1 can lose affinity to diverse antibodies by acquiring a few large-effect mutations and can reduce affinity to others through several small-effect mutations. However, our results also reveal alternative pathways to antibody escape that does not include every large-effect mutation. Moreover, epistatic interactions are shown to constrain affinity decline in S309 but only modestly shape the affinity landscapes of other antibodies. Together with previous work on the ACE2 affinity landscape, our results suggest that the escape of each antibody is mediated by distinct groups of mutations, whose deleterious effects on ACE2 affinity are compensated by another distinct group of mutations (most notably Q498R and N501Y).

Molecular fate-mapping of serum antibody responses to repeat immunization.

The protective efficacy of serum antibodies results from the interplay of antigen-specific B cell clones of different affinities and specificities. These cellular dynamics underlie serum-level phenomena such as original antigenic sin (OAS)-a proposed propensity of the immune system to rely repeatedly on the first cohort of B cells engaged by an antigenic stimulus when encountering related antigens, in detriment to the induction of de novo responses1-5. OAS-type suppression of new, variant-specific antibodies may pose a barrier to vaccination against rapidly evolving viruses such as influenza and SARS-CoV-26,7. Precise measurement of OAS-type suppression is challenging because cellular and temporal origins cannot readily be ascribed to antibodies in circulation; its effect on subsequent antibody responses therefore remains unclear5,8. Here we introduce a molecular fate-mapping approach with which serum antibodies derived from specific cohorts of B cells can be differentially detected. We show that serum responses to sequential homologous boosting derive overwhelmingly from primary cohort B cells, while later induction of new antibody responses from naive B cells is strongly suppressed. Such 'primary addiction' decreases sharply as a function of antigenic distance, allowing reimmunization with divergent viral glycoproteins to produce de novo antibody responses targeting epitopes that are absent from the priming variant. Our findings have implications for the understanding of OAS and for the design and testing of vaccines against evolving pathogens.

Identification of broad, potent antibodies to functionally constrained regions of SARS-CoV-2 spike following a breakthrough infection.

The antiviral benefit of antibodies can be compromised by viral escape especially for rapidly evolving viruses. Therefore, durable, effective antibodies must be both broad and potent to counter newly emerging, diverse strains. Discovery of such antibodies is critically important for SARS-CoV-2 as the global emergence of new variants of concern (VOC) has compromised the efficacy of therapeutic antibodies and vaccines. We describe a collection of broad and potent neutralizing monoclonal antibodies (mAbs) isolated from an individual who experienced a breakthrough infection with the Delta VOC. Four mAbs potently neutralize the Wuhan-Hu-1 vaccine strain, the Delta VOC, and also retain potency against the Omicron VOCs through BA.4/BA.5 in both pseudovirus-based and authentic virus assays. Three mAbs also retain potency to recently circulating VOCs XBB.1.5 and BQ.1.1 and one also potently neutralizes SARS-CoV-1. The potency of these mAbs was greater against Omicron VOCs than all but one of the mAbs that had been approved for therapeutic applications. The mAbs target distinct epitopes on the spike glycoprotein, three in the receptor binding domain (RBD) and one in an invariant region downstream of the RBD in subdomain 1 (SD1). The escape pathways we defined at single amino acid resolution with deep mutational scanning show they target conserved, functionally constrained regions of the glycoprotein, suggesting escape could incur a fitness cost. Overall, these mAbs are novel in their breadth across VOCs, their epitope specificity, and include a highly potent mAb targeting a rare epitope outside of the RBD in SD1.

Compensatory epistasis maintains ACE2 affinity in SARS-CoV-2 Omicron BA.1.

The Omicron BA.1 variant emerged in late 2021 and quickly spread across the world. Compared to the earlier SARS-CoV-2 variants, BA.1 has many mutations, some of which are known to enable antibody escape. Many of these antibody-escape mutations individually decrease the spike receptor-binding domain (RBD) affinity for ACE2, but BA.1 still binds ACE2 with high affinity. The fitness and evolution of the BA.1 lineage is therefore driven by the combined effects of numerous mutations. Here, we systematically map the epistatic interactions between the 15 mutations in the RBD of BA.1 relative to the Wuhan Hu-1 strain. Specifically, we measure the ACE2 affinity of all possible combinations of these 15 mutations (215 = 32,768 genotypes), spanning all possible evolutionary intermediates from the ancestral Wuhan Hu-1 strain to BA.1. We find that immune escape mutations in BA.1 individually reduce ACE2 affinity but are compensated by epistatic interactions with other affinity-enhancing mutations, including Q498R and N501Y. Thus, the ability of BA.1 to evade immunity while maintaining ACE2 affinity is contingent on acquiring multiple interacting mutations. Our results implicate compensatory epistasis as a key factor driving substantial evolutionary change for SARS-CoV-2 and are consistent with Omicron BA.1 arising from a chronic infection.

Deep mutational scans for ACE2 binding, RBD expression, and antibody escape in the SARS-CoV-2 Omicron BA.1 and BA.2 receptor-binding domains.

SARS-CoV-2 continues to acquire mutations in the spike receptor-binding domain (RBD) that impact ACE2 receptor binding, folding stability, and antibody recognition. Deep mutational scanning prospectively characterizes the impacts of mutations on these biochemical properties, enabling rapid assessment of new mutations seen during viral surveillance. However, the effects of mutations can change as the virus evolves, requiring updated deep mutational scans. We determined the impacts of all single amino acid mutations in the Omicron BA.1 and BA.2 RBDs on ACE2-binding affinity, RBD folding, and escape from binding by the LY-CoV1404 (bebtelovimab) monoclonal antibody. The effects of some mutations in Omicron RBDs differ from those measured in the ancestral Wuhan-Hu-1 background. These epistatic shifts largely resemble those previously seen in the Alpha variant due to the convergent epistatically modifying N501Y substitution. However, Omicron variants show additional lineage-specific shifts, including examples of the epistatic phenomenon of entrenchment that causes the Q498R and N501Y substitutions present in Omicron to be more favorable in that background than in earlier viral strains. In contrast, the Omicron substitution Q493R exhibits no sign of entrenchment, with the derived state, R493, being as unfavorable for ACE2 binding in Omicron RBDs as in Wuhan-Hu-1. Likely for this reason, the R493Q reversion has occurred in Omicron sub-variants including BA.4/BA.5 and BA.2.75, where the affinity buffer from R493Q reversion may potentiate concurrent antigenic change. Consistent with prior studies, we find that Omicron RBDs have reduced expression, and identify candidate stabilizing mutations that ameliorate this deficit. Last, our maps highlight a broadening of the sites of escape from LY-CoV1404 antibody binding in BA.1 and BA.2 compared to the ancestral Wuhan-Hu-1 background. These BA.1 and BA.2 deep mutational scanning datasets identify shifts in the RBD mutational landscape and inform ongoing efforts in viral surveillance.

Molecular fate-mapping of serum antibodies reveals the effects of antigenic imprinting on repeated immunization.

The ability of serum antibody to protect against pathogens arises from the interplay of antigen-specific B cell clones of different affinities and fine specificities. These cellular dynamics are ultimately responsible for serum-level phenomena such as antibody imprinting or "Original Antigenic Sin" (OAS), a proposed propensity of the immune system to rely repeatedly on the first cohort of B cells that responded to a stimulus upon exposure to related antigens. Imprinting/OAS is thought to pose a barrier to vaccination against rapidly evolving viruses such as influenza and SARS-CoV-2. Precise measurement of the extent to which imprinting/OAS inhibits the recruitment of new B cell clones by boosting is challenging because cellular and temporal origins cannot readily be assigned to antibodies in circulation. Thus, the extent to which imprinting/OAS impacts the induction of new responses in various settings remains unclear. To address this, we developed a "molecular fate-mapping" approach in which serum antibodies derived from specific cohorts of B cells can be differentially detected. We show that, upon sequential homologous boosting, the serum antibody response strongly favors reuse of the first cohort of B cell clones over the recruitment of new, naÏve-derived B cells. This "primary addiction" decreases as a function of antigenic distance, allowing secondary immunization with divergent influenza virus or SARS-CoV-2 glycoproteins to overcome imprinting/OAS by targeting novel epitopes absent from the priming variant. Our findings have implications for the understanding of imprinting/OAS, and for the design and testing of vaccines aimed at eliciting antibodies to evolving antigens.

Evolution of antibody immunity following Omicron BA.1 breakthrough infection.

Understanding the evolution of antibody immunity following heterologous SAR-CoV-2 breakthrough infection will inform the development of next-generation vaccines. Here, we tracked SARS-CoV-2 receptor binding domain (RBD)-specific antibody responses up to six months following Omicron BA.1 breakthrough infection in mRNA-vaccinated individuals. Cross-reactive serum neutralizing antibody and memory B cell (MBC) responses declined by two- to four-fold through the study period. Breakthrough infection elicited minimal de novo Omicron-specific B cell responses but drove affinity maturation of pre-existing cross-reactive MBCs toward BA.1. Public clones dominated the neutralizing antibody response at both early and late time points, and their escape mutation profiles predicted newly emergent Omicron sublineages. The results demonstrate that heterologous SARS-CoV-2 variant exposure drives the evolution of B cell memory and suggest that convergent neutralizing antibody responses continue to shape viral evolution.

Mosaic RBD nanoparticles protect against challenge by diverse sarbecoviruses in animal models.

To combat future severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants and spillovers of SARS-like betacoronaviruses (sarbecoviruses) threatening global health, we designed mosaic nanoparticles that present randomly arranged sarbecovirus spike receptor-binding domains (RBDs) to elicit antibodies against epitopes that are conserved and relatively occluded rather than variable, immunodominant, and exposed. We compared immune responses elicited by mosaic-8 (SARS-CoV-2 and seven animal sarbecoviruses) and homotypic (only SARS-CoV-2) RBD nanoparticles in mice and macaques and observed stronger responses elicited by mosaic-8 to mismatched (not on nanoparticles) strains, including SARS-CoV and animal sarbecoviruses. Mosaic-8 immunization showed equivalent neutralization of SARS-CoV-2 variants, including Omicrons, and protected from SARS-CoV-2 and SARS-CoV challenges, whereas homotypic SARS-CoV-2 immunization protected only from SARS-CoV-2 challenge. Epitope mapping demonstrated increased targeting of conserved epitopes after mosaic-8 immunization. Together, these results suggest that mosaic-8 RBD nanoparticles could protect against SARS-CoV-2 variants and future sarbecovirus spillovers.

The SARS-CoV-2 Delta variant induces an antibody response largely focused on class 1 and 2 antibody epitopes.

Exposure histories to SARS-CoV-2 variants and vaccinations will shape the specificity of antibody responses. To understand the specificity of Delta-elicited antibody immunity, we characterize the polyclonal antibody response elicited by primary or mRNA vaccine-breakthrough Delta infections. Both types of infection elicit a neutralizing antibody response focused heavily on the receptor-binding domain (RBD). We use deep mutational scanning to show that mutations to the RBD's class 1 and class 2 epitopes, including sites 417, 478, and 484-486 often reduce binding of these Delta-elicited antibodies. The anti-Delta antibody response is more similar to that elicited by early 2020 viruses than the Beta variant, with mutations to the class 1 and 2, but not class 3 epitopes, having the largest effects on polyclonal antibody binding. In addition, mutations to the class 1 epitope (e.g., K417N) tend to have larger effects on antibody binding and neutralization in the Delta spike than in the D614G spike, both for vaccine- and Delta-infection-elicited antibodies. These results help elucidate how the antigenic impacts of SARS-CoV-2 mutations depend on exposure history.

Shifting mutational constraints in the SARS-CoV-2 receptor-binding domain during viral evolution.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has evolved variants with substitutions in the spike receptor-binding domain (RBD) that affect its affinity for angiotensin-converting enzyme 2 (ACE2) receptor and recognition by antibodies. These substitutions could also shape future evolution by modulating the effects of mutations at other sites-a phenomenon called epistasis. To investigate this possibility, we performed deep mutational scans to measure the effects on ACE2 binding of all single-amino acid mutations in the Wuhan-Hu-1, Alpha, Beta, Delta, and Eta variant RBDs. Some substitutions, most prominently Asn501→Tyr (N501Y), cause epistatic shifts in the effects of mutations at other sites. These epistatic shifts shape subsequent evolutionary change-for example, enabling many of the antibody-escape substitutions in the Omicron RBD. These epistatic shifts occur despite high conservation of the overall RBD structure. Our data shed light on RBD sequence-function relationships and facilitate interpretation of ongoing SARS-CoV-2 evolution.

An antibody-escape estimator for mutations to the SARS-CoV-2 receptor-binding domain.

A key goal of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) surveillance is to rapidly identify viral variants with mutations that reduce neutralization by polyclonal antibodies elicited by vaccination or infection. Unfortunately, direct experimental characterization of new viral variants lags their sequence-based identification. Here we help address this challenge by aggregating deep mutational scanning data into an 'escape estimator' that estimates the antigenic effects of arbitrary combinations of mutations to the virus's spike receptor-binding domain. The estimator can be used to intuitively visualize how mutations impact polyclonal antibody recognition and score the expected antigenic effect of combinations of mutations. These scores correlate with neutralization assays performed on SARS-CoV-2 variants and emphasize the ominous antigenic properties of the recently described Omicron variant. An interactive version of the estimator is at https://jbloomlab.github.io/SARS2_RBD_Ab_escape_maps/escape-calc/ (last accessed 11 March 2022), and we provide a Python module for batch processing. Currently the calculator uses primarily data for antibodies elicited by Wuhan-Hu-1-like vaccination or infection and so is expected to work best for calculating escape from such immunity for mutations relative to early SARS-CoV-2 strains.