Department of Bioinformatics and Biological Statistics, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai, 200240, PR China.
Rashid Latif Medical College, Lahore, Punjab, Pakistan.
Comput Biol Med. 2022 Jul;146:105574. doi: 10.1016/j.compbiomed.2022.105574. Epub 2022 Apr 30.
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.
随着德尔塔和奥密克戎变异株的出现,包括 A.30 在内的导致 2019 冠状病毒病的 SARS-CoV-2 的许多其他重要变异株,被认为增加了全球大流行带来的担忧。在坦桑尼亚和其他国家报告的 A.30 变异株,其刺突基因发生突变,有助于该毒株更牢固地结合并逃避中和抗体。本研究使用分子建模和基于模拟的方法来研究导致这种毒株具有更高传染性的关键特征。刺突蛋白的蛋白-蛋白对接结果表明,额外的相互作用,特别是由突变残基 Lys484 形成的两个盐桥,增加了结合亲和力,而 N 端结构域 (NTD) 关键残基的缺失导致与单克隆抗体结合构象的改变,从而逃避其中和作用。此外,我们通过分子模拟深入研究了这些结合复合物的原子特征,发现与野生型相比,它们具有不同的动力学。使用 MM/GBSA 分析结合自由能表明,野生型受体结合结构域 (RBD) 复合物的总结合自由能 (TBE) 为-58.25 kcal/mol,而 A.30 RBD 复合物报告的 TBE 为-65.59 kcal/mol。A.30 RBD 复合物更高的 TBE 表明 A.30 变异株 RBD 与 ACE2 之间的相互作用更加强烈,使该变异株能够更快地结合和传播。野生型 NTD 复合物的 BFE 计算为-65.76 kcal/mol,而 A.30 NTD 复合物估计为-49.35 kcal/mol。这表明 A.30 变异株 NTD 中报道的取代和缺失的影响,从而降低了 mAb 的结合,使其能够逃避宿主的免疫反应。报告的结果将有助于开发针对 SARS-CoV-2 及其变异株的交叉保护药物。