Abstract
The purpose of this work is to provide an in silico molecular rationale of the role eventually played by currently circulating mutations in the receptor binding domain of the SARS-CoV-2 spike protein (S-RBD(CoV‑2)) in evading the immune surveillance effects elicited by the two Eli Lilly LY-CoV555/bamlanivimab and LY-CoV016/etesevimab monoclonal antibodies. The main findings from this study show that, compared to the wild-type SARS-CoV-2 spike protein, mutations E484A/G/K/Q/R/V, Q493K/L/R, S494A/P/R, L452R and F490S are predicted to be markedly resistant to neutralization by LY-CoV555, while mutations K417E/N/T, D420A/G/N, N460I/K/S/T, T415P, and Y489C/S are predicted to confer LY-CoV016 escaping advantage to the viral protein. A challenge of our global in silico results against relevant experimental data resulted in an overall 90% agreement. Thus, the results presented provide a molecular-based rationale for all relative experimental findings, constitute a fast and reliable tool for identifying and prioritizing all present and newly reported circulating spike SARS-CoV-2 variants with respect to antibody neutralization, and yield substantial structural information for the development of next-generation vaccines and monoclonal antibodies more resilient to viral evolution.