Spike Protein of SARS-CoV-2: Impact of Single Amino Acid Mutation and Effect of Drug Binding to the Variant-in Silico Analysis
Novel SARS-CoV-2, a bat based virus originated in Wuhan, China that caused a global pandemic in December, 2019 belongs to the Betacorona virus family and contains single stranded genome of ~29Kbp. The host cell invasion of SARS-CoV-2 is facilitated by interaction of C-Terminal Domain (CTD) of Spike (S) protein of virus and host ACE2 receptor in the presence of TMPRSS seine protease secreted by the host cell. In this study the mutation hotspots of S-protein will be identified and the impact of such mutation in the binding affinity will be studied. Additionally, the lead molecule which can bind to the mutated protein also will be identified. Multiple sequence alignment of the spike protein sequence of SARS-CoV-2 shows the number of single amino acid mutation hotspots such as L5F, R214L, R408I, G476S, V483A, H519Q, A520S, T572I, D614G and H655Y. Among these mutations D614G has 57.5% occurrence and G476S, V483A has 7.5% occurrence. The mutated proteins were modelled based on wild type homolog and docked to ACE2 receptor. When the mutated S protein is docked, the ∆G (binding free energy) value is very minimal in mutated protein showed the stability of variants. By the drug repurposing method, 1000 FDA approved drugs were virtually screened for its binding to RBD of S1 domain. Among these drugs Digitoxin, Gliquidone and Zorubicin Hcl binds to spike proteins with higher docking score (lesser than -8.5 Kcal/mol) to both wild type and mutants.