Mutations in the SARS-CoV-2 spike protein modulate the virus affinity to the human ACE2 receptor, an in silico analysis
Open Access
- 8 March 2021
- journal article
- research article
- Published by Forschungsgesellschaft für Arbeitsphysiologie und Arbeitsschutz e.V. in EXCLI Journal
- Vol. 20, 585-600
Abstract
The increasing number of SARS-CoV-2 variants associated with highly transmissible phenotypes is a health-public concern in the current pandemic scenario. Herein, we developed a comprehensive in silico analysis of the changes occurring upon mutations in the viral spike. We focused on mutants located in the receptor-binding domain of the viral spike protein and analyzed whether these mutants modulate the interaction with the human host receptor angiotensin-converting enzyme II (ACE2). Thirty-two highly prevalent mutants were retrieved from the GISAID database, and their structural models were built using the SWISS-Model server. The stabilization effect for each mutation was assessed by the DUET and DeepDGG software. By applying molecular docking using both Z-Dock and Haddock software we found that multiple mutations, including A475V, V455E, V445L, and V445I, resulted in the higher binding free energy as compared to the wild type (WT) spike protein, thus had a destabilizing effect on the binding to ACE2. On the other hand, several mutants, including the most prevalent N501Y and B.1.1.7 variants, as well as the K444R, L455F, Q493R, and Y505W variants exhibited lower binding free energy as compared to the WT spike. These mutants showed an increased number of electrostatic interactions with ACE2 than the WT spike protein, and they changed the interaction pattern of the neighboring residues. Together, the results presented in this study contribute to a better understanding of the changes in the interaction between SARS-CoV-2 and the human host ACE2 receptor associated with point mutations in the viral spike protein.Keywords
This publication has 30 references indexed in Scilit:
- DUET: a server for predicting effects of mutations on protein stability using an integrated computational approachNucleic Acids Research, 2014
- ZDOCK server: interactive docking prediction of protein–protein complexes and symmetric multimersBioinformatics, 2014
- Atomic-Level Protein Structure Refinement Using Fragment-Guided Molecular Dynamics Conformation SamplingStructure, 2011
- CHARMM general force field: A force field for drug‐like molecules compatible with the CHARMM all‐atom additive biological force fieldsJournal of Computational Chemistry, 2009
- ProSA-web: interactive web service for the recognition of errors in three-dimensional structures of proteinsNucleic Acids Research, 2007
- Monoclonal Antibodies Targeting the HR2 Domain and the Region Immediately Upstream of the HR2 of the S Protein Neutralize In Vitro Infection of Severe Acute Respiratory Syndrome CoronavirusJournal of Virology, 2006
- The SWISS-MODEL workspace: a web-based environment for protein structure homology modellingBioinformatics, 2005
- Scalable molecular dynamics with NAMDJournal of Computational Chemistry, 2005
- A structural and thermodynamic escape mechanism from a drug resistant mutation of the HIV‐1 proteaseProteins: Structure, Function, and Bioinformatics, 2004
- Antibody neutralization and escape by HIV-1Nature, 2003