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

Abstract

A key goal of 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 calculator" that estimates the antigenic effects of arbitrary combinations of mutations to the virus's spike receptor-binding domain (RBD). The calculator 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 calculator is at https://jbloomlab.github.io/SARS2_RBD_Ab_escape_maps/escape-calc/ , and we provide a Python module for batch processing.

Figure 1

Escape map for a hypothetical polyclonal mix consisting of an equipotent mixture of three monoclonal antibodies targeting distinct epitopes on the SARS-CoV-2 RBD. (A) Experimentally measured escape maps for three antibodies, and the mean of these maps (thick black line). Each point on the x-axis represents a site in the RBD, and the y-axis represents the total measured escape by all mutations at that site scaled so the maximum for each antibody is one. (B) Escape map if the contribution of antibody LY-CoV555 is ablated. (C) Escape map if the contributions of antibodies LY-CoV555 and LY-CoV016 are ablated. An interactive version of this figure is at https://jbloomlab.github.io/SARS2_RBD_Ab_escape_maps/mini-example-escape-calc/.

Figure 2

An escape calculator generated by aggregating deep mutational scanning for 33 neutralizing antibodies targeting the SARS-CoV-2 RBD. (A) The blue line shows the extent of escape mediated by mutations at each site, as estimated by simply averaging the data for all the individual antibodies. (B) The blue line shows escape map after a mutation to site 484 (red point) ablates recognition by antibodies strongly targeting that site, while the gray line shows the original escape map in the absence of any mutations. (C) The escape map after mutating sites 417, 484, and 501 (the three RBD sites mutated in the Beta variant). An interactive version of this figure is at https://jbloomlab.github.io/SARS2_RBD_Ab_escape_maps/escape-calc/.

Figure 3

Correlation of calculated binding scores with experimentally measured fold changes in neutralization for SARS-CoV-2 variants and mutants (smaller values indicate worse neutralization). The data of Lucas et al. (2021) was generated using authentic SARS-CoV-2 and sera from vaccinated individuals who were (top) or were not (bottom) previously infected with SARS-CoV-2. The data of Uriu et al. (2021) and Wang et al. (2021) was generated using pseudovirus against convalescent (top) or vaccine (bottom) sera, with vaccine sera from Pfizer BNT162b2 or Moderna mRNA-1273 vaccines, respectively. The fold changes are geometric means over all subjects in each cohort. An interactive version of this figure that allows mousing over points to see details is at https://jbloomlab.github.io/RBD_escape_calculator_paper/neut_studies.html.

Figure 4

Escape calculations for the Omicron variant. (A) The calculated binding scores for SARS-CoV-2 variants and the artificial polymutant spike (PMS20) generated by Schmidt et al. (2021). Scores of one indicate no mutations affect binding, and scores of zero indicate no antibody binding remains. An interactive version of this plot that allows mousing over points to see details is at https://jbloomlab.github.io/RBD_escape_calculator_paper/variants.html. (B) The calculated escape map for the Omicron variant’s RBD (blue) compared to an unmutated RBD (gray), with sites of mutations in the Omicron variant in red. The mutated RBD sites for each variant are in Table 1.