Repurposing of FDA‐Approved Drugs for Treating Iatrogenic Botulism: A Paired 3D‐QSAR/Docking Approach†

Botulinum neurotoxin (BoNT) is widely used for the treatment of spasticity, focal dystonia, chronic migraine, facial hemispasm, and facial aesthetic treatments. Generally, treatment with botulinum toxin is a safe procedure when conducted by clinicians with expertise, and local side effects are rare and transient. However, occasionally adverse effects can occur due to the spread of the drug to nontargeted muscles and organs, producing dry mouth, fatigue, and flu‐like symptoms, up to signs of systemic botulism, which appears to be more frequent in children treated for spasticity than in adults. In silico 3D‐QSAR and molecular docking studies were performed to build a structure‐based model on selected potent known botulinum neurotoxin type A inhibitors; this was used to screen the US Food and Drug Administration (FDA) database. Thirty molecules were identified as possible light‐chain BoNT/A inhibitors. In this study, we applied a well‐established ligand‐ and structure‐based methodology for the identification of hit compounds among a database of FDA‐approved drugs. The identification of budesonide, protirelin, and ciclesonide followed by other compounds can be considered a starting point for investigations of selected compounds that could bypass much of the time and costs involved in the drug approval process.

[1]  J. Sobel,et al.  Making the Best of the Evidence: Toward National Clinical Guidelines for Botulism. , 2018, Clinical infectious diseases : an official publication of the Infectious Diseases Society of America.

[2]  J. Gracies,et al.  Effects of repeated abobotulinumtoxinA injections in upper limb spasticity , 2017, Muscle & nerve.

[3]  V. Abbate,et al.  3D-QSAR assisted identification of FABP4 inhibitors: An effective scaffold hopping analysis/QSAR evaluation. , 2019, Bioorganic chemistry.

[4]  S. D. Jong SIMPLS: an alternative approach to partial least squares regression , 1993 .

[5]  Tim Cheeseright,et al.  Molecular field technology applied to virtual screening and finding the bioactive conformation , 2007, Expert opinion on drug discovery.

[6]  Kim D Janda,et al.  Identification of a potent botulinum neurotoxin a protease inhibitor using in situ lead identification chemistry. , 2006, Organic letters.

[7]  Gert Vriend,et al.  YASARA View—molecular graphics for all devices—from smartphones to workstations , 2014, Bioinform..

[8]  P. Roy,et al.  Exploring the impact of size of training sets for the development of predictive QSAR models , 2008 .

[9]  B. Jabbari Botulinum Toxin Treatment in Clinical Medicine , 2018, Springer International Publishing.

[10]  V. Sorrenti,et al.  Targeting heme Oxygenase-1 with hybrid compounds to overcome Imatinib resistance in chronic myeloid leukemia cell lines. , 2018, European journal of medicinal chemistry.

[11]  Márcia M. C. Ferreira,et al.  Basic validation procedures for regression models in QSAR and QSPR studies: theory and application , 2009 .

[12]  Andy Vinter,et al.  Molecular Field Extrema as Descriptors of Biological Activity: Definition and Validation , 2006, J. Chem. Inf. Model..

[13]  S. Bavari,et al.  Recent advances in botulinum neurotoxin inhibitor development. , 2014, Current topics in medicinal chemistry.

[14]  Christopher I. Bayly,et al.  Fast, efficient generation of high‐quality atomic charges. AM1‐BCC model: II. Parameterization and validation , 2002, J. Comput. Chem..

[15]  Gert Vriend,et al.  Increasing the precision of comparative models with YASARA NOVA—a self‐parameterizing force field , 2002, Proteins.

[16]  Carmela Parenti,et al.  Development of a Sigma‐2 Receptor affinity filter through a Monte Carlo based QSAR analysis , 2017, European journal of pharmaceutical sciences : official journal of the European Federation for Pharmaceutical Sciences.

[17]  H. Nagashima Studies on the Different Modes of Action of the Anticoagulant Protease Inhibitors DX-9065a and Argatroban , 2002, The Journal of Biological Chemistry.

[18]  S. Wold,et al.  PLS-regression: a basic tool of chemometrics , 2001 .

[19]  Wei Zhang,et al.  A point‐charge force field for molecular mechanics simulations of proteins based on condensed‐phase quantum mechanical calculations , 2003, J. Comput. Chem..

[20]  Carmela Parenti,et al.  Sigma-2 receptor ligands QSAR model dataset , 2017, Data in brief.

[21]  S. Long BabyBIG has BIG advantages for treatment of infant botulism. , 2018, The Journal of pediatrics.

[22]  V. Sorrenti,et al.  Novel Structural Insight into Inhibitors of Heme Oxygenase-1 (HO-1) by New Imidazole-Based Compounds: Biochemical and In Vitro Anticancer Activity Evaluation , 2018, Molecules.

[23]  James J P Stewart,et al.  Optimization of parameters for semiempirical methods IV: extension of MNDO, AM1, and PM3 to more main group elements , 2004, Journal of molecular modeling.

[24]  Carlos Aleman,et al.  Suitability of the PM3‐derived molecular electrostatic potentials , 1993, J. Comput. Chem..

[25]  Giuseppe Floresta,et al.  Hyphenated 3D-QSAR statistical model-scaffold hopping analysis for the identification of potentially potent and selective sigma-2 receptor ligands. , 2017, European journal of medicinal chemistry.

[26]  Anil K. Saxena,et al.  Consensus Superiority of the Pharmacophore-Based Alignment, Over Maximum Common Substructure (MCS): 3D-QSAR Studies on Carbamates as Acetylcholinesterase Inhibitors , 2009, J. Chem. Inf. Model..

[27]  G. Vriend,et al.  Fast empirical pKa prediction by Ewald summation. , 2006, Journal of molecular graphics & modelling.

[28]  R. Turnaturi,et al.  A Structure- and Ligand-Based Virtual Screening of a Database of “Small” Marine Natural Products for the Identification of “Blue” Sigma-2 Receptor Ligands , 2018, Marine drugs.

[29]  V. Abbate,et al.  FABP4 inhibitors 3D-QSAR model and isosteric replacement of BMS309403 datasets , 2018, Data in brief.

[30]  Roland L. Dunbrack,et al.  Assignment of protonation states in proteins and ligands: combining pKa prediction with hydrogen bonding network optimization. , 2012, Methods in molecular biology.

[31]  P. Wipf,et al.  Pharmacophore-guided lead optimization: the rational design of a non-zinc coordinating, sub-micromolar inhibitor of the botulinum neurotoxin serotype a metalloprotease. , 2009, Bioorganic & medicinal chemistry letters.

[32]  F. Galvano,et al.  Musician’s dystonia , 2015, Practical Neurology.

[33]  L. Salerno,et al.  Identification of Potentially Potent Heme Oxygenase 1 Inhibitors through 3D‐QSAR Coupled to Scaffold‐Hopping Analysis , 2018, ChemMedChem.

[34]  V. Sorrenti,et al.  Potholing of the hydrophobic heme oxygenase-1 western region for the search of potent and selective imidazole-based inhibitors. , 2018, European journal of medicinal chemistry.