Dynamic Allostery Highlights the Evolutionary Differences between the CoV-1 and CoV-2 Main Proteases

The SARS-CoV-2 coronavirus has become one of the most immediate and widely-studied systems since its identification and subsequent global outbreak from 2019-2020. In an effort to understand the biophysical changes as a result of mutations, the mechanics of multiple different proteins within the SARS-CoV-2 virus have been studied and compared with SARS-CoV-1. Focusing on the main protease (mPro), we first explored the long range dynamic-relationship, particularly in cross-chain dynamics, using the Dynamic Coupling Index (DCI) to investigate the dynamic coupling between the catalytic site residues and the rest of the protein, both inter and intra chain for the CoV-1 and CoV-2 mPro. We found that there is significant cross-chain coupling between these active sites and distal residues in the CoV-2 mPro but it was missing in CoV-1. The enhanced long distance interactions, particularly between the two chains, suggest subsequently enhanced cooperativity for CoV-2. A further comparative analysis of the dynamic flexibility using the Dynamic Flexibility Index (DFI) between the CoV-1 and CoV-2 mPros shows that the inhibitor binding near active sites induces change in flexibility to a distal region of the protein, opposite in behavior between the two systems; this region becomes more flexible upon inhibitor binding in CoV-1 while it becomes less flexible in the CoV-2 mPro. Upon inspection, we show that, on average, the dynamic flexibility of the sites substituted from CoV-1 to CoV-2 changes significantly less than the average calculated across all residues within the structure, indicating that the differences in behaviors between the two systems is likely the result of allosteric influence, where the new substitutions in COV-2 induce flexibility and dynamical changes elsewhere in the structure.