A computational evaluation of structural stability of omicron and delta mutations of SARS-CoV-2 spike proteins and human ACE-2 interactions

Several more infectious SARS-CoV-2 variants have emerged globally since SARS-CoV-2 pandemic and the discovery of the first D614G variant of SARS-CoV-2 spike proteins in 2020. Delta (B.1.617.2) and Omicron (B.1.1.529) variants have proven to be of major concern out of all the reported variants, considering their influence on the virus' transmissibility and severity. This study aimed at evaluating the impact of mutations on these two variants on stability and molecular interactions between the viral Spike protein and human angiotensin converting enzyme-2 (hACE-2). The spike proteins receptor binding domain (RBD) was docked with the hACE-2 using HADDOCK servers. To understand and establish the effects of the mutations on the structural stability and flexibility of the RBD-hACE-2 complex, molecular dynamic (MD) simulation of the docked complex was performed and evaluated. The findings from both molecular docking analysis and binding free energy showed that the Omicron (OM) variant has high receptiveness towards hACE-2 versus Delta variant (DT), thereby, responsible for its increase in transmission. The structural stability and flexibility evaluation of variants’ systems showed that mutations on DT and OM variants disturbed the stability of either the spike protein or the RBD-hACE-2 complex, with DT variant having greater instability impact. This study, therefore, assumed this obvious instability observed in DT variant might be associated or responsible for the reported severity in DT variant disease over the OM variant disease. This study provides molecular insight into the effects of OM and DT variants on stability and interactions between SARS-CoV-2 protein and hACE-2.