Differential dynamic behavior of prefusion spike proteins of SARS coronaviruses 1 and 2

by Kumar, Vivek Govind; Ogden, Dylan S.; Isu, Ugochi H.; Polasa, Adithya; Losey, James; Moradi, Mahmoud

Within the last two decades, severe acute respiratory syndrome (SARS) coronaviruses 1 and 2 (SARS-CoV-1 and SARS-CoV-2) have caused two major outbreaks. For reasons yet to be fully understood the COVID-19 outbreak caused by SARS-CoV-2 has been significantly more widespread than the 2003 SARS epidemic caused by SARS-CoV-1, despite striking similarities between the two viruses. One of the most variable genes differentiating SARS-CoV-1 and SARS-CoV-2 is the S gene that encodes the spike glycoprotein. This protein mediates a crucial step in the infection, i.e., host cell recognition and viral entry, which starts with binding to the host cell angiotensin converting enzyme 2 (ACE2) protein for both viruses. Recent structural and functional studies have shed light on the differential binding behavior of the SARS-CoV-1 and SARS-CoV-2 spike proteins. In particular, cryogenic electron microscopy (cryo-EM) studies show that ACE2 binding is preceded by a large-scale conformational change in the spike protein to expose the receptor binding domain (RBD) to its binding partner. Unfortunately, these studies do not provide detailed information on the dynamics of this activation process. Here, we have used an extensive set of unbiased and biased microsecond-timescale all-atom mol. dynamics (MD) simulations of SARS-CoV-1 and SARS-CoV-2 spike protein ectodomains in explicit solvent to determine the differential behavior of spike protein activation in the two viruses. Our results based on nearly 50µs of equilibrium and nonequilibrium MD simulations indicate that the active form of the SARS-CoV-2 spike protein is considerably more stable than the active SARS-CoV-1 spike protein. Unlike the active SARS-CoV-2 spike, the active SARS-CoV-1 spike spontaneously undergoes a large-scale conformational transition to a pseudo-inactive state, which occurs in part due to interactions between the N-terminal domain (NTD) and RBD that are absent in the SARS-CoV-2 spike protein. Steered MD (SMD) simulations indicate that the energy barriers between active and inactive states are comparatively lower for the SARS-CoV-1 spike protein. Based on these results, we hypothesize that the greater propensity of the SARS-CoV-2 spike protein to remain in the active conformation contributes to the higher transmissibility of SARS-CoV-2 in comparison to SARS-CoV-1. These results strongly suggest that the differential binding behavior of the active SARS-CoV-1 and 2 spike proteins is not merely due to differences in their RBDs as other domains of the spike protein such as the NTD could play a crucial role in the effective binding process, which involves the prebinding activation. Therefore, our hypothesis predicts that mutations in regions such as the NTD, which is not directly involved in binding, may lead to a change in the effective binding behavior of the coronavirus.

Journal
bioRxiv
Year
2020
Start Page
1-42
URL
https://dx.doi.org/10.1101/2020.12.25.424008
DOI
10.1101/2020.12.25.424008