Biomolecules (Nov 2023)

Pharmacophore-Based Screening, Molecular Docking, and Dynamic Simulation of Fungal Metabolites as Inhibitors of Multi-Targets in Neurodegenerative Disorders

  • Danish Iqbal,
  • Mohammed Alsaweed,
  • Qazi Mohammad Sajid Jamal,
  • Mohammad Rehan Asad,
  • Syed Mohd Danish Rizvi,
  • Moattar Raza Rizvi,
  • Hind Muteb Albadrani,
  • Munerah Hamed,
  • Sadaf Jahan,
  • Hadeel Alyenbaawi

DOI
https://doi.org/10.3390/biom13111613
Journal volume & issue
Vol. 13, no. 11
p. 1613

Abstract

Read online

Neurodegenerative disorders, such as Alzheimer’s disease (AD), negatively affect the economic and psychological system. For AD, there is still a lack of disease-altering treatments and promising cures due to its complex pathophysiology. In this study, we computationally screened the natural database of fungal metabolites against three known therapeutic target proteins of AD. Initially, a pharmacophore-based, drug-likeness category was employed for screening, and it filtered the 14 (A–N) best hits out of 17,544 fungal metabolites. The 14 best hits were docked individually against GSK-3β, the NMDA receptor, and BACE-1 to investigate the potential of finding a multitarget inhibitor. We found that compounds B, F, and L were immuno-toxic, whereas E, H, I, and J had a higher LD50 dose (5000 mg/kg). Among the examined metabolites, the Bisacremine-C (compound I) was found to be the most active molecule against GSK-3β (ΔG: −8.7 ± 0.2 Kcal/mol, Ki: 2.4 × 106 M−1), NMDA (ΔG: −9.5 ± 0.1 Kcal/mol, Ki: 9.2 × 106 M−1), and BACE-1 (ΔG: −9.1 ± 0.2 Kcal/mol, Ki: 4.7 × 106 M−1). It showed a 25-fold higher affinity with GSK-3β, 6.3-fold higher affinity with NMDA, and 9.04-fold higher affinity with BACE-1 than their native ligands, respectively. Molecular dynamic simulation parameters, such as RMSD, RMSF, Rg, and SASA, all confirmed that the overall structures of the targeted enzymes did not change significantly after binding with Bisacremine-C, and the ligand remained inside the binding cavity in a stable conformation for most of the simulation time. The most significant hydrophobic contacts for the GSK-3β-Bisacremine-C complex are with ILE62, VAL70, ALA83, and LEU188, whereas GLN185 is significant for H-bonds. In terms of hydrophobic contacts, TYR184 and PHE246 are the most important, while SER180 is vital for H-bonds in NMDA-Bisacremine-C. THR232 is the most crucial for H-bonds in BACE-1-Bisacremine-C and ILE110-produced hydrophobic contacts. This study laid a foundation for further experimental validation and clinical trials regarding the biopotency of Bisacremine-C.

Keywords