From atomic to global connectivity in the structure of the SARS-CoV2-Human ACE2 receptor complex
Pith reviewed 2026-05-24 04:21 UTC · model grok-4.3
The pith
SARS-CoV2 forms a more robust connection to the ACE2 receptor than SARS-CoV1 through stronger interface networks.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
SARS-CoV2 forms a more dominant, robust connection with the ACE2-receptor as compared to the less virulent SARS-CoV1. This is reflected by percolation analysis where the interface cluster persists when restricted to stronger and stronger bonds, and by richer atomic-level clique structure.
What carries the argument
Protein side chain-based network method that averages clusters and cliques over MD simulation snapshots to map from atomic non-covalent interactions to global amino acid connectivity.
If this is right
- The interface cluster in SARS-CoV2 persists under stronger bond restrictions unlike in SARS-CoV1.
- Key functional residues in SARS-CoV2 are pinpointed for their role in higher connectivity.
- Many mutations in variants of concern occur at the ACE2 interface.
- The method provides an objective way to map spatial connectivity across scales.
Where Pith is reading between the lines
- If the connectivity difference holds, it could guide design of receptor-like peptides to block escape variants.
- Similar network analysis might apply to other virus-receptor pairs to predict virulence.
- Variations in MD snapshots highlight conformational landscapes that could be targeted.
Load-bearing premise
The side-chain contact definition and chosen MD snapshots from simulations are free of artifacts and their average represents the dominant binding mode.
What would settle it
A recalculation of the percolation threshold showing that the SARS-CoV2 interface cluster breaks at the same bond strength as the SARS-CoV1 cluster, or fewer cliques in SARS-CoV2.
Figures
read the original abstract
We investigate connectivity properties of the SARS-CoV2 spike protein-human ACE2-receptor complex employing a protein side chain-based network method that allows us to span a range from atomic to global protein scales. We analyze network topology in terms of clusters and cliques obtained from averaging over snapshots of MD simulations (from D.E. Shaw Research). We demonstrate that SARS-CoV2 forms a more dominant, robust connection with the ACE2-receptor as compared to the less virulent SARS-CoV1. Globally, this stronger connectivity is reflected by our percolation analysis where the interface cluster for the SARS-CoV2-ACE2 complex persists when restricted to stronger and stronger bonds, as compared to the SARSCoV1- ACE2 complex. At the atomic level, interface clique structure reflects a stronger connectivity in the former complex. We pinpoint key functional residues in SARSCoV2 that play important roles in establishing this higher connectivity. Thus, our studies provide an objective method to map spatial connectivity of atomic level non-covalent interactions to global connectivity between any two amino acids in the complex. We also analyze specific snapshots of the MD simulations to highlight prominent variations in network topology that explore diverse conformational landscapes. Finally, we demonstrate that a majority of mutations that occur in the SARSCoV2 spike protein in variants of concern/interest (including the currently circulating JN.1) have been observed at the interface with the ACE2 receptor. Our analyses highlight the importance of interface interactions and provide a rationale for designing receptor-like peptides and proteins to combat immunity-escaping variants.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript applies a side-chain contact network analysis to MD simulation snapshots of the SARS-CoV2–ACE2 and SARS-CoV1–ACE2 complexes. It reports that the SARS-CoV2 interface exhibits stronger global connectivity, demonstrated by a percolation threshold that survives stricter bond-strength cutoffs, together with richer atomic-scale clique structure; key functional residues are identified and the locations of mutations in variants of concern are mapped onto the interface.
Significance. If the percolation and clique results are robust to the contact definition and snapshot ensemble, the work supplies a concrete, multi-scale mapping from atomic non-covalent interactions to global interface connectivity that could inform both mechanistic understanding of virulence differences and the design of receptor-mimetic inhibitors. The approach is parameter-light once the contact criterion is fixed and directly uses publicly available trajectories, which are positive features.
major comments (3)
- [Methods] Methods (percolation analysis paragraph): the claim that the SARS-CoV2 interface cluster “persists when restricted to stronger and stronger bonds” is load-bearing for the central robustness conclusion, yet no numerical bond-strength cutoff values, occupancy thresholds, or energy criteria are stated, nor is any sensitivity test to modest changes in these criteria reported.
- [Methods] Methods (snapshot selection): the percolation and clique results rest on averaging over an unspecified subset of D.E. Shaw trajectories; no justification is given for the choice of frames, no error bars or bootstrap estimates on the percolation threshold are supplied, and no comparison against the corresponding crystal structures is presented to rule out force-field or sampling artifacts.
- [Results] Results (clique enumeration): the atomic-level claim of “richer clique structure” requires an explicit description of how cliques are enumerated (e.g., Bron–Kerbosch parameters, minimum size, weighting) and whether the reported difference survives variation in those parameters; this information is absent.
minor comments (2)
- [Abstract] Abstract: inconsistent hyphenation (“SARSCoV1” vs. “SARS-CoV1”) and missing space before “ACE2” in one instance should be standardized.
- [Figures] Figure captions: the percolation plots should explicitly label the bond-strength axis with the numerical cutoff values used in each panel.
Simulated Author's Rebuttal
We thank the referee for the constructive comments that highlight areas where methodological details can be strengthened for clarity and reproducibility. We address each point below and will incorporate the requested information in a revised manuscript.
read point-by-point responses
-
Referee: [Methods] Methods (percolation analysis paragraph): the claim that the SARS-CoV2 interface cluster “persists when restricted to stronger and stronger bonds” is load-bearing for the central robustness conclusion, yet no numerical bond-strength cutoff values, occupancy thresholds, or energy criteria are stated, nor is any sensitivity test to modest changes in these criteria reported.
Authors: We agree that explicit numerical values and sensitivity tests are essential. In the revised manuscript we will state the precise bond-strength cutoffs, occupancy thresholds, and any energy criteria used to define contacts in the percolation analysis. We will also add a sensitivity analysis demonstrating how the percolation threshold for the SARS-CoV2 interface responds to modest changes in these parameters, thereby supporting the robustness claim with quantitative evidence. revision: yes
-
Referee: [Methods] Methods (snapshot selection): the percolation and clique results rest on averaging over an unspecified subset of D.E. Shaw trajectories; no justification is given for the choice of frames, no error bars or bootstrap estimates on the percolation threshold are supplied, and no comparison against the corresponding crystal structures is presented to rule out force-field or sampling artifacts.
Authors: We acknowledge the need for greater transparency on trajectory handling. The revised Methods section will specify the exact frames selected from the D.E. Shaw trajectories, the rationale for that selection (e.g., equilibration and convergence checks), and will report error bars or bootstrap estimates on the percolation threshold. We will also include a direct comparison of key network metrics between the MD snapshots and the corresponding crystal structures to address potential force-field or sampling concerns. revision: yes
-
Referee: [Results] Results (clique enumeration): the atomic-level claim of “richer clique structure” requires an explicit description of how cliques are enumerated (e.g., Bron–Kerbosch parameters, minimum size, weighting) and whether the reported difference survives variation in those parameters; this information is absent.
Authors: We agree that the clique-analysis procedure must be fully specified. In revision we will describe the enumeration algorithm (including Bron–Kerbosch parameters if used), the minimum clique size, any weighting scheme, and will demonstrate that the observed difference in clique richness between the SARS-CoV2–ACE2 and SARS-CoV1–ACE2 interfaces remains robust under reasonable variations of these parameters. These details will be added to the Results section. revision: yes
Circularity Check
No significant circularity; standard network analysis on external MD data
full rationale
The paper constructs side-chain contact networks from external D.E. Shaw MD trajectories, then applies standard graph operations (clustering, clique detection, percolation under bond-strength cutoffs) to compare SARS-CoV2 vs SARS-CoV1 interfaces. No equations, fitted parameters, or self-definitional reductions appear; the central robustness claim is a direct output of these operations on the input contact graphs rather than a quantity forced by construction or by a self-citation chain. The method itself is cited as prior work but functions as an independent, externally applicable tool whose results on this dataset are falsifiable against other contact definitions or trajectories.
Axiom & Free-Parameter Ledger
Reference graph
Works this paper leans on
-
[1]
This feature is consistent with the largest cluster transition profiles as function ofImin at 50% dynamic stability as seen in Figure 1(a). The stability of the largest interface cluster at higher interaction strengths in SARS-CoV2 complex as compared to SARS-CoV1 complex indicates that information regarding higher binding affinity can be extracted from t...
work page 2019
-
[2]
D Mannar, et al., Sars-cov-2 omicron variant: Antibody evasion and cryo-em structure of spike protein–ace2 complex. Science 375, 760–764 (2022)
work page 2022
-
[3]
PH Guzzi, et al., Computational analysis of the sequence-structure relation in sars-cov-2 spike protein using protein contact networks. Sci. Reports 13, 2837 (2023)
work page 2023
-
[4]
A Zabiegala, Y Kim, KO Chang, Roles of host proteases in the entry of sars-cov-2. Animal Dis. 3, 12 (2023)
work page 2023
-
[5]
DA Jamison, et al., A comprehensive sars-cov-2 and covid-19 review, part 1: Intracellular overdrive for sars-cov-2 infection. Eur. J. Hum. Genet . 30, 889–898 (2022)
work page 2022
-
[6]
SA Narayanan, et al., A comprehensive sars-cov-2 and covid-19 review, part 2: host extracel- lular to systemic effects of sars-cov-2 infection. Eur. J. Hum. Genet . 32, 10–20 (2024)
work page 2024
-
[7]
N Maroli, Riding the wave: Unveiling the conformational waves from rbd of sars-cov-2 spike protein to ace2. The J. Phys. Chem. B 127, 8525–8536 (2023)
work page 2023
-
[8]
BT Kaynak, et al., Sampling of protein conformational space using hybrid simulations: A critical assessment of recent methods. Front. Mol. Biosci. 9 (2022)
work page 2022
-
[9]
B Wingert, J Krieger, H Li, I Bahar, Adaptability and specificity: how do proteins balance opposing needs to achieve function? Curr Opin Struct Biol 67, 25–32 (2020)
work page 2020
-
[10]
D Scaramozzino, RL Jernigan, Special issue on computational approaches for protein dynam- ics and function. Appl. Sci. 13 (2023)
work page 2023
-
[11]
B Jawad, P Adhikari, R Podgornik, WY Ching, Key interacting residues between rbd of sars- cov-2 and ace2 receptor: Combination of molecular dynamics simulation and density functional calculation. J. Chem. Inf. Model . 61, 4425–4441 (2021)
work page 2021
-
[12]
AS Pillai, GK Hochberg, JW Thornton, Simple mechanisms for the evolution of protein complexity. Protein Sci. 31, e4449 (2022)
work page 2022
-
[13]
A Ali, R Vijayan, Dynamics of the ace2–sars-cov-2/sars-cov spike protein interface reveal unique mechanisms. Sci. Reports 10, 14214 (2020)
work page 2020
-
[14]
C Bai, A Warshel, Critical differences between the binding features of the spike proteins of sars-cov-2 and sars-cov. The J. Phys. Chem. B 124, 5907–5912 (2020)
work page 2020
-
[15]
Science 369, 1586– 1592 (2020)
Y Cai, et al., Distinct conformational states of sars-cov-2 spike protein. Science 369, 1586– 1592 (2020)
work page 2020
-
[16]
C Xu, et al., Conformational dynamics of sars-cov-2 trimeric spike glycoprotein in complex with receptor ace2 revealed by cryo-em. Sci. Adv. 7, eabe5575 (2021)
work page 2021
-
[17]
M Noval Rivas, RA Porritt, MH Cheng, I Bahar, M Arditi, Multisystem inflammatory syndrome in children and long covid: The sars-cov-2 viral superantigen hypothesis. Front. Immunol. 13 (2022)
work page 2022
-
[18]
G Verkhivker, S Agajanian, R Kassab, K Krishnan, Computer simulations and network-based profiling of binding and allosteric interactions of sars-cov-2 spike variant complexes and the host receptor: Dissecting the mechanistic effects of the delta and omicron mutations. Int. J. Mol. Sci. 23 (2022). 7
work page 2022
-
[19]
R Nussinov, M Zhang, Y Liu, H Jang, Alphafold, artificial intelligence (ai), and allostery. The J. Phys. Chem. B 126, 6372–6383 (2022)
work page 2022
-
[20]
JPA Ioannidis, F Zonta, M Levitt, Variability in excess deaths across countries with different vulnerability during 2020–2023. Proc. Natl. Acad. Sci . 120, e2309557120 (2023)
work page 2020
-
[21]
DES Research, Molecular dynamics simulations related to sars-cov-2 (2020) D. E. Shaw Research Technical Data, https://www.deshawresearch.com/downloads/download_trajectory_ sarscov2.cgi/
work page 2020
-
[22]
A Halder, et al., Surveying the side-chain network approach to protein structure and dynamics: The sars-cov-2 spike protein as an illustrative case. Front. Mol. Biosci. 7 (2020)
work page 2020
-
[23]
V Gadiyaram, S Vishveshwara, S Vishveshwara, From quantum chemistry to networks in biology: A graph spectral approach to protein structure analyses. J. Chem. Inf. Model . 59, 1715–1727 (2019) PMID: 30912941
work page 2019
-
[24]
Curr Protein Pept Sci 17, 4–25 (2016)
M Bhattacharyya, S Ghosh, S Vishveshwara, Protein structure and function: Looking through the network of Side-Chain interactions. Curr Protein Pept Sci 17, 4–25 (2016)
work page 2016
-
[25]
M Bhattacharyya, CR Bhat, S Vishveshwara, An automated approach to network features of protein structure ensembles. Protein Sci. 22, 1399–1416 (2013)
work page 2013
-
[26]
N Kannan, S Vishveshwara, Identification of side-chain clusters in protein structures by a graph spectral method11edited by j. m. thornton. J. Mol. Biol . 292, 441–464 (1999)
work page 1999
-
[27]
K Guruprasad, Mutations in human sars-cov-2 spike proteins, potential drug binding and epitope sites for covid-19 therapeutics development. Curr. Res. Struct. Biol . 4, 41–50 (2022)
work page 2022
-
[28]
B Cosar, et al., Sars-cov-2 mutations and their viral variants. Cytokine & Growth Factor Rev. 63, 10–22 (2022)
work page 2022
-
[29]
C Roemer, et al., Sars-cov-2 evolution, post-omicron (2022) https://virological.org/t/ sars-cov-2-evolution-post-omicron/911
work page 2022
-
[30]
Y Cao, et al., Imprinted sars-cov-2 humoral immunity induces convergent omicron rbd evolution. Nature 614, 521–529 (2023)
work page 2023
-
[31]
JD Bloom, RA Neher, Fitness effects of mutations to sars-cov-2 proteins. Virus Evol . 9, vead055 (2023)
work page 2023
-
[32]
Y Kaku, et al., Virological characteristics of the sars-cov-2 jn.1 variant. bioRxiv (2023)
work page 2023
-
[33]
MC Chan, KK Chan, E Procko, D Shukla, Machine learning guided design of high-affinity ace2 decoys for sars-cov-2 neutralization. The J. Phys. Chem. B 127, 1995–2001 (2023)
work page 1995
-
[34]
Y Huang, C Y ang, Xf Xu, W Xu, Sw Liu, Structural and functional properties of sars-cov-2 spike protein: potential antivirus drug development for covid-19. Acta Pharmacol. Sinica 41, 1141–1149 (2020)
work page 2020
-
[35]
C Chen, et al., Computational prediction of the effect of amino acid changes on the binding affin- ity between sars-cov-2 spike rbd and human ace2. Proc. Natl. Acad. Sci . 118, e2106480118 (2021)
work page 2021
-
[36]
T rends Pharmacol Sci 43, 838–851 (2022)
T Arimori, et al., Engineering ace2 decoy receptors to combat viral escapability. T rends Pharmacol Sci 43, 838–851 (2022)
work page 2022
-
[37]
J Lan, et al., Structure of the sars-cov-2 spike receptor-binding domain bound to the ace2 receptor. Nature 581, 215–220 (2020)
work page 2020
-
[38]
protein folding: Differing questions, different challenges
SJ Chen, et al., Protein folds vs. protein folding: Differing questions, different challenges. Proc. Natl. Acad. Sci . 120, e2214423119 (2023)
work page 2023
-
[39]
J Jumper, et al., Highly accurate protein structure prediction with alphafold. Nature 596, 583–589 (2021)
work page 2021
-
[40]
R Sathyapriya, MS Vijayabaskar, S Vishveshwara, Insights into protein–dna interactions through structure network analysis. PLOS Comput. Biol . 4, 1–15 (2008)
work page 2008
-
[41]
K Brinda, S Vishveshwara, S Vishveshwara, Random network behaviour of protein structures. Mol. BioSystems 6, 391–398 (2010) Copyright: Copyright 2010 Elsevier B.V., All rights reserved
work page 2010
-
[42]
?a" for residues that contribute to aromaticity,
D Deb, S Vishveshwara, S Vishveshwara, Understanding protein structure from a percolation perspective. Biophys. journal 97, 1787–1794 (2009). Fig. 1. Mesoscopic view of global connectivity maps of spike protein – ACE2 receptor complexes reveal stronger percolation behavior in SARS-CoV2 compared to SARS-CoV1 : (a-c) Percolation analyses of the largest clus...
work page 2009
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.