AbstractAt the time of this writing, August 2021, potential emergence of vaccine escape variants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a grave global concern. The interface between the receptor-binding domain (RBD) of SARS-CoV-2 spike (S) protein and the host receptor (ACE2) overlap with the binding site of principal neutralizing antibodies (NAb), limiting the repertoire of viable mutations. Nonetheless, variants with multiple mutations in the RBD have rose to dominance. Non-additive, epistatic relationships among RBD mutations are apparent, and assessing the impact of such epistasis on the mutational landscape is crucial. Epistasis can substantially increase the risk of vaccine escape and cannot be completely characterized through the study of the wild type (WT) alone. We employed protein structure modeling using Rosetta to compare the effects of all single mutants at the RBD-NAb and RBD-ACE2 interfaces for the WT, Gamma (417T, 484K, 501Y), and Delta variants (452R, 478K). Overall, epistasis at the RBD surface appears to be limited and the effects of most multiple mutations are additive. Epistasis at the Delta variant interface weakly stabilizes NAb interaction relative to ACE2, whereas in the Gamma variant, epistasis more substantially destabilizes NAb interaction. These results suggest that the repertoire of potential escape mutations for the Delta variant is not substantially different from that of the WT, whereas Gamma poses a moderately greater risk for enhanced vaccine escape. Thus, the modest ensemble of mutations relative to the WT shown to reduce vaccine efficacy might constitute the majority of all possible escape mutations.SignificancePotential emergence of vaccine escape variants of SARS-CoV-2 is arguably the most pressing problem during the COVID-19 pandemic as vaccines are distributed worldwide. We employed a computational approach to assess the risk of antibody escape resulting from mutations in the receptor-binding domain of the spike protein of the wild type SARS-CoV-2 virus as well as the Gamma and Delta variants. The results indicate that emergence of escape mutants is somewhat less likely for the Delta variant than for the wild type and moderately more likely for the Gamma variant. We conclude that the small set of escape-enhancing mutations already identified for the wild type is likely to include the majority of all possible mutations with this effect, a welcome finding.