==== Front Sci Rep Sci Rep Scientific Reports 2045-2322 Nature Publishing Group UK London 37386052 36342 10.1038/s41598-023-36342-7 Article Repurposing and computational design of PARP inhibitors as SARS-CoV-2 inhibitors Rampogu Shailima shailima.rampogu@gmail.com 1 Jung Tae Sung jungts@gnu.ac.kr 2 Ha Min Woo minuha@jejunu.ac.kr 3 Lee Keun Woo kwlee@gnu.ac.kr 1 1 grid.256681.e 0000 0001 0661 1492 Department of Bio and Medical Big Data (BK4 Program), Division of Life Sciences, Research Institute of Natural Science (RINS), Gyeongsang National University (GNU), Jinju, Republic of Korea 2 grid.256681.e 0000 0001 0661 1492 Laboratory of Aquatic Animal Diseases, College of Veterinary Medicine, Research Institute of Natural Science, Gyeongsang National University, Jinju, 52828 Republic of Korea 3 grid.411277.6 0000 0001 0725 5207 Interdisciplinary Graduate Program in Advanced Convergence Technology and Science, Jeju National University, 102 Jejudaehak-ro, Jeju, 63243 Republic of Korea 29 6 2023 29 6 2023 2023 13 105835 10 2022 1 6 2023 © The Author(s) 2023 https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. Coronavirus disease 2019 (COVID-19) is a recent pandemic that caused serious global emergency. To identify new and effective therapeutics, we employed a drug repurposing approach. The poly (ADP ribose) polymerase inhibitors were used for this purpose and were repurposed against the main protease (Mpro) target of severe acute respiratory syndrome Coronavirus 2 (SARS-CoV-2). The results from these studies were used to design compounds using the ‘Grow Scaffold’ modules available on Discovery Studio v2018. The three designed compounds, olaparib 1826 and olaparib 1885, and rucaparib 184 demonstrated better CDOCKER docking scores for Mpro than their parent compounds. Moreover, the compounds adhered to Lipinski’s rule of five and demonstrated a synthetic accessibility score of 3.55, 3.63, and 4.30 for olaparib 1826, olaparib 1885, and rucaparib 184, respectively. The short-range Coulombic and Lennard-Jones potentials also support the potential binding of the modified compounds to Mpro. Therefore, we propose these three compounds as novel SARS-CoV-2 inhibitors. Subject terms Virtual drug screening Drug discovery and development issue-copyright-statement© Springer Nature Limited 2023 ==== Body pmcIntroduction Recent Coronavirus disease 2019 (COVID-19) events have led the world into unprecedented circumstances1,2. COVID-19 is caused by the severe acute respiratory syndrome Coronavirus 2 (SARS-CoV-2). It belongs to the family Coronaviridae and subfamily Orthocoronavirinae and consists of four genera3. The genome of SARS-CoV-2 is 30 kb in size, and codes for a large non-structural polyprotein that further generates 15/16 proteins, five accessory proteins, and four structural proteins after proteolytic cleavage3. The identification of new inhibitors with effective and quick therapeutics is still in progress for this disease. An approach to discovering novel inhibitors is to repurpose approved drugs to identify new indications for older drugs4. This approach is one of the most effective methods for drug discovery5. Drug repurposing (DR) also called drug reprofiling, drug recycling, drug rescue, therapeutic switching, drug retasking, and drug redirection5. Many studies have adopted this approach to discover COVID-19 therapeutics, which has resulted in the discovery of plausible inhibitors6–9. Several studies have also reported the intervention of computational approaches in detecting SARS-CoV-2 inhibitors using predominantly molecular docking, molecular dynamics simulation, and pharmacophore modelling studies10–24. In this study, we applied the DR approach to a SARS-CoV-2 main protease (Mpro) target using PARP inhibitors. PARP enzymes contribute greatly to genome stability, and their inhibitors occupy the catalytic domain of PARP, hindering poly ADP-ribosylation (PARylation) of target proteins25. The present study is focused on Mpro, as it has the largest active site volume (Supplementary Fig. 1). The SARS-CoV-2 Mpro processes the proteolytic step during the replication of the virus. This 33.8-kDa protein, also known as 3-chymotrypsin-like protease (3C-like protease), is a promising target for developing and identifying new inhibitors26–29. The homodimer Mpro consists of two protomers30. Mpro consists of 306 amino acid residues divided into three domains31. The residues 8–101 form domain I; residues 102–184 form domain II; and residues 201–303 form domain III26,30,32. Domain III is connected to domain II by residues 185–200. A Cys-His catalytic dyad is present between domains I and domain II26,30–32. The S1 subsite residues were Phe140, Leu141, Asn142, His163, Glu166, and His172. The residues in S1 were Thr25, Thr26, and Leu2733. The hydrophobic S2 subsite consists of His41, Met49, Tyr54, Met165, and Asp187 residues. The S4 binding subsite contains Met165, Leu167, Phe185, Gln189, and Gln19233,34. The residues Cys145-His41 form the catalytic dyad33,34. All these residues bind to the ligand and are termed key residues. In the present study, we used approved PARP inhibitors to target COVID-19 Mpro. PARP inhibitors exhibit anticancer ability35–37. It is reported that the PARP inhibitors can reduce the levels of interleukin 1 (IL‐1), interleukin 6 (IL6), and tumour necrosis factor alpha (TNF‐α); and decrease lung fibrosis38. These inhibitors act against cell death promoted by inflammation, thereby favouring cell survival38. In another study, computational and experimental analyses revealed that the PARP1 inhibitor CVL218 is a potential COVID-19 inhibitor39. Another study reported that the in vitro PARP inhibitor tenoparib prevented the replication of SARS-CoV-2 and human coronavirus NL63 (HCoV-NL63)40. The PARP inhibitor mefuparib has been reported to bind to the nucleocapsid (N) protein39 and inhibit SARS-CoV-2 replication with a half maximal inhibitory concentration (IC50) of 5.12 µM41. These findings highlight the potential of novel PARP inhibitors in treating COVID-19. Based on these reports, we attempted to computationally design new PARP compounds with high affinity for SARS-CoV-2 targets by using Grow Scaffold, a drug design approach. This approach locates spaces within the active site after molecular docking to grow the scaffold. Since we selected PARP inhibitors for this study, we diligently sought compounds that exhibited greater affinity towards the COVID-19 Mpro target. The infographic workflow is illustrated in Fig. 1.Figure 1 Illustration presenting the study methodology. (a) Identification of the unoccupied space in the active site of Mpro. (b) Molecular docking of the known PARP inhibitors at the active site of Mpro and applying Grow Scaffold to design the compounds (pink oval). (c) Database creation with compounds obtained from Grow Scaffold. (d) Target ready for molecular docking. (e) Obtaining a protein- ligand complex with modified compounds. The icons are taken from BioRender.com and are modified accordingly. Results Step 1: Molecular docking based binding affinity studies to identify the empty spaces in the binding pocket Molecular docking was performed using the validated SARS-CoV-2 target, Mpro. All the chosen compounds, except Veliparib and Pamiparib, demonstrated a high docking score compared to the reference compound (co-crystallised ligand). Based on the docking scores, the high-scoring compounds, rucaparib and olaparib, were considered for further studies (Supplementary Table 1). Here, molecular docking was performed to discover compounds with better docking scores towards the COVID-19 target and to select empty space for designing new compounds. Strategy to design new compounds To design a novel compound, the target-ligand active site of Mpro was scrupulously examined, as it had a large volume. After careful visual analysis using DS, we identified two empty sites with olaparib (site 1 and site 2) and one empty site with rucaparib (site 1) and called them R1 (Supplementary Fig. 2). At the identified R1 position, we applied the “Grow Scaffold” approach, which has inbuilt libraries that facilitate the addition of groups to the small molecule (Supplementary Fig. 3). The olaparib site 1 generated a total of 3,201 compounds, site 2 generated 304 compounds; and site 1 of rucaparib generated 437 compounds, accounting for 3942 compounds after these identified sites were subjected to Grow Scaffold modules (Fig. 2a). The simplified molecular-input line-entry system (SMILES) representations of these compounds are provided in (Supplementary Table 2).Figure 2 Generation of novel compounds and drug-like database creation. (a) Marking the site for growing the scaffold in the protein–ligand complex. (b) Screening the compound for molecular docking. The 3942 compounds were subjected to Lipinski’s rule of five, enabling the ‘Filter by Lipinski and Veber Rules’ protocol. According to the “Rule of Five”, a drug-like molecule should have maximum five hydrogen bond donors, 10 hydrogen bond acceptors, a molecular weight less than 500 Dalton, and a LogP value less than five. This process resulted in the production of 19 compounds. These compounds were subsequently docked to Mpro to delineate their binding affinities (Fig. 2b). Step 2: Binding affinity studies and clustering of small molecules In this step, all 19 modified compounds were docked into the Mpro-binding pocket. The results revealed that three compounds displayed higher binding affinities than their parent structures (Table 1 and Supplementary Table 1). Manual clustering of the compounds was performed to understand the prospective binding modes of the small molecules (Fig. 3). The pose with the best docking score from the largest cluster was further examined using molecular dynamics simulation (MDS).Table 1 Various intermolecular interactions of the modified compounds and the target Mpro. -CDOCKER interaction energy (kcal/mol) Various intermolecular interactions Compounds Mpro (6LU7) Hydrogen bond interactions Alkyl/π-alkyl interactions Van der Waals interactions Olaparib 1826 58.89 Gly143 and Gln189 His41, Cys145, His164, and Met165 Met49, Phe140, Leu141, Asn142, Ser144, His163, Glu166, His172, Asp187, and Arg188 Olaparib 1885 57.51 Gly143, Ser144, and Cys145 His41, Cys145, and Met164 Met49, Phe140, Leu141, Glu166, His172, Asp187, Arg188, and Glu189 Rucaparib 184 58.80 Arg40, His41, and Gln192 Met49, His 41, and Met165 Thr45, Ser46, Tyr54, Cys85, Glu166, Leu167, Pro168, Phe181, Thr190, and Ala191 Figure 3 Clustering of small molecules at the active site of Mpro. (a) Olaparib1826; (b) Olaparib1885; (c) Rucaparib 184. The two Olaparib-modified compounds were well aligned within the binding pocket of the target. Olaparib 1826 formed 27 poses in a cluster (Fig. 3a), whereas olaparib 1885 formed 18 poses in a cluster (Fig. 3b). Whereas, rucaparib184 formed a cluster with 7 poses (Fig. 3c). These results indicate that olaparib 1826 may have a higher affinity for Mpro, followed by olaparib 1885 and rucaparib184. The three compounds demonstrated satisfactory results, with higher docking scores than the parent compounds (Table 1 and Supplementary Table 1). The parent compounds, olaparib and rucaparib, generated dock scores of 50.38 kcal/mol and 58.80 kcal/mol, respectively. The modified olaparib compounds, olaparib 1826 and olaparib 1885, demonstrated a dock score of 58.89 kcal/mol and 57.51 kcal/mol, respectively. Similarly, the modified rucaparib compound, rucaparib 184, had a dock score of 58.80 kcal/mol. These findings suggest that the modified compounds demonstrate a stronger affinity for Mpro than the parent structure (Table 1 and Supplementary Table 1). From the largest cluster, the compounds with the best docking scores were upgraded to molecular dynamics simulation (MDS) to understand the binding potential of the modified compounds to the binding pocket of the target. MDS analysis MDS was performed for 100 ns to analyse the binding potential of the ligands at the binding pocket of the targets. The MDS analysis was based on the root mean square deviation (RMSD), radius of gyration (Rg), root mean square fluctuation (RMSF), and hydrogen bond number. Stability analysis by RMSD The protein backbone demonstrated stability42 throughout the simulation. All the systems converged well and exhibited RMSD below 0.3 nm. The system with olaparib 1826 was largely stable during MDS evolution. At ~ 11,000 ps, a minute surge in the RMSD profiles was noted, and thereafter the system remained steady with an average of 0.16 nm (Fig. 4a). The RMSD profiles of olaparib 1885 and rucaparib 184 were stable throughout the MDS evolution without any noticeable variations, with an average of 0.17 nm and 0.2 nm, respectively (Fig. 4a).Figure 4 MDS results for the three systems. (a) Deviation analysis of the three systems according to RMSD. (b) Compactness analysis according to Rg. (c) Fluctuation analysis of the three systems according to RMSF. (d) Analysis to determine the number of hydrogen bonds for the three systems. Compactness analysis by Rg Rg is “defined as the ratio of the protein accessible surface area to that of an ideal sphere of the same volume”43. This principle governs the protein-folding mechanism43. Our analysis revealed that Rg defines the distance between each atom of a protein and its centroid, thereby defining its compactness44. These three systems showed that the protein backbone was highly stable and compact. The readings exhibit a range of 2.20 nm to 2.28 nm with an average of 2.21 nm, 2.22 nm, and 2.22 nm for olaparib 1826, olaparib 1885, and rucaparib 184, respectively (Fig. 4b). This finding demonstrates that the systems were remarkably compact. Fluctuation analysis by RMSF The RMSF plots reveal fluctuations during the MDS run. All the systems were stable, with no major fluctuations. The systems were stable with RMSF below 0.3 nm (Fig. 4c). The average RMSF was 0.17 nm, 0.12 nm, and 0.13 nm for olaparib 1826, olaparib 1885, and rucaparib 184, respectively (Fig. 4c). This result indicates that each residue in the system was stable. Number of hydrogen bonds The number of hydrogen bonds between the protein and ligand was assessed during the entire simulation run. The presence of stable interactions suggested that ligands were held within the binding pocket throughout the simulation run. The three systems demonstrated hydrogen bond interactions during the entire simulation run, with an average of 2.0 for olaparib 1826, 1.3 for olaparib 1885, and 4.4 for rucaparib 184 (Fig. 4d). Rucaparib184 demonstrated a higher number of hydrogen bonds than the olaparib derivatives (Fig. 4d). This result implied that the compounds were firmly presented in the binding pocket of the protein. Binding mode analysis From the stable RMSD of the last 10 ns (~ 90,000–100,000 ps), representative complex structures were extracted and superimposed onto the X-ray crystal structure of Mpro. The results demonstrated that the ligands occupied a binding pocket similar to that of the co-crystallised ligand (Fig. 5), held by various intermolecular interactions. The modified compounds settled at different subsites of the target-binding pocket held by various interactions.Figure 5 Superimposition of MDS derived structures against the X-ray structure. (a) Compounds present at the Mpro binding pocket. (b) The zoomed version of the accommodation of the compounds at the binding pocket. Key residue interaction of the modified compounds with the target residues Olaparib1826 The modified compound olaparib 1826, formed hydrogen bonds with Gly143 and Gln189 (Fig. 9a). Interestingly, oxygen atoms were involved in these interactions. The distances between the atoms involved in hydrogen bonding were calculated. The Gly143: HN-O3 interaction was highly stable, with an average of 0.21 nm throughout the simulation run (Fig. 6a). Similarly, the interaction between Gln189:HE21-O2 was also stable with an average of 0.28 nm, which indicates that the ligand tightly adheres to the binding pocket of the target (Fig. 6b). The key residue, His164, formed carbon hydrogen bonds to hold the ligand in the binding pocket. The residue Cys145 formed the π-alkyl interaction with the modified compound. The residues His41, His164, and Met165 adhered to the compound via alkyl interactions. His164 interacted with the fluorine atom of the ligand, thereby settling the ligand in the binding pocket of the target. Other residues, such as Met49, Phe140, Leu141, Asn142, Ser144, His163, Glu166, His172, Asp187, and Arg188 (Table 1) firmly held the compound in the binding pocket via van der Waals interactions (Fig. 9d).Figure 6 Hydrogen bond distances between the interacting protein atoms and the ligand atoms. (a) Distance between Gly143: HN-O3. (b) Distance between Gln189:HE21-O2. Olaparib 1885 Olaparib 1885 formed hydrogen bonds with the Gly143, Ser144, and Cys145 residues (Fig. 9b). The hydrogen bond distance between Gly143: HN-O3 was stable, with no major aberrations. The average distance measured was found to be 0.2 nm (Fig. 7a). The hydrogen bond distance between Ser144: HN-O3 was analysed during the evolution of MDS (Fig. 7b). These findings reveal that the distance was largely stable, with a dip in the profile at 89,900 ps, which remained stable thereafter. Interestingly, the overall average distance was measured at 0.33 nm, while the last 10 ns were measured at 0.26 nm, suggesting that Ser144 strongly interacted with the ligand (Fig. 7b). Hydrogen bond interactions between Cys145: HN-O3 were stable throughout the simulation. There was a minor depression in the distance plot at 85,000 ps, which remained stable thereafter. While the overall distance was measured to be 0.31 nm, the average measurement for the last 10 ns was 0.23 nm (Fig. 7c). Generally, these three interactions are firm and hold the ligand in the binding pocket. Residues Asn142, His163, and His164 formed carbon–hydrogen bonds that firmly held the ligand at the active site. The key residues His41 and Cys145 have generated alkyl and π-alkyl interactions, thereby positioning the ligand at the binding site. The residue Met164 has prompted an amide π-stacked interaction. Residues Met49, Phe140, Leu141, Glu166, His172, Asp187, Arg188, and Glu189 (Table 1) promoted van der Waals interactions, thereby aiding the ligand to be seated at the binding pocket. The Ser144 residue interacts with the fluorine atom of the ligand (Fig. 9e).Figure 7 Hydrogen bond distance between the interacting protein atoms and the ligand atoms. (a) Distance between Gly143: HN-O3. (b) Distance between Ser144: HN-O3. (c) Distance between Cys145: HN-O3. Rucaparib184 Residues Arg40, His41, and Gln192 formed hydrogen bonds with rucaparib 184 (Fig. 9c). The hydrogen bond distances between the protein residue atoms and ligand atoms indicated that the interactions were stable during the simulation run. The Arg40: HE-O47 interaction has seen an elevation during the initial steps of the simulation run to 15,710 ps; however, it was soon stable thereafter with an overall average of 0.19 nm (Fig. 8a). Another interaction was observed between His41: HD1-O49. The distance profile displayed a stable interaction, although a few fluctuations were observed from 10 to 11 ns. The average distance was measured to be 0.25 nm during the evolution of the simulation (Fig. 8b). The hydrogen bond interaction distance between Gln192: HE21–N5 was also stable throughout the simulation. However, the average distance was projected to be 0.46 nm. This distance measured seems to be marginally weaker than the previous interactions that allowed the ligand to settle at the binding pocket (Fig. 8c). His164 and Glu189 residues adhered to the compound via carbon-hydrogen bonds. The key residues Met49, His41 and Met165 have prompted π-alkyl and π–π T shaped interactions with the ligand. The residue Cys44 formed an interaction with the fluorine atom, and Arg40 generated an attractive charge interaction, holding the ligand in the binding pocket of the compound. The residues Thr45, Ser46, Tyr54, Cys85, Glu166, Leu167, Pro168, Phe181, Thr190 and Ala191 (Table 1) have interacted with the ligand via the van der Waals interactions, thus positioning the ligand at the binding pocket (Fig. 9f). Analysis of the modified groups of the three compounds revealed that rucaparib 184 generated two hydrogen bonds with key residues.Figure 8 Hydrogen bond distances between the interacting protein atoms and the ligand atoms. (a) Distance between Arg40:HE-O47; (b) Distance between His41:HD1-O49; (c) Distance between Gln192:HE21-N5. Figure 9 Comprehensive intermolecular interactions. (a–c) represent the hydrogen bond interactions between protein and ligand atoms of olaparib 1826, olaparib 1885, and rucaparib184. (d–f) represent the overall 2D interactions of olaparib1826, olaparib1885, and rucaparib 184 with the residues of the target. Interaction energies between protein and ligands Additionally, protein-ligand (Protein-LIG) interactions were confirmed by retrieving the short-range Coulomb (Coul-SR) and short-range Lennard-Jones (LJ-SR) energy terms. During the simulation, the complexes projected negative values. The average Coul-SR: Protein-LIG values for olaparib 1826, olaparib 1885, and rucaparib 184 were − 93.14 kJ/mol, 86.77 kJ/mol, and − 322.26 kJ/mol, respectively (Fig. 10a). The LJ-SR: Protein-LIG values for olaparib 1826, olaparib 1885 and rucaparib 184 are − 78.40 kJ/mol, − 165.11 kJ/mol, and − 131.55 kJ/mol, respectively (Fig. 10b). The average results of LJ-SR: Protein-LIG were similar for all ligands. However, the results of Coul-SR: Protein-LIG revealed that olaparib 1826 and olaparib 1885 demonstrated similar results, but rucaparib 184 projected a better result than both olaparib 1826 and olaparib 1885. These findings suggest that the modified compounds formed thermodynamically firm bonds with Mpro, as reported earlier45.Figure 10 Interaction energy calculation. (a) Coul-SR: Protein-LIG. (b) LJ-SR: Protein-LIG. Discussion PARP proteins demonstrate structural similarity and function, with two riboses and two phosphates in a unit polymer36. PARP1 was discovered in 196336. They perform critical functions including apoptosis, DNA damage response, and transcription modulation46,47. Their inhibitors exert anticancer activities36,46,48. In the current study, PARP inhibitors were used against COVID-19 target to specifically design putative inhibitors with affinity for SARS-CoV-2 Mpro. To identify effective therapeutics for COVID-19, PARP inhibitors were employed as starting structures. Here, we computationally designed compounds that demonstrate higher affinity for Mpro than their parent structures. To design new molecules computationally, we first identified a target with a large active site volume. We then marked the unoccupied spaces within the binding pocket after molecular docking. Correspondingly, Mpro was identified as the target and PARP inhibitors were docked to it. Olaparib and rucaparib had higher docking scores for Mpro. Therefore, these compounds were used in subsequent experiments. The Mpro-ligand complex was selected to design new compounds by identifying the empty space located in the binding pocket (unoccupied by the ligand). Accordingly, with olaparib, two sites were identified (site 1 and site 2) and one site was identified with rucaparib (site 1) (Fig. 11, Supplementary Fig. 2). These sites were marked R1 and the Grow Scaffold module was initiated.Figure 11 Clustering of the modified compounds post docking. The modified group occupies the same designated position that is marked. Molecular docking was initiated, followed by clustering analysis, to evaluate whether the new compounds induced conformational changes. The results revealed that the modified groups occupied positions at the desired sites, which aligned with the hypothesis. Furthermore, the binding affinity results showed that two compounds (olaparib 1826 and olaparib 1885) from olaparib site1 showed stronger affinity than site2 and occupied a defined position (Fig. 11), as did one compound from site3 (rucaparib 184). The olaparib-modified compounds showed better results against Mpro than the rucaparib-modified compounds. However, the modified compounds generated higher docking scores than the parent compounds. Olaparib 1826 was obtained from the fragment library of organosilanes, reaction name yama coupling, and fragment name 2-(triethoxysilyl)ethylamine. Olaparib 1885 was obtained from the fragment library GrignardReagents, the reaction name Kumada Coupling, and the fragment name cyclobutylmagnesium chloride. Rucaparib 184 was synthesized using the fragment library acids with the fragment name monosodium L-aspartate dehydrate #2 and a reaction called esterification. In the MDS studies, olaparib 1826 and olaparib 1885 maintained a similar binding mode to that seen in molecular docking, while rucaparib 184 changed dramatically from its initial position. However, all the compounds were accommodated at the active site of the target while adhering to key residue interactions (Fig. 12).Figure 12 Binding mode examination for the initial pose and MDS pose. The molecular docked pose is represented in grey and the final MDS pose is represented in colour code. Molecular docking interactions revealed that with Mpro, olaparib 1826 had hydrogen bond interactions with the key residues Gly143, Ser144, Gln189, and Thr190. MDS retained hydrogen bonds with Gly143 and Gln189. Interactions with Gly143 have been reported previously28,30,49. Similarly, a hydrogen bond was also observed with Gln189 in the molecular docking and MDS pose, as reported previously50–52. Hydrogen bonds between these residues were also observed in the X-ray structures. In addition, several other interactions originating from the binding pocket of the target firmly hold the ligand at its active site. The distance measured for these residues was also within the acceptable length below 0.3 nm during the progression of MDS (Fig. 6). Moreover, targeting the catalytic dyad residues His41 and Cys145 is important for the developing of Mpro inhibitors53. Our results show that olaparib 1826 formed van der Waals interactions with His41, and alkyl interactions with Cys145. No change in the binding mode was observed in the molecular docking or MDS poses (Fig. 12a). The compound olaparib 1885 formed hydrogen bonds with Gly143, Ser144, and Cys145. Similar bonds have been observed in a previous study54–56. The residue Glu189 that formed a hydrogen bond with the docked pose, also formed van der Waals interactions after MDS. However, the remaining two hydrogen bonds with Gly143 and Ser144 were retained. Hydrogen bond interactions with Gly143 and Ser144 were also observed in the X-ray structure; and Cys145 formed van der Waals interactions. The binding modes of olaparib 1826 and olaparib 1885 did not significantly change (Fig. 12b). Rucaparib 184 formed hydrogen bonds with Arg40, His41, and Gln192 after MDS which was different from the molecular docked pose. Among the three compounds, rucaparib underwent drastic changes from its initial structure and was buried deep in the active site (Figs. 11 and 12c). The interactions between Arg40 and His41 were stronger, with an acceptable distance measured below 0.3 nm. Various interactions with Gln192 have been reported earlier57,58. The key residue Met165 had formed π-alkyl interactions with olaparib 1885 and olaparib 1826 while with rucaparib it has formed an alkyl interaction as was seen in the X-ray structure. However, in the molecular docking pose, this residue generated a carbon-hydrogen bond (Supplementary Fig. 4). Several other key residues originating from different subsites59,60 also prompted van der Waals interactions. In particular, Glu166 promoted van der Waals interactions with all the ligands (Table 1). We also calculated the binding energies for the three complexes along with the co-crystallised compound. This was done by initiating the ‘Calculate Binding Energies’ tool available on the DS. This calculation allows for the assessment of the binding energy between the receptor and the ligand. The binding energy for the co-crystallised compound was estimated to be − 136.724 kcal/mol. The binding energy for olaparib 1826, olaparib 1885, and rucaparib 184 were calculated to be − 44.2554 kcal/mol, − 63.7637 kcal/mol, and − 86.8231 kcal/mol, respectively. The higher results for the co-crystallised ligand with respect to the molecular docking score and binding energy may be due to its larger size than the modified compounds. The identified compounds also demonstrated an acceptable Lipinski’s rule of five which was calculated using the Filter by Lipinski’s and Veber. These compounds showed acceptable results (Table 2). Since the compounds were novel, we estimated their ability to be synthesised by adapting the SWISS ADME61 web tool and read it according to the synthetic accessibility score. The results ranged from 1(very easy) to 10 (very difficult)61. The synthetic accessibility score for olaparib 1826, olaparib 1885, and rucaparib 184 were 3.55, 3.63, and 4.30, respectively. Scores close to 5, indicated that these compounds may be easily synthesised (Table 2). These elegant findings show the compounds olaparib 1826, olaparib 1885, and rucaparib 184 are possible Mpro inhibitors.Table 2 Synthetic accessibility score and drug-like properties of the modified compounds. Compound Name SMILES Synthetic accessibility Hydrogen bond Acceptor Hydrogen bond Donor MW ALogP Rotatable bonds Polar surface area (Å2) Olaparib 1826 NCCc1cccc2C(=NNC(=O)c12)Cc3ccc(F)c(c3)C(=O)N4CCN(CC4)C(=O)C5CC5 3.55 8 3 477.531 1.554 6 108.1 Olaparib 1885 Fc1ccc(CC2=NNC(=O)c3c(cccc23)C4CCC4)cc1C(=O)N5CCN(CC5)C(=O)C6CC6 3.63 7 1 488.553 3.406 5 82.08 Rucaparib184 CNCc1ccc(cc1)c2c3CCNC(=O)c4cc(F)cc(c34)n2OC(=O)C[C@H](N)C(=O)[O-] 4.30 9 4 453.443 − 2.051 8 138.51 The novelty of the compounds was verified using SMILES as an input to PubChem62. These results indicate that these compounds were not yet used. We speculated that these designed compounds (Fig. 13) have not yet been synthesised and could be new compounds for the treatment of COVID-19.Figure 13 The 2D structures of the modified compounds in comparison with the parent structures. Materials and methods Selection of the ligands PARP inhibitors were sketched using Biovia Draw v2017 and saved as a molfile (. Mol) format (Fig. 14) and were exported to the DS. The ligands were minimised using the “minimize ligands” protocol available with the DS. The CHARMM force field, was adapted using a Smart Minimizer algorithm. This was executed using 1000 steps of steepest descent with a RMS gradient tolerance of 3; subsequently, conjugate gradient minimisation was applied.Figure 14 2D structures of the PARP inhibitors that are used to design new compounds. Selection of the targets The target selected for this study was Mpro (PDB ID:6LU7) from SARS-CoV-226. The protein was prepared by enabling the ‘Prepare Protein’ protocol available with the DS. This protocol prepares and checks the given protein and performs actions that include modelling the missing loop regions, standardising the atom names, dislodging the water molecules, inserting missing atoms in the incomplete residues63,64, deleting the alternate conformations, and protonating titratable residues using the predicted pKs. The active site for molecular docking was selected around the co-crystallised ligand for all atoms and residues at 13.82 Ǻ. This was done by enabling the tool ‘Define and Edit Binding Site’ that correspondingly creates a binding sphere. It has been reported that the active site contains subsites of origin of different residues65. Prior to the initiation of molecular docking, the co-crystallised ligand N3 was dislodged and redocked into the selected binding pocket to ensure that the binding mode was reproduced. The results have shown that the ligand generated a similar binding mode as that of the co-crystallised ligand with an acceptable RMSD of 0.9 Ǻ (Supplementary Fig. 5). Binding affinity studies Binding affinity studies were performed between Mpro and small molecules using the CDOCKER program available in the DS. Prior to the initiation of docking studies, the binding pockets of the targets were examined to determine the active sites. We aimed to identify a target with a large active site to facilitate the designing the SARS-CoV-2 specific inhibitors (Supplementary Fig. 1). Accordingly, Mpro was selected, and PARP inhibitors were docked into the active site. Correspondingly, the active site was chosen to be around the co-crystallised inhibitor N3. The selected ligands were used to generate 30 conformations. The best pose was selected after clustering the conformers to determine the best binding mode. From the largest cluster, the compound with the best docking score (binding affinity) revealing interactions with the key residues, was chosen. Strategy to design the new compounds To design new compounds computationally, we examined the target ligand complex, after molecular docking, to locate the empty spaces. Subsequently, certain points are identified and marked as R1. The “Grow Scaffold” module available with the DS was enabled which resulted in a total of 3942 compounds. These compounds were docked into Mpro to estimate their binding affinities after a drug-like assessment. The ‘Grow Scaffold’ module permits the user to accomplish reaction-based ligand listing inside the protein’s active site by performing lead optimisation. It begins with the positioning of the scaffold at the receptor-binding site. The user can select a position that can act as a reaction vector, followed by the reactions and reagents to be used. Furthermore, the ligand “novelty” is computed by sorting and ranking the ligands by number of chain assemblies, number of double and aromatic bonds; and N, S, O atom count. Molecular dynamics simulation (MDS) analysis The protein-LIG complex structures obtained from molecular docking were escalated to MDS to understand their binding potency at the binding pocket of the targets. GROningen MAchine for Chemical Simulations (GROMACS) v2016.6 was used to study the MDS66,67. A CHARMM27 all-atom force field was utilized68. The topologies of the ligands were acquired from SwissParam69 and the topol.top file was updated accordingly. Subsequently, the systems were solvated in a dodecahedral water box using the TIP3P water model. The system was neutralised using counterions. Energy minimisation was performed in 50,000 steps to remove bad contacts and clashes. After successful energy minimisation, the protein and ligand were coupled using the gmx make_ndx command to progress through the two-step equilibration process. The first step of equilibration was (constant number of particles, volume, and temperature) NVT equilibration, which was performed for 100 ps using a V-rescale thermostat at 300 K. The second equilibration was conducted with a constant number of particles, pressure, and temperature (NPT) for 100 ps using a Parrinello-Rahman barostat to monitor the pressure at 1 bar. During equilibration, the protein backbone was restrained and the non-protein was permitted to wobble. The long-range electrostatic interactions were evaluated by the Particle Mesh Ewald (PME) method and the short-range interactions and interactions by van der Waals were calculated after applying upper limit of 9 Å and 14 Å, correspondingly. MDS progressed under periodic boundary conditions for 100 ns. The corresponding analysis was conducted using various GROMACS tools and visual molecular dynamics (VMD)42,70. Analysis of the trajectory The different tools available to analyse the MDS results, and GROMACS were utilized71. RMSD was assessed using the gmx rms. These calculations provide knowledge of the deviation present, if any, in the protein from the initial to final conformation during the MDS run. Logically, the lower the RMSD, the greater the stability of the protein72. Rg, which determines the compactness of the protein, was studied using gmx gyrate73. The RMSF of the protein residues was determined using the gmx rmsf tool. Here, the fluctuations in each residue of the protein were examined. The RMSD, Rg and RMSF were computed for the protein backbone. Using the gmx hbond, the number of hydrogen bond interactions between protein and ligand atoms were evaluated during MDS progression. Furthermore, the interaction energy was calculated using the gmx energy, which computes both the short-ranged Coulombic and Lennard–Jones energy interactions45. Conclusion PARP inhibitors are well-known for their anticancer properties. Recently, in vitro studies also demonstrated their anti-COVID-19 properties. In the present study, we aimed to design new PARP inhibitors with high affinity for the COVID-19 Mpro. Correspondingly, two compounds originating from olaparib and one compound originating from rucaparib displayed high binding affinities for Mpro. This indicates their potential use as SARS-CoV-2 inhibitors. These compounds also demonstrated favourable synthetic accessibility scores. Additionally, our method could aid the scientific community to design new compounds and discover new horizons for drug discovery against SARS-CoV-2. Supplementary Information Supplementary Figures. Supplementary Table 1. Supplementary Table 2. Supplementary Information The online version contains supplementary material available at 10.1038/s41598-023-36342-7. Acknowledgements This study was supported by a Korea Research Foundation grant (NRF 2021R1A2B5B02002220). Author contributions S.R. and K.W.L. initiated the idea of the project. S.R. conducted the experiment and analyzed the M.D.S. S.R. and K.W.L. analyzed the results. S.R., T.S.H. and M.W.H. performed the M.D. simulation studies and analysis. S.R wrote the manuscript. All the authors have read and agreed to the final version of the submitted manuscript. Data availability The generated dataset resulted from the ‘Grow Scaffold’ are provided as the Supplementary Table 2. A SMILES representation of the compounds is provided. Competing interests The authors declare no competing interests. Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. ==== Refs References 1. Major C Unprecedented times and innovation Innov. High. Educ. 2020 45 435 436 10.1007/s10755-020-09528-4 33132491 2. Hitchings E Maclean M Unprecedented times: some thoughts on the consequences of the COVID-19 pandemic from a family and social welfare law perspective J. Soc. Welf. Fam. Law 2020 42 277 280 10.1080/09649069.2020.1796219 3. Li H Liu SM Yu XH Tang SL Tang CK Coronavirus disease 2019 (COVID-19): Current status and future perspectives Int. J. Antimicrob. Agents 2020 10.1016/j.ijantimicag.2020.105951 33352235 4. Pushpakom S Drug repurposing: Progress, challenges and recommendations Nat. Rev. Drug Discov. 2018 10.1038/nrd.2018.168 30310233 5. Rudrapal M Khairnar SJ Drug repurposing (DR): An emerging approach in drug discovery Drug Repurposing - Hypothesis, Molecular Aspects and Therapeutic Applications 2020 IntechOpen 6. Dae JW Sangeun J Seungtaek K Yup LS Drugs repurposed for COVID-19 by virtual screening of 6218 drugs and cell-based assay Proc. Natl. Acad. Sci. 2021 118 e2024302118 10.1073/pnas.2024302118 34234012 7. Guy RK DiPaola RS Romanelli F Dutch RE Rapid repurposing of drugs for COVID-19 Science 2020 368 829 830 10.1126/science.abb9332 32385101 8. Gupta RK Nwachuku EL Zusman BE Jha RM Puccio AM Drug repurposing for COVID-19 based on an integrative meta-analysis of SARS-CoV-2 induced gene signature in human airway epithelium PLoS ONE 2021 16 e0257784 10.1371/journal.pone.0257784 34582497 9. Bakowski MA Drug repurposing screens identify chemical entities for the development of COVID-19 interventions Nat. Commun. 2021 12 3309 10.1038/s41467-021-23328-0 34083527 10. Omer SE Drug repurposing for SARS-CoV-2 main protease: Molecular docking and molecular dynamics investigations Biochem. Biophys. Rep. 2022 29 101225 35128086 11. Rampogu S Lee KW Pharmacophore modelling-based drug repurposing approaches for SARS-CoV-2 therapeutics Front. Chem. 2021 9 38 10.3389/fchem.2021.636362 12. Rampogu S Gajula RG Lee G Kim MO Lee KW Unravelling the therapeutic potential of marine drugs as SARS-CoV-2 inhibitors: An insight from essential dynamics and free energy landscape Comput. Biol. Med. 2021 135 104525 10.1016/j.compbiomed.2021.104525 34252682 13. Singh R Bhardwaj VK Sharma J Purohit R Kumar S In-silico evaluation of bioactive compounds from tea as potential SARS-CoV-2 nonstructural protein 16 inhibitors J. Tradit. Complement. Med. 2022 12 35 43 10.1016/j.jtcme.2021.05.005 34099976 14. Aljindan RY Investigation of nonsynonymous mutations in the spike protein of SARS-CoV-2 and its interaction with the ACE2 receptor by molecular docking and MM/GBSA approach Comput. Biol. Med. 2021 135 104654 10.1016/j.compbiomed.2021.104654 34346317 15. Borgio JF State-of-the-art tools unveil potent drug targets amongst clinically approved drugs to inhibit helicase in SARS-CoV-2 Arch. Med. Sci. 2020 16 508 518 10.5114/aoms.2020.94567 32399096 16. Singh R Bhardwaj VK Purohit R Inhibition of nonstructural protein 15 of SARS-CoV-2 by golden spice: A computational insight Cell Biochem. Funct. 2022 40 926 934 10.1002/cbf.3753 36203381 17. Fadlalla M Ahmed M Ali M Elshiekh AA Yousef BA Molecular docking as a potential approach in repurposing drugs against COVID-19: A systematic review and novel pharmacophore models Curr. Pharmacol. Reports 2022 8 3 212 226 10.1007/s40495-022-00285-w 18. Lazniewski M Drug repurposing for identification of potential spike inhibitors for SARS-CoV-2 using molecular docking and molecular dynamics simulations Methods 2022 203 498 510 10.1016/j.ymeth.2022.02.004 35167916 19. Al-Karmalawy AA Molecular docking and dynamics simulation revealed the potential inhibitory activity of aceis against SARS-CoV-2 targeting the hACE2 receptor Front. Chem. 2021 10.3389/fchem.2021.661230 34631658 20. Eweas AF Alhossary AA Abdel-Moneim AS Molecular docking reveals ivermectin and remdesivir as potential repurposed drugs against SARS-CoV-2 Front. Microbiol. 2021 10.3389/fmicb.2020.592908 33746908 21. Rampogu S Lee KW Old drugs for new purpose—fast pace therapeutic identification for SARS-CoV-2 infections by pharmacophore guided drug repositioning approach Bull. Korean Chem. Soc. 2021 10.1002/bkcs.12171 22. Hosseini M Chen W Xiao D Wang C Computational molecular docking and virtual screening revealed promising SARS-CoV-2 drugs Precis. Clin. Med. 2021 4 1 16 10.1093/pcmedi/pbab001 33842834 23. Mohammed AO Abo-Idrees MI Makki AA Ibraheem W Alzain AA Drug repurposing against main protease and RNA-dependent RNA polymerase of SARS-CoV-2 using molecular docking, MM-GBSA calculations and molecular dynamics Struct. Chem. 2022 33 1553 1567 10.1007/s11224-022-01999-9 35789829 24. Matondo A In silico drug repurposing of anticancer drug 5-FU and analogues against SARS-CoV-2 main protease: Molecular docking, molecular dynamics simulation, pharmacokinetics and chemical reactivity studies Adv. Appl. Bioinform. Chem. AABC 2022 15 59 35996620 25. Capoluongo E PARP-inhibitors in a non-oncological indication as COVID-19: Are we aware about its potential role as anti-thrombotic drugs? Discuss. Open. Biomed. Pharmacother. 2020 130 110536 10.1016/j.biopha.2020.110536 32688139 26. Jin Z Structure of Mpro from SARS-CoV-2 and discovery of its inhibitors Nature 2020 582 289 293 10.1038/s41586-020-2223-y 32272481 27. Ullrich S Nitsche C The SARS-CoV-2 main protease as drug target Bioorg. Med. Chem. Lett. 2020 30 127377 10.1016/j.bmcl.2020.127377 32738988 28. Zhang L Crystal structure of SARS-CoV-2 main protease provides a basis for design of improved a-ketoamide inhibitors Science 2020 10.1126/science.abb3405 33335060 29. Hu Q The SARS-CoV-2 main protease (Mpro): Structure, function, and emerging therapies for COVID-19 MedComm 2022 3 e151 10.1002/mco2.151 35845352 30. Mengist HM Dilnessa T Jin T Structural basis of potential inhibitors targeting SARS-CoV-2 main protease Front. Chem. 2021 9 622898 10.3389/fchem.2021.622898 33889562 31. Katre SG Review on development of potential inhibitors of SARS-CoV-2 main protease (MPro) Future J. Pharm. Sci. 2022 8 36 10.1186/s43094-022-00423-7 32. Citarella A Scala A Piperno A Micale N SARS-CoV-2 M(pro): A potential target for peptidomimetics and small-molecule inhibitors Biomolecules 2021 11 607 10.3390/biom11040607 33921886 33. Kanhed AM Identification of potential Mpro inhibitors for the treatment of COVID-19 by using systematic virtual screening approach Mol. Divers. 2020 10.1007/s11030-020-10130-1 32737681 34. Cui W Yang K Yang H recent progress in the drug development targeting SARS-CoV-2 main protease as treatment for COVID-19 Front. Mol. Biosci. 2020 10.3389/fmolb.2020.616341 33385011 35. Wang Y PARP inhibitors in gastric cancer: Beacon of hope J. Exp. Clin. Cancer Res. 2021 40 211 10.1186/s13046-021-02005-6 34167572 36. Chen A PARP inhibitors: its role in treatment of cancer Chin. J. Cancer 2011 30 463 471 10.5732/cjc.011.10111 21718592 37. Du Y Yamaguchi H Hsu JL Hung M-C PARP inhibitors as precision medicine for cancer treatment Natl. Sci. Rev. 2017 4 576 592 10.1093/nsr/nwx027 38. Curtin N Repositioning PARP inhibitors for SARS-CoV-2 infection(COVID-19); A new multi-pronged therapy for acute respiratory distress syndrome? Br. J. Pharmacol. 2020 177 3635 3645 10.1111/bph.15137 32441764 39. Ge Y An integrative drug repositioning framework discovered a potential therapeutic agent targeting COVID-19 Signal Transduct. Target. Ther. 2021 6 165 10.1038/s41392-021-00568-6 33895786 40. Stone NE Stenoparib, an inhibitor of cellular Poly(ADP-Ribose) polymerase, blocks replication of the SARS-CoV-2 and HCoV-NL63 human coronaviruses in vitro MBio 2021 12 e03495-20 10.1128/mBio.03495-20 33468703 41. Badawy AA-B Immunotherapy of COVID-19 with poly (ADP-ribose) polymerase inhibitors: Starting with nicotinamide Biosci. Rep. 2020 40 BSR20202856 10.1042/BSR20202856 33063092 42. Rampogu S Explicit molecular dynamics simulation studies to discover novel natural compound analogues as Mycobacterium tuberculosis inhibitors Heliyon 2023 9 e13324 10.1016/j.heliyon.2023.e13324 36816262 43. Galzitskaya OV Bogatyreva NS Ivankov DN Compactness determines protein folding type J. Bioinform. Comput. Biol. 2008 6 667 680 10.1142/S0219720008003618 18763735 44. Muhammad S Exploring the inhibitory potential of novel bioactive compounds from mangrove actinomycetes against nsp10 the major activator of SARS-CoV-2 replication Chem. Pap. 2022 76 3051 3064 10.1007/s11696-021-01997-x 45. Bepari A Reza H Identification of a novel inhibitor of SARS-CoV-2 3CL-PRO through virtual screening and molecular dynamics simulation PeerJ 2021 9 e11261 10.7717/peerj.11261 33954055 46. Rose M Burgess JT O’Byrne K Richard DJ Bolderson E PARP inhibitors: clinical relevance, mechanisms of action and tumor resistance Front. Cell Dev. Biol. 2020 8 564601 10.3389/fcell.2020.564601 33015058 47. Morales J Review of poly (ADP-ribose) polymerase (PARP) mechanisms of action and rationale for targeting in cancer and other diseases Crit. Rev. Eukaryot. Gene Expr. 2014 24 15 28 10.1615/CritRevEukaryotGeneExpr.2013006875 24579667 48. Chan CY Tan KV Cornelissen B PARP inhibitors in cancer diagnosis and therapy PARP imaging and therapy Clin. Cancer Res. 2021 27 1585 1594 10.1158/1078-0432.CCR-20-2766 33082213 49. Huff S Discovery and mechanism of SARS-CoV-2 main protease inhibitors J. Med. Chem. 2022 65 2866 2879 10.1021/acs.jmedchem.1c00566 34570513 50. Hattori S A small molecule compound with an indole moiety inhibits the main protease of SARS-CoV-2 and blocks virus replication Nat. Commun. 2021 12 668 10.1038/s41467-021-20900-6 33510133 51. Mahmud S Plant-derived compounds effectively inhibit the main protease of SARS-CoV-2: An in silico approach PLoS ONE 2022 17 e0273341 10.1371/journal.pone.0273341 35998194 52. Bahun M Inhibition of the SARS-CoV-2 3CLpro main protease by plant polyphenols Food Chem. 2022 373 131594 10.1016/j.foodchem.2021.131594 34838409 53. Ferreira JC Fadl S Villanueva AJ Rabeh WM Catalytic dyad residues His41 and Cys145 impact the catalytic activity and overall conformational fold of the main SARS-CoV-2 protease 3-chymotrypsin-like protease Front. Chem. 2021 9 491 10.3389/fchem.2021.692168 54. Sacco MD Structure and inhibition of the SARS-CoV-2 main protease reveal strategy for developing dual inhibitors against Mpro and cathepsin L Sci. Adv. 2020 10.1126/sciadv.abe0751 33158912 55. Su H Identification of pyrogallol as a warhead in design of covalent inhibitors for the SARS-CoV-2 3CL protease Nat. Commun. 2021 12 3623 10.1038/s41467-021-23751-3 34131140 56. Singh R Benchmarking the ability of novel compounds to inhibit SARS-CoV-2 main protease using steered molecular dynamics simulations Comput. Biol. Med. 2022 146 105572 10.1016/j.compbiomed.2022.105572 35551011 57. Suárez D Díaz N SARS-CoV-2 main protease: A molecular dynamics study J. Chem. Inf. Model. 2020 60 5815 5831 10.1021/acs.jcim.0c00575 32678588 58. Verma S Patel CN Chandra M Identification of novel inhibitors of SARS-CoV-2 main protease (M(pro) ) from Withania sp. by molecular docking and molecular dynamics simulation J. Comput. Chem. 2021 42 1861 1872 10.1002/jcc.26717 34287986 59. Rampogu S Lee KW Old drugs for new purpose—Fast pace therapeutic identification for SARS-CoV-2 infections by pharmacophore guided drug repositioning approach Bull. Korean Chem. Soc. 2021 42 2 212 226 10.1002/bkcs.12171 60. Rampogu S Gajula RG Lee G Kim MO Lee KW Unravelling the therapeutic potential of marine drugs as SARS-CoV-2 inhibitors: An insight from essential dynamics and free energy landscape Comput. Biol. Med. 2021 10.1016/j.compbiomed.2021.104525 34252682 61. Daina A Michielin O Zoete V SwissADME: A free web tool to evaluate pharmacokinetics, drug-likeness and medicinal chemistry friendliness of small molecules Sci. Rep. 2017 10.1038/srep42717 28256516 62. Kim S PubChem substance and compound databases Nucleic Acids Res. 2016 10.1093/nar/gkv951 27907895 63. Spassov VZ Flook PK Yan L LOOPER: A molecular mechanics-based algorithm for protein loop prediction Protein Eng. Des. Sel. 2008 21 91 100 10.1093/protein/gzm083 18194981 64. Spassov VZ Yan L A fast and accurate computational approach to protein ionization Protein Sci. 2008 17 1955 1970 10.1110/ps.036335.108 18714088 65. Kneller DW Structural plasticity of SARS-CoV-2 3CL Mpro active site cavity revealed by room temperature X-ray crystallography Nat. Commun. 2020 11 3202 10.1038/s41467-020-16954-7 32581217 66. Van Der Spoel D GROMACS: Fast, flexible, and free J. Comput. Chem. 2005 26 1701 1718 10.1002/jcc.20291 16211538 67. Lemkul J From proteins to perturbed hamiltonians: A suite of tutorials for the GROMACS-2018 molecular simulation package [Article v1.0] Living J. Comput. Mol. Sci. 2018 10.33011/livecoms.1.1.5068 68. Bjelkmar P Larsson P Cuendet MA Hess B Lindahl E Implementation of the CHARMM force field in GROMACS: Analysis of protein stability effects from correction maps, virtual interaction sites, and water models J. Chem. Theory Comput. 2010 6 459 466 10.1021/ct900549r 26617301 69. Zoete V Cuendet MA Grosdidier A Michielin O SwissParam: A fast force field generation tool for smallorganic molecules J. Comput. Chem. 2011 32 2359 2368 10.1002/jcc.21816 21541964 70. Humphrey W Dalke A Schulten KVMD Visual molecular dynamics J. Mol. Graph. 1996 14 33 38 10.1016/0263-7855(96)00018-5 8744570 71. Kushwaha PP Identification of natural inhibitors against SARS-CoV-2 drugable targets using molecular docking, molecular dynamics simulation, and MM-PBSA approach Front. Cell. Infect. Microbiol. 2021 10.3389/fcimb.2021.730288 34458164 72. Sharma J An in-silico evaluation of different bioactive molecules of tea for their inhibition potency against non structural protein-15 of SARS-CoV-2 Food Chem. 2021 346 128933 10.1016/j.foodchem.2020.128933 33418408 73. Lobanov MY Bogatyreva NS Galzitskaya OV Radius of gyration as an indicator of protein structure compactness Mol. Biol. 2008 10.1134/S0026893308040195