Molecular characterization of interactions between the D614G variant of SARS-CoV-2 S-protein and neutralizing antibodies: A computational approach

The D614G variant of SARS-CoV-2 S-protein emerged in early 2020 and quickly became the dominant circulating strain in Europe and its environs. The variant was characterized by the higher viral load, which is not associated with disease severity, higher incorporation into the virion, and high cell entry via ACE-2 and TMPRSS2. Previous strains of the coronavirus and the current SARS-CoV-2 have demonstrated the selection of mutations as a mechanism of escaping immune responses. In this study, we used molecular dynamics simulation and MM-PBSA binding energy analysis to provide insights into the behaviour of the D614G S-protein at the molecular level and describe the neutralization mechanism of this variant. Our results show that the D614G S-protein adopts distinct conformational dynamics which is skewed towards the open-state conformation more than the closed-state conformation of the wild-type S-protein. Residue-specific variation of amino acid flexibility and domain-specific RMSD suggest that the mutation causes an allosteric conformational change in the RBD. Evaluation of the interaction energies between the S-protein and neutralizing antibodies show that the mutation may enhance, reduce or not affect the neutralizing interactions depending on the neutralizing antibody, especially if it targets the RBD. The results of this study have shed insights into the behaviour of the D614G S-protein at the molecular level and provided a glimpse of the neutralization mechanism of this variant.