Abstract
SARS-CoV-2 infection is initiated by binding of the receptor-binding domain (RBD) of its spike glycoprotein to the peptidase domain (PD) of angiotensin-converting enzyme 2 (ACE2) receptors in host cells. Recently detected Omicron variant of SARS-CoV-2 (B.1.1.529) is heavily mutated on RBD. First the BA.1 and later the BA.2 variant became the most dominant strains of the Omicron variant. To investigate how the mutations of these strains affect RBD-PD interactions, we performed all-atom molecular dynamics simulations of the BA.1 and BA.2 RBD-PD in the presence of full-length glycans, explicit water, and ions. Simulations revealed that RBDs of BA.1 and BA.2 variants exhibit a more dispersed interaction network and make an increased number of salt bridges and hydrophobic interactions with PD compared to wild-type RBD. Although BA.1 and BA.2 differ in two residues at the RBD-ACE2 interface, no major difference in RBD-PD interactions and binding strengths were observed between these variants. Using the conformations sampled in each trajectory, the Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) method estimated ∼34% and ∼51% stronger binding free energies to PD for BA.1 and BA.2 RBD, respectively, than wild-type RBD, which may result in higher binding efficiency of the Omicron variant to infect host cells.