Molecular dynamics simulation, 3D-pharmacophore and scaffold hopping analysis in the design of multi-target drugs to inhibit potential targets of COVID-19

• Published in: Journal of biomolecular structure & dynamics 2022-01, Vol.40 (22), p.11787-11808
• Main Authors: Fayyazi, Neda, Mostashari-Rad, Tahereh, Ghasemi, Jahan B., Ardakani, Mehran Mirabzadeh, Kobarfard, Farzad
• Format: Article
• Language: English
• Subjects: ADMET prediction; COVID-19; Humans; MD simulations; Molecular Docking Simulation; Molecular Dynamics Simulation; Pharmacophore; Pharmacophore modeling; Quantitative Structure-Activity Relationship; repurposing approach; RNA-Dependent RNA Polymerase; SARS-CoV-2; scaffold hopping
• ISSN: 0739-1102; 1538-0254
• DOI: 10.1080/07391102.2021.1965914

SARS-CoV-2 has posed serious threat to the health and has inflicted huge costs in the world. Discovering potent compounds is a critical step to inhibit coronavirus. 3CL pro and RdRp are the most conserved targets associated with COVID-19. In this study, three-dimensional pharmacophore modeling, scaffold hopping, molecular docking, structure-based virtual screening, QSAR-based ADMET predictions and molecular dynamics analysis were used to identify inhibitors for these targets. Binding free energies estimated by molecular docking for each ligand in different binding sites of RdRp were used to predict the active site. Previously reported active 3CL pro and RdRp inhibitors were used to build a pharmacophore model to develop different scaffolds. Structure-based simulations and pharmacophore modeling based on Hip Hop algorithm converged in a state that suggest hydrogen bond acceptor and donor features have a critical role in the two binding sites. Further validations indicated that the best pharmacophore model has fairly good correlation values compared with approved inhibitors. Structure-based simulation results approved that GLu166 and Gln189 in 3CL pro and Lys551 and Glu811 in RdRp, are critical residues for dual activities. Ten compounds were extracted from pharmacophore-based virtual screening in six databases. The results, gained by repurposing approach, suggest the effectiveness of these ten compounds with different scaffolds as possible inhibitors of the two targets. Some quinoline-based hybrid derivatives also were designed. QSAR descriptors plot predicted that the scaffolds have had accepted pharmacokinetic profiles. Multiple molecular dynamics simulations in 100 ns and MM/PBSA studies of some reference inhibitors and the novel compounds in complex with both targets demonstrated stable complexes and confirmed the interaction modes. Based on different computational methods, COVID-19 multi-target inhibitors are proposed. Communicated by Ramaswamy H. Sarma Highlights Novel multitarget compounds against COVID-19 is predicted . Pharmacophore model is built based on approved inhibitors of RdRPp and 3CL pro . Structure based molecular modeling is used in six large databases to scaffold hopping. Replicated molecular dynamics simulation, MM/PBSA and Pharmacokinetic profile are analyzed to new scaffold inhibitors against COVID-19.