Microsecond simulations and CD spectroscopy reveals the intrinsically disordered nature of SARS-CoV-2 Spike-C-terminal cytoplasmic tail (residues 1242-1273) in isolation

All available SARS-CoV-2 spike protein crystal and cryo-EM structures have shown missing electron densities for cytosolic C-terminal regions. Generally, the missing electron densities point towards the intrinsically disordered nature of the protein region. This curiosity has led us to investigate the C terminal cytosolic region of the spike glycoprotein of SARS-CoV-2 in isolation. The cytosolic regions is supposed to be from 1235-1273 residues or 1242-1273 residues depending on what prediction tool we use. Therefore, we have demonstrated the structural conformation of cytosolic region and its dynamics through computer simulations up to microsecond timescale using OPLS and CHARMM forcefields. The simulations have revealed the unstructured conformation of cytosolic region. Also, in temperature dependent replica-exchange molecular dynamics simulations it has shown its conformational dynamics. Further, we have validated our computational observations with circular dichroism (CD) spectroscopy-based experiments and found its signature spectra at 198 nm which is also adding the analysis as its intrinsically disordered nature. We believe that our findings will surely help us understand the structure-function relationship of the spike protein’s cytosolic region.