Theoretical modeling and design of some pyrazolopyrimidine derivatives as Wolbachia inhibitors, targeting lymphatic filariasis and onchocerciasis

Virtual screening, MMGBSA, and molecular dynamics approaches for identification of natural products from South African biodiversity as potential Onchocerca volvulus pi-class glutathione S-transferase inhibitors

We investigated 1012 molecules from natural products previously isolated from the South African biodiversity (SANCDB, https://sancdb.rubi.ru.ac.za/), for putative inhibition of Onchocerca volvulus pi-class glutathione S-transferase (Ov-GST2) by virtual screening, MMGBSA, and molecular dynamics approaches. ADMET, docking, and MMGBSA shortlisted 12 selected homoisoflavanones-type hit molecules, among which two namely SANC00569, and SANC00689 displayed high binding affinities of -46.09 and -46.26 kcal mol-1, respectively towards π-class Ov-GST2, respectively. The molecular dynamics results of SANC00569 showed the presence of intermolecular H-bonding, hydrophobic interactions between the ligand and key amino acids of Ov-GST2, throughout the simulation period. This hit molecule had a stable binding pose and occupied the binding pockets throughout the 200 ns simulation. To the best of our knowledge, there is no report of any alleged anti-onchocerciasis activity referring to homoisoflavanones or flavonoids. Nevertheless, homoisoflavanones, which are a subclass of flavonoids, exhibit a plethora of biological activities. All these results led to the conclusion that SANC00569 is the most hypothetical Ov-GST2, which could lead the development of new drugs against Onchocerca volvulus pi-class glutathione S-transferase. Further validation of these findings through in vitro and in vivo studies is required.

Computational design, molecular properties, ADME, and toxicological analysis of substituted 2,6-diarylidene cyclohexanone analogs as potent pyridoxal kinase inhibitors

Leishmaniasis is one of the tropical diseases which affects over 12 million people mainly in the tropical regions of the world and is caused by the leishmanial parasites transmitted by the female sand fly. The lack of vaccines to prevent leishmaniasis, as well as limitations of existing therapies necessitated this study which was focused on a combined virtual docking screening and 3-D QSAR modeling approach to design some diarylidene cyclohexanone analogs, while also performing pharmacokinetic analysis and Molecular Dynamic (MD) simulation to ascertain their drug-ability. As a result, the built 3-D QSAR model was found to satisfy the requirement of a good model with R2 = 0.9777, SDEC = 0.0593, F-test = 105.028, and Q2 LOO = 0.6592. The template (compound 9, MolDock score =  - 161.064) and all seven newly designed analogs were found to possess higher docking scores than the reference drug (Pentamidine, Moldock score = - 137.827). The results of the pharmacokinetic analysis suggest 9 and the new molecules (9a, b, c, e, and f) as orally bioavailable with good ADME and safe toxicological profiles. These molecules also showed good binding interactions with the receptor (pyridoxal kinase). Additionally, the MD simulation result confirmed the stability of the tested protein-ligand complexes, with an estimated ∆G binding (MM/GBSA) of - 65.2177 kcal/mol and - 58.433 kcal/mol for 9_6K91 and 9a_6K91 respectively. Hence, the new compounds, especially 9a could be considered potential anti-leishmanial inhibitors.

Molecular docking-based virtual screening, molecular dynamic simulation, and 3-D QSAR modeling of some pyrazolopyrimidine analogs as potent anti-filarial agents

Lymphatic filariasis and onchocerciasis are common filarial diseases caused by filarial worms, which co-habit symbiotically with the Wolbachia organism. One good treatment method seeks Wolbachia as a drug target. Here, a computer-aided molecular docking screening and 3-D QSAR modeling were conducted on a series of Fifty-two (52) pyrazolopyrimidine derivatives against four Wolbachia receptors, including a pharmacokinetics study and Molecular Dynamic (MD) investigation, to find a more potent anti-filarial drug. The DFT approach (B3LYP with 6-31G** option) was used for the structural optimization. Five ligand-protein interaction pairs with the highest binding affinities were identified in the order; 23_7ESX (-10.2 kcal/mol) > 14_6EEZ (- 9.0) > 29_3F4R (- 8.0) > 26_6W9O (- 7.7) ≈ doxycycline_7ESX (- 7.7), with good pharmacological interaction profiles. The built 3-D QSAR model satisfied the requirement of a good model with R2 = 0.9425, Q2 LOO = 0.5019, SDEC = 0.1446, and F test = 98.282. The selected molecules (14, 23, 26, and 29) perfectly obeyed Lipinski's RO5 for oral bio-availability, and showed excellent ADMET properties, except 14 with positive AMES toxicity. The result of the MD simulation showed the great stability associated with the binding of 23 onto 7ESX's binding pocket with an estimated binding free energy (MM/GBSA) of - 60.6552 kcal/mol. Therefore, 23 could be recommended as a potential anti-filarial drug molecule, and/or template for the design of more prominent inhibitors.

Supplementary Information

The online version contains supplementary material available at 10.1007/s40203-022-00136-y.