High-Throughput Screening Uncovers Novel Botulinum Neurotoxin

Diversity Center, 3805 Old Easton Road, Doylestown, Pennsylvania 18902, United States. ACS Comb. Sci. , 2016, 18 (8), pp 461–474. DOI: 10.1021/a...
0 downloads 3 Views 2MB Size
Subscriber access provided by University of Sussex Library

Article

High-Throughput Screening Uncovers Novel Botulinum Neurotoxin Inhibitor Chemotypes Kristin M. Bompiani, Dejan Caglic, Michelle C. Krutein, Galit Benoni, Morgan Hrones, Luke L. Lairson, Haiyan Bian, Garry R. Smith, and Tobin J Dickerson ACS Comb. Sci., Just Accepted Manuscript • DOI: 10.1021/acscombsci.6b00033 • Publication Date (Web): 17 Jun 2016 Downloaded from http://pubs.acs.org on July 12, 2016

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

ACS Combinatorial Science is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 51

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Combinatorial Science

High-Throughput Screening Uncovers Novel Botulinum Neurotoxin Inhibitor Chemotypes Kristin M. Bompiani,¶ ‡♯ Dejan Caglič,¶ ‡♮ Michelle C. Krutein,¶ Galit Benoni,¶ Morgan Hrones,¶ Luke L. Lairson,¶ Haiyan Bian,§ Garry R. Smith,§ and Tobin J. Dickerson¶*



Department of Chemistry, The Scripps Research Institute, 10550 North Torrey Pines Road, La

Jolla, CA, 92037 §

Fox Chase Chemical Diversity Center, 3805 Old Easton Road, Doylestown, PA, 18902



Current address UCSD Moores Cancer Center, 3855 Health Sciences Drive, La Jolla, CA

92093 ♮Current address BASF Corporation, 3550 John Hopkins Ct., San Diego, CA, 92121 ‡

These authors contributed equally to this work

*Corresponding author

KEYWORDS: high-throughput screening, botulinum neurotoxin, small molecule, inhibitor, activator

ACS Paragon Plus Environment

1

ACS Combinatorial Science

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 51

ABSTRACT: Botulism is caused by potent and specific bacterial neurotoxins that infect host neurons and block neurotransmitter release. Treatment for botulism is limited to administration of an antitoxin within a short time window, before the toxin enters neurons. Alternatively, current botulism drug development targets the toxin light chain, which is a zinc-dependent metalloprotease that is delivered into neurons and mediates long-term pathology. Several groups have identified inhibitory small molecules, peptides, or aptamers, although no molecule has advanced to the clinic due to a lack of efficacy in advanced models. Here we used a homogenous high-throughput enzyme assay to screen three libraries of drug-like small molecules for new chemotypes that modulate recombinant botulinum neurotoxin light chain activity. Highthroughput screening of 97,088 compounds identified numerous small molecules that activate or inhibit metalloprotease activity. We describe four major classes of inhibitory compounds identified, detail their structure-activity relationships, and assess their relative inhibitory potency. A previously unreported chemotype in any context of enzyme inhibition is described with potent sub-micromolar inhibition (Ki 200-300 nM). Additional detailed kinetic analyses and cellular cytotoxicity assays indicate the best compound from this series is a competitive inhibitor with cytotoxicity values around 4-5 µM. Given the potency and drug-like character of these lead compounds, further studies, including cellular activity assays and DMPK analysis, are justified.

ACS Paragon Plus Environment

2

Page 3 of 51

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Combinatorial Science

Introduction: Botulinum neurotoxins (BoNTs), which cause the disease botulism, are proteins produced by the Gram-positive bacteria Clostridium botulinum. Among the most potent toxins known, botulinum neurotoxins (BoNTs) bind to and enter host neurons, cleave intra-neuronal proteins, and impede acetylcholine neurotransmitter release [1]. Disruption of acetylcholine release by motor neurons results in flaccid muscle paralysis and, in severe cases, may result in death [2]. BoNTs have garnered considerable interest to date because of their potential use in widespread bioterrorist attacks, as reports show that several countries including Iran, Iraq, North Korea, and Syria have previously or are attempting to purify/weaponize the toxin [3, 4]. A 1995 United Nations report stated that Iraq purified 19,000 liters of toxin and weaponized approximately half of the material in delivery systems [3]. Correspondingly, researchers predict that widespread exposure of BoNTs due to a terrorist attack would cost billions of dollars in intervention and treatment, which could be significantly decreased if treatments are widely available and immediately implemented [5]. Eight immunologically distinct serotypes of BoNTs have been identified (A-H), where types A, B, and E typically cause disease in humans. Serotype A (BoNT/A) is the most studied form due to its extreme potency (LD50 in humans of approximately 1 ng/kg) [6] and long-term activity in neurons (up to 3-6 months) [7]. The mature BoNT/A toxin is a 150-kDa heterodimer comprised of a heavy and light chain linked via a disulfide bond, where both chains are required for toxicity. The heavy chain (HC) contains binding and internalization domains that facilitate toxin binding and translocation into host motor neurons, while the light chain (LC) is a zinc dependent metalloprotease that persists in neurons and hydrolyzes host proteins involved in neurotransmitter signaling [8]. HC binding to neuronal surface gangliosides and proteins triggers

ACS Paragon Plus Environment

3

ACS Combinatorial Science

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 51

receptor-mediated endocytosis, and the translocation domain subsequently chaperones the LC from the endosome into the cytosol [9]. The translocated LC/A localizes to the inner side of the cytoplasmic membrane where it cleaves its membrane-bound substrate, SNAP-25 (synaptosomal-associated protein of 25 kDa), a member of the SNARE (soluble NSF attachment protein receptor) protein family [10]. SNAP-25 cleavage impedes acetylcholine release into the neuromuscular junction, thereby impairing muscle contraction and resulting in flaccid paralysis. Current botulism treatment is limited and includes the administration of a neutralizing equine or human anti-toxin for adult or infant cases, respectively [11, 12]. However, the antitoxins are only effective if administered within a finite time window, that is, before the toxin enters neurons and becomes unavailable to the antibodies. After cellular internalization, patients who have been exposed to large amounts of toxin can only be managed in a clinical setting and placed on a ventilator to prevent respiratory arrest [13]. Moreover, limited stockpiles of these anti-toxins are available as they are not cost effective to produce, may have batch-to batch variation, and often involve immunizing animals or collecting sera from voluntarily immunized humans. Thus, new botulism post-exposure treatments are urgently needed in the event of widespread outbreak, such as a bioterrorist attack. One current therapeutic design approach is to identify compounds that inhibit intracellular LC activity, as these can be administered postexposure to aid in long-term treatment. While peptides [14-19] and aptamers [20] have been described that inhibit LC/A activity, no compound has advanced clinically. These compounds are challenging to develop as therapeutics given that they do not readily enter cells. Thus, many researchers have focused on developing small molecule therapeutics to inhibit LC/A activity. Several small molecule HTS campaigns with recombinant LC/A (rLC/A) have been described, although the majority of inhibitory compounds initially identified possess poor (i.e.,

ACS Paragon Plus Environment

4

Page 5 of 51

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Combinatorial Science

micromolar) inhibitory constants (Ki) [21-24]. More recently, inhibitory compounds with submicromolar Ki values have been reported, including a betulin-derived triterpenoid with a Ki of 800 nM [25], a quinolinol inhibitor with an estimated IC50 of 40 nM in a rat synaptosomal model [26], and a “hybrid” small molecule with a Ki of 600 nM [27]. Several potent hydroxamate inhibitors have been reported, including 2,4-dichlorocinnamic hydroxamate (Ki = 300 nM) [28], a related benzothiophene hydroxamic acid (Ki = 77 nM) [29], and a recently reported adamantine hydroxamate (Ki = 27 nM) [30]. Although screening has successfully identified a number of inhibitory chemotypes, no compound has advanced to clinical testing due to a lack of efficacy in advanced cellular or in vivo models. Ideally, candidate compounds would have a high aqueous solubility, low cross-reactivity with host metalloproteases, high cellpermeability and efficacy, and low cytotoxicity. Because no BoNT small molecules have advanced to clinical trials, new small molecule LC/A inhibitors with improved efficacy ex vivo/in vivo or improved pharmacokinetic properties may provide a basis for therapeutic advancement. Herein, we describe the results of a homogenous high-throughput screening (HTS) campaign of three different drug-like small molecule libraries with the goal of identifying unrecognized BoNT inhibitor chemotypes. Screening of approximately 97,000 compounds identified numerous compounds that activate or inhibit rLC/A activity. Four inhibitory scaffolds were further studied: two previously known metalloprotease/LC scaffolds and two new scaffolds, one of which is previously unreported and among the most potent botulinum inhibitors to date.

ACS Paragon Plus Environment

5

ACS Combinatorial Science

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 51

Results and Discussion: HTS assay optimization Several in vitro and in silico high-throughput screens (HTS) with the botulinum neurotoxin light chain (BoNT LC) have been reported, where most involved screening enzymatic activity of soluble LC with a fluorescent substrate [21-24]. In our studies, we utilized a previously established FRET assay with recombinant LC/A (rLC/A) and a commercial SNAPtide® substrate to identify compounds that increase or decrease the rate of enzyme activity based on fluorescent product generation [31, 32]. Initially, we optimized and scaled down the reaction relative to previous studies [31] to facilitate cost-effective screening of a larger library. The assay was performed in 384-well low volume microtiter plates, with a total assay volume of 10 µL instead of the 100 µL total volume for 96-well format [31]. The concentration of SNAPtide substrate was decreased from 5 µM to 0.7 µM. This substrate concentration did not significantly reduce the assay signal-to-noise ratio, and pilot experiments performed under these conditions indicated that the reaction operated under steady-state conditions with a linear increase in fluorescent product formation over a period of 2 hours (data not shown). The final DMSO concentration was 2%, which had a negligible impact on the rate of substrate hydrolysis (data not shown). Collectively, the optimization and miniaturization described here led to a 70-fold decrease in the amount of commercial SNAPtide substrate used without affecting the quality of the assay readout. The Z’factor for the SNAPtide assay was 0.97 and the signal window was 171.9, suggesting the assay is well suited to screening with a high predictive value of identified hits. Additionally, we utilized automated screening equipment to significantly increase our screening capacity (Figure 1).

ACS Paragon Plus Environment

6

Page 7 of 51

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Combinatorial Science

Figure 1: Graphical overview of the high-throughput screening method. Primary high-throughput screening (HTS) involved single concentration screening with 97,088 compounds from the three libraries indicated. The compounds were pre-incubated with recombinant light chain, FRET substrate peptide was added, and the rate of substrate cleavage over two hours was monitored. Compounds that significantly decreased LC/A activity (greater than 3 standard deviations for ChemBridge, or greater than 90% inhibition for ChemDiv and Life) were identified as hits, where secondary screening was performed with new batches of the small molecules at a single concentration in the FRET enzyme assay. Lead compounds were identified by comparing the relative potencies (IC50s) of compounds in dose-response experiments.

Three different libraries of small molecules from separate commercial vendors, ChemBridge, ChemDiv, and Life Chemicals were chosen for screening. The libraries are commercially available diversity sets that were filtered based on the following properties: (a) the compounds have desirable drug-like properties, such as low molecular weight, favorable number of hydrogen bond acceptors and donors, and hydrophilicity according to the Lipinski’s Rule of 5 [33]; (b) the compounds were amenable to structural modifications in order to optimize activity while still retaining drug-like properties, and (c) the compounds were devoid of potentially reactive or known problematic moieties including aldehydes, hydrazines, ketones, charged compounds, disulfides, phosphates, carboxylic and sulphonic acids, phosphonates, and sulfates (Supplemental Table 1). To perform the filtering, first, a similarity (Tanimoto) coefficient of

ACS Paragon Plus Environment

7

ACS Combinatorial Science

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 51

3σ below mean ≥90% inhibition ≥90% inhibition n/a

No. hits identified 20 36 26 82

Hit Rate (%) 0.10 0.07 0.09 0.08

Hit compounds were identified and analyzed independently for each library by clustering based on substructure similarities. From the 82 compounds analyzed, four major scaffolds were identified: an 8-hydroxyquinoline scaffold, an 8-sulfonamidoquinoline scaffold, a spiro(indolthiadiazole) scaffold, and a benzimidazole pyrazolopiperidinone scaffold (Figure 3). The 8hydroxyquinoline scaffold, which is a well-known LC/A inhibitor [35, 36], was present in all three libraries and comprised 35% of the hit pool (5 compounds in ChemBridge, 14 in ChemDiv, and 10 in Life). In addition to identifying potentially new hydroxyquinoline iterations, classification of this scaffold as hits in all three libraries indicates that the HTS was robust and sensitive. Spiro(indol-thiadiazole) compounds were present in both ChemBridge (1 compound)

ACS Paragon Plus Environment

10

Page 11 of 51

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Combinatorial Science

and ChemDiv libraries (2 compounds); and the 8-sulfonamidoquinoline scaffold was only present in the ChemBridge library (4 compounds). Finally, the benzimidazole pyrazolopiperidinone scaffold, which represents a previously unreported chemotype for enzyme inhibition, was present only in the ChemDiv library (5 compounds).

Figure 3: Inhibitory small molecule scaffolds identified in the HTS campaign.

Secondary screening validates the four inhibitory scaffolds identified and indicates they have a range of potencies Substructure searches with the four inhibitory scaffolds were performed within each respective commercial database to identify related compounds; these additional compounds were purchased from the vendor and included in the secondary screen. The compounds were validated in a secondary screen at the same concentration used for primary screening (12.5 µM final concentration for ChemBridge and 20 µM for ChemDiv and Life compounds) in the SNAPtide

ACS Paragon Plus Environment

11

ACS Combinatorial Science

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 51

FRET assay [32]. To compare the relative potencies and identify lead compounds, we additionally performed dose-titrations and calculated IC50 values. The 8-hydroxyquinoline scaffold has been previously reported in the BoNT literature [37]; here we selected an additional 21 compounds from ChemBridge and ChemDiv hit lists to validate with secondary screening. All 19 compounds from ChemBridge and ChemDiv libraries that were primary hits were positive in the secondary screen (0% false positive rate). Five compounds that were not primary hits exhibited greater than 90% inhibition of LC/A (14.7% false negative rate). Dose-titrations with the 19 compounds indicated that all but one hydroxyquinoline compound had IC50 values below 10 µM, where the most potent compound (5) exhibited an IC50 of 1.6 µM (Table 2). These IC50 values are consistent with previous studies that reported low-micromolar IC50 values [26, 35, 38]. Using these compounds as a starting point, this series has been expanded and detailed structure-activity and ADME studies have been performed [37]. Table 2: Selected 8-hydroxyquinoline compounds assessed in secondary screening that showed inhibitory activity.

Compound ID

Vendor ID

R1

R2

1

CB5175086

H

H

2

CB6948044

R3

H

R4

IC50 (µM)

H

9.3

H

6.7

ACS Paragon Plus Environment

12

Page 13 of 51

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Combinatorial Science

3

CB7012006

CH3

H

2.7

4

CB7633607

CH3

H

1.7

5

CB7637950

CH3

H

1.6

6

7706-0972

H

Cl

4.0

7

8397-0126

H

H

5.6

8

8397-0127

H

H

1.9

9

8397-0140

H

H

7.1

10

8397-0166

H

H

2.9

11

8397-0180

H

H

4.9

12

8397-0181

H

H

4.2

ACS Paragon Plus Environment

13

ACS Combinatorial Science

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 51

13

8397-0271

H

H

4.2

14

8397-0490

H

H

5.5

15

8397-0558

H

H

5.1

16

8397-0559

H

H

4.4

17

8397-0572

H

H

5.2

18

8397-0573

H

H

2.3

19

8398-1033

H

H

10.9

Next, we determined the potency of 27 spiro(indol-thiadiazole) compounds (1 hit and 26 additional structures) available from the ChemBridge small molecule repository. Structurally, these compounds possess a distinct similarity to a previously reported compound that potently inhibits the matrix metalloprotease ADAMTS-5 [39]. The single compound that was a primary hit exhibited inhibitory activity towards LC/A in the secondary screen (0% false positive rate). All spiro(indol-thiadiazole) compounds in the initial ChemBridge library were identified as hits (0% false negative rate). Dose-titrations identified 9 compounds with IC50 values below 5 µM,

ACS Paragon Plus Environment

14

Page 15 of 51

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Combinatorial Science

and 3 compounds (25, 31 and 33) had submicromolar IC50 values (Table 3). The remaining 15 compounds had no impact on LC/A enzyme activity when tested at 100 µM (Supplemental Table 3). The most potent compound in the spiro inhibitor series, 33, displayed an IC50 value of 0.64 µM. Some SAR trends could be identified based on this series. Substituted (compounds 27 and 29) as well as unsubstituted (compounds 25, 26, 28, 30, 31, 32 and 33) benzyl rings at the R1 position were well tolerated. Halogenation at the R4 position was tolerated (compounds 28 and 30), and as was methylation (compound 26) or ethylation (compounds 25 and 27). N-methylacetamide substitution at the R1 position and/or N-methyl-acetamide at the R3 position were not tolerated, as these compounds were completely inactive (compounds 134-149; Supplemental Table 3). Interestingly, addition of a 5-chloro group on the existing 2-chloro-benzyl substitution at the R2 position of compound 31 rendered compound 132, inactive. Our data agree with the previously published data for this scaffold [39] and suggest that these compounds may be broad metalloprotease inhibitors with modest inhibitory constants. SAR analysis of the spiro(indolthiadiazole) LC/A inhibitor series described here implicates a relative plasticity of the LC/A, which agrees with previous observations with other inhibitory chemotypes [37, 40]. Table 3: Spiro(indol-thiadiazole) compounds assessed in secondary screening that showed inhibitory activity.

Compound ID

Vendor ID

25 26

R1

IC50 (µM)

R2

R3

R4

CB5668345

H

H

-CH2CH3

0.71

CB5668437

H

H

CH3

1.62

ACS Paragon Plus Environment

15

ACS Combinatorial Science

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 51

27

CB5672440

H

H

-CH2CH3

1.33

28

CB5677976

CH3

H

Br

2.90

29

CB5679472

CH3

H

H

2.82

30

CB5681294

CH3

H

H

3.44

31

CB5825244

H

H

0.95

32

CB5830995

H

Br

2.25

33

CB5833514

H

H

0.64

H

We expanded upon the four initial 8-sulfonamidoquinoline hit compounds and screened an additional 47 compounds from the ChemBridge repository. This core is a well-known zinc chelator [41], and related compounds inhibit metalloproteases [42]. IC50 values were determined for all 51 compounds (Table 4 and Supplemental Table 4). Three compounds identified as primary hits were positive in the secondary screen, whereas one hit from the primary screen showed only moderate inhibition of LC/A in the secondary screen and was a false positive (7.7% false positive rate). There were no 8-sulfonamidoquinoline compounds in the initial ChemBridge library that were unidentified as LC/A inhibitors (0% false negative rate). A majority of the additional compounds (84%) were inactive and had IC50 values greater than 25 µM. Closer analysis and correlation with the primary screen data indicated that 12 of the inactive compounds

ACS Paragon Plus Environment

16

Page 17 of 51

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Combinatorial Science

were present in the initial library and showed no inhibitory activity in the primary or secondary screens, validating our screening protocol. The sulfonamidoquinoline compounds tested consisted of compounds with quinoline and naphthalene backbones (Table 4 and Supplemental Table 4). The two most potent compounds, 35 and 37, exhibited low micromolar inhibitory activities with IC50 values of 3.5 and 3.2 µM, respectively. Interestingly, only these two hits and one sulfonamidoquinoline compound with moderate inhibitory activity (36; IC50 = 6.7 µM) have a methoxy substitution at the R4 position of the quinoline backbone, apparently increasing inhibitory activity. Substitution of the R4 thio-triazole group on the naphthalene backbone was tolerated, where these compounds have IC50 values between 9.7 and 25.5 µM (compounds 3841). Finally, a moderately active sulfonamidoquinoline compound, 34, with 2,4,6-trimethylbenzyl substitution at the R1 position and the quinoline backbone, had an IC50 value of 18.8 µM. Overall, the sulfonamidoquinoline compounds exhibited only moderate LC/A inhibition with IC50 values in the micromolar range, similar to what has been previously observed for zincdependent matrix metalloproteases [42]. This scaffold does not appear as potent as other compounds, such as the hydroxamates, which have sub- to low-micromolar IC50 values [29]. While these compounds have higher IC50 values, they may be easier to optimize for solubility or stability in future DMPK studies. These factors may prove to play a large a role in clinical success. Table 4: 8-sulfonamidoquinoline compounds assessed in secondary screening that showed inhibitory activity.

ACS Paragon Plus Environment

17

ACS Combinatorial Science

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Compound ID

Vendor ID

X

34

CB5469620

35

Page 18 of 51

R2

R3

R4

IC50 (µM)

N

CH3

H

H

18.8

CB6626965

N

H

H

-O-CH3

3.5

36

CB6633275

N

H

H

-O-CH3

6.7

37

CB6991782

N

H

H

-O-CH3

3.2

38

CB7788159

C

H

OH

15.3

39

CB7796127

C

H

OH

9.7

40

CB7803009

C

H

OH

25.5

41

CB7816522

C

H

OH

12.3

R1

Substructure searches with the benzimidazole pyrazolopiperidinone scaffold identified 80 additional commercial ChemDiv compounds, 54 of which were present in the initial library but inhibited less than 90% of enzyme activity at the screening concentration of 20 µM and were not identified as hits. All 85 compounds were evaluated in parallel with a secondary screen, where 21 compounds were validated (inhibited >80% of LC/A activity at 20 µM, data not shown). All five compounds that were present in the initial screen hit list (≥90% inhibition in the primary screen, compounds 69, 70, 94, 97 and 98) were positive in the secondary screen and had a

ACS Paragon Plus Environment

18

Page 19 of 51

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Combinatorial Science

statistically similar percentage of inhibition (90 ± 12%; 0% false positive rate). Four compounds that were present in the initial library but displayed poor (0-75%) inhibition in the primary screen inhibited greater than 90% of enzyme activity in the secondary screen (7.4% false negative rate; compounds 66, 67, 84 and 100; data not shown). In addition to the 85 commercial compounds, 67 additional benzimidazole pyrazolopiperidinone compounds or scaffold substructures were synthesized and evaluated in dose-titrations to further evaluate SAR and identify active moieties within the scaffold (Tables 5 and 6, Supplemental Tables 5 and 6). Similar to other previously reported LC/A inhibitory scaffolds, the SAR landscape for the benzimidazole compounds appears relatively flat, and a variety of substitutions were well tolerated. Compounds containing a phenyl R1 substituent were extremely potent and demonstrated submicromolar IC50 values (Table 5). Substitution with fluorine (43-45), chlorine (57-59), and bromine (54-56) at the 2, 3, or 4 position on the phenyl was well tolerated compared to phenyl alone (42; IC50 = 0.14 µM). Di-substitutions with the same halogen (46-51, 60-63) or different halogens (64-68), as well as tri-fluoro methyl substituents (51-53) were also well tolerated at each position of the phenyl ring. Bromo-ether (69, 70), ester (72), nitrile (73), nitro (74-76), and O-trifluoromethyl (77) substituents were tolerated, although carboxylation at the 4-position significantly decreased potency (71; IC50 = 23.5 µM) (Table 5; panel of representative IC50 curves Supplemental Figure 1). Collectively, these data indicate that the phenyl R1 substituent may play a steric, but not an electrostatic role in LC/A inhibition, as the presence and position of electronegative substitutions on the phenyl did not significantly increase or decrease potency.

ACS Paragon Plus Environment

19

ACS Combinatorial Science

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 20 of 51

Table 5: Commercial benzimidazole pyrazolopiperidinone compounds that showed inhibitory activity toward BoNT LC/A.

Compound ID

ChemDiv or FoxChase ID

42

R1

R2

IC50 (µM)

J094-0647

CH3

0.14

43

J094-0172

CH3

0.36

44

J094-0173

CH3

0.27

45

J094-0168

CH3

0.32

46

J094-0170

CH3

0.18

47

J094-0169

CH3

0.15

48

J094-0171

CH3

0.10

49

J094-0655

CH3

0.14

ACS Paragon Plus Environment

20

Page 21 of 51

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Combinatorial Science

50

J094-0184

CH3

0.18

51

J094-0182

CH3

0.15

52

J094-0197

CH3

0.20

53

J094-0183

CH3

0.12

54

J094-0177

CH3

0.09

55

J094-0176

CH3

0.37

56

J094-0178

CH3

0.16

57

J094-0179

CH3

0.12

58

J094-0180

CH3

0.12

59

J094-0181

CH3

0.28

60

J094-0648

CH3

0.06

ACS Paragon Plus Environment

21

ACS Combinatorial Science

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 22 of 51

61

J094-0199

CH3

0.19

62

J094-0649

CH3

0.33

63

J094-0198

CH3

0.48

64

J094-0167

CH3

0.27

65

J094-0174

CH3

0.19

66

J094-0175

CH3

0.51

67

J094-0195

CH3

0.26

68

J094-0196

CH3

0.11

69

J094-0185

CH3

0.88

70

J094-0186

CH3

0.28

ACS Paragon Plus Environment

22

Page 23 of 51

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Combinatorial Science

71

J094-0187

CH3

23.50

72

J094-0188

CH3

0.30

73

J094-0190

CH3

0.78

74

J094-0194

CH3

0.14

75

J094-0192

CH3

0.43

76

J094-0193

CH3

0.74

77

FC2538

CH3

0.43

78

J094-0652

CH3

8.89

79

FC2549

CH3

2.23

ACS Paragon Plus Environment

23

ACS Combinatorial Science

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 24 of 51

80

FC2556

CH3

0.75

81

FC2557

CH3

0.24

82

J094-0626

CH3

0.88

83

J094-0633

CH3

0.43

84

J094-0625

CH3

0.70

85

J094-0628

CH3

0.66

86

J094-0634

CH3

0.10

87

J094-0632

CH3

1.13

88

FC3027

CH3

0.11

89

J094-0630

CH3

0.30

ACS Paragon Plus Environment

24

Page 25 of 51

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Combinatorial Science

90

J094-0631

CH3

0.53

91

FC3025

CH3

0.17

92

FC3026

CH3

0.11

93

J094-0627

CH3

0.25

94

J094-0640

CH3

0.48

95

FC2537

CH3

0.42

96

FC3024

CH3

0.11

97

J094-0642

CH3

6.44

98

J094-0643

CH3

4.52

99

J094-0644

CH3

8.06

ACS Paragon Plus Environment

25

ACS Combinatorial Science

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 26 of 51

100

J094-0645

CH3

3.43

101

FC2592

CH3

>50

102

J094-0897

CH3

0.45

103

FC2539

CH3

14.60

104

FC2548

CH3

0.85

105

FC-3138

CH3

0.33

106

FC2838

0.13

107

FC2839

0.25

108

FC2840

0.44

109

FC2841

0.31

ACS Paragon Plus Environment

26

Page 27 of 51

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Combinatorial Science

Pyridyl substitution at the R1 position significantly decreased potency compared to phenyl substitution (compound 78 IC50 = 8.89 µM vs. compound 42 IC50 = 0.14 µM). However, conjugation of a pyrrolidine ring to the pyridine increased potency, particularly when added to the 4-position (80; IC50 = 0.75 µM). Similarly, conjugation of a piperidine to the 4-position further increased potency, resulting a compound with similar affinity to the R1-phenyl compounds (81; IC50 = 0.24 µM). Compounds with a thiophene at the R1 position were generally potent inhibitors with IC50s comparable to the phenyl compounds (82, IC50 = 0.88, 83 IC50 = 0.43). Placement of the sulfur within the ring had a minor (less than 2-fold) impact on potency, and methylation and ethylation of the thiophene also had a minor impact on potency (84, 85); in contrast, dimethylation appreciably increased potency (86; IC50 = 0.10 µM). Addition of a t-butyl to the thiophene significantly decreased activity (87; IC50 = 1.13 µM), while a nitro- (93), bromo(90-92), chloro- (88), or dichloro- (89) substituent increased thiophene compound potency and resulted in compounds with similar potencies to the phenyl compounds. Conjugation of additional aromatic rings to the thiophene R1 substituent was also well tolerated (94-96). Pyrazole substitution at R1 severely decreased inhibitory activity, and these compounds demonstrated micromolar IC50s (97-100). R2 substitution with a t-butyl was well tolerated, as compounds with this moiety (106-109) displayed sub-micromolar IC50 values similar to the methylated analogs (Table 5). Interestingly, a majority of the cognate benzimidazole compounds that contain an Nmethyl imidazole have severely decreased inhibitory, suggesting that imidazole methylation may result in an activity cliff. The R2 methylated analogs of compounds 64 (IC50 = 0.27 µM), 45 (IC50 = 0.32 µM), and 47 (IC50 = 0.15 µM), all have a greater than 300-fold decrease in activity (compounds 202, 190 and 192, IC50 > 100 µM; Supplemental Table 5). R2 methylation resulting

ACS Paragon Plus Environment

27

ACS Combinatorial Science

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 28 of 51

in an activity cliff holds true for 25 of 27 of the methylated compounds, except for compounds 197 (IC50 = 0.78 µM) and 211 (IC50 = 0.75 µM). We hypothesize that the amine hydrogen may play a critical role in LC/A binding and inhibition (e.g., potentially coordinate the active site zinc), or the methyl group may be too bulky for incorporation into the binding site. However, given the small size of this addition and the relatively known plasticity of the enzyme, such a dramatic drop in activity due to steric interference seems unlikely [37, 40]. Substructure analysis of the benzimidazole pyrazolopiperidinone scaffold was performed by synthesizing a series of compounds with various scaffold fragments from one of the most potent compounds (47). Strikingly, analogs of compound 47 that lack the benzimidazole and have a hydrogen or methyl as R3 (110, 115) were largely inactive, save compound 125. In this case, installation of a trifluoro-methyl group at R1 restored inhibitory activity (125; IC50 = 0.75 µM), indicating that the benzimidazole moiety is not absolutely required for potent inhibition. Benzimidazole pyrazolopiperidinone substructures bearing a hydrogen (121), cyclopropyl (122), isopropyl (123), or t-butyl (124) R3 group were inactive (Table 6). Table 6: Substructure analysis of the benzimidazole pyrazolopiperidinone scaffold.

Compound ID 110

R1

R2

R3

IC50 (µM)

CH3

H

>100

ACS Paragon Plus Environment

28

Page 29 of 51

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Combinatorial Science

111

H

>100

112

H

>100

113

H

>100

114

H

>100

CH3

>100

116

CH3

>100

117

CH3

>100

115

CH3

>100

118

ACS Paragon Plus Environment

29

ACS Combinatorial Science

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 30 of 51

>100

119

CH3

>100

CH3

>100

122

CH3

>50

123

CH3

>50

124

CH3

>50

125

CH3

0.75

120

121

H

Compounds where the benzimidazole was replaced with other aromatic substitutions were also assessed to determine if activity could be rescued. Initially we mined the original HTS library and filtered the screening data for compounds containing solely the pyrazolopiperidinone

ACS Paragon Plus Environment

30

Page 31 of 51

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Combinatorial Science

core. Analogs containing a benzathiazole, benzoxazole, pyrimidine, and triazole-pyridazine in place of the benzimidazole were present in the primary screen; all of these compounds had significantly decreased activity compared to their benzimidazole counterparts and were not identified as hits (Supplemental Table 6). Additionally, we synthesized analogs that possessed a pyridine (251-269) or pyrimidine (270-273), which had no activity, and a ketophenyl group (284, 285), which had decreased activity (284; IC50 = 10 µM) (Supplemental Table 6). Compounds 274 and 276, which have a methoxybenzothiazole substitution, did possess some weak inhibitory activity with micromolar IC50s. Collectively, these data indicate that the benzimidazole core is specifically required for high affinity LC/A inhibition, and that while conjunction with the pyrazolopiperidinone may not be required, the presence of this group may increase potency. Strikingly, the most potent benzimidazole compound (60) has an IC50 of 60 nM, which places it among the most potent botulinum inhibitors reported. While several groups have previously described LC/A inhibitors that contain a benzimidazole, presumably as a zinc chelator, the combined benzimidazole pyrazolopiperidinone scaffold has not yet been reported in the context of any small molecule enzyme inhibitor. Eichorn et al. performed in silico screening of the NCI library and identified a benzimidazole dendrimer (NSC348884) that had protective activity in BoNT-intoxicated mice, although its effect was modest [43]. Cardinale et al. screened libraries from ChemBridge, ChemDiv, and Microsource with a similar FRET-based enzyme assay and identified several benzimidazole acrylonitrile and vinylbenzimidazole compounds that are hypothesized to be covalent inhibitors. The most potent compound displayed an IC50 of approximately 10 µM. While this compound was not active in cells, a derived compound, MBX1533, prevented SNAP-25 cleavage in chick neuronal cells; however, both compounds

ACS Paragon Plus Environment

31

ACS Combinatorial Science

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 32 of 51

were cytotoxic to human cells [21]. Compared to the previously reported compounds, the benzimidazole scaffold described here is significantly more potent. Detailed kinetic and cytotoxic analyses of lead benzimidazole compounds Lead compounds 45, 47, and 90 were further evaluated in an alternative enzyme assay with SNAP-25 derived substrate peptide (SNAP141-206) [44]. Unlike the fluorophore SNAPtide substrate, which only engages the LC/A active site and utilizes a fluorescent substrate, this unmodified 66mer peptide represents the critical residues of endogenous SNAP-25 that interact with the LC/A (active site and exosites). Substrate cleavage under steady state conditions is monitored via semi-quantitative liquid chromatography/mass spectrometry (LC/MS). Dose titration with 45 indicated dose-dependent inhibition. Initially, a non-competitive curve fit was used to estimate a Ki value of 368 nM (Figure 4a), which is similar to the IC50 value from the SNAPtide assay (IC50 = 0.32 µM). Similarly, 45 and 47 titrations demonstrated dose-dependent inhibition with Ki = 213 and 211 nM, respectively (Figure 4b,c). These data support that these novel compounds are among the most potent botulinum LC inhibitors to date.

ACS Paragon Plus Environment

32

Page 33 of 51

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Combinatorial Science

Figure 4: LC/MS data from the LC/A enzyme assay with 66mer substrate in the presence of the three lead benzimidazole compounds. Dose-response Ki titrations with lead benzimidazole pyrazolopiperidinone small molecules. Dose titrations with A) compound 45; B) compound 47; and C) compound 90. The data represent the mean ± SEM of two independent experiments. The data were fit with a non-competitive model of enzyme inhibition to calculate inhibitory constants with GraphPad Prism software. Representative Lineweaver-Burk linear replots with various concentrations of D) compound 45; E) compound 47; and F) compound 90. The data represent the linear replot of the average velocities shown in panels A, B, and C. Note that compounds 45 and 47 show significant deviation from linearity at high substrate and high compound concentrations, indicating more than one active enzyme species is present.

Closer analysis of the non-competitive fit for the LC/MS data indicated the model did not fit well at high substrate and high inhibitor concentrations. Lineweaver-Burk linear replots of the LC/MS kinetic data for 45 and 47 indicated noticeable curvature and deviation from linearity at high substrate concentrations and high inhibitor concentrations (Figure 4d,e). In contrast, the data for 90 appear more linear (Figure 4f). The deviation from linearity with the fluorophenyl benzimidazoles at high substrate concentrations implies that there is a mixture of enzyme species present. Two enzyme forms that have similar Km values, but different Vmax values may be mistaken for a single enzyme species, as the Vmax observed is the sum of the individual enzyme species[45]. Such non-linear data is often initially mistaken for non-competitive enzyme inhibition unless carefully assessed at high substrate concentrations, as we have done here. The data indicate that multiple rLC/A species are present in the LC/MS assay, and we propose that all three compounds are true competitive inhibitors. Compounds 45 and 47 appear to preferentially

ACS Paragon Plus Environment

33

ACS Combinatorial Science

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 34 of 51

inhibit one form of the enzyme species, while 90 may either inhibit both (or all) enzyme forms present, or may bias the equilibrium between the species and stabilize one conformation. While several papers have discussed the plasticity of rLC/A, no reports have described the presence of multiple enzyme forms. Typically screens are performed with a truncated form of the rLC/A because of its increased solubility and stability [46]. However, it’s possible that truncated rLC/A may exist in multiple forms in this screening assay. Whether any of these forms is representative of native toxin on cell membranes remains unknown and should be assessed for future screening campaigns. Lead small molecules show specificity toward the LC/A over other zinc metalloproteases A majority of BoNT LC/A small molecule scaffolds reported in the literature are competitive inhibitors that presumably coordinate the catalytic zinc, although a few noncompetitive inhibitors also have been reported [47, 48]. Correspondingly, screening with the SNAPtide substrate, which is a small peptide that only engages the LC/A active site, may bias hit identification toward classical competitive inhibitors. Previous studies have shown that quinolinols chelate divalent cations [36], and consequentially this scaffold has been previously identified through HTS against BoNT/A LC [35]. Similarly, spiro(indol-thiadiazole) and 8sulfonamidoquinoline scaffolds have been reported to inhibit zinc metalloproteases. Crystal structures of 8-sulfonamidoquinoline compounds indicate that the amide nitrogens coordinate zinc [41]. While the benzimidazole pyrazolopiperidinone scaffold has not been previously described, reports have shown that compounds bearing a benzimidazole moiety can chelate zinc. A crystal structure of the well-known serine protease inhibitor BABIM (bis(5-amidino-2benzimidazolyl)methane) in complex with trypsin shows that the zinc is coordinated by the two

ACS Paragon Plus Environment

34

Page 35 of 51

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Combinatorial Science

imidazole nitrogens [49]. Therefore, we hypothesize that all four scaffolds inhibit LC/A activity through zinc chelation. To test potential off-target effects, we characterized the inhibitory activity of lead compounds 4 (IC50 = 1.7 µM; 242-fold molar excess compound to LC/A), 33 (0.64 µM; 91x), 37 (IC50 = 3.2 µM; 457x), 47 (0.15 µM; 21x), and 90 (0.53 µM; 75x) toward other zinc metalloproteases. Kinetic assays with recombinant matrix metalloproteinases (MMP) 1 and 12 show that the compounds have 2,000 molar excess) the compounds inhibited 9-33% activity (Supplemental Figure 2a, b). In contrast, the MMP inhibitor marimastat completely inhibited activity of both MMPs at 100 nM (not shown). Similarly, these five lead compounds showed no inhibition of thermolysin activity at 1 µM (14-fold molar excess) and only 10-30% inhibition at 30 µM (400-fold molar excess) (Supplemental Figure 2c), while a positive control compound (phosphoramidon) inhibited >90% of thermolysin activity at 5 µM (data not shown). Collectively, these data suggest that the lead compounds, particularly the benzimidazole compounds 47 and 90, may have increased specificity toward the LC/A. Additional kinetic studies or crystallography data are needed to determine the mechanism of enzyme inhibition and define the inhibitor binding site. Previous reports have described small molecules bearing a benzimidazole moiety that showed protection in cellular assays, yet were cytotoxic to human cells [43]. Dose titrations with 47 and 90 from 0-25 µM in an MTS viability assay yielded LD50 values of 4.4 ± 0.2 µM and 6.4 ± 0.2 µM, respectively (Supplemental Figure 3). While these compounds appear cytotoxic at micromolar concentrations, there may be an optimal dose that engenders LC/A inhibition without imparting cellular toxicity, as only nanomolar concentrations are required for robust

ACS Paragon Plus Environment

35

ACS Combinatorial Science

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 36 of 51

LC/A inhibition in vitro. Additional DMPK and cellular studies are needed to determine compound solubility and cell permeability. Preliminary SAR analysis of all four of the scaffolds suggests a relatively flat landscape, where a majority of substitutions are well tolerated. This phenomenon has been previously reported with several LC/A inhibitory scaffolds [37, 40], and it has been hypothesized that the LC/A adopts a flexible conformation in solution that can accommodate a variety of spatial combinations and chemical substituents. Correspondingly, several recent studies have suggested that LC/A can exist in multiple conformations in solution, which agrees with our kinetic data suggesting two enzyme forms are present [50-54]. It is currently unknown if any of the rLC/A forms present in solution are representative of the LC/A active on cell membranes in vivo. Consequentially, rLC/A enzyme plasticity may contribute to activity discrepancies between in vitro and ex vivo assays. Assays utilizing the full-length LC are challenging, given that the recombinant protein is difficult to purify and tends to precipitate rapidly. Future cell-based HTS assays may provide better information about compounds that can inhibit LC/A within cells. Although such systems are currently being developed [55], none has yet been applied to botulinum HTS. Typically, LC/A inhibitors identified through homogeneous HTS are subsequently tested in cell-based assays, and if functional, progress to in vivo studies with mice. However, given that most compounds identified to date have significant liabilities such as poor inhibitory constants [23, 24], poor specificity [56], instability [57], poor solubility [25], and/or cytotoxicity [21, 22, 58] in vitro and ex vivo, only a few compounds have been tested in vivo [22, 26, 43, 59]. No botulinum therapeutic has advanced to phase I clinical trials due to lack of efficacy in animal models. In this context, new scaffolds such as those reported here are critical to the long-term

ACS Paragon Plus Environment

36

Page 37 of 51

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Combinatorial Science

development of a viable botulism therapeutic that could treat exposure to this potent toxin. Early integration of drug metabolism and pharmacokinetic (DMPK) studies during the discovery process may further select for better lead compounds and hasten progression into advanced models.

Experimental Methods: Compound libraries The 20,000-member ChemBridge small molecule library (ChemBridge, Inc., San Diego, CA), a representative subset of a much larger library from ChemBridge, was designed based on the following criteria: (a) the compounds have desirable drug-like properties, such as low molecular weight, a favorable number of hydrogen bond acceptors and donors, and hydrophilicity; (b) the compounds are amenable for structural modifications in order to optimize activity, while still retaining drug-like properties, and (c) the compounds are devoid of potentially reactive or known problematic moieties including aldehydes, hydrazines, ketones, charged compounds, disulfides, phosphates, carboxylic and sulphonic acids, phosphonates, and sulfates. The stock compounds were dissolved in DMSO to a final concentration of 2.5 mM and transferred into 63, 384-well plates. All compounds were subsequently screened at a final concentration of 12.5 µM. Selected individual compounds were ordered from ChemBridge for secondary screening and potency determination. The ChemDiv library (ChemDiv, Inc., San Diego, CA) used in this study is derived from a larger commercial diversity set and contained 48,928 compounds stored in 139, 384-well

ACS Paragon Plus Environment

37

ACS Combinatorial Science

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 38 of 51

plates; the stock compounds were dissolved in DMSO to a final concentration of 2 mM and screened at a final concentration of 20 µM. Library compounds were initially filtered based on the criteria previously described. Selected individual compound hits or related compounds were ordered from ChemDiv for secondary screening and additional studies. The Life library (Life Chemicals, Inc., Ontario, Canada) is a subset of a commercial diversity set and contains 28,160 compounds in 80, 384-well plates; stock compounds were dissolved in DMSO to a final concentration of 2 mM and screened at a final concentration of 20 µM. Library compounds were initially filtered based on the previously described guidelines.

High-throughput screening Library compounds were evaluated for their impact on recombinant BoNT/A LC enzymatic activity towards SNAPtide® FITC/DABCYL fluorogenic substrate (List Biological Laboratories, Campbell, CA). The assay was performed in 40 mM Hepes, pH 7.4, with 0.01% Tween-20 at room temperature (22 °C). Recombinant BoNT/A LC (expressed in Escherichia coli and purified as described elsewhere[46]) was diluted in Hepes assay buffer to 8.75 nM and 8 µL of the enzyme solution was dispensed into a low-volume 384-well black assay plate (Greiner) using a MultiFlo microplate dispenser (BioTek, Winooski, VT). ChemBridge small molecule compound libraries (2.5 mM stock concentration; 50 nL), as well as ChemDiv and Life small molecule compound libraries (2.0 mM stock concentration; 100 nL) were transferred from 384well source plates using a 50- or 100-nl 384-slot pin tool (V&P Scientific, San Diego, CA) and a Biomek FXP liquid handler (Beckman Coulter, Brea, CA). Compounds were preincubated with the enzyme for at least 15 min at room temperature. Finally, 2 µL of 3.5 µM SNAPtide solution

ACS Paragon Plus Environment

38

Page 39 of 51

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Combinatorial Science

in Hepes buffer with 10% DMSO was dispensed using a MultiFlo microplate dispenser (BioTek). The final DMSO concentration in the screening assay was 2%. Stacks of 20 assay plates were simultaneously transferred to a stacked EnVision multilabel plate reader (PerkinElmer, Waltham, MA); rapid fluorescence reads with λex = 490 nm, λem = 523 nm were utilized to sequentially read all 20 plates within time points. Fluorescence was read every 30 minutes for 3 hours, and DMSO negative controls showed linear fluorescent product generation over the entire 3 hours. Slopes were calculated for each compound and compared to the slopes of the uninhibited enzyme (DMSO negative control). Normalized rates were calculated based on the average DMSO negative control rate for each respective individual library of compounds.

Statistical analyses Assay performance measures of the low-volume 384-well SNAPtide assay were determined as follows. Z’-factor for this assay was determined by comparing relative velocities of BoNT LC/A cleavage of SNAPtide in the presence of the positive (2,4-dichloro-cynnamyl hydroxamate[28]) and negative (DMSO) controls in low-volume 384-well plates. It was calculated as Z’ = 1 – 3(st + sb)/(mt - mb), where st and sb are standard deviations of the negative and positive controls, respectively, and mt and mb mean velocities of negative and positive controls, respectively[60]. Signal window was determined using the same data set as for determining the Z’-factor and was subsequently calculated as SW = ((mt – mb) – 3(st + sb)) / st [61].

Determination of compound potencies (IC50 values)

ACS Paragon Plus Environment

39

ACS Combinatorial Science

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 40 of 51

Serial dilutions of compounds were prepared in DMSO, and 1 µL of each compound dilution was preincubated with 10 µl of 70 nM BoNT LC/A and 79 µl of 40 mM Hepes assay buffer in a 96-well black plate (Greiner, Monroe, NC) for 5 min at room temperature. Subsequently, 10 µL of 7 µM SNAPtide® was added to initiate the reaction. The final assay concentrations were 7 nM BoNT LC/A and 0.7 µM SNAPtide® in 2% DMSO. The well fluorescence was recorded every 5 minutes for a total of 1 hr and 45 min at room temperature on a Synergy MX plate reader (BioTek; λex = 490 nm, λem = 523 nm). Enzyme velocities used to determine the IC50 values were calculated from the linear portion of the fluorescence versus time data. IC50 values were determined with a quadratic four-parameter fit with the BioTek Gen5 software.

66mer LC/MS kinetics Recombinant BoNT/A LC activity toward SNAP141-206 (66mer) substrate was assessed at ambient temperature with 2% final DMSO in 40 mM Hepes, pH 7.4 and analyzed via LC/MS as previously

described[44].

SNAP141-206

66mer

peptide

substrate

(H2N



ARENEMDENLEQVSGIIGNLRHMALDMGNEIDTQNRQIDRIMEKADSNKTRIDEANQRA TKMLGSG-COOH) and

13

C standard internal control peptide (H2N-R-(13C)A-TKM-(13C)L-

(13C)G-S-(13C)G-COOH) were custom peptides purchased from Bachem (Bubendorf, Switzerland). Briefly, LC/A (20 µL; 2 nM final) was added to assay buffer (40 mM HEPES, pH 7.4; 20 µL) and appropriate concentrations of diluted small molecule (5 µL in 10% DMSO and assay buffer). The mixture was incubated at ambient temperature for 10 min, and the assay was initiated with the addition of various concentrations of 66mer (0-30 µM; 5 µL in 10% DMSO and assay buffer). After 30 minutes, 40 µL aliquots were removed and quenched with 5 µL of

ACS Paragon Plus Environment

40

Page 41 of 51

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Combinatorial Science

15% TFA (1.5% final), and 13C-labeled internal control was added to a final concentration of 200 nM (5 µL) to yield 50 µL total sample. Ten microliters were injected onto a Zorbax 300SB-C8 column and analyzed with an Agilent 1100 LC/MS system. The product and internal control peak areas were integrated with ChemStation software (Agilent, Santa Clara, CA) and the concentration of product formed was calculated by normalizing the product peak area to the internal control peak area. Binding inhibitory constants (KI) were calculated using models of non-competitive enzyme inhibition with the Prism GraphPad software (GraphPad, La Jolla, CA).

Cytotoxicity and cell viability Murine neuroblastoma cells (Neuro-2as, ATCC #CCL-131) were seeded in 96-well tissue culture plates (Corning Life Sciences, Corning, NY) at a density of 10,000 cells per well (100 µL) in quadruplicate. Cells were grown in DMEM media (Life Technologies, Carlsbad, CA) supplemented with 10% FBS (Atlanta Biologicals, Flowery Branch, GA) for 24 hours at 37°C with 5% CO2. The media was removed and replaced with various concentrations of smallmolecules diluted in a serum-free neuronal media comprised of a 1:1 mixture of DMEM and F:12 media (Life Technologies) supplemented with 1x N-2 supplement (Life Technologies) and 1x B-27 supplement without Vitamin A (Life Technologies). The small molecules were diluted in 100 µL of serum-free neuronal media at a final DMSO concentration of 1%. Cells were incubated with compound or vehicle DMSO control for an additional 24 hours at 37°C with 5% CO2. Cell viability was quantified with a CellTiter 96 AQueous One Solution cell proliferation assay (Promega, Seattle, WA) according to the manufacturer’s instructions. The cell viability

ACS Paragon Plus Environment

41

ACS Combinatorial Science

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 42 of 51

data for the experimental absorbance at A495 nm was normalized to the average absorbance for the DMSO negative control to calculate the percent viability.

Matrix metalloproteinase and thermolysin enzyme assays Recombinant MMP-1 and MMP-12 were purchased from Anaspec (Fremont, CA). Enzymatic activity was quantified with an MMP activity kit (Anaspec) in a 96-well plate according to the manufacturer’s instructions. Briefly, rMMP (50 µL) was incubated with small molecule compound (1 µL) at room temperature for 10 minutes. MMP substrate (50 µL) was added and the fluorescence (λex = 490 nm, λem = 523 nm) was measured with a Synergy MX plate reader (BioTek) every 5 minutes for 60 minutes total. The final assay volume was 101 µL per well, which contained 20 ng rMMP, 0-30 µM small molecule compounds or DMSO vehicle control (1% DMSO final), and 1X substrate. Marimastat (100 nM; Sigma Aldrich, St. Louis, MO) was used as a positive control. Fluorescent values were corrected for background fluorescence and the enzyme velocity was calculated with BioTek Gen5 software. Velocities were normalized to the DMSO vehicle control. Recombinant thermolysin was purchased from Sigma Aldrich, and enzyme activity was quantified with a protease activity assay kit from Abcam (ab112153; Cambridge, MA) in a 96-well plate according to the manufacturer’s instructions. Briefly, thermolysin (10 µL) was incubated with small molecule compound (1 µL) and 80 µL 1X assay buffer at room temperature for 10 minutes. Fluorescent substrate (10 µL) was added and the fluorescence (λex = 540 nm, λem = 590 nm) was measured with a Synergy MX plate reader (BioTek) every 5 minutes for 45 minutes total. The final assay volume was 101 µL per well, which contained 70 nM thermolysin, 0-30 µM small molecule compounds or DMSO vehicle

ACS Paragon Plus Environment

42

Page 43 of 51

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Combinatorial Science

control (1% DMSO final), and 1X substrate. Phosphoramidon (5 µM; Sigma Aldrich, St. Louis, MO) was used as a positive control. Fluorescent values were corrected for background fluorescence and the enzyme velocity was calculated with BioTek Gen5 software. Velocities were normalized to the DMSO vehicle control.

ASSOCIATED CONTENT

Supporting Information. Supplemental Figure 1 shows representative dose-response curves for the 3 lead benzimidazole compounds 45, 47, and 90 in the SNAPtide FRET-based LC/A enzyme assay. Supplemental Figure 2 shows off target zinc metalloprotease activity for compounds 4, 33, 37, 47, and 90 with MMP-1, MMP-12 and thermolysin. Supplemental Figure 3 shows the murine neuroblastoma cytotoxicity data for compounds 45 and 90. Supplemental Table 1 lists the potentially reactive moieties and problematic structures that were filtered for during the initial library generation. Supplemental Table 2 contains the structures of inactive spiro(indolthiadiazole) compounds, Supplemental Table 3 contains the structures of inactive 8sulfonamidoquinoline compounds, and Supplemental Tables 4, 5 and 6 include the structures of inactive benzimidazole pyrazolopiperidinone compounds. This material is available free of charge via the Internet at http://pubs.acs.org.

ACS Paragon Plus Environment

43

ACS Combinatorial Science

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 44 of 51

AUTHOR INFORMATION Corresponding Author *Tobin J. Dickerson - phone: +1 (858) 784-2522; e-mail: [email protected] Author Contributions K. Bompiani and D. Caglic conceived and performed experiments, analyzed data, and co-wrote the manuscript; M. Krutein, M. Hrones, B. Haiyan, and G. Benoni performed experiments and analyzed data; L. Lairson, G. Smith, and T. Dickerson conceived experiments, co-analyzed data, and co-wrote the manuscript. All authors have given approval to the final version of the manuscript. Funding Sources This work was supported by a grant from the National Institutes of Health (AI082190 to T.J.D.) and the California Institute of Regenerative Medicine (TB1-01186 and CL1-00502). Conflict of Interest The authors have no conflicts of interest to declare.

ACKNOWLEDGMENTS The authors would like to thank Sujata Godbole for technical assistance.

ABBREVIATIONS BoNT: botulinum neurotoxin; LC: light chain; HC: heavy chain; HTS: high-throughput screen; FRET: Förster resonance energy transfer, SAR: structure activity relationship

ACS Paragon Plus Environment

44

Page 45 of 51

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Combinatorial Science

REFERENCES [1] Montal, M. Botulinum neurotoxin: a marvel of protein design. Annu. Rev. Biochem., 2010, 79, 591-617. [2] Simpson, L. L. Identification of the major steps in botulinum toxin action. Annu. Rev. Pharmacol. Toxicol., 2004, 44, 167-93. [3] Council, U. N. S., Tenth report of the executive committee of the special commission established by the secretary-general pursuant to paragraph 9(b)(I) of security council resolution 687 (1991), and paragraph 3 of resolution 699 (1991) on the activities of the Special Commission. Council, U. N. S., Ed. New York, NY, 1995. [4]

Zilinskas, R. A. Iraq's biological weapons. The past as future? JAMA, 1997, 278, 418-24.

[5] St John, R.; Finlay, B.; Blair, C. Bioterrorism in Canada: An economic assessment of prevention and postattack response. Can. J. Infect. Dis., 2001, 12, 275-84. [6] Schantz, E. J.; Johnson, E. A. Properties and use of botulinum toxin and other microbial neurotoxins in medicine. Microbiol. Rev., 1992, 56, 80-99. [7] Sloop, R. R.; Cole, B. A.; Escutin, R. O. Human response to botulinum toxin injection: type B compared with type A. Neurology, 1997, 49, 189-94. [8] Lacy, D. B.; Tepp, W.; Cohen, A. C.; DasGupta, B. R.; Stevens, R. C. Crystal structure of botulinum neurotoxin type A and implications for toxicity. Nat. Struct. Biol., 1998, 5, 898-902. [9] Binz, T.; Rummel, A. Cell entry strategy of clostridial neurotoxins. J. Neurochem., 2009, 109, 1584-95. [10] Fernandez-Salas, E.; Steward, L. E.; Ho, H.; Garay, P. E.; Sun, S. W.; Gilmore, M. A.; Ordas, J. V.; Wang, J.; Francis, J.; Aoki, K. R. Plasma membrane localization signals in the light chain of botulinum neurotoxin. Proc. Natl. Acad. Sci. U. S. A., 2004, 101, 3208-13. [11] Arnon, S. S. Creation and development of the public service orphan drug Human Botulism Immune Globulin. Pediatrics, 2007, 119, 785-9. [12] Franz, D. R.; Pitt, L. M.; Clayton, M. A.; Hanes, M. A.; Rose, K. J. In: Efficacy of prophylactic and therapeutic administration of antitoxin for inhalation botulism; DasGupta, B. R. Ed.; Plenum Press: New York, 1993; pp. 473-476. [13] Brook, I. Botulism: the challenge of diagnosis and treatment. Rev. Neurol. Dis., 2006, 3, 182-9.

ACS Paragon Plus Environment

45

ACS Combinatorial Science

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 46 of 51

[14] Schmidt, J. J.; Stafford, R. G.; Bostian, K. A. Type A botulinum neurotoxin proteolytic activity: development of competitive inhibitors and implications for substrate specificity at the S1' binding subsite. FEBS Lett., 1998, 435, 61-4. [15] Schmidt, J. J.; Stafford, R. G. A high-affinity competitive inhibitor of type A botulinum neurotoxin protease activity. FEBS Lett., 2002, 532, 423-6. [16] Kumar, G.; Kumaran, D.; Ahmed, S. A.; Swaminathan, S. Peptide inhibitors of botulinum neurotoxin serotype A: design, inhibition, cocrystal structures, structure-activity relationship and pharmacophore modeling. Acta Crystallogr. Sect. D. Biol. Crystallogr., 2012, 68, 511-20. [17] Yiadom, K. P.; Muhie, S.; Yang, D. C. Peptide inhibitors of botulinum neurotoxin by mRNA display. Biochem. Biophys. Res. Commun., 2005, 335, 1247-53. [18] Zdanovskaia, M. V.; Los, G.; Zdanovsky, A. G. Recombinant derivatives of clostridial neurotoxins as delivery vehicles for proteins and small organic molecules. J. Protein Chem., 2000, 19, 699-707. [19] Moore, G. J.; Moore, D. M.; Roy, S. S.; Hayden, L. J.; Hamilton, M. G.; Chan, N. W.; Lee, W. E. Hinge peptide combinatorial libraries for inhilbitors of botulinum neurotoxins and saxitoxin: deconvolution strategy. Mol. Divers., 2006, 10, 9-16. [20] Chang, T. W.; Blank, M.; Janardhanan, P.; Singh, B. R.; Mello, C.; Blind, M.; Cai, S. In vitro selection of RNA aptamers that inhibit the activity of type A botulinum neurotoxin. Biochem. Biophys. Res. Commun., 2010, 396, 854-60. [21] Cardinale, S. C.; Butler, M. M.; Ruthel, G.; Nuss, J. E.; Wanner, L. M.; Li, B.; Pai, R.; Peet, N. P.; Bavari, S.; Bowlin, T. L. Novel Benzimidazole Inhibitors of Botulinum Neurotoxin/A Display Enzyme and Cell-Based Potency. Botulinum J., 2011, 2, 16-29. [22] Eubanks, L. M.; Hixon, M. S.; Jin, W.; Hong, S.; Clancy, C. M.; Tepp, W. H.; Baldwin, M. R.; Malizio, C. J.; Goodnough, M. C.; Barbieri, J. T.; Johnson, E. A.; Boger, D. L.; Dickerson, T. J.; Janda, K. D. An in vitro and in vivo disconnect uncovered through highthroughput identification of botulinum neurotoxin A antagonists. Proc. Natl. Acad. Sci. U. S. A., 2007, 104, 2602-7. [23] Eubanks, L. M.; Silhar, P.; Salzameda, N. T.; Zakhari, J. S.; Xiaochuan, F.; Barbieri, J. T.; Shoemaker, C. B.; Hixon, M. S.; Janda, K. D. Identification of a Natural Product Antagonist against the Botulinum Neurotoxin Light Chain Protease. ACS Med. Chem. Lett., 2010, 1, 268272.

ACS Paragon Plus Environment

46

Page 47 of 51

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Combinatorial Science

[24] Cai, S.; Lindo, P.; Park, J. B.; Vasa, K.; Singh, B. R. The identification and biochemical characterization of drug-like compounds that inhibit botulinum neurotoxin serotype A endopeptidase activity. Toxicon, 2010, 55, 818-26. [25] Silhar, P.; Alakurtti, S.; Capkova, K.; Xiaochuan, F.; Shoemaker, C. B.; Yli-Kauhaluoma, J.; Janda, K. D. Synthesis and evaluation of library of betulin derivatives against the botulinum neurotoxin A protease. Bioorg. Med. Chem. Lett., 2011, 21, 2229-31. [26] Singh, P.; Singh, M. K.; Chaudhary, D.; Chauhan, V.; Bharadwaj, P.; Pandey, A.; Upadhyay, N.; Dhaked, R. K. Small-molecule quinolinol inhibitor identified provides protection against BoNT/A in mice. PLoS One, 2012, 7, e47110. [27] Burnett, J. C.; Wang, C.; Nuss, J. E.; Nguyen, T. L.; Hermone, A. R.; Schmidt, J. J.; Gussio, R.; Wipf, P.; Bavari, S. Pharmacophore-guided lead optimization: the rational design of a non-zinc coordinating, sub-micromolar inhibitor of the botulinum neurotoxin serotype a metalloprotease. Bioorg. Med. Chem. Lett., 2009, 19, 5811-3. [28] Boldt, G. E.; Kennedy, J. P.; Janda, K. D. Identification of a potent botulinum neurotoxin a protease inhibitor using in situ lead identification chemistry. Org. Lett., 2006, 8, 1729-32. [29] Capek, P.; Zhang, Y.; Barlow, D. J.; Houseknecht, K. L.; Smith, G. R.; Dickerson, T. J. Enhancing the Pharmacokinetic Properties of Botulinum Neurotoxin Serotype A Protease Inhibitors Through Rational Design. ACS Chem Neurosci, 2011, 2, 288-293. [30] Silhar, P.; Silvaggi, N. R.; Pellett, S.; Capkova, K.; Johnson, E. A.; Allen, K. N.; Janda, K. D. Evaluation of adamantane hydroxamates as botulinum neurotoxin inhibitors: synthesis, crystallography, modeling, kinetic and cellular based studies. Bioorg. Med. Chem., 2013, 21, 1344-8. [31] Boldt, G. E.; Kennedy, J. P.; Hixon, M. S.; McAllister, L. A.; Barbieri, J. T.; Tzipori, S.; Janda, K. D. Synthesis, characterization and development of a high-throughput methodology for the discovery of botulinum neurotoxin a inhibitors. J. Comb. Chem., 2006, 8, 513-21. [32] Caglic, D.; Bompiani, K. M.; Krutein, M. C.; Capek, P.; Dickerson, T. J. A Highthroughput-compatible FRET-based Platform for Identification and Characterization of Botulinum Neurotoxin Light Chain Modulators. J. Vis. Exp., 2013, 50908. [33] Lipinski, C. A.; Lombardo, F.; Dominy, B. W.; Feeney, P. J. Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings. Adv. Drug Del. Rev., 2001, 46, 3-26.

ACS Paragon Plus Environment

47

ACS Combinatorial Science

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 48 of 51

[34] McAllister, L. A.; Hixon, M. S.; Kennedy, J. P.; Dickerson, T. J.; Janda, K. D. Superactivation of the botulinum neurotoxin serotype A light chain metalloprotease: a new wrinkle in botulinum neurotoxin. J. Am. Chem. Soc., 2006, 128, 4176-7. [35] Roxas-Duncan, V.; Enyedy, I.; Montgomery, V. A.; Eccard, V. S.; Carrington, M. A.; Lai, H.; Gul, N.; Yang, D. C.; Smith, L. A. Identification and biochemical characterization of small-molecule inhibitors of Clostridium botulinum neurotoxin serotype A. Antimicrob. Agents Chemother., 2009, 53, 3478-86. [36] Jones, J. G.; Poole, J. B.; Tomkinson, J. C.; Williams, R. J. P. The Relationship between Proton Dissociation Constants and the Stability Constants of Complex Ions. J. Chem. Soc., 1958, 2001-2009. [37] Caglic, D.; Krutein, M. C.; Bompiani, K. M.; Barlow, D. J.; Benoni, G.; Pelletier, J.; Reitz, A. B.; Lairson, L. L.; Houseknecht, K. L.; Smith, G. R.; Dickerson, T. J. Identification of clinically viable quinolinol inhibitors of botulinum neurotoxin A light chain. J. Med. Chem., 2014, 57, 669-76. [38] Lai, H.; Feng, M.; Roxas-Duncan, V.; Dakshanamurthy, S.; Smith, L. A.; Yang, D. C. Quinolinol and peptide inhibitors of zinc protease in botulinum neurotoxin A: effects of zinc ion and peptides on inhibition. Arch. Biochem. Biophys., 2009, 491, 75-84. [39] Bursavich, M. G.; Gilbert, A. M.; Lombardi, S.; Georgiadis, K. E.; Reifenberg, E.; Flannery, C. R.; Morris, E. A. 5'-Phenyl-3'H-spiro[indoline-3,2'-[1,3,4]thiadiazol]-2-one inhibitors of ADAMTS-5 (aggrecanase-2). Bioorg. Med. Chem. Lett., 2007, 17, 5630-3. [40] Smith, G. R.; Caglic, D.; Capek, P.; Zhang, Y.; Godbole, S.; Reitz, A. B.; Dickerson, T. J. Reexamining hydroxamate inhibitors of botulinum neurotoxin serotype A: Extending towards the beta-exosite. Bioorg. Med. Chem. Lett., 2012, 22, 3754-7. [41] Nasir, M. S.; Fahrni, C. J.; Suhy, D. A.; Kolodsick, K. J.; Singer, C. P.; O'Halloran, T. V. The chemical cell biology of zinc: structure and intracellular fluorescence of a zincquinolinesulfonamide complex. J. Biol. Inorg. Chem., 1999, 4, 775-83. [42] Rouffet, M.; de Oliveira, C. A.; Udi, Y.; Agrawal, A.; Sagi, I.; McCammon, J. A.; Cohen, S. M. From sensors to silencers: quinoline- and benzimidazole-sulfonamides as inhibitors for zinc proteases. J. Am. Chem. Soc., 2010, 132, 8232-3. [43] Eichhorn, T.; Dolimbek, B. Z.; Deeg, K.; Efferth, T.; Atassi, M. Z. Inhibition in vivo of the activity of botulinum neurotoxin A by small molecules selected by virtual screening. Toxicon, 2012, 60, 1180-90.

ACS Paragon Plus Environment

48

Page 49 of 51

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Combinatorial Science

[44] Capkova, K.; Hixon, M. S.; McAllister, L. A.; Janda, K. D. Toward the discovery of potent inhibitors of botulinum neurotoxin A: development of a robust LC MS based assay operational from low to subnanomolar enzyme concentrations. Chem. Commun. (Camb.), 2008, 3525-7. [45] Segel, I. H. In: Enzyme kinetics : behavior and analysis of rapid equilibrium and steadystate enzyme systems; Wiley: New York, 1975. [46] Baldwin, M. R.; Bradshaw, M.; Johnson, E. A.; Barbieri, J. T. The C-terminus of botulinum neurotoxin type A light chain contributes to solubility, catalysis, and stability. Protein Expr. Purif., 2004, 37, 187-95. [47] Salzameda, N. T.; Eubanks, L. M.; Zakhari, J. S.; Tsuchikama, K.; DeNunzio, N. J.; Allen, K. N.; Hixon, M. S.; Janda, K. D. A cross-over inhibitor of the botulinum neurotoxin light chain B: a natural product implicating an exosite mechanism of action. Chem Commun (Camb), 2011, 47, 1713-5. [48] Silhar, P.; Capkova, K.; Salzameda, N. T.; Barbieri, J. T.; Hixon, M. S.; Janda, K. D. Botulinum neurotoxin A protease: discovery of natural product exosite inhibitors. J. Am. Chem. Soc., 2010, 132, 2868-9. [49] Katz, B. A.; Clark, J. M.; Finer-Moore, J. S.; Jenkins, T. E.; Johnson, C. R.; Ross, M. J.; Luong, C.; Moore, W. R.; Stroud, R. M. Design of potent selective zinc-mediated serine protease inhibitors. Nature, 1998, 391, 608-12. [50] Mizanur, R. M.; Frasca, V.; Swaminathan, S.; Bavari, S.; Webb, R.; Smith, L. A.; Ahmed, S. A. The C terminus of the catalytic domain of type a botulinum neurotoxin may facilitate product release from the active site. J. Biol. Chem., 2013, 288, 24223-33. [51] Silhar, P.; Lardy, M. A.; Hixon, M. S.; Shoemaker, C. B.; Barbieri, J. T.; Struss, A. K.; Lively, J. M.; Javor, S.; Janda, K. D. The C-terminus of Botulinum A Protease Has Profound and Unanticipated Kinetic Consequences Upon the Catalytic Cleft. ACS Med. Chem. Lett., 2013, 4, 283-287. [52] Silvaggi, N. R.; Boldt, G. E.; Hixon, M. S.; Kennedy, J. P.; Tzipori, S.; Janda, K. D.; Allen, K. N. Structures of Clostridium botulinum Neurotoxin Serotype A Light Chain complexed with small-molecule inhibitors highlight active-site flexibility. Chem. Biol., 2007, 14, 533-42. [53] Thompson, A. A.; Jiao, G. S.; Kim, S.; Thai, A.; Cregar-Hernandez, L.; Margosiak, S. A.; Johnson, A. T.; Han, G. W.; O'Malley, S.; Stevens, R. C. Structural characterization of three novel hydroxamate-based zinc chelating inhibitors of the Clostridium botulinum serotype A neurotoxin light chain metalloprotease reveals a compact binding site resulting from 60/70 loop flexibility. Biochemistry (Mosc). 2011, 50, 4019-28.

ACS Paragon Plus Environment

49

ACS Combinatorial Science

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 50 of 51

[54] Burnett, J. C.; Schmidt, J. J.; McGrath, C. F.; Nguyen, T. L.; Hermone, A. R.; Panchal, R. G.; Vennerstrom, J. L.; Kodukula, K.; Zaharevitz, D. W.; Gussio, R.; Bavari, S. Conformational sampling of the botulinum neurotoxin serotype A light chain: implications for inhibitor binding. Bioorg. Med. Chem., 2005, 13, 333-41. [55] Dong, M.; Tepp, W. H.; Johnson, E. A.; Chapman, E. R. Using fluorescent sensors to detect botulinum neurotoxin activity in vitro and in living cells. Proc. Natl. Acad. Sci. U. S. A., 2004, 101, 14701-6. [56] Sheridan, R. E.; Deshpande, S. S. Interactions between heavy metal chelators and botulinum neurotoxins at the mouse neuromuscular junction. Toxicon, 1995, 33, 539-49. [57] Capkova, K.; Hixon, M. S.; Pellett, S.; Barbieri, J. T.; Johnson, E. A.; Janda, K. D. Benzylidene cyclopentenediones: First irreversible inhibitors against botulinum neurotoxin A's zinc endopeptidase. Bioorg. Med. Chem. Lett., 2010, 20, 206-8. [58] Burnett, J. C.; Ruthel, G.; Stegmann, C. M.; Panchal, R. G.; Nguyen, T. L.; Hermone, A. R.; Stafford, R. G.; Lane, D. J.; Kenny, T. A.; McGrath, C. F.; Wipf, P.; Stahl, A. M.; Schmidt, J. J.; Gussio, R.; Brunger, A. T.; Bavari, S. Inhibition of metalloprotease botulinum serotype A from a pseudo-peptide binding mode to a small molecule that is active in primary neurons. J. Biol. Chem., 2007, 282, 5004-14. [59] Pang, Y. P.; Davis, J.; Wang, S.; Park, J. G.; Nambiar, M. P.; Schmidt, J. J.; Millard, C. B. Small molecules showing significant protection of mice against botulinum neurotoxin serotype A. PLoS One, 2010, 5, e10129. [60] Zhang, J. H.; Chung, T. D.; Oldenburg, K. R. A Simple Statistical Parameter for Use in Evaluation and Validation of High Throughput Screening Assays. J. Biomol. Screen., 1999, 4, 67-73. [61] Iversen, P. W.; Eastwood, B. J.; Sittampalam, G. S.; Cox, K. L. A comparison of assay performance measures in screening assays: signal window, Z' factor, and assay variability ratio. J. Biomol. Screen., 2006, 11, 247-52.

ACS Paragon Plus Environment

50

Page 51 of 51

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Combinatorial Science

TOC Graphic:

ACS Paragon Plus Environment

51