key: cord-0744839-kidzs3kz
authors: Ortega, Joseph Thomas; Jastrzebska, Beata; Rangel, Hector Rafael
title: Omicron SARS-CoV-2 Variant Spike Protein Shows an Increased Affinity to the Human ACE2 Receptor: An In Silico Analysis
date: 2021-12-31
journal: Pathogens
DOI: 10.3390/pathogens11010045
sha: aba6e772b0fbb72a30c4f91b1ac9fff9c3258fda
doc_id: 744839
cord_uid: kidzs3kz

The rise of SARS-CoV-2 variants, with changes that could be related to an increased virus pathogenicity, have received the interest of the scientific and medical community. In this study, we evaluated the changes that occurred in the viral spike of the SARS-CoV-2 Omicron variant and whether these changes modulate the interactions with the angiotensin-converting enzyme 2 (ACE2) host receptor. The mutations associated with the Omicron variant were retrieved from the GISAID and covariants.org databases, and a structural model was built using the SWISS-Model server. The interaction between the spike and the human ACE2 was evaluated using two different docking software, Zdock and Haddock. We found that the binding free energy was lower for the Omicron variant as compared to the WT spike. In addition, the Omicron spike protein showed an increased number of electrostatic interactions with ACE2 than the WT spike, especially the interactions related to charged residues. This study contributes to a better understanding of the changes in the interaction between the Omicron spike and the human host ACE2 receptor.

The viral infection related to SARS-CoV-2 in humans is a recent evolutionary event. However, the occurrence of genetic mutations is increasing continuously while the virus is disseminated worldwide, allowing the identification of several viral variants that are actively circulating in the population [1] . These mutations, depending on the location and the type, can produce changes in the viral fitness and/or pathogenicity [2, 3] . The spike protein, which plays a pivotal role in the early events in viral replication, is a hot spot for such changes [4, 5] . The viral spike participates in the interaction with the host receptor angiotensin-converting enzyme 2 (ACE2), which is highly expressed in human airway cells [6] . The wild-type (WT) SARS-CoV-2 spike already showed changes that likely increased the efficiency of the virus-ACE2 receptor interactions as compared to SARS-CoV [7] . Major changes were noted in the region located within the residues 493 and 505 [7] . Thus, the SARS-CoV-2 viral variants could have a different transmissibility and response to neutralizing antibodies [8] [9] [10] . The mutations found in the different viral variants can occur in any part of the genome. However, some mutations that arose within the receptor-binding domain (RBD) received the most interest from researchers and medical groups due to its importance in host recognition and possibly changes in viral fitness [8, 11] . One of the first variants reported and widely distributed worldwide was the Alpha variant (B.1.1.7 linage) [12] . Later, the other variants, including Beta, Delta, and Mu, were found in the population (https://www.cdc.gov/coronavirus/2019-ncov/science/science-briefs/ scientific-brief-emerging-variants.html (accessed on 2 December 2021)). These variants contain several mutations in their genome that could modulate viral infectivity. In this study, we focused our interest on the mutations that occurred in the spike region interacting with ACE2, which could play a direct role in the modulation of the spike-ACE2 binding pattern. The Alpha variant has a non-synonymous substitution of N501Y, while the Beta variant, as well as Gamma and Mu variants additionally contain the E484K mutation (Table 1) . Previous reports showed that these mutations most likely increase the binding affinity of the spike to ACE2 as compared to the WT spike, without affecting its overall structure [13] [14] [15] . These mutations are located primarily in the protein loops, which due to their high flexibility could likely tolerate the conformational changes associated with the residue substitution. Interestingly, most of the single residue mutants exhibit a lower binding free energy, thus higher binding affinity to the host receptor in comparison to the WT variant [15] . The increased affinity of these viral variants to the host receptor could be explained by an increase in the number of hydrogen bonds and pi-pi interactions, as well as the change in charged residue interactions. The information gained from the analysis of the single residue mutants allowed us to hypothesize that variants carrying multiple mutations in the RBD will exhibit an even higher stabilizing effect on the interaction between the viral spike and the host ACE2 receptor [15] . The recently emerged Omicron variant could prove the correctness of our hypotheses. The Omicron variant carries over 30 mutations in the viral spike, among which eight mutations are located directly in the receptor-binding region of the RBD [16] . These mutations possibly could increase the binding stability of the spike to the ACE2. This variant shows a significant increase in the number of mutations in comparison to the Alpha variant that has one residue mutated in the RBD as well as the Beta and Gamma variants that contain three residues mutated in the RBD (Table 1 and Figure S1 ). 

To gain a deeper understanding on how these changes in the Omicron spike variant impact its binding affinity to the human host ACE2, we performed a molecular docking analysis using two different approaches, Zdock/Prodigy and Haddock (see Supplementary Materials for details). The obtained results indicated that the mutations that occurred in the RBD of the Omicron variant resulted in an increased number of interactions between the spike and ACE2 receptor, especially, the number of interactions associated with the charged residues was increased ( Table 2 ). Importantly, it has been previously described that the interactions between the charged residues of the viral spike and the host ACE2 play a major role in the complex stabilization and host selectivity [7, 18, 19] . These changes in the number and type of interactions were reflected in a decreased binding free energy (Table 2 ) and a decrease in the electrostatic, van der Waals interactions, and desolvation energies with an overall lower haddock score obtained for the Omicron variant as compared to the WT variant (Table 3) . The main interactions between the viral spike and the host ACE2 for the Omicron variant are shown in Figure 1 . The mutated residues in the RBD that directly interact with ACE2 are shown in orange. Interestingly, most of these residues are localized in the region between 493 and 505, which was previously described as important for the evolution between SARS-CoV and SARS-CoV-2 [7] . Furthermore, evaluation of the individual SARS-CoV-2 spike mutations in residue 505 revealed that the change from tyrosine to tryptophan resulted in the lower binding free energy interaction complex. In the Omicron variant, Y505 is substituted by the histidine residue, which resulted in fewer pi-pi interactions than observed for the double aromatic ring present in tryptophan [15] . However, histidine flexibility produces an increase in the number of interactions of the viral spike with other residues of ACE2, including E35 and D355, which could further stabilize the interaction between the spike and the host receptor (Table 4 ). In addition, other mutated residues, such as R493, Y501, and R498, showed changes in the type and number of interactions with ACE2 that are most likely related to conformational changes occurring in the spike protein, enhancing the overall stability of the spike-ACE2 complex (Table 4) .

Thus, this cluster located in the region between residues 498 to 505 could represent an evolutionary region that could be evaluated to predict the binding affinity of the emerging viral variants toward ACE2 in the human host.

To validate further the results obtained by the molecular docking, the stability of the SARS-CoV-2 Omicron spike-ACE2 complex was analyzed using MD simulations and compared to the WT spike-ACE2 complex. The results showed a similar trend of change in the residue fluctuation in the RBD region for both the WT and Omicron spike protein, confirming our results obtained by molecular docking (Figure 2 ). The median average root-mean-square fluctuation (RMSF) obtained with the CASB-flex software for the WT and Omicron variant was 1.05 ± 0.37 Å and 1.09 ± 0.25 Å, respectively. However, some differences in the fluctuations were noted for the residues 493, 501, and 505 due to the presence of the mutations. [20] . ProSA-web and PROCHECK programs were used to validate the quality of the structure [21, 22] . Hydrogen atoms were added, and partial charges were assigned for the energy refinement, as described in [23, 24] . The crystal structure of the SARS-CoV-2 spike protein bound to the human ACE2 receptor (PDB code: 6M0J) and the structure of the human ACE2 receptor (PDB code: 1R42) were downloaded from the Protein Data Bank. The Omicron variant and the WT spike proteins were evaluated. First, the docking between a ligand (RBD of the WT or Omicron spike protein) and a receptor (ACE2) was performed with Z-Dock software [25] , and the obtained complexes were processed and analyzed by using the tools in the PRODIGY software [26] . Furthermore, docking analysis for the Omicron spike protein was assayed using the Haddock server [27] .

To validate further the results obtained by the molecular docking, the stability of the SARS-CoV-2 Omicron spike-ACE2 complex was analyzed using MD simulations and compared to the WT spike-ACE2 complex. The results showed a similar trend of change in the residue fluctuation in the RBD region for both the WT and Omicron spike protein, The results obtained from molecular docking analysis for the interaction between the SARS-CoV-2 spike (blue) and ACE2 (pink) are shown. The residues mutated in the Omicron variant (R493, S496, R498, Y501, H505) that interact directly with the ACE2 are shown in orange. The sequence for the SARS-CoV-2 viral spike protein was retrieved from the Uniprot server (sequence number: P0DTC2) and the homology structural model for the Omicron variant was built by using the tools of the SWISS-MODEL modeling server and the DeepView/Swiss-PdbViewer 4.01 software [20] . ProSA-web and PROCHECK programs were used to validate the quality of the structure [21, 22] . Hydrogen atoms were added, and partial charges were assigned for the energy refinement, as described in [23, 24] . The crystal structure of the SARS-CoV-2 spike protein bound to the human ACE2 receptor (PDB code: 6M0J) and the structure of the human ACE2 receptor (PDB code: 1R42) were downloaded from the Protein Data Bank. The Omicron variant and the WT spike proteins were evaluated. First, the docking between a ligand (RBD of the WT or Omicron spike protein) and a receptor (ACE2) was performed with Z-Dock software [25] , and the obtained complexes were processed and analyzed by using the tools in the PRODIGY software [26] . Furthermore, docking analysis for the Omicron spike protein was assayed using the Haddock server [27] . confirming our results obtained by molecular docking (Figure 2 ). The median average root-mean-square fluctuation (RMSF) obtained with the CASB-flex software for the WT and Omicron variant was 1.05 ± 0.37 Å and 1.09 ± 0.25 Å, respectively. However, some differences in the fluctuations were noted for the residues 493, 501, and 505 due to the presence of the mutations. Changes in the pattern of flexibility in the interaction between the spike protein-ACE2 receptor complexes for Omicron and WT. The root-mean-square fluctuation (RMSF) of the protein Cα atoms with the respect to the initial structure for the evaluated spike protein-ACE2 complex was obtained by molecular docking. The spike protein-ACE2 receptor complexes were evaluated using the CABS-flex software [28] . The PDB files were submitted to the CABS-flex server to further assess the stability of the spike protein mutants-ACE2 complexes, and the parameters were adjusted to the default. The MD simulations output data obtained with both software were analyzed according to RMSF. This software enables an efficient modeling procedure for short simulations, being able to produce an analysis of the protein dynamics consistent with the dynamics obtained from 10nanoseconds MD simulations with the most popular force fields.

Herein, we discussed the importance of the viral spike for the entry into the cell and focused on the effects of the point mutations on the binding affinity to the host receptor ACE2. However, it is important to stress that the viral spike, despite targeting the host for virus replication, also plays a major role in the infection resolution, by being blocked with the antibodies either generated by the human host defense system upon infection or after vaccination, or by therapeutic monoclonal antibodies [29, 30] . Mutations occurring in the Figure 2 . Changes in the pattern of flexibility in the interaction between the spike protein-ACE2 receptor complexes for Omicron and WT. The root-mean-square fluctuation (RMSF) of the protein Cα atoms with the respect to the initial structure for the evaluated spike protein-ACE2 complex was obtained by molecular docking. The spike protein-ACE2 receptor complexes were evaluated using the CABS-flex software [28] . The PDB files were submitted to the CABS-flex server to further assess the stability of the spike protein mutants-ACE2 complexes, and the parameters were adjusted to the default. The MD simulations output data obtained with both software were analyzed according to RMSF. This software enables an efficient modeling procedure for short simulations, being able to produce an analysis of the protein dynamics consistent with the dynamics obtained from 10nanoseconds MD simulations with the most popular force fields.

Herein, we discussed the importance of the viral spike for the entry into the cell and focused on the effects of the point mutations on the binding affinity to the host receptor ACE2. However, it is important to stress that the viral spike, despite targeting the host for virus replication, also plays a major role in the infection resolution, by being blocked with the antibodies either generated by the human host defense system upon infection or after vaccination, or by therapeutic monoclonal antibodies [29, 30] . Mutations occurring in the viral genome could result in the evolution of so-called 'escape variants'; some of these variants can show reduced susceptibility to antibodies [31, 32] . Though, the previous analysis showed that changes in the spike binding affinity to the human host receptor did not correlate with changes in the response to monoclonal antibodies [15, 33, 34] . Moreover, some reports elucidated that the host antibodies generated after vaccination with the currently available mRNA vaccine could neutralize the prior circulating variants such as Alpha and Delta [35] [36] [37] [38] . However, a change in neutralization efficiency was observed for variants carrying the E484K mutation [33] . This is related to the antibody recognition site within the S1 region of the spike protein located between the residues 301 and 430 in the main receptor-binding motif [33, 39, 40] . The Omicron variant carries several mutations situated in this region. Thus, further analyses to assess the efficiency of the immune response elicited by vaccines to this variant are necessary. Importantly, the vaccination efforts need to be boosted up, ensuring broad coverage across the world to prevent the selection of new viral variants. In the current pandemic scenario, the world has observed that new variants can change the SARS-CoV-2 virus distribution dynamics in terms of prevalence in the population. The Delta variant has been the most prevalent before the rise of Omicron. This variant carries several mutations that could affect viral fitness. However, within the RBD of the spike protein there are only two mutated residues, L452R and T478K. These residues do not interact directly with the ACE2 receptor. However, the residue substitution produced a conformational change in the viral spike, enhancing the binding affinity of the virus to the host receptor. On the other hand, the Omicron variant carries at least eight mutations in key residues that interact directly with the host ACE2. The substitutions in the directly interacting residues resulted in stabilization of the viral spike-ACE2 interaction, evidenced by a decrease in the binding free energy found for Omicron as compared to WT by our molecular docking analysis. Thus, sequence analysis evaluating the mutations present in the viral variants within the region interacting with the host receptor ACE2 could enhance our understanding of the changes in the prevalence dynamics in the SARS-CoV-2 infection. Together, the results presented in this study indicate that mutations that emerged in the new SARS-CoV-2 Omicron variant enhance the binding affinity of this variant to the host ACE2 receptor, which likely is an important driver of its increased transmissibility. However, it is important to note that SARS-CoV-2 besides ACE2 could also use other cellular receptors as a secondary binding way to target the host cells [19, [41] [42] [43] . Thus, it cannot be excluded that the changes that emerged in the Omicron viral spike could have an impact on the virus interaction with these secondary cellular receptors, additionally increasing the virus's ability to reach its target. Nevertheless, further analyses employing the structural, computational, and biochemical approaches are required to gain a deeper understanding of the structural details of the residue rearrangements and changes in the specific virus-host interactions triggering higher transmissibility of the SARS-CoV-2 Omicron variant. Funding: This research was support in part by "Desarrollo de capacidades para el estudio y control del coronavirus causante de la COVID-19" Codecyt Caracas-Venezuela (To HR).

Informed Consent Statement: Not applicable.

Data Availability Statement: Data on SARS-CoV-2 variants was retrieved from the GISAID [17] and covariants.org webservers, on 26 November 2021. The sequence for the SARS-CoV-2 viral spike protein was retrieved from the Uniprot server (sequence number: P0DTC2). The crystal structure of the SARS-CoV-2 spike protein bound to the human ACE2 receptor (PDB code: 6M0J) and the structure of the human ACE2 receptor (PDB code: 1R42) were downloaded from the Protein Data Bank.

COVID-19 and SARS-CoV-2 Variants: Current Challenges and Health Concern

SARS-CoV-2 Quasispecies Provides an Advantage Mutation Pool for the Epidemic Variants

The biological and clinical significance of emerging SARS-CoV-2 variants

Recombination events are concentrated in the spike protein region of Betacoronaviruses

Evolutionary and structural analysis elucidates mutations on SARS-CoV2 spike protein with altered human ACE2 binding affinity

Structure of SARS-CoV-2 spike protein

Role of changes in SARS-CoV-2 spike protein in the interaction with the human ACE2 receptor: An in silico analysis

Mutations of SARS-CoV-2 spike protein: Implications on immune evasion and vaccine-induced immunity

Emerging SARS-CoV-2 variants can potentially break set epidemiological barriers in COVID-19

Replication Kinetics of B.1.351 and B.1.1.7 SARS-CoV-2 Variants of Concern Including Assessment of a B.1.1.7 Mutant Carrying a Defective ORF7a Gene

Implication of SARS-CoV-2 Immune Escape Spike Variants on Secondary and Vaccine Breakthrough Infections

Emergence of SARS-CoV-2 B.1.1.7 Lineage-United States

In silico binding profile characterization of SARS-CoV-2 spike protein and its mutants bound to human ACE2 receptor

Bioinformatics analysis of the differences in the binding profile of the wild-type and mutants of the SARS-CoV-2 spike protein variants with the ACE2 receptor

Mutations in the SARS-CoV-2 spike protein modulate the virus affinity to the human ACE2 receptor, an in silico analysis

CoVariance: SARS-CoV-2 Mutations and Variants of Interest

Data, disease and diplomacy: GISAID's innovative contribution to global health

Structure, Function, and Antigenicity of the SARS-CoV-2 Spike Glycoprotein

SARS-CoV-2 host tropism: An in silico analysis of the main cellular factors

The SWISS-MODEL workspace: A web-based environment for protein structure homology modelling

ProSA-web: Interactive web service for the recognition of errors in three-dimensional structures of proteins

Main-Chain Bond Lengths and Bond Angles in Protein Structures

Unrevealing sequence and structural features of novel coronavirus using in silico approaches: The main protease as molecular target

Class A G Protein-Coupled Receptor Antagonist Famotidine as a Therapeutic Alternative against SARS-CoV2: An In Silico Analysis

Interactive docking prediction of proteinprotein complexes and symmetric multimers

PRODIGY: A web server for predicting the binding affinity of protein-protein complexes

Bonvin, A. The HADDOCK2.2 Web Server: User-Friendly Integrative Modeling of Biomolecular Complexes

CABS-flex 2.0: A web server for fast simulations of flexibility of protein structures

SARS-CoV-2 Virus-Host Interaction: Currently Available Structures and Implications of Variant Emergence on Infectivity and Immune Response

Potent SARS-CoV-2 neutralizing antibodies with protective efficacy against newly emerged mutational variants

Antibody resistance of SARS-CoV-2 variants B.1.351 and B.1.1.7

Molecular insights into receptor binding energetics and neutralization of SARS-CoV-2 variants

SARS-CoV-2 variants, spike mutations and immune escape

Computational prediction of the effect of amino acid changes on the binding affinity between SARS-CoV-2 spike RBD and human ACE2

SARS-COV-2 Delta variant displays moderate resistance to neutralizing antibodies and spike protein properties of higher soluble ACE2 sensitivity, enhanced cleavage and fusogenic activity

Immune Evasion of SARS-CoV-2 Emerging Variants: What Have We Learnt So Far? Viruses

Reduced sensitivity of SARS-CoV-2 variant Delta to antibody neutralization

Neutralization of SARS-CoV-2 spike 69/70 deletion, E484K and N501Y variants by BNT162b2 vaccine-elicited sera

T cell and antibody responses induced by a single dose of ChAdOx1 nCoV-19 (AZD1222) vaccine in a phase 1/2 clinical trial

Antibody signature induced by SARS-CoV-2 spike protein immunogens in rabbits

Receptome profiling identifies KREMEN1 and ASGR1 as alternative functional receptors of SARS-CoV-2

On the Selectivity of Heparan Sulfate Recognition by SARS-CoV-2 Spike Glycoprotein

Structural-Dynamics and Binding Analysis of RBD from SARS-CoV-2 Variants of Concern (VOCs) and GRP78 Receptor Revealed Basis for Higher Infectivity

The authors thanks Flor Pujol for her comments on the manuscript.

The authors declare no conflict of interest.