Computational evaluation and benchmark study of 342 crystallographic holo-structures of SARS-CoV-2 Mpro enzyme

Hamlet Khachatryan,Mher Matevosyan,Vardan Harutyunyan,Smbat Gevorgyan,Anastasiya Shavina,Irina Tirosyan,Yeva Gabrielyan,Marusya Ayvazyan,Marine Bozdaganyan,Zeynab Fakhar,Sajjad Gharaghani,Hovakim Zakaryan

DOI: https://doi.org/10.1038/s41598-024-65228-5

IF: 4.6

2024-06-22

Scientific Reports

Abstract:The coronavirus disease 19 pandemic, caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has led to a global health crisis with millions of confirmed cases and related deaths. The main protease (Mpro) of SARS-CoV-2 is crucial for viral replication and presents an attractive target for drug development. Despite the approval of some drugs, the search for effective treatments continues. In this study, we systematically evaluated 342 holo-crystal structures of Mpro to identify optimal conformations for structure-based virtual screening (SBVS). Our analysis revealed limited structural flexibility among the structures. Three docking programs, AutoDock Vina, rDock, and Glide were employed to assess the efficiency of virtual screening, revealing diverse performances across selected Mpro structures. We found that the structures 5RHE, 7DDC, and 7DPU (PDB Ids) consistently displayed the lowest EF, AUC, and BEDROCK scores. Furthermore, these structures demonstrated the worst pose prediction results in all docking programs. Two structural differences contribute to variations in docking performance: the absence of the S1 subsite in 7DDC and 7DPU, and the presence of a subpocket in the S2 subsite of 7DDC, 7DPU, and 5RHE. These findings underscore the importance of selecting appropriate Mpro conformations for SBVS, providing valuable insights for advancing drug discovery efforts.

multidisciplinary sciences

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to identify the conformation most suitable for structure - based virtual screening (SBVS) by systematically evaluating the crystal structures of 342 SARS - CoV - 2 main proteases (Mpro). Specifically, the researchers aim to:

1. **Evaluate structural flexibility**: Analyze the crystal structures of these Mpro to evaluate their structural flexibility. This helps to understand the differences in the performance of different structures in virtual screening.

2. **Select the best conformation**: Determine which Mpro structures are most suitable for SBVS, thereby improving the efficiency of drug discovery. Three different docking software (AutoDock Vina, rDock, and Glide) were used in the study to evaluate the efficiency of virtual screening, and the performance of different structures was measured by multiple performance indicators (such as enrichment factor EF, area under the ROC curve AUC, BEDROCK score).

3. **Analyze structural differences**: Through structural analysis, especially the analysis of binding sites, identify the key factors that cause differences in the docking performance of different structures. The study found that the absence of the S1 subsite in the 7DDC and 7DPU structures and the presence of the secondary pocket of the S2 subsite in the 7DDC, 7DPU, and 5RHE structures are the main reasons for affecting the docking performance.

4. **Provide guidance and suggestions**: Based on the above analysis, provide valuable insights for future drug discovery work, especially when selecting appropriate Mpro conformations for SBVS, emphasizing the importance of selecting appropriate conformations.

Through these research objectives, the paper hopes to provide more effective strategies and methods for the drug development of SARS - CoV - 2 Mpro.