Department of Pharmacoinformatics, National Institute of Pharmaceutical Education and Research (NIPER), S.A.S. Nagar 160062, Punjab, India.
Department of Clinical Microbiology and Laboratory for Molecular Infection Medicine Sweden (MIMS), Umeå University, SE-90185 Umeå, Sweden.
J Chem Inf Model. 2020 Dec 28;60(12):5781-5793. doi: 10.1021/acs.jcim.0c00546. Epub 2020 Aug 4.
The COVID-19 disease is caused by a new strain of the coronavirus family (SARS-CoV-2), and it has affected at present millions of people all over the world. The indispensable role of the main protease (M) in viral replication and gene expression makes this enzyme an attractive drug target. Therefore, inhibition of SARS-CoV-2 M as a proposition to halt virus ingression is being pursued by scientists globally. Here we carried out a study with two objectives: the first being to perform comparative protein sequence and 3D structural analysis to understand the effect of 12 point mutations on the active site. Among these, two mutations, viz., Ser46 and Phe134, were found to cause a significant change at the active sites of SARS-CoV-2. The Ser46 mutation present at the entrance of the S5 subpocket of SARS-CoV-2 increases the contribution of other two hydrophilic residues, while the Phe134 mutation, present in the catalytic cysteine loop, can cause an increase in catalytic efficiency of M by facilitating fast proton transfer from the Cys145 to His41 residue. It was observed that active site remained conserved among M of both SARS-CoVs, except at the entrance of the S5 subpocket, suggesting sustenance of substrate specificity. The second objective was to screen the inhibitory effects of three different data sets (natural products, coronaviruses main protease inhibitors, and FDA-approved drugs) using a structure-based virtual screening approach. A total of 73 hits had a combo score >2.0. Eight different structural scaffold classes were identified, such as one/two tetrahydropyran ring(s), dipeptide/tripeptide/oligopeptide, large (approximately 20 atoms) cyclic peptide, and miscellaneous. The screened hits showed key interactions with subpockets of the active site. Further, molecular dynamics studies of selected screened compounds confirmed their perfect fitting into the subpockets of the active site. This study suggests promising structures that can fit into the SARS-CoV-2 M active site and also offers direction for further lead optimization and rational drug design.
新型冠状病毒(SARS-CoV-2)引发的 COVID-19 疫情目前已在全球范围内影响了数百万人。主蛋白酶(M)在病毒复制和基因表达过程中发挥着不可或缺的作用,这使其成为极具吸引力的药物靶标。因此,全球科学家都在致力于研究抑制 SARS-CoV-2 M 以阻止病毒入侵的方法。在此,我们开展了一项研究,旨在实现两个目标:其一,通过比较蛋白序列和三维结构分析,了解 12 个点突变对活性位点的影响。其中,两个突变,即丝氨酸 46 位(Ser46)和苯丙氨酸 134 位(Phe134),被发现对 SARS-CoV-2 的活性位点产生了显著影响。SARS-CoV-2 中 S5 亚口袋入口处的 Ser46 突变增加了其他两个亲水残基的贡献,而催化半胱氨酸环中的 Phe134 突变可以通过促进 Cys145 到 His41 残基的快速质子转移,提高 M 的催化效率。结果表明,除了 S5 亚口袋入口处,M 中的活性位点在两种 SARS-CoV 之间均保持保守,这表明其维持了底物特异性。其二,我们采用基于结构的虚拟筛选方法,筛选了三个不同数据集(天然产物、冠状病毒主蛋白酶抑制剂和 FDA 批准的药物)的抑制效果。共有 73 个命中物的组合得分>2.0。确定了 8 种不同的结构骨架类别,如一个/两个四氢吡喃环(s)、二肽/三肽/寡肽、大环(约 20 个原子)环肽和杂环。筛选出的命中物与活性位点的亚口袋存在关键相互作用。进一步的分子动力学研究证实,所选筛选化合物能够完美地适配活性位点的亚口袋。本研究为进一步的先导优化和合理药物设计提供了方向,表明有希望的结构可以与 SARS-CoV-2 M 的活性位点结合。