Abstract
Heparan sulfate proteoglycans serve as initial attachment sites for several viruses and bacteria. Recent studies suggest that SARS-CoV-2 similarly exploits these glycosaminoglycans, facilitating conformational changes in the spike protein that promote the interaction between the receptor-binding domain (S1-RBD) and the cellular angiotensin-converting enzyme 2 receptor (ACE2), thereby triggering the virus internalization process. The molecular details that drive this process, particularly the co-receptor role of heparan sulfate (HS), remain incompletely understood. The interaction between an HS hexasaccharide (hexa) and the N343 glycosylated S1-RBD of the wild-type (WT) and Omicron variant of SARS-CoV-2 was investigated. The conformational properties of hexa with these S1-RBDs in unbound and bound states are explored using multiple independent MD simulations; the protein binding epitope of hexa, as well as the details of its interaction with S1-RBD of the Omicron variant, are characterized by comparing experimental and theoretical (1)H STD NMR signals. This investigation identifies the role played by the glycosyl moiety at N343 in potentially affecting this interaction in both WT and Omicron S1-RBD, explaining the observed low specificity and multi-modal nature of the interaction between HS oligosaccharides and these S1-RBDs.