Abstract
Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike glycoprotein is the prime target for vaccines, diagnostics, and therapeutic antibodies against the virus. While anchored in the viral envelope, for effective virulence, the spike needs to maintain structural flexibility to recognize the host cell surface receptors and bind to them, a property that can heavily depend upon the dynamics of the unresolved domains, most prominently the stalk. Construction of the complete, membrane-bound spike model and the description of its dynamics are critical steps in understanding the inner working of this key element of the viral infection by SARS-CoV-2. Combining homology modeling, protein-protein docking, and molecular dynamics (MD) simulations, we have developed a full spike structure in a native membrane. Multimicrosecond MD simulations of this model, the longest known single trajectory of the full spike, reveal conformational dynamics employed by the protein to explore the surface of the host cell. In agreement with cryogenic electron microscopy (cryo-EM), three flexible hinges in the stalk allow for global conformational heterogeneity of spike in the fully glycosylated system mediated by glycan-glycan and glycan-lipid interactions. The dynamical range of the spike is considerably reduced in its nonglycosylated form, confining the area explored by the spike on the host cell surface. Furthermore, palmitoylation of the membrane domain amplifies the local curvature that may prime the fusion. We show that the identified hinge regions are highly conserved in SARS coronaviruses, highlighting their functional importance in enhancing viral infection, and thereby, provide points for discovery of alternative therapeutics against the virus.