Ligand-based pharmacophore modeling, virtual screening, and molecular dynamics simulations of Pfhsp90 fingerprints in Plasmodium malaria treatment

Abstract

The World Health Organization (WHO) documents malaria as one of the leading causes of high morbidity and mortality worldwide. The disease affects millions and kills thousands of people annually. Efforts to reduce the global burden of malaria have prompted the WHO to recommend prevention strategies such as using antimalarial drugs. However, these strategies have been ineffective because of antimalarial drug resistance. The only efficacious malaria treatment is Artemisinin-based Combination Therapy (ACT). However, the extended ACT clearance times, linked to the emergence of artemisinin monotherapy resistance recorded most recently in Africa and the Great Mekong region, pose a danger to its efficacy. Therefore, better efficacious antimalarial drugs are required. Since P. falciparum heat shock protein 90 (PfHsp90) is a well-characterized malaria drug target, this study uses it to discover more efficacious antimalarial drugs. An in silico approach was used to discover PfHsp90 inhibitors with pharmacological properties against Plasmodium malaria using a molecular dynamics simulation (MDS) and hierarchical virtual screening. Geldanamycin (GDM), a well-known anti-PfHsp90 compound, was used to identify PfHsp90 inhibitors with pharmacological properties against Plasmodium malaria by screening it against the ZINC20 database via the ZINCPHARMER web server. This virtual screening process resulted in 17 hits. These ZINCPHARMER hits were subjected to drug-likeness and pharmacokinetics properties analysis in the SwissADME web server, and nine of them satisfied the requirements. The nine ZINC20 compounds were docked with PfHsp90 using the PyRx software version 0.8 to understand their interactions. From the molecular docking results, ZINC09060002 (−8.2 kcal/mol), ZINC72133064 (−7.8 kcal/mol), ZINC72163401 (−7.7 kcal/mol), ZINC72358537 (−8.1 kcal/mol), and ZINC72358557 (−7.6 kcal/mol) had better binding affinities to PfHsp90 than GDM (−7.5 kcal/mol). The stability of these molecularly docked protein–inhibitor complexes was assessed through MDS using GROMACS 2022. ZINC72163401, ZINC72358537, and ZINC72358557 formed stable complexes with PfHsp90. The lead compounds were subjected to in vitro validation for their inhibitory capability. They showed promising inhibition of parasite growth with IC50 values ranging between 200 and 400 ng/mL. In this regard, the three PfHsp90 inhibitors can be further developed as potential antimalarial drugs. However, further structural optimization studies and clinical (in vivo) tests are necessary to ascertain the antimalarial activity of these compounds in humans.