Abstract
The pandemic of COVID-19, caused by SARS-CoV-2, has globally affected the human health and economy. Since the emergence of the novel coronavirus SARS-CoV-2, the life-threatening virus continues to mutate and evolve. Irrespective of acquired natural immunity and vaccine-induced immunity, the emerging multiple variants are growing exponentially, crossing the territorial barriers of the modern world. The rapid emergence of SARS-CoV-2 multiple variants challenges global researchers regarding the efficacy of available vaccines and variant transmissibility. SARS-CoV-2 surface-anchored S-protein recognizes and interacts with the host-cell ACE2, facilitating viral adherence and entrance into the cell. Understanding the interfacial interactions between the spike protein of SARS-CoV-2 variants and human ACE2 receptor is important for the design and development of antiviral therapeutics against SARS-CoV-2 emerging variants. Despite extensive research, the crucial determinants related to the molecular interactions between the spike protein of SARS-CoV-2 variants and host receptors are poorly understood. Thus, in this study, we explore the comparative interfacial binding pattern of SARS-CoV-2 spike RBD of wild type, Delta, and Omicron with the human ACE2 receptor to determine the crucial determinants at the atomistic level, using MD simulation and MM/GBSA energy calculations. Based on our findings, the substitution of Q493R, G496S, Q498R, and Y505H induced internal conformational changes in Omicron spike RBD, which leads to higher binding affinity than Delta spike RBD with the human ACE2 receptor, eventually contributing to higher transmission and infectivity. Taken together, these results could be used for the structure-based design of effective antiviral therapeutics against SARS-CoV-2 variants.