Increasing the frequency of omicron variant mutations boosts the immune response and may reduce the virus virulence.

The coronavirus has posed a serious threat to the world since its discovery in Wuhan in 2019. Beta, gamma, delta, and the final omicron variants have emerged as a result of several mutations in the virion structure. The Australian Omicron S protein variant contains 37 mutations out of a total of 67 mutations. According to preliminary data from South Africa, Omicron variant infection is not associated with any particular symptoms. The purpose of this research was to determine how changes in the structure of the S protein alter the protein's interaction with the ACE2 receptor. The Omicron variant stimulates the immune response more than the wild strain.

Introduction

The globe is still on high COVID-19 alert over two years after the onset of the SARS-CoV-2 epidemic, which has claimed the lives of over 5 million people. Weekly, if not daily, reports of new SARS-CoV-2 mutations have been made [1,2]. Poor public health infrastructure, poor vaccination rates, and the existence of a significant number of immunocompromised persons with weakened immune defenses and a higher susceptibility to infection might provide an ideal ground for the development of novel variations [3,4]. In addition, RNA viruses are highly susceptible to mutations due to the inexistence repairing mechanisms or very rudimentary repairing during RNA synthesis. Travel restrictions have been imposed by a number of nations for visitors coming in or departing from southern African countries. South Africa informed the WHO on November 24th, 2021, that a new SARS-CoV-2 strain, had been discovered. B.1.1.529 was initially discovered in Botswana on November 11th, 2021, and in South Africa on November 14th, 2021. This variant contains many mutations that may open up a great deal of speculation about upcoming problems [5,6]. The Omicron variant's spike protein has at least 30 amino acid mutations, three tiny deletions, and one short insertion. The receptor binding domain contains 15 of the 30 amino acid mutations [7,8]. This variant has been found at higher rates than earlier outbreaks, suggesting that it may have a growth advantage. It is yet unknown if infection with the Omicron form causes more severe illness. Assessment of illness severity is challenging due to the minimal number of cases linked to the Omicron variants [9,10]. According to preliminary data from South Africa, Omicron variant infection is not associated with any particular symptoms, and some individuals are asymptomatic, as is the case with other variations [11,12].

Materials and methods

From the GISAID website, 16 strains of the SARS-CoV-2 carrying the code (B.1.1.529) of the Omicron variant were selected. Information of cases have been registered in Botswana, South Africa, Hong Kong and Australia. Spike protein was chosen because it binds to the angiotensin-converting enzyme 2 (ACE2) which is essential for viral entry into the host cell [13,14].

For the process of evaluating, identifying mutations, and examining the structural changes in the protein, the registered strain in Australia was chosen. The identified strain in Wuhan, ID (6VXX), was used as a template reference from PDB. The Australian Omicron S protein variant contains 37 mutations. Table 1 shows description of wild and mutant amino acid residues.

Table 1

Alteration, deletion and insertion of amino acids in omicron variant according to the Wuhan strain. del: deletion amino acid. ins: insertion amino acid. AA: amino acid.

Type	Residue	Wild AA	Omicron AA	Type	Residue	Wild AA	Omicron AA
H69del	69	Histidine	–	T478K	478	Threonine	Lysine
A67V	67	Alanine	Valine	E484A	484	Glutamate	Alanine
V70del	70	Valine	–	Q493R	493	Glutamine	Arginine
T95I	95	Threonine	Isoleucine	G496S	496	Glycine	Serine
G142D	142	Glycine	Aspartate	Q498R	498	Glutamine	Arginine
V143del	143	Valine	–	N501Y	501	Asparagine	Tyrosine
Y144del	144	Tyrosine	–	Y505H	505	Tyrosine	Histidine
Y145del	145	Tyrosine	–	T547K	547	Threonine	Lysine
N211del	211	Asparagine	–	D614G	614	Aspartate	Glycine
L212I	212	Leucine	Isoleucine	H655Y	655	Histidine	Tyrosine
ins214EPE	214	–	Glutamate, Proline Glutamate	N679K	679	Asparagine	Lysine
G339D	339	Glycine	Aspartate	P681H	681	Proline	Histidine
S371L	371	Serine	Leucine	N764K	764	Asparagine	Lysine
S373P	373	Serine	Proline	D796Y	796	Aspartate	Tyrosine
S375F	375	Serine	Phenylalanine	N856K	856	Asparagine	Lysine
K417 N	417	Lysine	Asparagine	Q954H	954	Glutamine	Histidine
N440K	440	Asparagine	Lysine	N969K	969	Asparagine	Lysine
G446S	446	Glycine	Serine	L981F	981	Leucine	Phenylalanine
S477 N	477	Serine	Asparagine

Constructed 3D structure of the Australian Omicron S protein

In the current study, we used two approaches to construct missing loops in the S protein (6VXX): the PyMOL Builder tool and the ModLoop Queue online site. According to the reference S protein isolated from Wuhan state, the Australian Omicron S protein 3D structure was constructed using PyMol.

Molecular docking between chain A of S protein and ACE2

HDOCK and PatchDock servers were used to investigate the interaction of S protein with the ACE2 (1R42) receptor. Based on the RMSDs and docking scores, the top five models were chosen.

Immunity system simulation

Using in silico immunological simulations, the immunogenicity and immune response profile were further defined. In this investigation, immunological simulations were performed using the C-ImmSim server, which is an agent-based model with a position-specific scoring matrix (PSSM) [15,16].

Applied RBD mutations

RBD-ACE2 (PDB: 6M0J) was subjected to more scrutiny. Since the RBD is the primary section of S protein due to its relationship with the (ACE2) receptor, we assigned mutations to it and noted the effect of those changes in the form of interaction. 15 of the 37 protein S mutations were carried out inside the RBD residues (333–526). To compare the distances between the wild and the mutated strains, we assessed the distance between atoms in the docking residues (RBD-ACE2). The immunogenicity of RBD-ACE2 wild and mutant strains was determined using the C-ImmSim server.

Results

The residues created using Builder PyMOL and ModLoop in the current investigation had no discernible differences. There are obvious conformational changes between wild type and Omicron S protein variants, fig. 1 .

Fig. 1

Cartoon representation of S protein. A: Wild strain. B: Omicron variant strain shows the 37 mutation residues in different color dots. C: Chain A of wild strain shows wild residues. D: Mutated chain A of Omicron variant.

Predict molecular docking of chain A and ACE2 receptor

The output models were evaluated after the ACE2 receptor was uploaded to chain A of S protein. The best 5 models were chosen after examining the docking outcomes and calculating the RMSD. Fig. 2 depicts the ACE2 receptor's ability to bind with the A chain of wild-type and Omicron S proteins. The protein model quality was measured between 3.0 ≤ LGscore <5.0. The LGscore for wild S protein is (4.615) while in the omicron variant was higher (4.450). On the other hand, the MaxSub score were (0.355) and (0.333) respectively. The energy docking score of the omicron variant is higher than the wild type of all predicted models. On the contrary, for the RMSD, a discrepancy was observed in the proposed models for the two types, as shown in Table 2 .

Fig. 2

The best 5 docking models of S protein and ACE2 receptor. 3D green structure is the wild type of chain A. 3D magenta structure is the omicron chain A.

Table 2

The docking energy score and RMSD in all predict models of wild and omicron variant's S protein.

Rank	Model 1	Model 2	Model 3	Model 4	Model 5
Wild S protein (Docking score)	−245.91	−245.69	−242.93	−236.36	−235.79
Omicron S variant (Docking score)	−275.32	−249.02	−247.10	−246.39	−241.10
Wild RMSD	268.66	263.14	243.94	285.37	249.59
Omicron RMSD	266.54	243.85	323.89	213.96	300.68

Immunogenicity analysis

The immunogenicity and the immune response of the predicted molecular docking of (chain A-ACE2) were determined using immunoinformatics for both wild (6VXX) and omicron variant. Despite having equal immunological characteristics, there was a distinct difference between the wild type and the omicron when the outcomes of the immune response were analyzed. The duration of the immune response has been shorter for the omicron variant. As a result of the immunogenicity, the antibody titers are greater for the omicron variant including (IgM + IgG), (IgG1+IgG2), and IgG1 alone, Fig. 3 . There is an elevation in the count of active B-lymphocytes between the wild and the mutant type, yet they are comparable in terms of presenting on class-II, internalized Ag, replicating, and being anergic, Fig. 4 . Count of T-regulatory memory lymphocyte is higher in the S protein wild type and on the contrary for the count of T-regulatory not memory lymphocytes. The activity and resting of omicron reaches 97 (cells per mm3) while it reaches 87(cells per mm3) for wild type. The regulatory T-cell (not memory) is higher for omicron, Fig. 5 . The number of cytokines and interleukins stimulated by the S protein is similar in the both types. IFN-g and TGF-b are the highest types stimulated by the immune system for both types of S protein. IFN-γ of the wild type is 620,000 ng/ml and 630,000 ng/ml for omicron, Fig. 6 .

Fig. 3

Immunogenicity and immune response against wild and omicron S protein. Ag count per ml and all Ab titers through 30 days of infection. Peaks of (IgM + IgG), (IgG1+IgG2), and IgG1 of omicron is higher than of wild type.

Fig. 4

B-lymphocytes population per entity-state, showing counts for active, presenting on class-II, internalized the Ag, duplicating and anergic. The number of presenting class II in omicron is higher but the count activity is less than wild type.

Fig. 5

Immunogenicity and immune response against wild and omicron S protein. CD4 T-regulatory lymphocytes count. Both total, memory and per entity-state counts are plotted. The activity and resting of omicron reaches 97 (cells per mm3) while it reaches 87(cells per mm3) for the wild type. Regulatory T-cell not memory is higher for omicron.

Fig. 6

Concentration of cytokines and interleukins. D in the inset plot is danger signal. IFN-γ of the wild type is 620,000 ng/ml and 630,000 ng/ml for omicron.

Performed RBD-ACE2 mutations and immunity response

PDB (6M0J) is the molecular interaction between RBD and ACE2. 15 mutations were constructed on the wild type to predict new omicron variant (B.1.1.529). These mutations are: E484A, G339D, G446S, G496S, K417 N, N440K, N501Y, Q493R, Q498R, S371L, S373P, S375F, S477 N, T478K, and Y505H. 10 out of 15 of these mutations are closer to the docking region. As a result of these changes, the shape of the micron variant differs significantly from the original type, which might lead to a change in the strength of the fusion between RBD and ACE2, Fig. 7 . In terms of the immune response, our research found that various characteristics, such as Ab titers amount, IgM, IgG1 and IgG2 arbitrary scale, CD8 T-cytotoxic lymphocytes count per entity-state. There is no significant difference in account of antibodies produced and stimulated by the antigens per ml between wild type and omicron variants, Fig. 8 . Anergic CD8 T-cell peak reaches 1180 cells/mm3 while it is 1390 cells/mm3 for omicron variant, Fig. 9 . T-regulatory lymphocyte count (CD4), memory and per entity-state (non-memory) counts reach 88 and 86 cells/mm3, respectively for wild type. TR memory and non-memory cells reach 91 and 89 cells/mm3, respectively for omicron variant, Fig. 10 . The levels of cytokines and interleukins, including as IFN-, IL-12, TGF-b, and IL-10, are similar in wild and omicron variants, Fig. 11 .

Fig. 7

Applied 15 mutations on the RBD-ACE2 on the wild (6M0J). RBD: light brown color. ACE2: cyan color. A: RBD-ACE2 wild (6MOJ) shows 15 residues in res dots. B: RBD-ACE2 omicron variant shows 15 mutated resides in different colors.

Fig. 8

Immunogenicity and immune response against wild and omicron of RBD protein. Ag count per ml and all Ab titers through 30 days of infection. No significant difference in account of Ab produced and stimulated by the Ag/ml between wild type and omicron variants.

Fig. 9

CD8 T-cytotoxic lymphocytes count per entity-state, showing counts for active, presenting on class-II, internalized the Ag, duplicating and anergic. Anergic CD8 T-cell peak reaches 1180 cells/mm3 while it is 1390 cells/mm3 for omicron variant.

Fig. 10

Immunogenicity and immune response against wild and omicron of RBD protein. CD4 T-regulatory lymphocytes count. Both total, memory and per entity-state counts are plotted. TR regulatory memory and non-memory cells reache 88 and 86 cells/mm3, respectively for wild type. TR memory and non-memory cells reach 91 and 89 cells/mm3, respectively for omicron.

Fig. 11

Concentration of cytokines and interleukins. D in the inset plot is danger signal. IFN-γ, IL-12, TGF-b and IL-10 are similar between wild and omicron variants.

Discussion

The rapid appearance of a new coronavirus strain alarmed the globe, particularly because it included several mutations in all of its proteins [17]. The virus's efficacy and the severity of illness symptoms are both severely impacted by the abrupt shift brought on by mutations. So far, no unusual trends in the number of deaths or viral efficacy have been seen with the new Omicron type [18,19]. The fact that this strain affects youngsters more than adults is noteworthy, but happily, the number of deaths is modest, suggesting that this strain is the weakest, as future studies will reveal. According to the WHO, omicron has a greater propensity to spread than prior strains [20].

There is one significant issue: the number of mutations in Omicron varies, starting at 62 in South Africa and increasing to 67 in the strain isolated from Australia that we are studying. The S protein in the omicron variant has 37 mutations out of a total of 67 mutations. In addition to the replacement mutations that arise in each strain, we've seen a loss of certain residues and an insertion of new base pairs. Because the S protein is involved in the insertion of viral RNA into the host cell via the ACE2 receptor, it was chosen as the subject of our research.

In the current study, we presented 37 out of 67 total mutations in the viral proteins. Based on the reference strain isolated from Wuhan, the protein S model of the Omicron variation was created. Mutations result in a distinct alteration in the structure of the protein. To determine the effect of the degree of fusion on the mutant and wild type, the docking stage between S protein and the ACE2 receptor was required. In all of the predicted docking models, the Omicron variant had a higher fusion energy than the wild strain. Table 1 shows that the degree of RMSD differed between the two types.

We demonstrated in this study that the novel strain activates the immune system and produces a distinct immunological response, as opposed to the wild strain, which has higher stimulation for both IgG and IgM.

To examine the influence of the Omicron variation on vaccination efficacy and outbreak infections, especially in those who have received booster doses, more laboratory and epidemiological research is needed. Vaccination, on the other hand, is expected to continue to provide protection against hospitalization and mortality, and vaccines will continue to play an important part in the COVID-19 pandemic's containment.

15 of the 37 mutations are found in RBD region of the S protein. This region is well recognized for its involvement in docking with the ACE2 receptor. We demonstrated a noticeable difference in the form of the mutant strain's RBD compared to the original strain in this study. The degree of immunological response in the Omicron strain differs significantly from that in the original strain.

Individuals over 60 years of age who are not completely vaccinated or do not have confirmation of past SARS-COV-2 infection, as well as those with underlying health issues, should postpone travel, according to the WHO [21,22].

Because Omicron has been designated as a Variant of Concern [23], WHO recommends that countries improve surveillance and sequencing of cases, share genome sequences on publicly available databases like GISAID [24], report initial cases or clusters to WHO [25,26], and conduct field investigations and laboratory assessments to better understand if Omicron has different transmission or disease characteristics, or has an impact on vaccine effectiveness.

Despite the Omicron strain's increased dissemination, we believe it is the weakest of the ancient strains, and that it will exhaust the virus, causing its symptoms to revert to those of a common cold, as did the ancestors of the Human coronavirus 229E (HCoV-229E).

Conclusions

Omicron variant has an increase in the number of mutations daily. One of the strain of the omicron has 67 mutations in all proteins. 37 out of 67 presented in the S protein 15 out of 37 of S mutation protein found in the RBD region. This study revealed that mutations have conformational changes in the S protein shape. The Omicron variant stimulates the immune response more than the wild strain. To effectively contain this pandemic, a global concerted effort including government agencies, pharmaceutical, and biotech firms, as well as academic and healthcare institutions, is necessary.

Author statement

Ali A. Dawood: Data acquisition; Data analysis / interpretation; Drafting manuscript; Critical revision of manuscript; Statistical analysis; Securing funding; Admin; technical; or material; support; Supervision; Final approval.