1
|
Tariq SS, Zia K, Nur-E-Alam M, Nerukh D, Farafonov VS, Ul-Haq Z. Impact of mutations in SARS-CoV-2 recombinant sub-variant XBB.1.16 on the binding affinity with human ACE2 receptor. J Mol Graph Model 2024; 131:108813. [PMID: 38885553 DOI: 10.1016/j.jmgm.2024.108813] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/27/2023] [Revised: 06/12/2024] [Accepted: 06/12/2024] [Indexed: 06/20/2024]
Abstract
Despite the waning threat of the COVID-19 pandemic, its detrimental impact on global health persists. Regardless of natural immunity or immunity obtained through vaccination, emerging variants of the virus continue to undergo mutations and propagate globally. The persistent mutations in SARS-CoV-2, along with the subsequent formation of recombinant sub-variants has become a challenge for researchers and health professionals, raising concerns about the efficacy of current vaccines. Gaining a better understanding of the biochemical interactions between the Spike Protein (RBD) of SARS-CoV-2 variants and the human ACE2 receptor can prove to be beneficial in designing and developing antiviral therapeutics that are equally effective against all strains and emerging variants. Our objective in this study was to investigate the interfacial binding pattern of the SARS-CoV-2 RBD-ACE2 complex of the Wild Type (WT), Omicron, and the Omicron recombinant sub-variant XBB.1.16. We aimed to examine the atomic level factors and observe how mutations influence the interaction between the virus and its host using Molecular Dynamics simulation, MM/GBSA energy calculations, and Principal Component Analysis. Our findings reveal a higher degree of structural deviation and flexibility in XBB.1.16 compared to WT and Omicron. PCA indicated a wider cluster and significant flexibility in the movements of XBB.1.16 which can also be observed in free energy landscapes, while the normal mode analysis revealed converging motions within the RBD-ACE2 complexes which can facilitate the interaction between them. A pattern of decreased binding affinity was observed in case of XBB.1.16 when compared to the WT and Omicron. These observed deviations in XBB.1.16 when compared to its parent lineage Omicron, and WT can be attributed to the mutations specific to it. Collectively, these results enhance our understanding of the impact of mutations on the interaction between this strain and the host, taking us one step closer to designing effective antiviral therapeutics against the continually mutating strains.
Collapse
Affiliation(s)
- Syeda Sumayya Tariq
- Dr. Panjwani Center for Molecular Medicine and Drug Research, International Center for Chemical and Biological Sciences, University of Karachi, Karachi, 75270, Pakistan
| | - Komal Zia
- Dr. Panjwani Center for Molecular Medicine and Drug Research, International Center for Chemical and Biological Sciences, University of Karachi, Karachi, 75270, Pakistan
| | - Mohammad Nur-E-Alam
- Department of Pharmacognosy, College of Pharmacy, King Saud University, P.O. Box. 2457, Riyadh, 11451, Saudi Arabia
| | - Dmitry Nerukh
- Department of Mathematics, Aston University, Birmingham, B4 7ET, United Kingdom
| | - Vladimir S Farafonov
- Department of Mathematics, Aston University, Birmingham, B4 7ET, United Kingdom; Department of Physical Chemistry, V.N.Karazin Kharkiv National University, Kharkiv, 61022, Ukraine
| | - Zaheer Ul-Haq
- Dr. Panjwani Center for Molecular Medicine and Drug Research, International Center for Chemical and Biological Sciences, University of Karachi, Karachi, 75270, Pakistan.
| |
Collapse
|
2
|
Abdel-Naim AB, Kumar P, Bazuhair MA, Rizg WY, Niyazi HA, Alkuwaity K, Niyazi HA, Alharthy SA, Harakeh S, Haque S, Prakash A, Kumar V. Computational insights into dynamics and conformational stability of N-acetylmannosamine kinase mutations. J Biomol Struct Dyn 2024:1-11. [PMID: 38502682 DOI: 10.1080/07391102.2024.2323702] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/02/2023] [Accepted: 02/21/2024] [Indexed: 03/21/2024]
Abstract
The activity of UDP-GlcNAc 2-epimerase/ManNAc kinase (GNE) is essential for the biosynthesis of sialic acid, which is involved in cellular processes in health and diseases. GNE contains an N-terminal epimerase domain and a C-terminal kinase domain (N-acetylmannosamine kinase, MNK). Mutations of the GNE protein led to hypoactivity of the enzyme and cause sialurea or autosomal recessive inclusion body myopathy/Nonaka myopathy. Here, we used all-atom molecular dynamics (MD) simulations to comprehend the folding, dynamics and conformational stability of MNK variants, including the wild type (WT) and three mutants (H677R, V696M and H677R/V696M). The deleterious and destabilizing nature of MNK mutants were predicted using different prediction tools. Results predicted that mutations modulate the stability, flexibility and function of MNK. The effect of mutations on the conformational stability and dynamics of MNK was next studied through the free-energy landscape (FEL), hydrogen-bonds and secondary structure changes. The FEL results show that the mutations interfere with various conformational transitions in both WT and mutants, exposing the structural underpinnings of protein destabilization and unfolding brought on by mutation. We discover that, when compared to the other two mutations, V696M and H677R/V696M, H677R has the most harmful effects. These findings have a strong correlation with published experimental studies that demonstrate how these mutations disrupt MNK activity. Hence, this computational study describes the structural details to unravel the mutant effects at the atomistic resolution and has implications for understanding the GNE's physiological and pathological role.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Ashraf B Abdel-Naim
- Department of Pharmacology and Toxicology, Faculty of Pharmacy, King Abdulaziz University, Jeddah, Saudi Arabia
- Mohamed Saeed Tamer Chair for Pharmaceutical Industries, Faculty of Pharmacy, King Abdulaziz University, Jeddah, Saudi Arabia
- Center of Excellence for Drug Research and Pharmaceutical Industries, King Abdulaziz University, Jeddah, Saudi Arabia
| | - Pawan Kumar
- School of Computational & Integrative Sciences, Jawaharlal Nehru University, New Delhi, India
| | - Mohammed A Bazuhair
- Center of Excellence for Drug Research and Pharmaceutical Industries, King Abdulaziz University, Jeddah, Saudi Arabia
- Department of Clinical Pharmacology, Faculty of Medicine, King Abdulaziz University, Jeddah, Saudi Arabia
| | - Waleed Y Rizg
- Mohamed Saeed Tamer Chair for Pharmaceutical Industries, Faculty of Pharmacy, King Abdulaziz University, Jeddah, Saudi Arabia
- Center of Excellence for Drug Research and Pharmaceutical Industries, King Abdulaziz University, Jeddah, Saudi Arabia
- Department of Pharmaceutics, Faculty of Pharmacy, King Abdulaziz University, Jeddah, Saudi Arabia
| | - Hatoon A Niyazi
- Department of Pharmacology and Toxicology, Faculty of Pharmacy, King Abdulaziz University, Jeddah, Saudi Arabia
| | - Khalil Alkuwaity
- Department of Pharmacology and Toxicology, Faculty of Pharmacy, King Abdulaziz University, Jeddah, Saudi Arabia
- Vaccines and Immunotherapy Unit, King Fahd Medical Research Center, King Abdulaziz University, Jeddah, Saudi Arabia
| | - Hanouf A Niyazi
- Department of Pharmacology and Toxicology, Faculty of Pharmacy, King Abdulaziz University, Jeddah, Saudi Arabia
| | - Saif A Alharthy
- Department of Medical Laboratory Sciences, Faculty of Applied Medical Sciences, King Abdulaziz University, Jeddah, Saudi Arabia
- Toxicology and Forensic Sciences Unit, King Fahd Medical Research Center, King Abdulaziz University, Jeddah, Saudi Arabia
| | - Steve Harakeh
- King Fahd Medical Research Center, King Abdulaziz University, Jeddah, Saudi Arabia
- Yousef Abdul Latif Jameel Scientific Chair of Prophetic Medicine Application, Faculty of Medicine, King Abdulaziz University, Jeddah, Saudi Arabia
| | - Shafiul Haque
- Research and Scientific Studies Unit, College of Nursing and Health Sciences, Jazan University, Jazan, Saudi Arabia
- Gilbert and Rose-Marie Chagoury School of Medicine, Lebanese American University, Beirut, Lebanon
- Centre of Medical and Bio-Allied Health Sciences Research, Ajman University, Ajman, United Arab Emirates
| | - Amresh Prakash
- Amity Institute of Integrative Sciences and Health, Amity University Haryana, Gurgaon, India
| | - Vijay Kumar
- Amity Institute of Neuropsychology & Neurosciences, Amity University, Noida, India
| |
Collapse
|
3
|
Haque S, Mathkor DM, Alkhanani MF, Bantun F, Momenah AM, Faidah H, Jalal NA, Kumar V. Comprehensive deep mutational scanning reveals the pH induced stability and binding differences between SARS-CoV-2 spike RBD and human ACE2. J Biomol Struct Dyn 2023; 41:15207-15218. [PMID: 36995177 DOI: 10.1080/07391102.2023.2194007] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/14/2022] [Accepted: 02/25/2023] [Indexed: 03/31/2023]
Abstract
The SARS-CoV-2 spike (S) glycoprotein with its mobile receptor-binding domain (RBD), binds to the human ACE2 receptor and thus facilitates virus entry through low-pH-endosomal pathways. The high degree of SARS-CoV-2 mutability has raised concern among scientists and medical professionals because it created doubt about the effectiveness of drugs and vaccinations designed specifically for COVID-19. In this study, we used computational saturation mutagenesis approach, including structure-based free energy calculations to analyse the effects of the missense mutations on the SARS-CoV-2 S-RBD stability and the S-RBD binding affinity with ACE2 at three different pH (pH 4.5, pH 6.5, and pH 7.4). A total of 3705 mutations in the S-RBD protein were analyzed, and we discovered that most of these mutations destabilize the RBD protein. Specifically, residues G404, G431, G447, A475, and G526 were important for RBD protein stability. In addition, RBD residues Y449, Y489, Y495, Q498, and N487 were critical for the RBD-ACE2 interaction. Next, we found that the distribution of the mean stability changes and mean binding energy changes of RBD due to mutations at both serological and endosomal pH correlated well, indicating the similar effects of mutations. Overall, this computational analysis is useful for understanding the effects of missense mutations in SARS-CoV-2 pathogenesis at different pH.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Shafiul Haque
- Research and Scientific Studies Unit, College of Nursing and Allied Health Sciences, Jazan University, Jazan-45142, Saudi Arabia
- Gilbert and Rose-Marie Chagoury School of Medicine, Lebanese American University, Beirut, Lebanon
- Centre of Medical and Bio-Allied Health Sciences Research, Ajman University, Ajman, United Arab Emirates
| | - Darin Mansor Mathkor
- Research and Scientific Studies Unit, College of Nursing and Allied Health Sciences, Jazan University, Jazan-45142, Saudi Arabia
| | - Mustfa Faisal Alkhanani
- Biology Department, College of Sciences, University of Hafr Al Batin, Hafr Al Batin, Saudi Arabia
| | - Farkad Bantun
- Department of Microbiology, Faculty of Medicine, Umm Al-Qura University, Makkah, Saudi Arabia
| | - Aiman M Momenah
- Department of Microbiology, Faculty of Medicine, Umm Al-Qura University, Makkah, Saudi Arabia
| | - Hani Faidah
- Department of Microbiology, Faculty of Medicine, Umm Al-Qura University, Makkah, Saudi Arabia
| | - Naif A Jalal
- Department of Microbiology, Faculty of Medicine, Umm Al-Qura University, Makkah, Saudi Arabia
| | - Vijay Kumar
- Amity Institute of Neuropsychology & Neurosciences, Amity University, Noida, Uttar Pradesh, India
| |
Collapse
|
4
|
Mathkor DM, Faidah H, Jalal NA, Qashqari FSI, Haque S, Bantun F. In silico identification of microRNAs targeting the PPARα/γ: promising therapeutics for SARS-CoV‑2 infection. Biotechnol Genet Eng Rev 2023:1-12. [PMID: 36708330 DOI: 10.1080/02648725.2022.2163867] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/06/2022] [Accepted: 12/26/2022] [Indexed: 01/29/2023]
Abstract
The SARS-CoV-2 lifecycle is dependent on the host metabolism machinery. It upregulates the PPARα and PPARγ genes in lipid metabolism, which supports the essential viral replication complex including lipid rafts and palmitoylation of viral protein. The use of PPAR ligands in SARS-CoV-2 infection may have positive effects by preventing cytokine storm and the ensuing inflammatory cascade. The inhibition of PPARα and PPARγ genes may alter the metabolism and may disrupt the lifecycle of SARS-CoV-2 and COVID-19 progression. In the present work, we have identified possible miRNAs targeting PPARα and PPARγ in search of modulators of PPARα and PPARγ genes expression. The identified miRNAs could possibly be viewed as new therapeutic targets against COVID-19 infection.
Collapse
Affiliation(s)
- Darin Mansor Mathkor
- Research and Scientific Studies Unit, College of Nursing and Allied Health Sciences, Jazan University, Jazan, Saudi Arabia
| | - Hani Faidah
- Department of Microbiology, Faculty of Medicine, Umm Al-Qura University, Makkah, Saudi Arabia
| | - Naif A Jalal
- Department of Microbiology, Faculty of Medicine, Umm Al-Qura University, Makkah, Saudi Arabia
| | - Fadi S I Qashqari
- Department of Microbiology, Faculty of Medicine, Umm Al-Qura University, Makkah, Saudi Arabia
| | - Shafiul Haque
- Research and Scientific Studies Unit, College of Nursing and Allied Health Sciences, Jazan University, Jazan, Saudi Arabia
- Gilbert and Rose-Marie Chagoury School of Medicine, Lebanese American University, Beirut, Lebanon
- Centre of Medical and Bio-Allied Health Sciences Research, Ajman University, Ajman, United Arab Emirates
| | - Farkad Bantun
- Department of Microbiology, Faculty of Medicine, Umm Al-Qura University, Makkah, Saudi Arabia
| |
Collapse
|
5
|
Chakraborty C, Bhattacharya M, Sharma AR, Dhama K, Lee SS. The rapid emergence of multiple sublineages of Omicron (B.1.1.529) variant: Dynamic profiling via molecular phylogenetics and mutational landscape studies. J Infect Public Health 2022; 15:1234-1258. [PMID: 36270226 DOI: 10.1016/j.jiph.2022.10.004] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/29/2022] [Revised: 09/27/2022] [Accepted: 10/05/2022] [Indexed: 11/28/2022] Open
Abstract
PURPOSE The recent Omicron (B.1.1.529) variant poses a significant threat to global health. This variant has spread worldwide, and several sublineages have rapidly emerged. Study tried to analyze the microevolution of this variant. METHODS We studied the molecular phylogenetics, divergence, geographical distributions, frequencies, risk mutations for antibody affinity, and mutational landscape for Omicron sublineages using in silico analysis and statistical models. The risk mutation of spike for nAb affinity was analyzed and illustrated by statistical plots. Finally, the mutational properties of the spike mutations and their stability were predicted and demonstrated. RESULTS First, we studied the microevolutionary Omicron sublineages using molecular phylogenetics. Simultaneously, we revealed divergence events of the Omicron sublineages and observed the lowest minimum divergence of 51 in clade 21K and the highest maximum divergence of 90 in clade 21L. We have demonstrated cluster analyses, geographical distributions, frequencies of Omicron and its sublineages. Finally, we evaluated the mutational landscape of the Omicron sublineages. In this mutational study, we performed a genome-wide analysis of general mutations, mutations in the non-spike genome, and spike mutations of Omicron sublineages. The risk mutation of S-glycoprotein for nAb affinity has been analyzed for Omicron sublineages. Here, we found that some sublineages have all four significant highly destabilizing mutations. Such sublineages are BA.1 (G446S, E484A, T95I, and D614G), BA.2 (H655Y, Q493R, G493S, and D614G), BA.4 (N501Y, Y505H, N969K, and D614G), and BA.2.75 (Q454H, T547K, N764K, D614G and G446S). Finally, from the mutation stability prediction through ΔΔG, we observed that BA.1 and BA.4 had two destabilizing and two stabilizing mutations. Similarly, BA.2, BA.5, and BA.2.12.1 have one destabilizing and three stabilizing mutations. However, all four mutations in BA.2.75 are stabilizing mutations. CONCLUSIONS Our molecular phylogenetic studies provided a deeper understanding of the microevolution of sublineages and the creation of Omicron. Similarly, this study might help scientists develop pan-coronavirus vaccines that consider their mutational properties.
Collapse
Affiliation(s)
- Chiranjib Chakraborty
- Department of Biotechnology, School of Life Science and Biotechnology, Adamas University, Kolkata, West Bengal 700126, India.
| | - Manojit Bhattacharya
- Department of Zoology, Fakir Mohan University, Vyasa Vihar, Balasore 756020, Odisha, India
| | - Ashish Ranjan Sharma
- Institute for Skeletal Aging &Orthopedic Surgery, Hallym University-Chuncheon Sacred Heart Hospital, Chuncheon, Gangwon-Do 24252, Republic of Korea
| | - Kuldeep Dhama
- Division of Pathology, ICAR-Indian Veterinary Research Institute, Izatnagar, Bareilly 243122, Uttar Pradesh, India
| | - Sang-Soo Lee
- Institute for Skeletal Aging &Orthopedic Surgery, Hallym University-Chuncheon Sacred Heart Hospital, Chuncheon, Gangwon-Do 24252, Republic of Korea.
| |
Collapse
|
6
|
Isolation and In Silico Inhibitory Potential against SARS-CoV-2 RNA Polymerase of the Rare Kaempferol 3-O-(6″-O-acetyl)-Glucoside from Calligonum tetrapterum. PLANTS 2022; 11:plants11152072. [PMID: 35956550 PMCID: PMC9370365 DOI: 10.3390/plants11152072] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/17/2022] [Revised: 08/04/2022] [Accepted: 08/05/2022] [Indexed: 11/17/2022]
Abstract
The phytochemical constituents of Calligonum tetrapterum Jaub. & Spach (Family Polygonaceae) were studied for the first time. The study resulted in the isolation of the rare flavonol glycoside, kaempferol 3-O-(6″-O-acetyl)-glucoside,(K3G-A). The potential inhibitive activity of K3G-A toward SARS-CoV-2 was investigated utilizing several in silico approaches. First, molecular fingerprints and structural similarity experiments were carried out for K3G-A against nine co-crystallized ligands of nine proteins of SARS-CoV-2 to reveal if there is a structural similarity with any of them. The conducted studies showed the high similarity of K3G-A and remdesivir, the co-crystallized ligand of SARS-CoV-2 RNA-dependent RNA polymerase (PDB ID: 7BV2), RdRp. To validate these findings, a DFT study was conducted and confirmed the proposed similarity on the electronic and orbital levels. The binding of K3G-A against RdRp was confirmed through molecular docking studies exhibiting a binding energy of −27.43 kcal/mol, which was higher than that of remdesivir. Moreover, the RdRp-K3G-A complex was subjected to several MD studies at 100 ns that authenticated the accurate mode of binding and the correct dynamic behavior. Finally, in silico ADMET and toxicity evaluation of K3G-A was conducted and denoted the safety and the drug-likeness of K3G-A. In addition to K3G-A, two other metabolites were isolated and identified to be kaempferol (K) and β-sitosterol (β-S).
Collapse
|