High-throughput molecular dynamics-based alchemical free energy calculations for predicting the binding free energy change associated with the common mutations in the spike receptor-binding domain of SARS-CoV-2

The ongoing pandemic caused by SARS-CoV-2 has gone through various phases. From the initial outbreak the virus has mutated several times, with some lineages showing even stronger infectivity and faster spread than the original virus. Among all the variants, beta, gamma, delta and the latest (omicron) are currently classified as variants of concern (VOC) while the remaining are labelled either as variants of interest (VOI) or variants under monitoring (VUM). In this work, we have focused on the mutations observed in important variants, particularly at the receptor-binding domain (RBD) of the spike protein that is responsible for the interactions with the host ACE2 receptor and binding of antibodies. Studying these mutations is particularly important for understanding the viral infectivity, spread of the disease and for tracking the escape routes of this virus from antibodies. Molecular dynamics (MD) based alchemical free energy calculations have been shown to be very accurate in predicting the free energy change due to a mutation that could have a deleterious or a stabilising effect on the protein itself or its binding affinity to another protein. Here, we investigated the significance of six commonly observed spike RBD mutations on the stability of the spike protein binding to ACE2 by free energy calculations using high throughput MD simulations. For comparison, we also used other (rigorous and non-rigorous) binding free energy prediction methods and compared our results with the experimental data if available. The alchemical free energy-based method consistently predicted the free-energy changes with an accuracy close to ±1.0 kcal/mol when compared with the available experimental values. As per our simulation data the most significant mutations responsible for stabilising the spike RBD interactions with human ACE2 are N501Y and L452R.