Discovery of novel choline acetyltransferase inhibitors using structure-based virtual screening

Alzheimer disease and related dementias are major challenges, demanding urgent needs for earliest possible diagnosis to optimize the success rate in finding effective therapeutic interventions. Mounting solid scientific premises point at the core acetylcholine-biosynthesizing cholinergic enzyme, ChAT as a legitimate in vivo target for developing positron emission tomography biomarker for early diagnosis and/or monitoring therapeutic responses in the neurodegenerative dementias. Up-to-date, no PET tracer ligands for ChAT are available. Here we report for the first time a novel hierarchical virtual screening approach on a commercial library of ~300,000 compounds, followed by in vitro screening of the hits by a new High-Throughput ChAT assay. We report detailed pharmacodynamic data for three identified selective novel ChAT ligands with IC50 and Ki values ranging from ~7 to 26 µM. In addition, several novel selective inhibitors of the acetylcholine-degrading enzymes, AChE and BuChE were identified, with one of the compounds showing an IC50-value of ~6 µM for AChE. In conclusion, this report provides an excellent starting platform for designing and optimizing potent and selective ChAT ligands, with high potential as PET-imaging probe for early diagnosis of AD, and related dementias, such as Down’s syndrome and Lewy body disorders.

At present, dementia is the major cause of death affecting approximately 47.5 million people worldwide and this figure is projected to be double by 2030 1 . Alzheimer's disease (AD) type dementia alone prevalent in nearly 60-70% cases and designated as a "major killer" 2 . Other forms include dementia with Lewy bodies (DLB), frontotemporal dementia and vascular dementia. Beside almost a century of research in this field, there is no treatment available to cure the disease and only symptomatic treatments are available mainly indicating the use of acetylcholinesterase inhibitors to increase the availability of acetylcholine (ACh) in the diseased brain.
In vitro positron emission tomography (PET) imaging is gaining immense clinical impact and is an invaluable scientific tool for understanding the early pathological events in neurodegenerative disorders. It is also essential for effective monitoring of novel therapies and early diagnosis of neurodegeneration in AD 3 . In last few decades, increased number of labeled amyloid beta (Aβ) imaging agents based on conjugated Aβ specific dyes such as Congo red, thioflavin-T and PIB were developed and successfully tested for clinical diagnosis of AD 4 . However, as many as 30% of healthy elderly subjects with no clinical signs of dementia show PIB-retention in the brain. Whilst, some patients with no PIB-retention in the brain show severe cognitive deficits 5 . Aβ deposition is also a feature of DLB brain. Thereby, new more suitable in vivo PET biomarkers for a better disease prognostic and therapeutic evaluation are desirable.
Choline acetyltransferase (ChAT) (EC: 2.3.1.6; Choline O-acetyltransferase) is an important enzyme catalyzing the transfer of acetyl group from Acetyl-CoA to choline for synthesis of acetylcholine (ACh), which is a major neurotransmitter in the brain.
The neurons expressing ChAT are called cholinergic neurons and their communication with target tissues such as muscles depends on the functional ChAT. It has been observed that there is a decreased ChAT expression and activity in AD 6 . Therefore, ChAT has been proposed as a legitimate biomarker for early detection of AD and other neurodegenerative dementia disorders. Thus a PET tracer that can specifically bind to ChAT and help to monitor the health of cholinergic neurons will provide an important tool for early prognosis of AD.
In vitro screening of the 35 hit compounds for ChAT inhibitory activities. At first, 35 identified hits from the virtual screening protocol, were procured and screened in vitro for their inhibition of human rChAT enzyme activity using our new in-house high throughput fluorometric assay, at a single concentration of 100 µM. As a positive control for the compounds, the only commercially available ChAT inhibitor α-NETA 7 was included which resulted in 93% inhibition of the enzyme. This initial screening assay identified a total of three compounds: B1, B4, and E1 for human rChAT with more than 50% inhibition of enzyme activity at the tested concentrations ( Fig. 3; Table 2). All these compounds were selected and further assayed for their dose-dependent inhibitory activity and mode of inhibition determination.

In vitro estimation of inhibition constant and mode of inhibition for human rChAT hits.
In the initial assay, compound B1, showed 84% inhibition of human rChAT activity as compared to only 9% and 15% inhibition for AChE and BuChE activity, respectively, indicating a high selectivity for the human rChAT enzyme ( Fig. 3A; Table 2). Next, Ki and IC 50 values were calculated from the dose-response curves at different inhibitor and substrate concentrations as shown in Fig. 4B and D. The analysis yielded a Ki value of ~9.4 µM and IC 50 value of ~7.0 µM at substrate concentrations ranging from 20 to 80 µM. In order to determine the mode of inhibition and kinetic parameters (K m and V max ), Lineweaver-Burk (double reciprocal) plot 13 , was generated using substrate-velocity curves. In Lineweaver-Burk plots, reciprocal of rate (1/v) is plotted as a function of reciprocal of substrate concentration (1/S) for various inhibitor concentrations. The resultant plot is a straight line, with X-and Y-axis intercepts represent −1/ K m and 1/V max , respectively and the slope is K m /V max 14 . The Lineweaver-Burk plot showed an increase in K m value but a decrease in V max with increasing inhibitor concentrations, indicating mixed-model inhibition of human rChAT by the compound B1 ( Fig. 4C; Table 3). In addition, K m and V max values for each Scientific RepoRts | 7: 16287 | DOI:10.1038/s41598-017-16033-w inhibitor concentration (Table 3) are also determined from the substrate-velocity curves (Figs 4B, 5B and 6B). The mixed-model inhibition mode of action of compound B1 was confirmed by nonlinear regression fit of the data with GraphPad Prism 7 15 .
The compound B4 also showed high selectivity for human rChAT with 71% activity inhibition in the initial assay ( Fig. 3A; Table 2). In contrast, it only showed 22% activity inhibition for BuChE and almost no inhibition of AChE activity ( Fig. 3B and C; Table 2). Analogous to compound B1, the Ki and IC 50 values were calculated from the dose-response curves at different inhibitor and substrate concentrations as shown in Fig. 5B and D. The Ki value of ~11.9 µM and IC 50 value of ~16.5 µM at substrate concentrations ranging from 20 to 80 µM were obtained for B4 from the analysis. Further, the compound B4 affected both K m and V max values with increasing concentration (Table 3) and exhibited mixed-model type inhibition similar to B1 in the Lineweaver-Burk plot (Fig. 5C) and also in nonlinear regression fit analysis.
The compound E1, similar to B1 and B4 showed high selectivity for human rChAT with 55% activity inhibition in comparison to merely 12% and 17% inhibition of AChE and BuChE activity, respectively, at the initial tested concentration ( Fig. 3B and C; Table 2). Thus in the initial screening assay, the percentage inhibition of rChAT by E1 was lower as compared to B1 and B4 (Fig. 3A). In addition, calculation of Ki and IC 50 values from the dose-response curves at different inhibitor and substrate concentrations ( Fig. 6B and D) resulted in ~15.2 µM and ~25.4 µM, respectively, further indicated weaker inhibition by E1 as compared B1 and B4. Surprisingly, the Lineweaver-Burk plots for E1 (Fig. 6C) showed an exact similar pattern like B1 and B4; increased K m and decreased V max value with increasing inhibitor concentrations (Table 3), indicating mixed-model inhibition of human rChAT, which was further confirmed by nonlinear regression fit analysis.
Overall, all three hits (B1, B4, and E1) showed similar behavior in the substrate-velocity curves (Figs 4B, 5B and 6B) and a dose-dependent decrease in the enzymatic rate of the human ChAT-substrate catalytic reaction (Figs 4D, 5D and 6D; Table 3). More interestingly, in spite of their different basic structural scaffold and inhibition efficiency, all three hit compounds showed a similar inhibition mechanism for human rChAT, indicating that they can bind to either enzyme or enzyme-substrate (ES) complex upon interaction, thereby interfering with the enzymatic reaction. Thus in the context of being a potential in vivo tracer this property may be favorable as it suggests the ligand may have access to overall available ChAT bindings sites in the brain.
In vitro screening of the 35 hit compounds for selectivity over anti-targets, AChE and BuChE. As our main objective was to discover highly specific novel human rChAT ligands using our structure-based virtual screening method, we further screened the top 35 hit compounds for AChE and BuChE inhibitory activities in order to see whether any of the identified human rChAT hits can also inhibit either both Overall workflow of the structure based virtual screening based identification of novel ChAT inhibitors. The protein structure was prepared followed by preparation of Asinex Gold/Platinum library comprising 296537 compounds. After an initial filtering of compounds for optimal physicochemical properties, the library was screened using Surflex-Dock docking algorithm. The compounds with a cutoff score of 7 were further subjected to flexible docking with Surflex-GeomX approach. The top scoring compounds were inspected visually for their molecular interaction and were procured and tested in vitro using our in house developed nonradioactive assay method. Additionally, in vitro cellular cytotoxicity studies were also performed. or all of the three enzymes. We used a modified version of Ellman's colorimetric assay 16,17 , at a single concentration of 200 µM for the inhibition assay. This initial inhibition assay identified five hit compounds for AChE (B7, C1, C5, D3 and D4) and four for BuChE (A2, C2, C4, and C5), with more than 50% inhibition of enzyme activity at the tested concentration ( Fig. 3B and C; Table 2). Like ChAT inhibitors, Ki and IC 50 values of all nine hit compounds were calculated from the dose-response curves at different inhibitor and substrate concentrations ( Table 2). The dose-response curves at different inhibitor and substrate concentrations are shown in Figs 7 and 8 for only two top hits; C5 (for AChE) and C2 (for BuChE). The Ki and IC 50 values for all these hits were in the higher micro-molar range, except for compound C5 (for AChE, Table 2), indicating weak inhibition of both AChE and BuChE enzymes by these compounds. Moreover, out of these nine hits, only compound C5 showed non-selective inhibition of both AChE (~82%) and BuChE (~52%) with almost no inhibition of human rChAT in the initial screening assay ( Fig. 3; Table 2). The dose-response curves analysis yielded a Ki and IC 50 value of ~23.7 µM and ~6.3 µM for AChE, respectively ( Fig. 7B and D; Table 3). The IC 50 value of C5 for BuChE could not be calculated due to lack of complete inhibition even at the highest tested concentration of 300 µM. This clearly indicates that C5 is a highly selective and potent inhibitor of AChE as compared to BuChE. Additionally, we also determined the C5 mode of inhibition using Lineweaver-Burk plot, which showed a decrease in V max value while K m value remained almost unaffected with increasing inhibitor concentrations, indicating the non-competitive mode of inhibition of AChE by C5 (Fig. 7C). The nonlinear regression fit of the data also suggested non-competitive inhibition (with >99% probability), further complementing the double reciprocal plot data.
On the contrary, compound C2 showed selective inhibition of BuChE activity (78%) with no inhibitory effect on both AChE and human rChAT activity in the initial assay ( Fig. 3; Table 2). The calculation of Ki and IC 50 values from the dose-response curves at different inhibitor and substrate concentrations ( Fig. 8B and D) Table 1. Surflex-Dock score of virtual screening based identified compounds against ChAT, AChE and BuChE along with cLogP and cLogD. a Total Score is the total Surflex-Dock score expressed as -log(Kd) 36 . b Crash is the degree of inappropriate penetration by the ligand into the protein and of interpenetration (self-clash) between ligand atoms that are separated by rotatable bonds. Crash scores close to 0 are favorable. c Polar is the contribution of the polar interactions to the total score. # Marvin was used for calculating pKa/Log D @ pH 7. 59.6 µM and ~364.5 µM, respectively, suggested C2 as a weak inhibitor as compared to other inhibitors of human rChAT and AChE identified in our studies. Further, Lineweaver-Burk plot ( Fig. 8C) showed increased K m value with little change in the V max value with increasing inhibitor concentrations signifying the possibility of a competitive mode of inhibition, which was further confirmed by nonlinear regression analysis of the data.
Importantly, hit compounds identified as AChE and BuChE inhibitors were completely different and weaker in potency as compared to identified ChAT inhibitors and also showed a very diverse mode of inhibition. Overall, the data suggests that our virtual screening method is highly effective in discovering novel ChAT ligands with high specificity and potency.
Cell viability assay. Finally, to assess the prospective of the identified compounds for further development into ChAT ligands, we assessed the compounds active against ChAT, AChE and BuChE, for their cellular toxicity on human embryonic kidney (HEK 293) cells using the cell viability 3-(4,5-dimethylthiazol-2-yl)-2,5-diph enyltetrazolium bromide (MTT) assay 18 . All compounds were evaluated at 10 and 50 µM concentrations. The cell viability data is given in Table 4 (also represented as a bar chart in Fig. 9). Compounds B1, B7 and C1 reduced cell viability to less than fifty percent at 10 µM and 50 µM concentrations. Among the compounds active against ChAT, compounds B4 and E1 tends to be feasible for further development as they have shown lower cellular toxicity (Table 4).

Structural analysis of active hits. Binding mode analysis of ChAT hits. The active compounds against
ChAT were analyzed for their binding mode into the active site of ChAT. The docked poses of the compounds were extracted from the final step of virtual screening. The docking was performed using the ChAT crystal structure (PDB ID: 2FY3). The binding site of the enzyme was defined using co-crystallized choline coordinates. The docking scores (Total Score), which represent -logK d , for ChAT, AChE and BuChE is given in Table 1. The compounds B1, B4, and E1 were identified as most active compounds against ChAT with the IC 50 values of 7.0, 16.5 and 25.4, respectively. The three-dimensional docked pose and two-dimensional ligand interaction diagram of Human rChAT was pre-incubated with 100 μM of different Asinex compounds for 10-30 minutes at room temperature. The activity was measured using our in-house developed fluorometric assay. (B,C) rAChE & BuChE were pre-incubated with 200 μM of different Asinex compounds for 10-30 minutes at room temperature. The activity was measured using modified version of Ellman's colorimetric assay. The effect of these compounds on rChAT, AChE and BuChE was compared to the activity of a control sample that was pre-incubated with the buffer alone. The data is presented as mean ± SD of three individual experiments performed in triplicate. The dashed lines are representing 50% inhibition rate (IR) value; the compound with over 50% IR were selected for further kinetic analysis. rChAT = recombinant choline acetyltransferase; rAChE = recombinant acetylcholinesterase; BuChE = butyrylcholinesterase. B1, B4, and E1 are shown in Fig. 10(A-C). As given in the Fig. 10A, the amide carbonyl of compound B1 formed a hydrogen bond with the SER 540 amino acid residue of ChAT with a distance of 2.15 Å. Further, the 5H- [1,2,4] triazino [5,6-b]indole nucleus forms a hydrogen bond with the GLY 329 amino acid residue of ChAT active site. Interestingly, the HIS 324 , which is the responsible for the catalytic transfer of acetyl group from acetyl-coenzyme A to choline for synthesis of acetylcholine is not accessible after binding of the compound B1 (Fig. 10A), which might explain its mechanism of action.
The compound B4 has shown a different binding mode as compared to that of B1. As depicted in Fig. 10B, the terminal amide carbonyl forms a hydrogen bond with SER 438 at a distance of 1.88 Å, while the thiazole nucleus forms π -π interactions with the TYR 436 amino acid residue of the active site of ChAT.  The compound E1 exhibited a binding mode similar to the compound B1, making the HIS 324 inaccessible for the catalysis. The triazine nitrogen and the furan oxygen form hydrogen bonds with ILE 330 at a distance of 2.36 Å and 2.11 Å, respectively, while the benzene nucleus forms π -π interaction with TYR 552 amino acid residue (Fig. 10C).

Binding mode analysis of AChE and BuChE hits.
Of all the 35 compounds, which were evaluated in vitro against the anti-targets, AChE and BuChE, the compound C5 and C2 showed the most favorable enzyme inhibition kinetics for AChE and BuChE, respectively. We hence evaluated further their binding modes by docking them into the active site of AChE and BuCHE. The three-dimensional binding pose and the two-dimensional ligand interaction diagram for both of the compounds are given in Fig. 11A and B.
The compound C5 bound well to the binding pocket of AChE. The oxadiazole nucleus and its adjacent benzene ring formed π -π interactions with the TRP 286 and the piperidine nitrogen formed π -π interactions with the indole nucleus of TRP 86 amino acid residue (Fig. 11A). Docking of compound C2 into the active site of BuChE indicated that the indole nucleus of the compound formed π -π interactions with the indole nucleus of TRP 82 the amide carbonyl formed a hydrogen bond with TYR 128 at a distance of 2.24 Å, while, the oxygen atom adjacent to the benzene ring formed a hydrogen bond with the GLY 115 at a distance of 2.01 Å.

Discussion
Even though ChAT was discovered about a century ago 19 , surprisingly only few ChAT ligands (mostly inhibitors) are known today and they all lack the desired stability as well as brain permeability and thus have very limited applicability 9 . Here we report for the first time three new and structurally diverse ligands of ChAT. Detailed in vitro pharmacodynamic analyses indicated that these three ligands have considerable selectivity and potency as ChAT inhibitors. In this study, we employed a structure based hierarchical virtual screening approach for identifying the potential novel and specific ChAT ligands with optimal BBB permeability from Asinex Gold/Platinum small molecule database. Based on our initial analysis, the top 35 compounds were selected, purchased and tested in vitro for their inhibitory potential against ChAT using our newly developed high-throughput non-radiometric fluorescence assay. The three identified hit compounds for ChAT had Ki and IC 50 values in the sub-micromolar range (~7 to 25 µM).
Alzheimer's disease is the dominating form of dementia, with a prevalent of about 50-70% of the dementia cases. Up-to-date, only four FDA-approved drugs are available on the market, which however provides symptomatic relief. Three out of these drugs are targeted at the cholinergic system which is primarily and selectively affected in AD, DLB and Down's syndrome disorders 6,20,21 . Furthermore, most of the new drugs for AD are failing in the clinical trials phases, suggesting that early diagnosis and prevention may be pivotal for successful therapeutic intervention. Keeping in mind this urging need for novel biomarkers for early diagnosis of AD, the present work aimed at the discovery of novel ChAT ligands with the potential of radio-labeling, which can be used as an in vivo PET tracer to map regional brain distribution of ChAT, and thereby study the health of cholinergic network in CNS and neurodegeneration for diagnostic purposes. Given that cholinergic system is affected early in major dementia disorders, a ChAT-PET tracer may additionally have the potential for monitoring therapeutic response since an effective treatment is expected to prevent and/or to regress the functional deficit of the cholinergic network in the brain.
Amongst the identified ChAT ligands, two hits exhibited very low toxicity in the MTT assay, indicating good therapeutic window for these compounds. Although ChAT inhibitors are not suitable as AD therapeutic, they might be useful in other disorders. Intriguingly in this context, recent studies have shown a very critical role of ChAT in promoting cancer cell proliferation in both colon and lung cancer. The studies indicate that the cancer cells by upregulating ChAT may increase the synthesis of non-neuronal ACh, which act as an autocrine agent [22][23][24] . Their findings further suggested that any agents/therapy, which could target and stop the non-neuronal synthesis and release of ACh can be valuable in combatting disease progression. It is also possible that high ACh production and release may allow the cancer cells to avoid the immune surveillance since ACh is shown to be an effective agent for dampening the immune responses 25 . In light of these findings, the potent ChAT inhibitors identified in our study might also be developed as potential therapeutic candidates for the prevention and treatment of different forms of cancer.
The success of our hierarchical virtual screening approach was demonstrated by the identification of three ChAT ligands with high specificity, inhibition pattern, structural diversity, low toxicity, and novelty. Moreover, all these three ChAT ligands also showed a ClogP value of less than 5 indicating a good lipophilicity and CNS penetration. Importantly, the chemical scaffold of these hits can also be used as a template for designing highly selective analogs that can bind ChAT with much higher affinity and possibly be used as in vivo PET tracer biomarker for monitoring the possible dysfunctionality in the central cholinergic system. Such PET tracer may also be useful for diagnosis of the certain type of cancer tumors, where ChAT is overexpressed [22][23][24] .  These new ChAT ligands could be further improved in term of higher selectivity and potency using the in silico strategy for lead optimization of ChAT ligands that is reported in extensive details in our previous paper 26 . In addition, our newly developed ChAT assay can also be easily adapted as a high throughput in vitro screening method (based on 384-wells microtiter plates), with proved aptness for identification of novel and effective ChAT ligands.
Moreover, the screening of our top 35 hits for selectivity against anti-targets, AChE and BuChE, identified in total nine hits; five for AChE (B7, C1, C5, D3, and D4) and four (A2, C2, C4, and C5) for BuChE. Interestingly, only one hit, C5 (out of total nine top hits) showed non-selective inhibition of both AChE and BuChE in the initial assay. On further kinetic analysis, C5 showed a very high selectivity and an IC 50 value of ~6.26 µM for AChE. More interestingly, all these hits, except B7, exhibited a high therapeutic index in the MTT toxicity assay.
Even though AChE is the main enzyme target for the current AD drug therapies, recent findings have also implicated the importance of BuChE in the later stages of disease advancement after the progressive loss of AChE in AD patients [27][28][29] . Thereby opening the possibility of further exploring and optimizing the highly selective AChE and BuChE inhibitors discovered in our analysis as experimental drugs for the treatment of AD.
In summary, we report for the first time, three new and potent ChAT ligands identified using our novel hierarchical virtual screening strategy, with an aim of further optimizing and developing these novel ligands into in vivo PET tracer biomarkers for the early diagnosis and prevention of AD and other related neurogenerative disorders. Protein structure preparation. The 3D X-ray crystallographic structure of hChAT in complex with choline (PDB ID: 2FY3, UniProt ID: P28329) with a resolution of 2.27 Å, R-value of 0.211 and R-free value of 0.233 was downloaded from Protein Data Bank (PDB) and used for the study. It consists of residues from 120 to 733 out of a total 748 amino acid residues and has two binding pockets, one for CoA and other for choline molecule and the catalytic residue HIS 324 lies in the middle of these two binding pockets 6 . The enzyme structure was prepared for docking prior to defining docking site by addition of any missing hydrogen, assigning a bond order, repairing side chain, treating termini, fixing protonation state and side chain amides, removal of water molecules and other ligands as implemented in prepare structure wizard of SYBYL-X suite. The co-crystallized choline substructure was extracted and finally, the energy of protein structure was minimized utilizing Powell method 30 with Tripos force field. Similarly, the structure of AChE (PDB ID: 4EY7 31 ) and BuChE (PDB ID: 4BDS 32 ) were prepared and used for docking to identify selective compounds against ChAT.
Asinex database preparation and virtual screening. Asinex Gold/Platinum collection library consisting of 296537 compounds (structure data file; SDF) was downloaded from www.asinex.com and used for the study. The library was prepared by addition of explicit hydrogens, neutralizing charged structures and Figure 6. Inhibitory studies of compound E1 against human rChAT. The specified concentrations of compound E1 were pre-incubated with human rChAT and different substrate concentration at room temperature for 10-30 min. Afterwards; the activity was measured using our in-house non-radioactive fluorescence assay as described in the methods section. The effect of these compounds on rChAT was compared to the activity of a control sample that was pre-incubated with the buffer alone. generation of most probable ionization and tautomer state at pH 7.0 ± 2.0 using ligand preparation module of the SYBYL x 2.1.1 molecular modeling suite. To reduce the computational cost and increase the chance of successful identification of potential pre-clinical drug candidates, the library was screened against modified Lipinski parameters for CNS drug with potential to penetrate BBB 33 . These includes molecular weight between 150-450 dalton; cLogP <5; cLogD <4; Total polar surface area between 35-90; Hydrogen bond donors <3; Hydrogen bond acceptors <7; and Rotatable bonds <8. Further, the compounds were filtered for undesirable groups potential PAINS which make them promiscuous false positives by reacting with the screening assay method 12 . The prepared ligands were further subjected to Concord to obtain accurate and efficient, 3D global minima conformation supposed to be close to the bioactive conformation of the dataset molecules. The final dataset of 99823 compounds was obtained after preparation and was further used in the virtual screening process. The overall workflow is schematically depicted in Fig. 1.
The virtual screening process was performed in two steps using Surflex-Dock module interfaced in SYBYL-X 2.1.1 which is a fully automatic flexible molecular docking algorithm with a combination of empirical scoring function and a surface-based molecular similarity-based search engine 34 . The 'protomol' docking site was defined using the co-crystallized choline ligand in the active site of the enzyme. The compounds were docked into the defined active site using Surflex-Dock Screen (SFXC) for initial docking. The compounds were ranked according to their Total score and the compounds with a Total score >7 were passed to the second phase of more flexible and exhaustive docking using Surflex-Dock GeomX (SFXC). The output ligands from this step were further docked to the active site of AChE and BuChE as an additional check to identify only selective ChAT ligands. Finally, 35 compounds were selected based on their low total score against AChE and BuChE and high total score    against ChAT, and importatnt interactions with the active site amino acid residues of ChAT. The selected compounds were purchased commercially for in vitro testing using our in-house developed non-radioactive enzymatic ChAT assay method.  Purification of recombinant ChAT. DYT media (16 g/l Tryptone, 10 g/l yeast extract, 5 g/l NaCl, 100 µg/ ml ampicillin, 34 µg/ml chloramphenicol) was inoculated with a preculture of E.Coli BL21 Rosetta2 transformed with pProEXHTa-ChAT (a gift from Dr. Brian Shilton, Department of Biochemistry, University of Western Ontario, London, Canada). The bacteria were grown in shaking incubator at 37 °C with 200 rpm until the optical density at 600 nm reached 0.5. After which, 0.5 mM IPTG was added and His 6 -ChAT was expressed for circa 16 h at 18 °C. The bacteria were harvested and stored at −80 °C. His 6 -ChAT was purified with "Ni-NTA fast start Kit" (Qiagen) following the manufacturer's instructions. The elution buffer was exchanged to storage buffer (10 mM Tris pH 7.4, 500 mM NaCl, 10% (v/v) glycerol) using 30 kDa molecular weight cutoff Amicon Ultra concentrators (Merck Millipore). The purified protein was aliquoted, and stored at −80 °C. The purity and molecular weight of the protein was determined using SDS-PAGE electrophoresis. The total protein concentration was measured with DC Protein Assay (BioRad).

In vitro AChE and BuChE activity inhibition assay.
A modified version of Ellman's colorimetric assay 16,17 was adapted to a high throughput assay for the enzymatic activity of BuChE and AChE. The reagent, butyrylthiocholine iodide (BTC), acetylthiocholine iodide (ATC), 5,5′-dithiobis (2-nitrobenzoic acid) (DTNB) were purchased from Sigma-Aldrich (St. Louis, MO, USA). The modified assay details are as previously described 34,35 . The main modification concerned with the high throughput adaptation of the assay for use in 384-well plates. Briefly, 25 μL/well of a 1:450 diluted solution of a pooled human plasma and 1:768 diluted (3.5 ng/ml final concentration) purified recombinant human AChE protein (Sigma, Cat no. C1682) was used for measurement of BuChE and AChE activity, respectively. In the initial in vitro screening step, the wells were pre-incubated with 25 μL/wells of different Asinex compounds for 10-30 minutes at room temperature. The concentration of the compounds working solution was 3X to give a final concentration of 200 μM in a final volume of 75 μL in each well. The stock concentration of the compounds were prepared in 100% DMSO. The stock concentration was chosen usually at 10 mM, allowing the final DMSO concentration in the wells to be less than 4%. Finally, 25 μL of a cocktail mix prepared in Na/K phosphate buffer, containing DTNB (final concentration 0.4 mM) and BTC, (final concentration 1 mM) or ATC (final concentration 0.5 mM) was added to each well and the changes in absorbance was monitored at 412 nm wavelength for 15-20 minutes with one-minute interval, using a microplate spectrophotometer reader (Infinite M1000, Tecan). The rate of the enzyme activity was determined from the linear part of the kinetic reaction curves as ∆OD/time.

In vitro
ChAT activity inhibition fluorometric assay. ChAT activity was measured using our newly developed fluorometric assay, using human recombinant ChAT (rChAT) protein. The reagents, choline chloride, Figure 11. Molecular interactions of the identified hits C5 (A) and C2 (B) with active site amino acid residues of AChE and BuChE repectively, determined using Surflex-Dock based docking protocol. The active site amino acid residues are rendered as stick and the protein backbone as ribbon. The compounds are rendered as ball and stick and the hydrogen bonds between the compound and the active site amino acid residues are shown as yellow dotted line. acetyl coenzyme-A (ACoA, A2181) and 7-Diethylamino-3-(4-maleimidophenyl)-4-methylcoumarin (CPM) were purchased from Sigma-Aldrich (St. Louis, MO, USA).
Each Asinex compounds were run in triplicates. On each 96-wells plate, several enzyme wells without inhibitor were also included during measurements as control and for estimating the inhibition level. Negative controls were wells without enzyme. The percentage inhibition for each Asinex compound was calculated based on the enzyme control value as a reference (100% activity). The compounds showing over 50% inhibition of enzyme activity were selected for further kinetic studies.
Kinetic studies or in vitro estimation of Ki, IC 50, and mode of action of hits. For kinetic studies, a similar protocol as inhibition assay was followed; a dilution series of five different concentrations ranging from 10 −6 to 10 −9 M were prepared for each selected Asinex compounds. For ChAT, the A-CoA was kept constant at 10 μM (final) but the concentration of choline chloride varied between 320-10 μM. For BuChE, BTC was used as the substrate in the concentrations ranging from 0.0625 to 2 mM. For AChE, ATC was used as the substrate in the concentrations ranging from 0.0312 to 1 mM. Each compound concentration was measured in duplicates. The rate of enzyme activity (as ∆OD/hr kinetic data) was calculated and processed using the GraphPad Prism 7 analysis software 15 . The inhibitory constant (Ki) values were determined from the dose-response curve and the half-maximal inhibitory concentration (IC 50 ) values were calculated by plotting the percentage enzyme activity vs. the log of the Asinex compound concentrations and fitting the data using the nonlinear regression Enzyme Kinetics-Inhibition function.
The Michaelis-Menten constant (K m ) and maximal velocity (V max ) values were calculated from substrate-velocity curve after fitting the data with non-linear regression Michaelis-Menten kinetic function. The K m and V max values were further used to plot the Lineweaver-Burk plots; the plots were fitted using linear regression function.
Cell viability assay. The in vitro cellular toxicity of the active compounds was evaluated using MTT assay 18 .
Briefly, Human embryonic kidney epithelial (HEK-293) cells cultured at 37 °C in a humidified environment with 5% CO 2 in Dulbecco's Modified Eagle's medium (DMEM) in 96 well plate to obtain 70-80% confluency. Further, the cells were incubated overnight at 37 °C with 10 µM and 50 µM compounds in triplicate. Cells treated with 1% DMSO were used as controls. At the end of the treatment (12 h), 10 μL of 12 mM MTT dissolved in Phosphate buffer saline pH 7.4, was added to each well and incubated for an additional 2 h at 37 °C. The media was then removed from the plate and 100 µL of DMSO was added to dissolve the violet formazan crystals. Finally, the absorbance was measured at 540 nm with a reference wavelength of 630 nm on Tecan spectrophotometer with shaking before the reading. The percentage viability was calculated taking 1% DMSO control as 100 percent and presented as mean ± SD.