6lu7-2point-25-Grel-0-0.25.xlsx from 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 no publicly-available vaccine or antiviral drugs at the time 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, targeting regions other than 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^pro homodimer protein and analyse its dynamics and thermodynamics. We found a rich and heterogeneous dynamical structure, 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, yet control 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 fragment-based binding screens. The results also provide new insights into the active biophysical research field of protein fluctuation allostery and its underpinning dynamical structure.

由新型冠状病毒（SARS-CoV-2）引发的COVID-19大流行在本文撰写之际，仍无公开可用的疫苗或抗病毒药物。冠状病毒的一类极具研究价值的药物靶点为主要蛋白酶（M^pro，又称3CL^pro），因其在病毒增殖周期中扮演着不可或缺的关键角色。目前已有大量研究聚焦于开发可结合并阻断主蛋白酶活性位点的抑制剂，但针对非竞争性抑制机制、靶向活性位点以外区域的相关研究却较为匮乏，部分原因在于此类变构调控的基础生物物理机制仍未得到充分阐释。本研究构建了SARS-CoV-2主蛋白酶同源二聚体蛋白的弹性网络模型（ENM），并对其动力学与热力学特性展开分析。研究团队发现该蛋白存在丰富且具有异质性的动力学结构，其中包括同源二聚体蛋白酶活性位点之间的变构关联运动。通过对该弹性网络模型开展全面的单点与两点突变扫描，并分析其对涨落自由能的影响，本研究不仅验证了此前实验已确认的生物活性残基，还发现了数个远离活性位点却可调控蛋白酶功能的全新候选区域。研究结果显示，这类以动态调控为核心的新型区域，可作为该蛋白酶非竞争性抑制结合位点的潜在候选靶点，有望为当前基于片段的结合筛选开发提供助力。本研究结果还为蛋白质涨落变构学这一活跃的生物物理研究领域及其支撑性动力学结构提供了全新的研究视角。