Computational studies on the interaction of Omicron subvariants (BA.1, BA.2, and BA.3) with ACE2 and polyphenols

The SARS-CoV-2 Omicron variant BA.2 is spreading widely across the globe. The World Health Organization (WHO) designated BA.2 as a variant of concern due to its high transmission rate and pathogenicity. To elucidate the structural changes caused by mutations, we conducted a comparative analysis of BA.2 with variants BA.1 and BA.3.

In the present study, we aimed to investigate the interactions of the spike glycoprotein receptor-binding domain (SGp RBD) of Omicron variants BA.1, BA.2, and BA.3 with the human receptor hACE2. Further, a library of 233 polyphenols was screened by molecular docking with the SGp RBDs of Omicron variants BA.1, BA.2, and BA.3.

Protein-protein and protein-ligand molecular docking simulations were performed with AutoDock Vina and the ClusPro 2.0 server, respectively. The protein-ligand interactions were evaluated by BIOVIA Discovery Studio and ChimeraX 1.4. The molecular dynamics simulations for 100 ns were performed using GROMACS 2021.

Compared to other variants of concern, the structural changes in Omicron caused by mutations at key positions improved its ability to cause infection. Despite multiple mutations, many important polyphenols bind effectively at the RBDs of Omicron variants, with the required pharmacokinetic and ADME features and obeying the Lipinski rule.

Even though Omicron variants have multiple mutations and their transmission rate is relatively high, the computed binding affinities of lead polyphenols like epigallocatechin-3-O-gallate (EGCG) and luteolin-7-O-glucuronide (L7G) indicate that traditional medicines and proper immunity booster diets may be useful in the long-term fight against SARS-CoV-2.