Computer-Aided Lead Optimization: Improved Small-Molecule Inhibitor of the Zinc Endopeptidase of Botulinum Neurotoxin Serotype A

Optimization of a serotype-selective, small-molecule inhibitor of botulinum neurotoxin serotype A (BoNTA) endopeptidase is a formidable challenge because the enzyme-substrate interface is unusually large and the endopeptidase itself is a large, zinc-binding protein with a complex fold that is difficult to simulate computationally. We conducted multiple molecular dynamics simulations of the endopeptidase in complex with a previously described inhibitor (Kiapp of 7±2.4 µM) using the cationic dummy atom approach. Based on our computational results, we hypothesized that introducing a hydroxyl group to the inhibitor could improve its potency. Synthesis and testing of the hydroxyl-containing analog as a BoNTA endopeptidase inhibitor showed a twofold improvement in inhibitory potency (Kiapp of 3.8±0.8 µM) with a relatively small increase in molecular weight (16 Da). The results offer an improved template for further optimization of BoNTA endopeptidase inhibitors and demonstrate the effectiveness of the cationic dummy atom approach in the design and optimization of zinc protease inhibitors

Potent New Small-Molecule Inhibitor of Botulinum Neurotoxin Serotype A Endopeptidase Developed by Synthesis-Based Computer-Aided Molecular Design

Botulinum neurotoxin serotype A (BoNTA) causes a life-threatening neuroparalytic disease known as botulism. Current treatment for post exposure of BoNTA uses antibodies that are effective in neutralizing the extracellular toxin to prevent further intoxication but generally cannot rescue already intoxicated neurons. Effective small-molecule inhibitors of BoNTA endopeptidase (BoNTAe) are desirable because such inhibitors potentially can neutralize the intracellular BoNTA and offer complementary treatment for botulism. Previously we reported a serotype-selective, small-molecule BoNTAe inhibitor with a Kiapp value of 3.8±0.8 µM. This inhibitor was developed by lead identification using virtual screening followed by computer-aided optimization of a lead with an IC50 value of 100 µM. However, it was difficult to further improve the lead from micromolar to even high nanomolar potency due to the unusually large enzyme-substrate interface of BoNTAe. The enzyme-substrate interface area of 4,840 Å2 for BoNTAe is about four times larger than the typical protein-protein interface area of 750–1,500 Å2. Inhibitors must carry several functional groups to block the unusually large interface of BoNTAe, and syntheses of such inhibitors are therefore time-consuming and expensive. Herein we report the development of a serotype-selective, small-molecule, and competitive inhibitor of BoNTAe with a Ki value of 760±170 nM using synthesis-based computer-aided molecular design (SBCAMD). This new approach accounts the practicality and efficiency of inhibitor synthesis in addition to binding affinity and selectivity. We also report a three-dimensional model of BoNTAe in complex with the new inhibitor and the dynamics of the complex predicted by multiple molecular dynamics simulations, and discuss further structural optimization to achieve better in vivo efficacy in neutralizing BoNTA than those of our early micromolar leads. This work provides new insight into structural modification of known small-molecule BoNTAe inhibitors. It also demonstrates that SBCAMD is capable of improving potency of an inhibitor lead by nearly one order of magnitude, even for BoNTAe as one of the most challenging protein targets. The results are insightful for developing effective small-molecule inhibitors of protein targets with large active sites