Computational Analysis of Dynamic Allostery and Control in the SARS-CoV-2 Main Protease

The COVID-19 pandemic caused by the novel coronavirus SARS-CoV-2 has generated a global pandemic and no vaccine or antiviral drugs exist at the moment of writing. An attractive coronavirus drug target is the main protease (M^pro, also known as 3CL^pro) because of its vital role in the viral cycle. A significant body of work has been focused on finding inhibitors which bind and block the active site of the main protease, but little has been done to address potential non-competitive inhibition which targets regions beyond the active site, partly because the fundamental biophysics of such allosteric control is still poorly understood. In this work, we construct an Elastic Network Model (ENM) of the SARS-CoV-2 M^prohomodimer protein and analyse the dynamics and thermodynamics of the main protease’s ENM. We found a rich and heterogeneous dynamical structure in the correlated motions, including allosterically correlated motions between the homodimeric protease’s active sites. Exhaustive 1-point and 2-point mutation scans of the ENM and their effect on fluctuation free energies confirm previously experimentally identified bioactive residues, but also suggest several new candidate regions that are distant from the active site for control of the protease function. Our results suggest new dynamically-driven control regions as possible candidates for non-competitive inhibiting binding sites in the protease, which may assist the development of current fragmentbased binding screens. The results also provide new insight into the protein physics of fluctuation allostery and its underpinning dynamical structure.