GRIND2-based 3D-QSAR and prediction of activity spectra for symmetrical bis-pyridinium salts with promastigote antileishmanial activity

Leishmaniasis is a major group of neglected tropical diseases caused by the protozoan parasite Leishmania. About 12 million people are affected in 98 countries and 350 million people worldwide are at risk of infection. Current leishmaniasis treatments rely on a relatively small arsenal of drugs, including amphotericin B, pentamidine and others, which in general have some type of inconvenience. Recently, we have synthesized antileishmanial bis-pyridinium derivatives and symmetrical bis-pyridinium cyclophanes. These compounds are considered structural analogues of pentamidine, where the amidino moiety, protonated at physiological pH, is replaced by a positively charged nitrogen atom as a pyridinium ring. In this work, a statistically significant GRIND2-based 3D-QSAR model was built and biological activity predictions were in silico carried out allowing rationalization of the different activities recently obtained against Leishmania donovani (in L. donovani promastigotes) for a data set of 19 bis-pyridinium compounds. We will emphasize the most important structural requirements to improve the biological activity and probable interactions with the biological receptor as a guide for lead and prototype optimization. In addition, since no information about the actual biological target for this series of active compounds is provided, we have used Prediction of Activity Spectra for Biologically Active Substances to propose our compounds as potential nicotinic α6β3β4α5 receptor antagonists. This proposal is reinforced by the high structural similarity observed between our compounds and several anthelmintic drugs in current clinical use, which have the same drug action mechanism here predicted. Such new findings would be confirmed with further and additional experimental assays.