Benchmarked molecular docking integrated molecular dynamics stability analysis for prediction of SARS-CoV-2 papain-like protease inhibition by olive secoiridoids.
J King Saud Univ Sci
; 35(1): 102402, 2023 Jan.
Article
em En
| MEDLINE
| ID: mdl-36338939
Objectives: We performed a virtual screening of olive secoiridoids of the OliveNetTM library to predict SARS-CoV-2 PLpro inhibition. Benchmarked molecular docking protocol that evaluated the performance of two docking programs was applied to execute virtual screening. Molecular dynamics stability analysis of the top-ranked olive secoiridoid docked to PLpro was also carried out. Methods: Benchmarking virtual screening used two freely available docking programs, AutoDock Vina 1.1.2. and AutoDock 4.2.1. for molecular docking of olive secoiridoids to a single PLpro structure. Screening also included benchmark structures of known active and decoy molecules from the DEKOIS 2.0 library. Based on the predicted binding energies, the docking programs ranked the screened molecules. We applied the usual performance evaluation metrices to evaluate the docking programs using the predicted ranks. Molecular dynamics of the top-ranked olive secoiridoid bound to PLpro and computation of MM-GBSA energy using three iterations during the last 50 ps of the analysis of the dynamics in Desmond supported the stability prediction. Results and discussions: Predictiveness curves suggested that AutoDock Vina has a better predictive ability than AutoDock, although there was a moderate correlation between the active molecules rankings (Kendall's correlation of rank (τ) = 0.581). Interestingly, two same molecules, Demethyloleuropein aglycone, and Oleuroside enriched the top 1 % ranked olive secoiridoids predicted by both programs. Demethyloleuropein aglycone bound to PLpro obtained by docking in AutoDock Vina when analyzed for stability by molecular dynamics simulation for 50 ns displayed an RMSD, RMSF<2 Å, and MM-GBSA energy of -94.54 ± 6.05 kcal/mol indicating good stability. Molecular dynamics also revealed the interactions of Demethyloleuropein aglycone with binding sites 2 and 3 of PLpro, suggesting a potent inhibition. In addition, for 98 % of the simulation time, two phenolic hydroxy groups of Demethyloleuropein aglycone maintained two hydrogen bonds with Asp302 of PLpro, specifying the significance of the groups in receptor binding. Conclusion: AutoDock Vina retrieved the active molecules accurately and predicted Demethyloleuropein aglycone as the best inhibitor of PLpro. The Arabian diet consisting of olive products rich in secoiridoids benefits from the PLpro inhibition property and reduces the risk of viral infection.
AD, AutoDock 4.2.1; ADV, AutoDock Vina 1.1.2; BEDROC, Boltzmann enhanced discrimination of ROC; Benchmarking docking; DEKOIS, Demanding evaluation kits for objective in-silico screening; EF, Enrichment factor; M, Moles; MD, Molecular dynamics; MM-GBSA, Molecular mechanics generalized Born surface area; MW, Molecular weight; Molecular docking; Molecular dynamics; OS, Olive secoiridoids; Olive secoiridoids; PC, Predictiveness curve; PLpro; PLpro, Papain-like protease; RIE, Robust initial enhancement; RMSD, Root mean square deviation; RMSF, Root mean square fluctuation; ROC, Receiver operating characteristic curve; ROC-AUC, Area under ROC; SARS-CoV-2; SARS-CoV-2, Severe acute respiratory syndrome coronavirus-2; TG, Total gain; g/mol, Grams/mole; kcal/mol, Kilocalorie/mole; ns, nanoseconds; pAUC, partial area under ROC; pTG, Partial total gain; ps, picoseconds
Texto completo:
1
Coleções:
01-internacional
Base de dados:
MEDLINE
Tipo de estudo:
Prognostic_studies
/
Risk_factors_studies
Idioma:
En
Ano de publicação:
2023
Tipo de documento:
Article