Structural insights into the effect of mutations in the spike protein of SARS-CoV-2 on the binding with human furin protein

Abstract

The SARS COV-2 and its variants are spreading around the world at an alarming speed, due to its higher transmissibility and the conformational changes caused by mutations. The resulting COVID-19 pandemic has imposed severe health consequences on human health. Several countries of the world including Pakistan have studied its genome extensively and provided productive findings. In the current study, the mCSM, DynaMut2, and I-Mutant servers were used to analyze the effect of identified mutations on the structural stability of spike protein however, the molecular docking and simulations approaches were used to evaluate the dynamics of the bonding network between the wild-type and mutant spike proteins with furin. We addressed the mutational modifications that have occurred in the spike protein of SARS-COV-2 that were found in 215 Pakistani's isolates of COVID-19 patients to study the influence of mutations on the stability of the protein and its interaction with the host cell. We found 7 single amino acid substitute mutations in various domains that reside in spike protein. The H49Y, N74K, G181V, and G446V were found in the S1 domain while the D614A, V622F, and Q677H mutations were found in the central helices of the spike protein. Based on the observation, G181V, G446V, D614A, and V622F mutants were found highly destabilizing and responsible for structural perturbation. Protein-protein docking and molecular simulation analysis with that of furin have predicted that all the mutants enhanced the binding efficiency however, the V622F mutant has greatly altered the binding capacity which is further verified by the KD value (7.1 E−14) and therefore may enhance the spike protein cleavage by Furin and increase the rate of infectivity by SARS-CoV-2. On the other hand, the total binding energy for each complex was calculated which revealed −50.57 kcal/mol for the wild type, for G181V −52.69 kcal/mol, for G446V −56.44 kcal/mol, for D614A −59.78 kcal/mol while for V622F the TBE was calculated to be −85.84 kcal/mol. Overall, the current finding shows that these mutations have increased the binding of Furin for spike protein and shows that D614A and V622F have significant effects on the binding and infectivity.