RESUMO
Spike glycoprotein has a significant role in the entry of SARS-CoV-2 to host cells, which makes it a potential drug target. Continued accumulation of non-synonymous mutations in the receptor binding domain of spike protein poses great challenges in identifying antiviral drugs targeting this protein. This study aims to identify potential entry inhibitors of SARS-CoV-2 using virtual screening and molecular dynamics (MD) simulations from three distinct chemical libraries including Pandemic Response Box, Drugbank and DrugCentral, comprising 6971 small molecules. The molecules were screened against a binding pocket identified in the receptor-binding domain (RBD) region of the spike protein which is known as the linoleic acid binding pocket, a highly conserved motif among several SARS-CoV-2 variants. Through virtual screening and binding free energy calculations, we identified four top-scoring compounds, MMV1579787 ([2-Oxo-2-[2-(3-phenoxyphenyl)ethylamino]ethyl]phosphonic acid), Tretinoin, MMV1633963 ((2E,4E)-5-[3-(3,5-dichlorophenoxy)phenyl]penta-2,4-dienoic acid) and Polydatin, which were previously reported to have antibacterial, antifungal or antiviral properties. These molecules showed stable binding on MD simulations over 100 ns and maintained stable interactions with TYR365, PHE338, PHE342, PHE377, TYR369, PHE374 and LEU368 of the spike protein RBD that are found to be conserved among SARS-CoV-2 variants. Our findings were further validated with free energy landscape, principal component analysis and dynamic cross-correlation analysis. Our in silico analysis of binding mode and MD simulation analyses suggest that the identified compounds may impede viral entrance by interacting with the linoleic acid binding site of the spike protein of SARS-CoV-2 regardless of its variants, and they thus demand for further in vitro and in vivo research.Communicated by Ramaswamy H. Sarma.
RESUMO
Since its advent in December 2019, SARS-CoV-2 has diverged into multiple variants with differing levels of virulence owing to the accumulation of mutations in its genome. The structural changes induced by non-synonymous mutations in major drug targets of the virus are known to alter the binding of potential antagonistic inhibitors. Here, we analyzed the effects of non-synonymous mutations in major targets of SARS-CoV-2 in response to potential peptide inhibitors. We screened 12 peptides reported to have anti-viral properties against RBD and 5 peptides against Mpro of SARS-CoV-2 variants using molecular docking and simulation approaches. The mutational landscape of RBD among SARS-CoV-2 variants had 21 non-synonymous mutations across 18 distinct sites. Among these, 14 mutations were present in the RBM region directly interacting with the hACE2 receptor. However, Only 3 non-synonymous mutations were observed in Mpro. We found that LCB1 - a de novo-synthesized peptide has the highest binding affinity to RBD despite non-synonymous mutations in variants and engages key residues of RBD-hACE2 interaction such as K417, E484, N487, and N501. Similarly, an antimicrobial peptide; 2JOS, was identified against Mpro with high binding affinity as it interacts with key residues in dimerization sites such as E166 and F140 crucial for viral replication. MD simulations affirm the stability of RBD-LCB1 and Mpro-2JOS complexes with an average RMSD of 1.902 and 2.476 respectively. We ascertain that LCB1 and 2JOS peptides are promising inhibitors to combat emerging variants of SARS-CoV-2 and thus warrant further investigations using in-vitro and in-vivo analysis.Communicated by Ramaswamy H. Sarma.