The spike protein receptor-binding domain (RBD) of SARS-CoV-2 is the molecular target for many vaccines and antibody-based prophylactics aimed at bringing COVID-19 under control. Such a narrow molecular focus raises the specter of viral immune evasion as a potential failure mode for these biomedical interventions. With the emergence of new strains of SARS-CoV-2 with altered transmissibility and immune evasion potential, a critical question is this: how easily can the virus escape neutralizing antibodies (nAbs) targeting the spike RBD? To answer this question, we combined an analysis of the RBD structure-function with an evolutionary modeling framework. Our structure-function analysis revealed that epitopes for RBD-targeting nAbs overlap one another substantially and can be evaded by escape mutants with ACE2 affinities comparable to the wild type, that are observed in sequence surveillance data and infect cells in vitro. This suggests that the fitness cost of nAb-evading mutations is low. We then used evolutionary modeling to predict the frequency of immune escape before and after the widespread presence of nAbs due to vaccines, passive immunization or natural immunity. Our modeling suggests that SARS-CoV-2 mutants with one or two mildly deleterious mutations are expected to exist in high numbers due to neutral genetic variation, and consequently resistance to vaccines or other prophylactics that rely on one or two antibodies for protection can develop quickly -and repeatedly- under positive selection. Predicted resistance timelines are comparable to those of the decay kinetics of nAbs raised against vaccinal or natural antigens, raising a second potential mechanism for loss of immunity in the population. Strategies for viral elimination should therefore be diversified across molecular targets and therapeutic modalities.

The rapid emergence of immune-evading viral variants of SARS-CoV-2 calls into question the practicality of a vaccine-only public-health strategy for managing the ongoing COVID-19 pandemic. It has been suggested that widespread vaccination is necessary to prevent the emergence of future immune-evading mutants. Here, we examined that proposition using stochastic computational models of viral transmission and mutation. Specifically, we looked at the likelihood of emergence of immune escape variants requiring multiple mutations and the impact of vaccination on this process. Our results suggest that the transmission rate of intermediate SARS-CoV-2 mutants will impact the rate at which novel immune-evading variants appear. While vaccination can lower the rate at which new variants appear, other interventions that reduce transmission can also have the same effect. Crucially, relying solely on widespread and repeated vaccination (vaccinating the entire population multiple times a year) is not sufficient to prevent the emergence of novel immune-evading strains, if transmission rates remain high within the population. Thus, vaccines alone are incapable of slowing the pace of evolution of immune evasion, and vaccinal protection against severe and fatal outcomes for COVID-19 patients is therefore not assured.

As the COVID-19 pandemic drags into its second year, there is hope on the horizon, in the form of SARS-CoV-2 vaccines which promise disease suppression and a return to pre-pandemic normalcy. In this study we critically examine the basis for that hope, using an epidemiological modeling framework to establish the link between vaccine characteristics and effectiveness in bringing an end to this unprecedented public health crisis. Our findings suggest that a return to pre-pandemic social and economic conditions without fully suppressing SARS-CoV-2 will lead to extensive viral spread, resulting in a high disease burden even in the presence of vaccines that reduce risk of infection and mortality. Our modeling points to the feasibility of complete SARS-CoV-2 suppression with high population-level compliance and vaccines that are highly effective at reducing SARS-CoV-2 infection. Notably, vaccine-mediated reduction of transmission is critical for viral suppression, and in order for partially-effective vaccines to play a positive role in SARS-CoV-2 suppression, complementary biomedical interventions and public health measures must be deployed simultaneously.

As many prophylactics targeting SARS-CoV-2 are aimed at the spike protein receptor-binding domain (RBD), we examined the risk of immune evasion from previously published RBD-targeting neutralizing antibodies (nAbs). Epitopes for RBD-targeting nAbs overlap one another substantially and can give rise to escape mutants with ACE2 affinities comparable to wild type that still infect cells in vitro . We used evolutionary modeling to predict the frequency of immune escape before and after the widespread presence of nAbs due to vaccines, passive immunization or natural immunity. Our modeling suggests that SARS-CoV-2 mutants with one or two mildly deleterious mutations are expected to exist in high numbers due to neutral genetic variation, and consequently resistance to single or double antibody combinations can develop quickly under positive selection. One Sentence Summary: SARS-CoV-2 will evolve quickly to evade widely deployed spike RBD-targeting monoclonal antibodies, requiring combinations that rely on at least three antibodies targeting distinct epitopes to suppress viral immune evasion.

The development and deployment of several SARS-CoV-2 vaccines in a little over a year is an unprecedented achievement of modern medicine. The high levels of efficacy against transmission for some of these vaccines makes it feasible to use them to suppress SARS-CoV-2 altogether in regions with high vaccine acceptance. However, viral variants with reduced susceptibility to vaccinal and natural immunity threaten the utility of vaccines, particularly in scenarios where a return to pre-pandemic conditions occurs before the suppression of SARS-CoV-2 transmission. In this work we model the situation in the United States in May-June 2021, to demonstrate how pre-existing variants of SARS-CoV-2 may cause a rebound wave of COVID-19 in a matter of months under a certain set of conditions. A high burden of morbidity (and likely mortality) remains possible, even if the vaccines are partially effective against new variants and widely accepted. Our modeling suggests that variants that are already present within the population may be capable of quickly defeating the vaccines as a public health intervention, a serious potential limitation for strategies that emphasize rapid reopening before achieving control of SARS-CoV-2.

The duration of natural immunity in response to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a matter of some debate in the literature at present. For example, in a recent publication characterizing SARS-CoV-2 immunity over time, the authors fit pooled longitudinal data, using fitted slopes to infer the duration of SARS-CoV-2 immunity. In fact, such approaches can lead to misleading conclusions as a result of statistical model-fitting artifacts. To exemplify this phenomenon, we reanalyzed one of the markers (pseudovirus neutralizing titer) in the publication, using mixed-effects modeling, a methodology better suited to longitudinal datasets like these. Our findings showed that the half-life was both longer and more variable than reported by the authors. The example selected by us here illustrates the utility of mixed-effects modeling in provide more accurate estimates of the duration and heterogeneity of half-lives of molecular and cellular biomarkers of SARS-CoV-2 immunity.