Back

Network analysis outlines strengths and weaknesses of emerging SARS-CoV-2 Spike variants

Manrique, P. D.; Chakraborty, S.; Nguyen, K.; Mansbach, R.; Korber, B.; Gnanakaran, S.

2021-09-04 molecular biology
10.1101/2021.09.03.458946 bioRxiv
Show abstract

The COVID-19 pandemic, caused by the SARS-CoV-2 virus, has triggered myriad efforts to dissect and understand the structure and dynamics of this complex pathogen. The Spike glycoprotein of SARS-CoV-2 has received special attention as it is the means by which the virus enters the human host cells. The N-terminal domain (NTD) is one of the targeted regions of the Spike protein for therapeutics and neutralizing antibodies against COVID-19. Though its function is not well-understood, the NTD is reported to acquire mutations and deletions that can accelerate the evolutionary adaptation of the virus driving antibody escape. Cellular processes are known to be regulated by complex interactions at the molecular level, which can be characterized by means of a graph representation facilitating the identification of key residues and critical communication pathways within the molecular complex. From extensive all-atom molecular dynamics simulations of the entire Spike for the wild-type and the dominant variant, we derive a weighted graph representation of the protein in two dominant conformations of the receptor-binding-domain; all-down and one-up. We implement graph theory techniques to characterize the relevance of specific residues at facilitating roles of communication and control, while uncovering key implications for fitness and adaptation. We find that many of the reported high-frequency mutations tend to occur away from the critical residues highlighted by our graph theory analysis, implying that these mutations tend to avoid targeting residues that are most critical for protein allosteric communication. We propose that these critical residues could be candidate targets for novel antibody therapeutics. In addition, our analysis provides quantitative insights of the critical role of the NTD and furin cleavage site and their wide-reaching influence over the protein at large. Many of our conclusions are supported by empirical evidence while others point the way towards crucial simulation-guided experiments.

Matching journals

The top 4 journals account for 50% of the predicted probability mass.

1
Journal of Molecular Biology
232 papers in training set
Top 0.1%
26.6%
2
Biomolecules
100 papers in training set
Top 0.1%
9.8%
3
PLOS Computational Biology
1863 papers in training set
Top 4%
9.7%
4
Computational and Structural Biotechnology Journal
242 papers in training set
Top 0.5%
5.5%
50% of probability mass above
5
Journal of Chemical Information and Modeling
238 papers in training set
Top 1%
4.3%
6
Viruses
332 papers in training set
Top 1%
4.3%
7
Biophysical Journal
631 papers in training set
Top 2%
4.0%
8
PLOS ONE
5266 papers in training set
Top 43%
2.4%
9
eLife
5828 papers in training set
Top 41%
2.4%
10
Journal of The Royal Society Interface
235 papers in training set
Top 2%
2.4%
11
Scientific Reports
3612 papers in training set
Top 44%
2.4%
12
Nature Communications
5641 papers in training set
Top 40%
2.4%
13
Proceedings of the National Academy of Sciences
2444 papers in training set
Top 24%
2.1%
14
International Journal of Molecular Sciences
494 papers in training set
Top 6%
2.0%
15
Communications Biology
993 papers in training set
Top 16%
1.5%
16
Proteins: Structure, Function, and Bioinformatics
88 papers in training set
Top 0.9%
1.3%
17
Entropy
21 papers in training set
Top 0.1%
1.3%
18
Physical Chemistry Chemical Physics
36 papers in training set
Top 0.4%
1.1%
19
PNAS Nexus
159 papers in training set
Top 2%
1.1%
20
Nucleic Acids Research
1281 papers in training set
Top 13%
0.8%
21
BMC Biology
265 papers in training set
Top 5%
0.8%
22
iScience
1154 papers in training set
Top 39%
0.6%