About

The interaction between the receptor-binding domain (RBD) of the SARS-CoV-2 spike glycoprotein and the ACE2 enzyme is believed to be the entry point of the virus into various cells in the body, including the lungs, heart, liver, and kidneys.

The current focus of several therapeutic design efforts explores attempts at affecting the binding potential between the two proteins to limit the activity of the virus and disease progression. In this work, we analyze the stability of the spike protein under all possible single-point mutations in the RBD and computationally explore mutations that can affect the binding with the ACE2 enzyme.

We unravel the mutation landscape of the receptor region and assess the toxicity potential of single and multi-point mutations, generating insights for future vaccine efforts on mutations that might further stabilize the spike protein and increase its infectivity.

We developed a tool, called SpikeMutator, to construct full atomic protein structures of the mutant spike proteins and shared a database of 3800 single-point mutant structures. We analyzed the recent 65,000 reported spike sequences across the globe and observed the emergence of stable multi-point mutant structures.

Using the landscape, we searched through 7.5 million possible 2-point mutation combinations and report that the (R355D K424E) mutation produces one of the strongest spike proteins that therapeutic efforts should investigate for the sake of developing effective vaccines.


Citation

When citing this tool in any scientific publication please refer to it as:

Smaoui, M.R., Yahyaoui, H. Unraveling the stability landscape of mutations in the SARS-CoV-2 receptor-binding domain. Sci Rep 11, 9166 (2021). https://doi.org/10.1038/s41598-021-88696-5

https://doi.org/10.1038/s41598-021-88696-5

Project information

Software In Action

Take a closer look at what we did with SpikeMutator! You can download all the COVID-19 mutant structures we generated and analyze them further from our repo.

Spike Protomer

3D Structure of a Spike Protomer

Spike Trimer

3D Structure of a Spike Trimer

Mutation Landscape

SARS-CoV-2 receptor-binding domain landscape

2D Heatmaps

2D Projection of the RBD Mutation Landscape

RMSD Plot

Full 50ns MD Simulation

RMSF Plot

Full 50ns MD Simulation

Non-mutant analysis

Superimposing non-mutant structure at time 50ns on starting confirmation

Mutant analysis

Superimposing mutant structure at time 50ns on starting confirmation

Contact Map

Residue contact map of the native and top mutant (R355D K424E) structures