Structural and Energetic Insights into SARS-CoV-2 Evolution: Analysis of hACE2–RBD Binding in Wild-Type, Delta, and Omicron Subvariants

The evolution of SARS-CoV-2, particularly the emergence of Omicron variants, has raised questions regarding changes in its binding affinity to the human angiotensin-converting enzyme 2 receptor (hACE2). Understanding the impact of mutations on the interaction between the receptor-binding domain (RBD) of the spike protein and hACE2 is critical for evaluating viral transmissibility, immune evasion, and the efficacy of therapeutic strategies. Here, we used molecular dynamics (MD) simulations and binding energy calculations to investigate the structural and energetic differences between the hACE2- RBD complexes of wild-type (WT), Delta, and Omicron subvariants. Our results indicate that the Delta and the first Omicron variants showed the highest and the second-highest binding energy among the variants studied. Furthermore, while Omicron variants exhibit increased structural stability and altered electrostatic potential at the hACE2-RBD interface when compared to the ancestral WT, their binding strength to hACE2 does not consistently increase with viral evolution. Moreover, newer Omicron subvariants like JN.1 exhibit a bimodal conformational strategy, alternating between a high-affinity state for hACE2 and a low-affinity state, which could potentially facilitate immune evasion. These findings suggest that, in addition to enhanced hACE2 binding affinity, other factors, such as immune evasion and structural adaptability, shape SARS-CoV-2 evolution.