Elucidating the Binding Mechanism of SARS-CoV-2 NSP6-TBK1 and Structure-based Designing of Phytocompounds Inhibitors for Instigating the Host Immune Response

Overview

Journal Front Chem

Specialty Chemistry

Date 2024 Jan 31

PMID 38293247

Authors

Muhammad Suleman

Iqra Ishaq

Haji Khan

Safir Ullah Khan

Rehana Masood

Norah A Albekairi

Abdulrahman Alshammari

Sergio Crovella

Affiliations

Soon will be listed here.

Abstract

SARS-CoV-2, also referred to as severe acute respiratory syndrome coronavirus 2, is the virus responsible for causing COVID-19, an infectious disease that emerged in Wuhan, China, in December 2019. Among its crucial functions, NSP6 plays a vital role in evading the human immune system by directly interacting with a receptor called TANK-binding kinase (TBK1), leading to the suppression of IFNβ production. Consequently, in the present study we used the structural and biophysical approaches to analyze the effect of newly emerged mutations on the binding of NSP6 and TBK1. Among the identified mutations, four (F35G, L37F, L125F, and I162T) were found to significantly destabilize the structure of NSP6. Furthermore, the molecular docking analysis highlighted that the mutant NSP6 displayed its highest binding affinity with TBK1, exhibiting docking scores of -1436.2 for the wildtype and -1723.2, -1788.6, -1510.2, and -1551.7 for the F35G, L37F, L125F, and I162T mutants, respectively. This suggests the potential for an enhanced immune system evasion capability of NSP6. Particularly, the F35G mutation exhibited the strongest binding affinity, supported by a calculated binding free energy of -172.19 kcal/mol. To disrupt the binding between NSP6 and TBK1, we conducted virtual drug screening to develop a novel inhibitor derived from natural products. From this screening, we identified the top 5 hit compounds as the most promising candidates with a docking score of -6.59 kcal/mol, -6.52 kcal/mol, -6.32 kcal/mol, -6.22 kcal/mol, and -6.21 kcal/mol. The molecular dynamic simulation of top 3 hits further verified the dynamic stability of drugs-NSP6 complexes. In conclusion, this study provides valuable insight into the higher infectivity of the SARS-CoV-2 new variants and a strong rationale for the development of novel drugs against NSP6.

References

Li G, Hilgenfeld R, Whitley R, De Clercq E . Therapeutic strategies for COVID-19: progress and lessons learned. Nat Rev Drug Discov. 2023; 22(6):449-475. PMC: 10113999. DOI: 10.1038/s41573-023-00672-y. View

Acharya D, Liu G, Gack M . Dysregulation of type I interferon responses in COVID-19. Nat Rev Immunol. 2020; 20(7):397-398. PMC: 7249038. DOI: 10.1038/s41577-020-0346-x. View

Jewkes R, Sikweyiya Y, Morrell R, Dunkle K . Gender inequitable masculinity and sexual entitlement in rape perpetration South Africa: findings of a cross-sectional study. PLoS One. 2012; 6(12):e29590. PMC: 3247272. DOI: 10.1371/journal.pone.0029590. View

Salomon-Ferrer R, Gotz A, Poole D, Le Grand S, Walker R . Routine Microsecond Molecular Dynamics Simulations with AMBER on GPUs. 2. Explicit Solvent Particle Mesh Ewald. J Chem Theory Comput. 2015; 9(9):3878-88. DOI: 10.1021/ct400314y. View

Roe D, Cheatham 3rd T . PTRAJ and CPPTRAJ: Software for Processing and Analysis of Molecular Dynamics Trajectory Data. J Chem Theory Comput. 2015; 9(7):3084-95. DOI: 10.1021/ct400341p. View

Sundar S, Thangamani L, Piramanayagam S, Rahul C, Aiswarya N, Sekar K . Screening of FDA-approved compound library identifies potential small-molecule inhibitors of SARS-CoV-2 non-structural proteins NSP1, NSP4, NSP6 and NSP13: molecular modeling and molecular dynamics studies. J Proteins Proteom. 2021; 12(3):161-175. PMC: 8188161. DOI: 10.1007/s42485-021-00067-w. View

Panahi Y, Mahdavi Gorabi A, Talaei S, Beiraghdar F, Akbarzadeh A, Tarhriz V . An overview on the treatments and prevention against COVID-19. Virol J. 2023; 20(1):23. PMC: 9906607. DOI: 10.1186/s12985-023-01973-9. View

Suleman M, Yousafi Q, Ali J, Ali S, Hussain Z, Ali S . 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. Comput Biol Med. 2021; 138:104936. PMC: 8501515. DOI: 10.1016/j.compbiomed.2021.104936. View

Wu Z, McGoogan J . Characteristics of and Important Lessons From the Coronavirus Disease 2019 (COVID-19) Outbreak in China: Summary of a Report of 72 314 Cases From the Chinese Center for Disease Control and Prevention. JAMA. 2020; 323(13):1239-1242. DOI: 10.1001/jama.2020.2648. View

10.

Gopalakrishnan K, Sowmiya G, Sheik S, Sekar K . Ramachandran plot on the web (2.0). Protein Pept Lett. 2007; 14(7):669-71. DOI: 10.2174/092986607781483912. View

11.

Kawai T, Akira S . Antiviral signaling through pattern recognition receptors. J Biochem. 2006; 141(2):137-45. DOI: 10.1093/jb/mvm032. View

12.

Webb B, Sali A . Comparative Protein Structure Modeling Using MODELLER. Curr Protoc Bioinformatics. 2016; 54:5.6.1-5.6.37. PMC: 5031415. DOI: 10.1002/cpbi.3. View

13.

Aziz A, Suleman M, Shah A, Ullah A, Rashid F, Khan S . Comparative mutational analysis of the Zika virus genome from different geographical locations and its effect on the efficacy of Zika virus-specific neutralizing antibodies. Front Microbiol. 2023; 14:1098323. PMC: 9992208. DOI: 10.3389/fmicb.2023.1098323. View

14.

Laha S, Chakraborty J, Das S, Manna S, Biswas S, Chatterjee R . Characterizations of SARS-CoV-2 mutational profile, spike protein stability and viral transmission. Infect Genet Evol. 2020; 85:104445. PMC: 7324922. DOI: 10.1016/j.meegid.2020.104445. View

15.

Zhu N, Zhang D, Wang W, Li X, Yang B, Song J . A Novel Coronavirus from Patients with Pneumonia in China, 2019. N Engl J Med. 2020; 382(8):727-733. PMC: 7092803. DOI: 10.1056/NEJMoa2001017. View

16.

Rodrigues C, Pires D, Ascher D . DynaMut2: Assessing changes in stability and flexibility upon single and multiple point missense mutations. Protein Sci. 2020; 30(1):60-69. PMC: 7737773. DOI: 10.1002/pro.3942. View

17.

Balachandran S, Roberts P, Brown L, Truong H, Pattnaik A, Archer D . Essential role for the dsRNA-dependent protein kinase PKR in innate immunity to viral infection. Immunity. 2000; 13(1):129-41. DOI: 10.1016/s1074-7613(00)00014-5. View

18.

Maier J, Martinez C, Kasavajhala K, Wickstrom L, Hauser K, Simmerling C . ff14SB: Improving the Accuracy of Protein Side Chain and Backbone Parameters from ff99SB. J Chem Theory Comput. 2015; 11(8):3696-713. PMC: 4821407. DOI: 10.1021/acs.jctc.5b00255. View

19.

Khan A, Hussain S, Ahmad S, Suleman M, Bukhari I, Khan T . Computational modelling of potentially emerging SARS-CoV-2 spike protein RBDs mutations with higher binding affinity towards ACE2: A structural modelling study. Comput Biol Med. 2022; 141:105163. PMC: 8716154. DOI: 10.1016/j.compbiomed.2021.105163. View

20.

Liu S, Cai X, Wu J, Cong Q, Chen X, Li T . Phosphorylation of innate immune adaptor proteins MAVS, STING, and TRIF induces IRF3 activation. Science. 2015; 347(6227):aaa2630. DOI: 10.1126/science.aaa2630. View