Mutations of Omicron Variant at the Interface of the Receptor Domain Motif and Human Angiotensin-Converting Enzyme-2

Increased preference for lysine over arginine in spike proteins of SARS-CoV-2 BA.2.86 variant and its daughter lineages

The COVID-19 pandemic offered an unprecedented glimpse into the evolution of its causative virus, SARS-CoV-2. It has been estimated that since its outbreak in late 2019, the virus has explored all possible alternatives in terms of missense mutations for all sites of its polypeptide chain. Spike protein of the virus exhibits the largest sequence variation in particular, with many individual mutations impacting target recognition, cellular entry, and endosomal escape of the virus. Moreover, recent studies unveiled a significant increase in the total charge on the spike protein during the evolution of the virus in the initial period of the pandemic. While this trend has recently come to a halt, we perform a sequence-based analysis of the spike protein of 2665 SARS-CoV-2 variants which shows that mutations in ionizable amino acids continue to occur with the newly emerging variants, with notable differences between lineages from different clades. What is more, we show that within mutations of amino acids which can acquire positive charge, the spike protein of SARS-CoV-2 exhibits a prominent preference for lysine residues over arginine residues. This lysine-to-arginine ratio increased at several points during spike protein evolution, most recently with BA.2.86 and its sublineages, including the recently dominant JN.1, KP.3, and XEC variants. The increased ratio is a consequence of mutations in different structural regions of the spike protein and is now among the highest among viral species in the Coronaviridae family. The impact of high lysine-to-arginine ratio in the spike proteins of BA.2.86 and its daughter lineages on viral fitness remains unclear; we discuss several potential mechanisms that could play a role and that can serve as a starting point for further studies.

Semi-Covariance Coefficient Analysis of Spike Proteins from SARS-CoV-2 and Its Variants Omicron, BA.5, EG.5, and JN.1 for Viral Infectivity, Virulence and Immune Escape

Semi-covariance has attracted significant attention in recent years and is increasingly employed to elucidate statistical phenomena exhibiting fluctuations, such as the similarity or difference in charge patterns of spike proteins among coronaviruses. In this study, by examining values above and below the average/mean based on the positive and negative charge patterns of amino acid residues in the spike proteins of SARS-CoV-2 and its current circulating variants, the proposed methods offer profound insights into the nonlinear evolving trends in those viral spike proteins. Our study indicates that the charge span value can predict the infectivity of the virus and the charge density can estimate the virulence of the virus, and both predicated infectivity and virulence appear to be associated with the capability of viral immune escape. This semi-covariance coefficient analysis may be used not only to predict the infectivity, virulence and capability of immune escape for coronaviruses but also to analyze the functionality of other viral proteins. This study improves our understanding of the trend of viral evolution in terms of viral infectivity, virulence or the capability of immune escape, which remains further validated by more future studies and statistical data.

Changes in total charge on spike protein of SARS-CoV-2 in emerging lineages

Motivation

Charged amino acid residues on the spike protein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) have been shown to influence its binding to different cell surface receptors, its non-specific electrostatic interactions with the environment, and its structural stability and conformation. It is therefore important to obtain a good understanding of amino acid mutations that affect the total charge on the spike protein which have arisen across different SARS-CoV-2 lineages during the course of the virus' evolution.

Results

We analyse the change in the number of ionizable amino acids and the corresponding total charge on the spike proteins of almost 2200 SARS-CoV-2 lineages that have emerged over the span of the pandemic. Our results show that the previously observed trend toward an increase in the positive charge on the spike protein of SARS-CoV-2 variants of concern has essentially stopped with the emergence of the early omicron variants. Furthermore, recently emerged lineages show a greater diversity in terms of their composition of ionizable amino acids. We also demonstrate that the patterns of change in the number of ionizable amino acids on the spike protein are characteristic of related lineages within the broader clade division of the SARS-CoV-2 phylogenetic tree. Due to the ubiquity of electrostatic interactions in the biological environment, our findings are relevant for a broad range of studies dealing with the structural stability of SARS-CoV-2 and its interactions with the environment.

Availability and implementation

The data underlying the article are available in the Supplementary material.

Predicting potential SARS-CoV-2 mutations of concern via full quantum mechanical modelling

Ab initio quantum mechanical models can characterize and predict intermolecular binding, but only recently have models including more than a few hundred atoms gained traction. Here, we simulate the electronic structure for approximately 13 000 atoms to predict and characterize binding of SARS-CoV-2 spike variants to the human ACE2 (hACE2) receptor using the quantum mechanics complexity reduction (QM-CR) approach. We compare four spike variants in our analysis: Wuhan, Omicron, and two Omicron-based variants. To assess binding, we mechanistically characterize the energetic contribution of each amino acid involved, and predict the effect of select single amino acid mutations. We validate our computational predictions experimentally by comparing the efficacy of spike variants binding to cells expressing hACE2. At the time we performed our simulations (December 2021), the mutation A484K which our model predicted to be highly beneficial to ACE2 binding had not been identified in epidemiological surveys; only recently (August 2023) has it appeared in variant BA.2.86. We argue that our computational model, QM-CR, can identify mutations critical for intermolecular interactions and inform the engineering of high-specificity interactors.

Impact of BA.1, BA.2, and BA.4/BA.5 Omicron mutations on therapeutic monoclonal antibodies

The emergence of Omicron SARS-CoV-2 subvariants (BA.1, BA.2, BA.4, and BA.5), with an unprecedented number of mutations in their receptor-binding domain (RBD) of the spike-protein, has fueled a resurgence of COVID-19 infections, posing a major challenge to the efficacy of existing vaccines and monoclonal antibody (mAb) therapeutics. We conducted a systematic molecular dynamics (MD) simulation to investigate how the RBD mutations of these subvariants affect the interactions with broad mAbs including AstraZeneca (COV2-2196 and COV2-2130), Brii Biosciences (BRII-196), Celltrion (CT-P59), Eli Lilly (LY-CoV555 and LY-CoV016), Regeneron (REGN10933 and REGN10987), Vir Biotechnology (S309), and S2X259. Our results show a complete loss of binding for COV2-2196, BRII-196, CT-P59, and LY-CoV555 with all Omicron RBDs. Additionally, REGN10987 totally loses its binding with BA.1, but retains a partial binding with BA.2 and BA.4/5. The binding reduction is significant for LY-CoV016 and REGN10933 but moderate for COV2-2130. S309 and S2X259 retain their binding with BA.1 but exhibit decreased binding with other subvariants. We introduce a mutational escape map for each mAb to identify the key RBD sites and the corresponding critical mutations. Overall, our findings suggest that the majority of therapeutic mAbs have diminished or missing activity against Omicron subvariants, indicating the urgent need for a new therapeutic mAb with a better design.

The anti-SARS-CoV-2 effect and mechanism of Chiehyuan herbal oral protection solution

The Chiehyuan herbal oral protection solution (GB-2) is a herbal mixture commonly utilized in Taiwan for combating severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) as per traditional Chinese medicine practices. This study assessed the clinical impact of GB-2 through prospective clinical trials. With twice-daily use for a week, GB-2 was shown to diminish the expression of angiotensin-converting enzyme 2 (ACE2) in oral mucosal cells. Moreover, after two weeks of use, it could reduce transmembrane protease, serine 2 (TMRPSS2) expression in these cells. Additionally, in vitro experiments demonstrated that GB-2 lessened the entry efficiency of the Omicron, L452R-D614G, T478K-D614G, and L452R-T478K-D614G variants of the SARS-CoV-2 pseudotyped lentivirus. It also impeded the interaction between ACE2 and the receptor-binding domain (RBD) presenting N501Y-K417N-E484A-G339D-Q493R-G496S-Q498R and L452R-T478K mutations. Glycyrrhizic acid, a major compound in GB-2, also hindered the entry of the Omicron variant (BA.1) of the SARS-CoV-2 pseudotyped lentivirus by obstructing the binding between ACE2 and the RBD presenting the N501Y-K417N-E484A-G339D-Q493R-G496S-Q498R mutation. To sum up, these findings suggest that GB-2 can decrease ACE2 and TMPRSS2 expression in oral mucosal cells. Both glycyrrhizic acid and GB-2 were found to reduce the entry efficiency of the Omicron variant (BA.1) of the SARS-CoV-2 pseudotyped lentivirus and block the binding between ACE2 and the RBD with the N501Y-K417N-E484A-G339D-Q493R-G496S-Q498R mutation. This evidence implies that GB-2 might be a potential candidate for further study as a preventative measure against SARS-CoV-2 infection.

Evolutionary changes in the number of dissociable amino acids on spike proteins and nucleoproteins of SARS-CoV-2 variants

The spike protein of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is responsible for target recognition, cellular entry, and endosomal escape of the virus. At the same time, it is the part of the virus that exhibits the greatest sequence variation across the many variants which have emerged during its evolution. Recent studies have indicated that with progressive lineage emergence, the positive charge on the spike protein has been increasing, with certain positively charged amino acids (AAs) improving the binding of the spike protein to cell receptors. We have performed a detailed analysis of dissociable AAs of more than 1400 different SARS-CoV-2 lineages, which confirms these observations while suggesting that this progression has reached a plateau with Omicron and its subvariants and that the positive charge is not increasing further. Analysis of the nucleocapsid protein shows no similar increase in positive charge with novel variants, which further indicates that positive charge of the spike protein is being evolutionarily selected for. Furthermore, comparison with the spike proteins of known coronaviruses shows that already the wild-type SARS-CoV-2 spike protein carries an unusually large amount of positively charged AAs when compared to most other betacoronaviruses. Our study sheds light on the evolutionary changes in the number of dissociable AAs on the spike protein of SARS-CoV-2, complementing existing studies and providing a stepping stone towards a better understanding of the relationship between the spike protein charge and viral infectivity and transmissibility.

Elucidation of novel compounds and epitope-based peptide vaccine design against C30 endopeptidase regions of SARS-CoV-2 using immunoinformatics approaches

There has been progressive improvement in immunoinformatics approaches for epitope-based peptide design. Computational-based immune-informatics approaches were applied to identify the epitopes of SARS-CoV-2 to develop vaccines. The accessibility of the SARS-CoV-2 protein surface was analyzed, and hexa-peptide sequences (KTPKYK) were observed having a maximum score of 8.254, located between amino acids 97 and 102, whereas the FSVLAC at amino acids 112 to 117 showed the lowest score of 0.114. The surface flexibility of the target protein ranged from 0.864 to 1.099 having amino acid ranges of 159 to 165 and 118 to 124, respectively, harboring the FCYMHHM and YNGSPSG hepta-peptide sequences. The surface flexibility was predicted, and a 0.864 score was observed from amino acids 159 to 165 with the hepta-peptide (FCYMHHM) sequence. Moreover, the highest score of 1.099 was observed between amino acids 118 and 124 against YNGSPSG. B-cell epitopes and cytotoxic T-lymphocyte (CTL) epitopes were also identified against SARS-CoV-2. In molecular docking analyses, -0.54 to -26.21 kcal/mol global energy was observed against the selected CTL epitopes, exhibiting binding solid energies of -3.33 to -26.36 kcal/mol. Based on optimization, eight epitopes (SEDMLNPNY, GSVGFNIDY, LLEDEFTPF, DYDCVSFCY, GTDLEGNFY, QTFSVLACY, TVNVLAWLY, and TANPKTPKY) showed reliable findings. The study calculated the associated HLA alleles with MHC-I and MHC-II and found that MHC-I epitopes had higher population coverage (0.9019% and 0.5639%) than MHC-II epitopes, which ranged from 58.49% to 34.71% in Italy and China, respectively. The CTL epitopes were docked with antigenic sites and analyzed with MHC-I HLA protein. In addition, virtual screening was conducted using the ZINC database library, which contained 3,447 compounds. The 10 top-ranked scrutinized molecules (ZINC222731806, ZINC077293241, ZINC014880001, ZINC003830427, ZINC030731133, ZINC003932831, ZINC003816514, ZINC004245650, ZINC000057255, and ZINC011592639) exhibited the least binding energy (-8.8 to -7.5 kcal/mol). The molecular dynamics (MD) and immune simulation data suggest that these epitopes could be used to design an effective SARS-CoV-2 vaccine in the form of a peptide-based vaccine. Our identified CTL epitopes have the potential to inhibit SARS-CoV-2 replication.

Electrostatic Interactions Are the Primary Determinant of the Binding Affinity of SARS-CoV-2 Spike RBD to ACE2: A Computational Case Study of Omicron Variants

To explore the mechanistic origin that determines the binding affinity of SARS-CoV-2 spike receptor binding domain (RBD) to human angiotensin converting enzyme 2 (ACE2), we constructed the homology models of RBD-ACE2 complexes of four Omicron subvariants (BA.1, BA.2, BA.3 and BA.4/5), and compared them with wild type complex (RBD_WT-ACE2) in terms of various structural dynamic properties by molecular dynamics (MD) simulations and binding free energy (BFE) calculations. The results of MD simulations suggest that the RBDs of all the Omicron subvariants (RBD_OMIs) feature increased global structural fluctuations when compared with RBD_WT. Detailed comparison of BFE components reveals that the enhanced electrostatic attractive interactions are the main determinant of the higher ACE2-binding affinity of RBD_OMIs than RBD_WT, while the weakened electrostatic attractive interactions determine RBD of BA.4/5 subvariant (RBD_BA.4/5) lowest ACE2-binding affinity among all Omicron subvariants. The per-residue BFE decompositions and the hydrogen bond (HB) networks analyses indicate that the enhanced electrostatic attractive interactions are mainly through gain/loss of the positively/negatively charged residues, and the formation or destruction of the interfacial HBs and salt bridges can also largely affect the ACE2-binding affinity of RBD. It is worth pointing out that since Q493R plays the most important positive contribution in enhancing binding affinity, the absence of this mutation in RBD_BA.4/5 results in a significantly weaker binding affinity to ACE2 than other Omicron subvariants. Our results provide insight into the role of electrostatic interactions in determining of the binding affinity of SARS-CoV-2 RBD to human ACE2.

Molecular Mechanisms of Pathogenesis, Prevention, and Therapy of COVID-19: Summarizing the Results of 2021

The purpose of this special issue is to highlight the main problems of the COVID-19 epidemic and to outline some ways to solve these problems, including research into the biology of the SARS-CoV-2 virus, general pathological and particular patterns of COVID-19 pathogenesis, acute and long-term complications of COVID-19, and evaluation of high-potential general and specific prevention methods and etiological and pathogenetic therapies for COVID-19 [...].

Quantum Chemical Computation of Omicron Mutations Near Cleavage Sites of the Spike Protein

The attachment of the spike protein in SARS-CoV-2 to host cells and the initiation of viral invasion are two critical processes in the viral infection and transmission in which the presence of unique furin (S1/S2) and TMPRSS2 (S2′) cleavage sites play a pivotal role. We provide a detailed analysis of the impact of the BA.1 Omicron mutations vicinal to these cleavage sites using a novel computational method based on the amino acid–amino acid bond pair unit (AABPU), a specific protein structural unit as a proxy for quantifying the atomic interaction. Our study is focused mainly on the spike region between subdomain 2 (SD2) and the central helix (CH), which contains both S1/S2 and S2’ cleavage sites. Based on ab initio quantum calculations, we have identified several key features related to the electronic structure and bonding of the Omicron mutations that significantly increase the size of the relevant AABPUs and the positive charge. These findings enable us to conjecture on the biological role of Omicron mutations and their specific effects on cleavage sites and identify the principles that can be of some value in analyzing new variants.

Effect of Delta and Omicron Mutations on the RBD-SD1 Domain of the Spike Protein in SARS-CoV-2 and the Omicron Mutations on RBD-ACE2 Interface Complex

The receptor-binding domain (RBD) is the essential part in the Spike-protein (S-protein) of SARS-CoV-2 virus that directly binds to the human ACE2 receptor, making it a key target for many vaccines and therapies. Therefore, any mutations at this domain could affect the efficacy of these treatments as well as the viral-cell entry mechanism. We introduce ab initio DFT-based computational study that mainly focuses on two parts: (1) Mutations effects of both Delta and Omicron variants in the RBD-SD1 domain. (2) Impact of Omicron RBD mutations on the structure and properties of the RBD-ACE2 interface system. The in-depth analysis is based on the novel concept of amino acid-amino acid bond pair units (AABPU) that reveal the differences between the Delta and/or Omicron mutations and its corresponding wild-type strain in terms of the role played by non-local amino acid interactions, their 3D shapes and sizes, as well as contribution to hydrogen bonding and partial charge distributions. Our results also show that the interaction of Omicron RBD with ACE2 significantly increased its bonding between amino acids at the interface providing information on the implications of penetration of S-protein into ACE2, and thus offering a possible explanation for its high infectivity. Our findings enable us to present, in more conspicuous atomic level detail, the effect of specific mutations that may help in predicting and/or mitigating the next variant of concern.

Structural aspects of SARS-CoV-2 mutations: Implications to plausible infectivity with ACE-2 using computational modeling approach

Some of the SARS-CoV-2 variants are said to be more infectious than the previous others and are causing panic around the globe. Cases related to Delta plus (δ+) and omicron (ο) variants are on the rise worldwide. This sudden surge warrants an investigation into the reasons for its binding with ACE-2. The present study attempts to find out the structural basis of binding interactions of SARS-CoV-2 mutants based on computational modeling and comparative analysis. In silico strategies including protein-protein docking, mutation analysis, molecular dynamics, and binding energy calculations were used to study the binding of the 'receptor binding domain' (RBD) of the seven 'variants of concern' which include Alpha (α), Beta (β), Gamma (γ), Kappa (κ), Delta (δ), Delta plus (δ+) and omicron (ο) with ACE-2 (human angiotensin-converting enzyme-2) and with antibodies. Among all the variants dealt with in this study, Delta plus and omicron were found to be binding more strongly to ACE-2 than others due to inherent mutations and the consequent change in the hydrophilic and hydrophobic environment of the binding site. Furthermore, molecular dynamic (MD) simulations and subsequent MM/PBSA calculations provided useful structural insights into key residues participating in the interaction. Infectivity of a virus could be dependent on the interplay of evading antibodies and simultaneously attaching strongly with the host receptor. A cross-correlation between mutant spike proteins' binding with ACE-2 and antibodies provides a holistic assessment of the binding nature of these mutants vis-à-vis native virus and offers opportunities for designing potential therapeutics against these new mutants.

MicroRNAs as Potential Tools for Predicting Cancer Patients’ Susceptibility to SARS-CoV-2 Infection and Vaccination Response

Coronavirus disease (COVID-19) is an infectious disease that is caused by a highly contagious and severe acute respiratory syndrome—coronavirus 2 (SARS-CoV-2). This infection started to spread across the world in 2019 and rapidly turned into a global pandemic, causing an urgent necessity for treatment strategies development. The mRNA vaccines against SARS-CoV-2 can trigger an immune response, providing genetic information that allows the production of spike glycoproteins. MiRNAs play a crucial role in diverse key cellular processes, including antiviral defense. Several miRNAs are described as key factors in SARS-CoV-2 human infection through the regulation of ACE2 levels and by the inhibition of SARS-CoV-2 replication and spike expression. Consequently, these molecules have been considered as highly promising biomarkers. In numerous human malignancies, it has been recognized that miRNAs expression is dysregulated. Since miRNAs can target SARS-CoV-2-associated mRNAs, in cancer patients, the deregulation of these molecules can impair the immune response to the vaccines. Therefore, in this review, we propose a miRNA profile of seven SARS-CoV-2-related miRNAs, namely miR-214, miR-98-5p, miR-7-5p, miR-24-3p, miR-145-5p, miR-223-3p and miR-15b-5p, that are deregulated in a high number of cancers and have the potential to be used as prognostic biomarkers to stratify cancer patients.

An Insight Based on Computational Analysis of the Interaction between the Receptor-Binding Domain of the Omicron Variants and Human Angiotensin-Converting Enzyme 2

Simple Summary

The Omicron variant has recently been divided into BA.1, BA.2, and BA.3 subvariants. In the present study, we focused on comparing the interaction between the receptor-binding domain (RBD) of BA.1 and BA.2 spike proteins with human angiotensin-converting enzyme 2 (hACE2) using a computational approach. The RBD BA.2 was modeled after the BA.1. The results from molecular docking and molecular dynamics studies showed that RBD BA.2 has a higher and more stable affinity for hACE2 compared to RBD BA.1.

Abstract

Concerns have been raised about the high number of mutations in the spike protein of the new emergence of the highly transmissible Omicron variant (B.1.1529 lineage) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). This variant’s extraordinary ability to evade antibodies would significantly impair the current vaccination program. This present study aimed to computationally analyze the interaction between the receptor-binding domain (RBD) in the spike protein of Omicron variants and human angiotensin-converting enzyme 2 (hACE2). The docking results indicated that Omicron BA.2 has exceptionally strong interactions with hACE2 in comparison to Omicron BA.1, Delta, and wild-type, as indicated by various parameters such as salt bridge, hydrogen bond, and non-bonded interactions. The results of the molecular dynamics simulation study corroborate these findings, indicating that Omicron BA.2 has a strong and stable interaction with hACE2. This study provides insight into the development of an effective intervention against this variant.