RESUMEN
Coronavirus disease-19 (COVID-19) caused by SARS-CoV-2 has already killed more than one million people worldwide. Since novel coronavirus is a new virus, mining its genome sequence is of crucial importance for drug/vaccine(s) development. Whole genome sequencing is a helpful tool in identifying genetic changes that occur in a virus when it spreads through the population. In this study, we performed complete genome sequencing of SARS-CoV-2 to unveil the genomic variation and indel, if present. We discovered thirteen (13) mutations in Orf1ab, S and N gene where seven (7) of them turned out to be novel mutations from our sequenced isolate. Besides, we found one (1) insertion and seven (7) deletions from the indel analysis among the 323 Bangladeshi isolates. However, the indel did not show any effect on proteins. Our energy minimization analysis showed both stabilizing and destabilizing impact on viral proteins depending on the mutation. Interestingly, all the variants were located in the binding site of the proteins. Furthermore, drug binding analysis revealed marked difference in interacting residues in mutants when compared to the wild type. Our analysis also suggested that eleven (11) mutations could exert damaging effects on their corresponding protein structures.
RESUMEN
The emerging variants of SARS coronavirus-2 (SARS-CoV-2) have been continuously spreading all over the world and have raised global health concerns. The B.1.1.7 (United Kingdom), P.1 (Brazil), B.1.351 (South Africa) and B.1.617 (India) variants, resulting from multiple mutations in the spike glycoprotein (SGp), are resistant to neutralizing antibodies and enable increased transmission. Hence, new drugs might be of great importance against the novel variants of SARS-CoV-2. The SGp and main protease (Mpro) of SARS-CoV-2 are important targets for designing and developing antiviral compounds for new drug discovery. In this study, we selected seventeen phytochemicals and later performed molecular docking to determine the binding interactions of the compounds with the two receptors and calculated several drug-likeliness properties for each compound. Luteolin, myricetin and quercetin demonstrated higher affinity for both the proteins and interacted efficiently. To obtain compounds with better properties, we designed three analogues from these compounds and showed their greater druggable properties compared to the parent compounds. Furthermore, we found that the analogues bind to the residues of both proteins, including the recently identified novel variants of SARS-CoV-2. The binding study was further verified by molecular dynamics (MD) simulation and molecular mechanics/Poisson Boltzmann surface area (MM/PBSA) approaches by assessing the stability of the complexes. MD simulations revealed that Arg457 of SGp and Met49 of Mpro are the most important residues that interacted with the designed inhibitors. Our analysis may provide some breakthroughs to develop new therapeutics to treat the proliferation of SARS-CoV-2 in vitro and in vivo.