An integrated pharmacophore/docking/3D-QSAR approach to screening a large library of products in search of future botulinum neurotoxin a inhibitors

Botulinum toxins are neurotoxins produced by Clostridium botulinum. This toxin can be lethal for humans as a cause of botulism; however, in small doses, the same toxin is used to treat different conditions. Even if the therapeutic doses are effective and safe, the adverse reactions could be local and could unmask a subclinical impairment of neuromuscular transmissions. There are not many cases of adverse events in the literature; however, it is possible that sometimes they do not occur as they are transient and, if they do occur, there is no possibility of a cure other than to wait for the pharmacological effect to end. Inhibition of botulinum neurotoxin type A (BoNT/A) effects is a strategy for treating botulism as it can provide an effective post-exposure remedy. In this paper, 13,592,287 compounds were screened through a pharmacophore filter, a 3D-QSAR model, and a virtual screening; then, the compounds with the best affinity were selected. Molecular dynamics simulation studies on the first four compounds predicted to be the most active were conducted to verify that the poses foreseen by the docking were stable. This approach allowed us to identify compounds with a calculated inhibitory activity in the range of 316-500 nM.