tetano
Editor, Senior Moderator
Comput Biol Med
. 2022 Apr 30;146:105574.
doi: 10.1016/j.compbiomed.2022.105574. Online ahead of print.
Investigation of the binding and dynamic features of A.30 variant revealed higher binding of RBD for hACE2 and escapes the neutralizing antibody: A molecular simulation approach
Athar Shafiq[SUP] 1 [/SUP], Farrukh Zubair[SUP] 2 [/SUP], Amna Ambreen[SUP] 3 [/SUP], Muhammad Suleman[SUP] 4 [/SUP], Qudsia Yousafi[SUP] 5 [/SUP], Zahid Rasul Niazi[SUP] 6 [/SUP], Zeeshan Anwar[SUP] 7 [/SUP], Abbas Khan[SUP] 1 [/SUP], Anwar Mohammad[SUP] 8 [/SUP], Dong-Qing Wei[SUP] 9 [/SUP]
Affiliations
Abstract
With the emergence of Delta and Omicron variants, many other important variants of SARS-CoV-2, which cause Coronavirus disease-2019, including A.30, are reported to increase the concern created by the global pandemic. The A.30 variant, reported in Tanzania and other countries, harbors spike gene mutations that help this strain to bind more robustly and to escape neutralizing antibodies. The present study uses molecular modelling and simulation-based approaches to investigate the key features of this strain that result in greater infectivity. The protein-protein docking results for the spike protein demonstrated that additional interactions, particularly two salt-bridges formed by the mutated residue Lys484, increase binding affinity, while the loss of key residues at the N terminal domain (NTD) result in a change to binding conformation with monoclonal antibodies, thus escaping their neutralizing effects. Moreover, we deeply studied the atomic features of these binding complexes through molecular simulation, which revealed differential dynamics when compared to wild type. Analysis of the binding free energy using MM/GBSA revealed that the total binding free energy (TBE) for the wild type receptor-binding domain (RBD) complex was -58.25 kcal/mol in contrast to the A.30 RBD complex, which reported -65.59 kcal/mol. The higher TBE for the A.30 RBD complex signifies a more robust interaction between A.30 variant RBD with ACE2 than the wild type, allowing the variant to bind and spread more promptly. The BFE for the wild type NTD complex was calculated to be -65.76 kcal/mol, while the A.30 NTD complex was estimated to be -49.35 kcal/mol. This shows the impact of the reported substitutions and deletions in the NTD of A.30 variant, which consequently reduce the binding of mAb, allowing it to evade the immune response of the host. The reported results will aid the development of cross-protective drugs against SARS-CoV-2 and its variants.
Keywords: A.30 variant; Dissociation constant; Free energy; Protein-protein docking; Simulation.
. 2022 Apr 30;146:105574.
doi: 10.1016/j.compbiomed.2022.105574. Online ahead of print.
Investigation of the binding and dynamic features of A.30 variant revealed higher binding of RBD for hACE2 and escapes the neutralizing antibody: A molecular simulation approach
Athar Shafiq[SUP] 1 [/SUP], Farrukh Zubair[SUP] 2 [/SUP], Amna Ambreen[SUP] 3 [/SUP], Muhammad Suleman[SUP] 4 [/SUP], Qudsia Yousafi[SUP] 5 [/SUP], Zahid Rasul Niazi[SUP] 6 [/SUP], Zeeshan Anwar[SUP] 7 [/SUP], Abbas Khan[SUP] 1 [/SUP], Anwar Mohammad[SUP] 8 [/SUP], Dong-Qing Wei[SUP] 9 [/SUP]
Affiliations
- PMID: 35533461
- DOI: 10.1016/j.compbiomed.2022.105574
Abstract
With the emergence of Delta and Omicron variants, many other important variants of SARS-CoV-2, which cause Coronavirus disease-2019, including A.30, are reported to increase the concern created by the global pandemic. The A.30 variant, reported in Tanzania and other countries, harbors spike gene mutations that help this strain to bind more robustly and to escape neutralizing antibodies. The present study uses molecular modelling and simulation-based approaches to investigate the key features of this strain that result in greater infectivity. The protein-protein docking results for the spike protein demonstrated that additional interactions, particularly two salt-bridges formed by the mutated residue Lys484, increase binding affinity, while the loss of key residues at the N terminal domain (NTD) result in a change to binding conformation with monoclonal antibodies, thus escaping their neutralizing effects. Moreover, we deeply studied the atomic features of these binding complexes through molecular simulation, which revealed differential dynamics when compared to wild type. Analysis of the binding free energy using MM/GBSA revealed that the total binding free energy (TBE) for the wild type receptor-binding domain (RBD) complex was -58.25 kcal/mol in contrast to the A.30 RBD complex, which reported -65.59 kcal/mol. The higher TBE for the A.30 RBD complex signifies a more robust interaction between A.30 variant RBD with ACE2 than the wild type, allowing the variant to bind and spread more promptly. The BFE for the wild type NTD complex was calculated to be -65.76 kcal/mol, while the A.30 NTD complex was estimated to be -49.35 kcal/mol. This shows the impact of the reported substitutions and deletions in the NTD of A.30 variant, which consequently reduce the binding of mAb, allowing it to evade the immune response of the host. The reported results will aid the development of cross-protective drugs against SARS-CoV-2 and its variants.
Keywords: A.30 variant; Dissociation constant; Free energy; Protein-protein docking; Simulation.