In Silico Genome Analysis Reveals the Evolution and Potential Impact of SARS-CoV-2 Omicron Structural Changes on Host Immune Evasion and Antiviral Therapeutics

Abstract

New variants of SARS-CoV-2 continue to evolve. The novel SARS-CoV-2 variant of concern (VOC) B.1.1.529 (Omicron) was particularly menacing due to the presence of numerous consequential mutations. In this study, we reviewed about 12 million SARS-CoV-2 genomic and associated metadata using extensive bioinformatic approaches to understand how evolutionary and mutational changes affect Omicron variant properties. Subsampled global data based analysis of molecular clock in the phylogenetic tree showed 29.56 substitutions per year as the evolutionary rate of five VOCs. We observed extensive mutational changes in the spike structural protein of the Omicron variant. A total of 20% of 7230 amino acid and structural changes exclusive to Omicron’s spike protein were detected in the receptor binding domain (RBD), suggesting differential selection pressures exerted during evolution. Analyzing key drug targets revealed mutation-derived differential binding affinities between Delta and Omicron variants. Nine single-RBD substitutions were detected within the binding site of approved therapeutic monoclonal antibodies. T-cell epitope prediction revealed eight immunologically important functional hotspots in three conserved non-structural proteins. A universal vaccine based on these regions may likely protect against all these SARS-CoV-2 variants. We observed key structural changes in the spike protein, which decreased binding affinities, indicating that these changes may help the virus escape host cellular immunity. These findings emphasize the need for continuous genomic surveillance of SARS-CoV-2 to better understand how novel mutations may impact viral spread and disease outcome.