D614G reshapes allosteric networks and opening mechanisms of SARS-CoV-2 spikes

The SARS-CoV-2 spike glycoprotein binds human epithelial cells and enables infection through a key conformational transition that exposes its receptor binding domain (RBD). Experimental evidence indicates that spike mutations, particularly the early D614G variant, alter the rate of this conformational shift, potentially increasing viral infectivity. To investigate how mutations reshape the conformational landscape, we conducted extensive weighted ensemble simulations of the Ancestral, Delta, and Omicron BA.1 spike strains along the RBD opening pathway. We observe that Ancestral, Delta, and Omicron BA.1 spike RBDs open differently, with Omicron BA.1 following a more direct opening profile until it reaches a “super-open” state wherein it begins to “peel”, suggesting increased S1 flexibility. Via dynamical network analysis, we identified two allosteric communication networks uniting all S1 domains: the established N2R linker and a newly discovered anti-parallel R2N linker. In Delta and Omicron BA.1 variant spikes, RBD opening is facilitated by both linkers, while the Ancestral strain relies predominantly on the N2R linker. In the ancestral spike, the D614-K854 salt bridge impedes allosteric communication through the R2N linker, whereas the loss of this salt bridge in all subsequent VOCs alleviates local frustration and, we believe, accelerates RBD opening. Hydrogen-deuterium mass spectrometry experiments validate these altered dynamics in the D614 region across Ancestral, D614G, and Omicron BA.1 spikes. This study unveils a ‘hidden’ allosteric network, connecting the NTD to the RBD via the 614-proximal region, and the D614G mutation reshapes the fitness landscape of these critical viral glycoproteins.