Structure-Based Virtual Screening to Discover Potential Lead Molecules for the SARS-CoV-2 Main Protease
Anuj Gahlawat,Navneet Kumar,Rajender Kumar,Hardeep Sandhu,Inder Pal Singh,Saranjit Singh,Anders Sjöstedt,Prabha Garg,Anders Sjöstedt
DOI: https://doi.org/10.1021/acs.jcim.0c00546
IF: 6.162
2020-07-20
Journal of Chemical Information and Modeling
Abstract: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<sup>pro</sup>) in viral replication and gene expression makes this enzyme an attractive drug target. Therefore, inhibition of SARS-CoV-2 M<sup>pro</sup> 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<sup>pro</sup> by facilitating fast proton transfer from the Cys145 to His41 residue. It was observed that active site remained conserved among M<sup>pro</sup> 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<sup>pro</sup> active site and also offers direction for further lead optimization and rational drug design.The Supporting Information is available free of charge at <a class="ext-link" href="/doi/10.1021/acs.jcim.0c00546?goto=supporting-info">https://pubs.acs.org/doi/10.1021/acs.jcim.0c00546</a>.Phylogenetic relationship of SARS-CoV-2 M<sup>pro</sup>; superimposed cocrystallized and redocked ligand; protein thermodynamic stability information upon point mutations; list of top M<sup>pro</sup> inhibitors with their structures, source and scaffold information; interaction fraction histogram obtained from MD simulations; protein-ligand contact timeline information; RMSF plot for protein in presence of different ligands; binding free energy plot; and analysis of impact of mutation on active site residues of M<sup>pro</sup> (<a class="ext-link" href="/doi/suppl/10.1021/acs.jcim.0c00546/suppl_file/ci0c00546_si_001.pdf">PDF</a>)This article has not yet been cited by other publications.
chemistry, multidisciplinary, medicinal,computer science, interdisciplinary applications, information systems