key: cord-0707338-sl0ybqfm authors: Kumar, Suresh; Thambiraja, Thiviya S.; Karuppanan, Kalimuthu; Subramaniam, Gunasekaran title: Omicron and Delta Variant of SARS-CoV-2: A Comparative Computational Study of Spike protein date: 2021-12-03 journal: bioRxiv DOI: 10.1101/2021.12.02.470946 sha: 012089e95e63df60ec5d3fd635ae3d1315ec09c1 doc_id: 707338 cord_uid: sl0ybqfm Emerging SARS-CoV-2 variants, especially those of concern, may have an impact on the virus’s transmissibility and pathogenicity, as well as diagnostic equipment performance and vaccine effectiveness. Even though the SARS-CoV-2 Delta variant (B.1.617.2) emerged during India’s second wave of infections, Delta variants have grown dominant internationally and are still evolving. On November 26, 2021, WHO identified the variant B.1.1.529 as a variant of concern, naming it Omicron, based on evidence that Omicron contains numerous mutations that may influence its behaviour. However, the mode of transmission and severity of the Omicron variant remains unknown. We used computational studies to examine the Delta and Omicron variants in this work and found that the Omicron variant had a higher affinity for human ACE2 than the Delta variant due to a significant number of mutations in the SARS-CoV-2 receptor binding domain, indicating a higher potential for transmission. Based on docking studies, the Q493R, N501Y, S371L, S373P, S375F, Q498R, and T478K mutations contribute significantly to high binding affinity with human ACE2. In comparison to the Delta variant, both the entire spike protein and the RBD in Omicron include a high proportion of hydrophobic amino acids such as leucine and phenylalanine. These amino acids are located within the protein’s core and are required for structural stability. Omicron has a higher percentage of alpha-helix structure than the Delta variant in both whole spike protein and RBD, indicating that it has a more stable structure. We observed a disorder-order transition in the Omicron variant between spike protein RBD regions 468-473, and it may be significant in the influence of disordered residues/regions on spike protein stability and binding to ACE2. A future study might investigate the epidemiological and biological consequences of the Omicron variant. SARS-CoV-2 (Severe acute respiratory syndrome coronavirus type 2) is a coronavirus that caused the Covid-19 disease outbreak in late 2019 in Wuhan China. By early 2020, the disease had rapidly spread across the world and was declared a global pandemic as a public health emergency of international concern. The virus spreads from person to person by respiratory droplets in close contact between sick and asymptomatic people (within 6 feet) 1 . Transmission by aerosols and maybe contact with fomites is also a possibility, although this is not considered to be the most probable route 2 . SARS-CoV-2 pathogenesis is dependent on the viral spike protein binding to angiotensin-converting enzyme 2 (ACE2) receptors, with cell entrance required ACE2 receptor cleavage by a type 2 transmembrane serine protease to activate the viral spike protein 3 . COVID-19 individuals have a wide range of clinical symptoms, from moderate to severe, fast progressive, and acute disease 4 . The diagnosis of COVID-19 is non-specific, and the virus may manifest itself in a variety of ways, ranging from no symptoms (asymptomatic) to severe pneumonia and death. The CoVID-19 pandemic response plan is based on the development of therapeutic alternatives and vaccination formulations [5] [6] [7] . The term "variant of concern" (VOC) for SARS-CoV-2 (which produces COVID- 19) refers to viral variants with mutations in their spike protein receptor-binding domain (RBD) that dramatically improve binding affinity in the RBD-hACE2 complex while also causing fast dissemination in human populations 8 . Increased viral replication increases the likelihood of SARS-CoV-2 mutations forming. Therefore, the only option to end the pandemic is for effective vaccinations against circulating variations to be extensively and fairly delivered globally. Because raising nations are rushing to vaccinate their people within months, they risk SARS-CoV-2 evolving into a new lineage that vaccines may not be able to protect against in other countries. To combat some emerging SARS-CoV-2 strains, new vaccinations may need to be developed regularly. With the introduction of extremely infectious SARS-CoV-2 variants, greater vaccine penetration will be required to build protective immunity, and children may also need to be vaccinated 9 . While the majority of SARS-CoV-2 sequence changes are projected to be detrimental and swiftly removed or to be neutral, a small number are predicted to influence functional characteristics, possibly modifying infection rate, disease severity, or interactions with the host immune system 10 . Nonetheless, beginning in late 2020, the development of SARS-CoV-2 has been marked by the introduction of 'variants of concern,' or changes in viral properties such as disease transmission and antigenicity, most likely because of the changing immunological composition of the human species. The delta variant (B.1.617.2) was discovered for the first time in India in late 2020. The Delta version may have invaded over 163 nations by August 24, 2021. The World Health Organization (WHO) stated in June 2021 that the Delta strain is on its way to becoming the most prevalent strain in the world 11 . Therefore, the Delta variant was changed from Variant of Interest (VOI) to Variant of Concern (VOC). According to present evidence, the SARS-CoV-2 Delta variant of concern (VOC) is 40-60% more transmissible than the Alpha (B.1.1.7) VOC and may be associated with an increased risk of hospitalisation. The Delta VOC mostly endangers those who are unvaccinated or just partially vaccinated 12 . Several SARS-CoV-2 variants have been identified, posing a long-term infection risk in immunocompromised individuals 26 . The "variant of concern" (VOC) SARS-CoV-2 (which generates refers to viral variants in which mutations in the spike protein receptor-binding domain (RBD) drastically increase binding affinity in the RBD-hACE2 complex while also being connected to rapid transmission in human populations 27 . A variety of computational approaches are utilised in this study to compare the currently categorised omicron variant to the delta variant to get insight into its characteristics and binding affinity for the ACE2 protein. The omicron variation includes 30 mutations in the Spike protein, half of which are in the receptor-binding domain, according to the multiple alignments (Figure 1) . From a previous study, it is observed that RBD T470-T478 loop and Y505 as viral determinants for specific recognition of SARS-CoV-2 RBD by ACE2 28 . T478 is a common mutation seen in Delta and Omicron variants (Figure 2 ). RBD has the potential to be developed into an efficient and safe subunit vaccine against SARS-CoV-2 due to its ability to produce very robust nAb responses. Numerous mutations in the spike protein's receptor-binding region in Omicron compared to the Delta variant ( Figure 3 ) suggests that the Omicron variant may be immunologically resistant to antibody-mediated protection ( Table 1) . While Wuhan-Hu-1 has 1273 amino acids, the delta variant has 1271 and the omicron variant has 1270; nevertheless, due to sequence loss, both the delta and omicron variants have a few fewer residues than the wild type. A protein's isoelectric point (pI) is the pH value at which its surface is completely charged but its net charge is zero. A pI value of more than 7 indicates that the protein is alkaline, whereas a value less than 7 indicates that it is acidic. The molecular weight of Wuhan-Hu-1 is 141178.47 with a theoretical pI of 6.24, the delta variant is 140986.31 with a theoretical pI of 6.78, and the omicron variant is 141328.11 with a theoretical pI of 7.14. Despite having three fewer amino acids than Wuhan-Hu-1, the omicron variant has a higher molecular weight and theoretical PI than the delta variant and Wuhan-Hu-1. In the current investigations, the omicron variant is expected to have an alkaline pI, while the delta and Wuhan-Hu-1 variants are expected to have an acidic pI. According to previous research 29 a stability score of less than 40 indicates that the protein structure is stable. A value of 40 or above suggests that the protein is structurally unstable. In our research, the range remained 32.81-34.69, indicating the great stability of all SARS-CoV-2 spike proteins. The average extinction coefficient is 11238.61, which indicates how much light the protein can absorb at 280 nm. The aliphatic index measures the volume of a protein that is filled by aliphatic amino acids on the side chain, such as alanine. A high aliphatic index of 84.50 to 84.95 indicates that the protein is temperature stable across a wide temperature range. The greater a protein's aliphatic index, the more thermostable it is. The degree to which amino acids in a protein sequence are hydrophobic or hydrophilic is referred to as hydropathicity. A protein with a low GRAVY (Grand average of hydrophobicity) value is nonpolar and has a stronger affinity for water, indicating that it is intrinsically hydrophilic. (Lys), Aspartic acid (Asp), and Glutamic acid (Glu), indicating that the omicron has more charged residues that contribute to salt bridge formation and that charged residues are exposed to a much greater degree. The higher amino acid composition of Phenylalanine (F), Isoleucine (I) in the omicron spike protein, when compared to the delta variant, suggests that the omicron spike protein includes more hydrophobic amino acids, which may be due to its positioning inside the protein core. When compared to the delta version, the omicron variant's amino acid composition is low in polar amino acids such as Asparagine (N), Glutamine (Q). Omicron RBD is high in non-polar amino acids such as Leucine (L), Phenylalanine (F), and Proline (P) 30 . These residues are located inside the protein core and are thus inaccessible to the solvent ( Table 2 ). Omicron has a higher fraction of alpha-helix structure (23.46%) than delta variant (22.03%), but less extended strand and random coil structure (Table 3) . Because it is largely made up of alpha-helices, the evidence suggests that the omicron variant protein is structurally highly stable. The Omicron form of RBD has a greater alpha helix composition than the Delta variation in secondary structure prediction, suggesting that the RBD has a stable structure, however, the random coil composition is slightly increased in Omicron 31 . Disordered areas of viral proteins are linked to viral pathogenicity and infectivity. PONDR® VLXT was used to predict the intrinsic disorder of Wuhan-Hu-1, Delta, and Omicron variants. Residues with anticipated disorder scores more than 0.5 are regarded inherently disordered, while residues with expected disorder values between 0.2 and 0.5 are considered flexible. According to the prediction, the Omicron variant has a less disordered area than the Delta variant and the wild type. We observed that disordered regions in entire spike protein as well as RBD in Omicron exhibit disorder-to-order transition when compared to Delta variant and wild type. According to prior research from the cryo-EM structure of Wuhan-Hu-1, the T470-F490 loop and Q498-Y505 within RBD are key contacting elements that interact with RBD and ACE2 28 . The disorder prediction ranges from 468-473 with residues ISTEIYQA in Wuhan-Hu-1-RBD, 469-471 with residues EIY in Delta variant-RBD, and there are no disorder residues predicted in this region in Omicron Variant-RBD. This implies that there is a chance of disorder-order transition between region 468-473 of spike protein, which could be important in the influence of disordered residues/regions on spike protein stability and binding to ACE2. (Table 4 ). An I-Mutant protein stability study predicted that all amino acid modifications in the delta variant reduce spike protein stability. Except for the N501Y mutation 32 , which is expected to improve the stability of the spike protein, all amino acid changes in the omicron variant result in a decrease in stability (Table 5 ). SIFT analysis revealed that, whereas the delta variant D950N impairs protein function, other mutations are tolerated. The N211I, Y505H, and N764K mutations in the omicron variant impair protein function, although other variants are tolerated (Table 5) . Although the RBD L452R and T478 delta mutations are tolerated, they reduce protein stability and increase disease risk. There are 15 Omicron variant mutations in RBD, the N501Y mutation being one of them. It is tolerated and enhances protein stability; however, it is disease-prone. Other mutations that decrease protein stability and increase disease risks, such as G339D, S371L, S373P, S375F, N440K, G446S, T478K, G496S, and Q498R, are tolerated. Tolerable mutations include K417N, S477N, E484A, Q493R, and Q498R, which decrease protein stability and increase disease vulnerability. Protein function is impaired by Y505H mutations, resulting in decreased protein stability and an increased risk of disease ( Table 6 ). The large type I transmembrane S glycoprotein on the viral envelope and the homologous receptor on the surface of host cells enable membrane fusion. The S glycoprotein's exposed surface not only allows membrane fusion but also drives host immune responses, making it a great target for neutralising antibodies 33 . Cleavage at the S1/S2 site results in the formation of a surface subunit S1, which attaches the virus to the host cell surface receptor, and a transmembrane component S2, which allows the viral and host cell membranes to merge. The S2 subunit of the transmembrane is made up of an N-terminal hydrophobic fusion peptide (FP), two heptad repeats (HR1 and HR2), a transmembrane domain (TM), and a cytoplasmic tail (CT), in the following order: FP-HR1-HR2-TM-CT 34 . The amino acid residues Y695, I923, S982, V1189, F1220, and I1221 in Wuhan-Hu-1 are extremely flexible. Residues I921, S980, V1187, F1218, and I1219 are extremely variable in the Delta variant. The I920, S979, V1186, F1217, and I1218 residues are particularly flexible in the Omicron variant, as predicted by PredyFlexy. Flexible prediction and local structure prediction from sequence show that the heptapeptide repeat sequence 1 (HR1) (912-984 residues), HR2 (1163-1213 residues), and TM domain (1213-1237 residues) of the S2 subunit are very flexible in both the Delta and Omicron variants. Understanding the SARS-CoV-2 virus's receptor recognition mechanism is critical since it governs the virus's infectivity, host range, and pathogenesis. The binding affinity of SARS-CoV-2 variants of RBD to ACE2 differs because of minor variations in ACE2 interactions. In this study, the PDB (6M0J) crystal structure of the SARS-CoV-2 spike receptor-binding domain associated with ACE2 was employed. The RBD of SARS-CoV-2 was isolated from ACE2 and used for protein-protein docking. Hex was used to dock SARS-CoV-2 RBD and ACE2, and its docking score (-500.37) is for Wuhan-Hu-1 (Wild type), which was utilised to compare the docking energies of Delta and Omicron. The pymol mutagenesis wizard was used to add mutations into the Delta and Omicron versions. Docking was performed using hACE2 between the Delta and Omicron Variant (Figure 4) . The docking score for the Omicron variation is the highest (-539.81), while the Delta variant is the lowest (-529.62). This suggests that the Omicron variant is more responsive to hACE2 than the Delta variant, indicating a higher potential for transmission (Table 7 ). In addition, the impact of each changed residue on hACE2 affinity was investigated. The highest binding affinity score of all 15 RBD mutations is Q493R (-581.53), followed by N501Y (-560.81), S371L (-54.34), S373P (-541.87), S375F (-530.07), Q498R (-527.38), and T478 (-517.03) ( Table 8) . For the Delta version, only two mutations were found in RBD, with L452R (-517.52) having the highest binding affinity, followed by T478 (-517.03). Point mutations at key residues have a significant impact on the interaction with ACE2. SARS-CoV-2 interacts with hACE2 through its C-terminal domain (SARS-CoV-2-CTD), indicating that it has a higher affinity for the receptor. In SARS-CoV-2-CTD, E484 forms ionic contacts with K31, increasing receptor affinity. Previous research has found that the single mutation E484 in the viral spike (S) protein (which is shared by the Beta and Gamma VOCs as well as the Mu VOI) may be critical in avoiding vaccination immunity; variants with the E484 mutation have demonstrated resistance to neutralising antibodies generated by prior infection 35 . The E484A mutation (478.49) has binding affinity in the Omicron variant may result in enhanced hACE2 binding. The Omicron form of ACE2 binds more strongly to SARS-CoV-2 than the Delta variant of hACE2. The existence of a high number of Omicron variant mutations is also a hallmark of the variants, indicating that viral evolution in immunocompromised persons may have played a significant role in their development. Because many people worldwide suffer from inherent or induced immunosuppression, the relationship between immunosuppression and the generation of highly transmissible or pathogenic SARS-CoV-2 variants must be investigated further and mitigation strategies devised. Both the Omicron and Delta variants were investigated in this study using several computational tools and a computational saturation mutagenesis model, examining structural, sequence-driven, and dynamic changes that effect overall protein stability were examined. According to the findings of this study, large changes in the RBD region of the Omicron variant contribute to higher binding with hACE2, which may result in a higher transmission rate when compared to the Delta variant. l e 3 : S e c o n d a r y s t r u c t u r e p r e d i c t i o n a n d c o m p a r i s o n o f D e l t a a n d O m i c r o n v a r i a n t w i t h r e f e r e n c e t o w i l d t y p e ( W u h 2 ) i s t o g r a m a s i n i n v o l v e s 2 b a s e c h a n g e s i n t h e c o r r e s p o n d i n g D N A c o d o n . G 3 3 9 D D e c r e a s e s t a b i l i t y -1 . 1 6 T O L E R A T E D ( 0 . 8 8 ) 1 . 9 3 S 3 7 1 L D e c r e a s e s t a b i l i t y -0 . The newly emerged COVID-19 disease: a systemic review Transmissibility and transmission of respiratory viruses Drug and Vaccine Design against Novel Coronavirus (2019-nCoV) Spike Protein through Computational Approach COVID-19: A Drug Repurposing and Biomarker Identification by Using Comprehensive Gene-Disease Associations through Protein-Protein Interaction Network Analysis Protein-Protein Interaction Network for the Identification of New Targets Against Novel Coronavirus Online Resource and Tools for the Development of Drugs Against Novel Coronavirus COVID-19 Vaccine Candidates by Identification of B and T Cell Multi-Epitopes Against SARS-COV-2 Emergence and spread of SARS-CoV-2 lineage B.1.620 with variant of concern-like mutations and deletions Vaccination strategies to combat novel corona virus SARS-CoV-2 Evaluation of the Effect of D614g, N501y and S477n Mutation in Sars-Cov-2 through Computational Approach Delta variant: What is happening with transmission, hospital admissions, and restrictions? Impact of the Delta variant on vaccine efficacy and response strategies Virus Pathogen Database and Analysis Resource (ViPR): A Comprehensive Bioinformatics Database and Analysis Resource for the Coronavirus Research Community Global initiative on sharing all influenza data -from vision to reality. Euro surveillance : bulletin Europeen sur les maladies transmissibles = European communicable disease bulletin ExPASy: the proteomics server for in-depth protein knowledge and analysis GOR V server for protein secondary structure prediction Clustal Omega, accurate alignment of very large numbers of sequences PONDR-FIT: a metapredictor of intrinsically disordered amino acids 0: predicting stability changes upon mutation from the protein sequence or structure PredyFlexy: flexibility and local structure prediction from sequence SIFT: Predicting amino acid changes that affect protein function VarSite: Disease variants and protein structure We acknowledge Faculty of Health and Life Sciences, Management and Science University, Malaysia; Division of Cardiovascular Medicine, Radcliffe Department of Medicine, Wellcome Centre for Human Genetics, University of Oxford, Oxford, UK and Department of Physiology, Anatomy and Genetics, University of Oxford, Oxford, UK. We also acknowledge GISAID (https://www.gisaid.org/) for facilitating open data sharing.