1
|
Oduselu GO, Afolabi R, Ademuwagun I, Vaughan A, Adebiyi E. Structure-based pharmacophore modeling, virtual screening, and molecular dynamics simulation studies for identification of Plasmodium falciparum 5-aminolevulinate synthase inhibitors. Front Med (Lausanne) 2023; 9:1022429. [PMID: 36714108 PMCID: PMC9877529 DOI: 10.3389/fmed.2022.1022429] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/18/2022] [Accepted: 12/23/2022] [Indexed: 01/13/2023] Open
Abstract
Plasmodium falciparum (Pf) 5-aminolevulinic acid synthase (5-ALAS) is an essential enzyme with high selectivity during liver stage development, signifying its potential as a prophylactic antimalarial drug target. The aim of this study was to identify important potential lead compounds which can serve as inhibitors of Pf 5-ALAS using pharmacophore modeling, virtual screening, qualitative structural assessment, in silico ADMET (Absorption, Distribution, Metabolism, Excretion and Toxicity) evaluation and molecular dynamics simulation. The best model of the tertiary structure of Pf 5-ALAS was obtained using MolProbity, while the following databases were explored for the pharmacophore-based virtual screening: CHEMBL, ChemDiv, ChemSpace, MCULE, MCULE-ULTIMATE, MolPort, NCI Open Chemical Repository, LabNetwork and ZINC databases. 2,621 compounds were screened against the modeled Pf 5-ALAS using AutoDock vina. The post-screening analysis was carried out using Discovery Studio while molecular dynamics simulation was performed on the best hits using NAMD-VMD and Galaxy Europe platform. Compound CSMS00081585868 was observed as the best hit with a binding affinity of -9.9 kcal/mol and predicted Ki of 52.10 nM, engaging in seven hydrogen bonds with the target's active site amino acid residues. The in silico ADMET prediction showed that all ten best hits possessed relatively good pharmacokinetic properties. The qualitative structural assessment of the best hit, CSMS00081585868, revealed that the presence of two pyridine scaffolds bearing hydroxy and fluorine groups linked by a pyrrolidine scaffold contributed significantly to its ability to have a strong binding affinity with the receptor. The best hit also showed stability in the active site of Pf 5-ALAS as confirmed from the RMSD obtained during the MD simulation.
Collapse
Affiliation(s)
- Gbolahan O. Oduselu
- Covenant University Bioinformatics Research (CUBRe), Covenant University, Ota, Ogun State, Nigeria
- Department of Chemistry, Covenant University, Ota, Ogun State, Nigeria
| | - Rufus Afolabi
- Covenant University Bioinformatics Research (CUBRe), Covenant University, Ota, Ogun State, Nigeria
- Department of Biochemistry, Covenant University, Ota, Ogun State, Nigeria
| | - Ibitayo Ademuwagun
- Covenant University Bioinformatics Research (CUBRe), Covenant University, Ota, Ogun State, Nigeria
- Department of Biochemistry, Covenant University, Ota, Ogun State, Nigeria
| | - Ashley Vaughan
- Center for Global Infectious Disease Research, Seattle Children’s Research Institute, Seattle, WA, United States
| | - Ezekiel Adebiyi
- Covenant University Bioinformatics Research (CUBRe), Covenant University, Ota, Ogun State, Nigeria
- Department of Computer and Information Science, Covenant University, Ota, Ogun State, Nigeria
- Covenant Applied Informatics and Communications ACE (CApIC-ACE), Covenant University, Ota, Ogun State, Nigeria
- Division of Applied Bioinformatics, German Cancer Research Center (DKFZ), Heidelberg, Germany
| |
Collapse
|