Exploring the immune evasion of SARS-CoV-2 variant harboring E484K by molecular dynamics simulations

Although the current coronavirus disease 2019 (COVID-19) vaccines have been used worldwide to halt spread of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the emergence of new SARS-CoV-2 variants with E484K mutation shows significant resistance to the neutralization of vaccine sera. To better understand the resistant mechanism, we calculated the binding affinities of 26 antibodies to wild-type (WT) spike protein and to the protein harboring E484K mutation, respectively. The results showed that most antibodies (~85%) have weaker binding affinities to the E484K mutated spike protein than to the WT, indicating the high risk of immune evasion of the mutated virus from most of current antibodies. Binding free energy decomposition revealed that the residue E484 forms attraction with most antibodies, while the K484 has repulsion from most antibodies, which should be the main reason of the weaker binding affinities of E484K mutant to most antibodies. Impressively, a monoclonal antibody (mAb) combination was found to have much stronger binding affinity with E484K mutant than WT, which may work well against the mutated virus. Based on binding free energy decomposition, we predicted that the mutation of four more residues on receptor-binding domain (RBD) of spike protein, viz., F490, V483, G485 and S494, may have high risk of immune evasion, which we should pay close attention on during the development of new mAb therapeutics.