Crystal Correlation Of Heterocyclic Imidazo[1,2-a]pyridine Analogues and Their Anticholinesterase Potential Evaluation

Imidazo[1,2-a]pyridine-based compounds are clinically important to the treatments of heart and circulatory failures, while many are under development for pharmaceutical uses. In this study, a series of imidazo[1,2-a]pyridine-based derivatives 2(a–o) were synthesized by reacting a-haloketones with 2-aminopyridines in a basic media at ambient temperature. Single crystal X-ray diffraction studies suggest that with low degree-of-freedom, the introduction of bulky adamantyl or electron-rich biphenyl moiety into the imidazopyridine derivatives will not affect its structural occupancy. Imidazo[1,2-a]pyridine-based derivatives with biphenyl side chain are potential AChE inhibitors. Compound 2h which bears a biphenyl side chain and methyl substituent at the position R4 of the imidazo[1,2-a]pyridine ring showed the strongest AChE inhibition with an IC50 value of 79 µM. However, imidazo[1,2-a]pyridine derivatives with phenyl side chain exhibit better BChE inhibition effect among the series. Compound 2j with 3,4-dichlorophenyl side chain and unsubstituted imidazo[1,2-a]pyridine ring appears to be the strongest BChE inhibitor with an IC50 value of 65 µM and good selectivity. The inhibitory effects of active compounds were further confirmed by computational molecular docking studies. The results unveiled that peripheral anionic sites of AChE and acyl pocket of BChE were the predominated binding sites for the subjected inhibitors.

(2-arylimidazo[1,2,-a]pyridinium salts) were subjected for the evaluation of electric eel acetylcholinesterase (AChE) inhibition. These compounds showed positive results with IC 50 values ranging from 0.2 to 50.0 μM, meanwhile several of them showed protective effects against the organophosphorus AChE inhibitor soman in mice 13 . Besides, N-(benzylidene)imidazo [1,2-a]pyridine derivatives with halogen as substituent on the attached phenyl ring, showed potential AChE inhibition 14 . Recently, the cholinesterase inhibition effects of imidazole analogues designed to be multi-targeted, were confirmed through computational studies, pharmacological evaluation using animal model and mechanistic in-vitro cholinesterase inhibition study 15 .
In view of the potential cholinesterase inhibition activities associated with the imidazo[1,2-a]pyridine ring system, it led to our present study on the synthesis of fifteen new imidazo [1,2-a]pyridine-based derivatives, 2(ao) and evaluation of their anti-cholinesterase activities. These new derivatives were characterized by spectroscopy methods and single crystal X-diffraction study. By using group comparison of single-crystal X-ray diffraction data, the structural conformation, structural occupancy and crystal packing similarity were studied. The AChE and BChE inhibitory activities of these derivatives were evaluated by Ellman's colorimetric test 16 . In addition, we also investigated in-silico binding mode of the biologically active ligands into AChE and BChE enzymes in comparison with tacrine as reference and protocol validation. The molecular docking procedure was validated by using tacrine as the native ligand.

Results and Discussion
Spectroscopy analysis. Cyclization of 2(a-o) was confirmed by FTIR, as the N-H stretching bonds were in the range of 3500-3300 cm −1 , which corresponded to the disappearance of primary amide. Furthermore, some peaks were observed at ~1370 and ~1200 cm −1 , which corresponded to the C-N stretching of imidazole. The presence of unsaturated C-H (pyridine) stretching was observed near 3100 cm −1 , and their v(C-C) and ω(C-H) were showed at ~1600 & 1450 cm −1 and ~750 cm −1 , respectively. Imidazo [1,2-a]pyridines-based derivatives contain adamantyl, biphenyl or phenyl moiety as their side chains. The signature absorption bands of these moieties can be identified from their FTIR spectra. Adamantyl's methylene and methine v(C-H) were observed as strong bands around 2900 and 2850 cm −1 . Meanwhile, the v(C-H) of biphenyl and phenyl moieties were revealed at ~3000 cm −1 . In the FTIR spectra of 2(a-o), the unexpected observations are the occurrences of O-H stretching in hydrated compounds 2b, 2c and 2e, as these compounds were co-crystallized with water molecules 17,18 .
In the 1 H-NMR spectra of 2(a-o), the singlet corresponding to the methylene (-CH 2 -) group in 1 was replaced by methine (-CH-) protons at a higher chemical shift (δ ≈ 7.1-8.2 ppm) which associated with the other required aromatic peaks. For compounds 2(a-e), adamantyl moiety was presented as broad peaks in the lower region, centering around δ ≈ 2.1, 2.0 and 1.8 ppm with the integration values of 3:6:6. For compounds 2(f-i), the first phenyl ring of the biphenyl moiety was shown as doublet and two sets of triplet near δ ≈ 7.7 − 7.4 ppm with the integration values of 2:2:1. Two well-resolved sets of doublet were centered around δ ≈ 8.2 and 7.7 ppm with the integration values of 2:2, ascribing to the -CH-protons of second phenyl ring. Identical J-coupling values were used to distinguish different methine protons from imidazo[1,2-a]pyridine, biphenyl and phenyl ring systems. In addition, protons of -CH 3 and -OCH 3 groups were revealed in the up-field region near δ ≈ 2.5 and 3.9 ppm. The numbers of protons were in agreement with the proposed values based on the integration values.
The 13 C NMR spectra of 2(a-o) showed four distinct sets of imidazole carbon, aromatic carbon, adamantyl carbon and saturated carbon signals. In the down-field region, both δ(C = N) and δ(C-N) imidazole carbon signals were centered around δ ≈ 146 and 145 ppm, respectively. The aromatic carbon signals of biphenyl and benzene groups were found in the range of δ ≈ 140 to 108 ppm. The adamantyl carbon signals were centered at δ ≈ 42, 37, 34 and 28 ppm, while the -CH 3 and -OCH 3 saturated carbon signals were observed in the up-field region, centering at δ ≈ 20 and 55 ppm, respectively [18][19][20] . All spectra are included in Supplementary Information. X-ray crystal structure description. The asymmetric unit (Z') of all studied compounds consist of a crystallographic independent molecule, except 2n and 2o, which consist of two molecules, while compound 2k consists of three molecules. Besides, the molecules of 2b and 2e lay on the mirror plane, therefore the asymmetric unit of 2b and 2e consist of half molecule. In addition, compounds 2b, 2c and 2e are crystallized with an extra water molecule as its solvate. Degree-of-freedom in the molecular conformations of the reported compounds can be characterized by the torsion angles between the imidazopyridine group with its attached adamantyl, biphenyl or phenyl groups C1-C7-C8-C9, τ1 (an extra torsion angle of C10-C11-C14-C15, τ2 for compounds containing biphenyl moiety, 2f-i).
For imidazo[1,2-a]pyridine derivatives with adamantyl side chain (2a-e), the torsion angles of C1-C7-C8-C9, τ1 are in the range of 0.0-7.1° (Fig. 1), while the imidazopyridine moiety is coplanar to the center of adamantyl (C9-C10-C13/14) with dihedral angle of 0.0-5.3°. Whereas, for compounds with biphenyl side chain (2f-i), the imidazopyridine and both phenyl rings are almost coplanar to each other with dihedral angles of 2.0-5.7° and 1.1-4.8°, respectively, as the torsion angle of τ1 and C10-C11-C14-C15, τ2 were nearly parallel with the range between 179.3-174.9° and 1.3-4.8°. As for imidazo [1,2-a]pyridine derivatives with phenyl side chain (2j-o), the phenyl ring had twisted away from the imidazopyridine ring, as the torsion angles of τ1 vary from 178.9° to 165.2°. Overall, the torsion angles in between imidazopyridine ring and its side chain (τ1) are always in periplanar conformation, although the dihedral angles have shown a discrepancy up to 25.7°. Selected torsional and dihedral angles for the reported compounds are tabulated in Table 1  were compared with 36 reported phenyl-imidazopyridine derivatives. In the aspect of molecular conformation, it is noteworthy to observe that no clinal conformation (τ1 between ± 30 °-± 150°) is adopted by current studied compounds, unlike the variations displayed in phenyl-imidazopyrine derivatives (Fig. 2). All present compounds adopt periplanar conformation (τ1 between 0-30° or 150°-180°) with C1-C7-C8-C9 torsion angles falling in the range from 0-7.1° degree and 165.2-179.3°, which is comparable to those previously reported phenyl-imidazolepyridine derivatives that adopt the same conformation (152.9-179.8°). The comparison of crystal structure occupancy between the search results and the present compounds are listed in Table 2. In the previously reported phenacyl benzoate derivatives, the replacement of phenyl moiety with adamantyl moiety had reduced the occurrence of π···π interaction, hence reduced the crystal packing coefficient 22 . Whereas, the introduction of biphenyl moiety as a replacement for phenyl moiety in phenacyl benzoate derivative had encouraged the formation of weak intermolecular π···π and C-H···π interaction, thus increased their crystal packing coefficient 23 . In contrast to the previously reported phenacyl benzoate derivatives, replacement of phenyl moiety with the adamantyl or biphenyl moiety does not arise a direct impact on the crystal packing coefficient in the present imidazopyridine derivatives.
As Fig. 3 shows the structure occupancy percentage, there is no obvious pattern as their structure occupancy were scattered within the range of the reported imidazopyridine derivatives. These might arise from the existence of imidazopyridine moiety and the low degree-of-freedom in the present compounds. As mentioned above, crystal packing coefficient were under direct influence of weak π···π and C-H···π interaction. While in the present compound, the existence of imidazopyridine ring system had play a vital role in those π···π and C-H···π interaction, as it is involves in 9 out of 14 crystal packing of those compound. In addition, low degree-of-freedom had limited the packing patterns of the present compound. Therefore, molecules in the present compound tend to  pack in a similar pattern with average crystal packing coefficient, which indicated by the occurrence of isostructural with 3D and 2D packing similarities.
In the investigation of the crystal structural similarity among the present compounds, two pairs of compounds (2g/2i and 2n/2o) were found to be crystallized in the same space group with similar lattice constants, these are the two main designations for 3D structural similarity. In fact, their crystal structures overlaid diagrams show that those two pairs of compounds as isostructure with similar packing pattern (Figs 4 and 5). Besides, 2D similarities are observed in between 2l/2n and 2l/2o, while compound 2f, 2g, 2n, 2i, 2k and 2m shared also 2D similarity (Fig. 6).  Cholinesterases Inhibitory Activities and Their Molecular Docking Studies. In the present study, all synthesized compounds were examined for their AChE and BChE inhibitory activities by Ellman's assay 16 , and their IC 50 values are tabulated in Table 3. The present imidazo[1,2-a]pyridine derivatives contain biologically active adamantyl, biphenyl or phenyl moiety as their side chains, which showed different levels of inhibitory activities towards AChE and BChE. Compounds 2a, 2f, 2j and 2l with unsubstituted imidazo[1,2-a]pyridine ring were synthesized as the reference compounds to accentuate the electronic effects of different substituents on the imidazo[1,2-a]pyridine ring, i.e. electron-donating (CH 3 ) and electron-withdrawing (Cl) groups, upon their cholinesterase inhibitory activities.
In contrast to the inactive activity in AChE inhibition, derivatives with phenyl side chain gave better BChE inhibition as compared to the adamantyl or biphenyl side chain if their phenyl side chains are substituted with strong electron-withdrawing Cl group, i.e. compounds 2j, 2k and 2n. Compounds 2j and 2k bearing with 3,4-dichlorophenyl side chain are two strongest BChE inhibitors with IC 50 values of 65 and 191 µM, respectively. There is a significant difference in the BChE inhibition activities of compounds 2n (IC 50 = 722 µM) and 2o (IC 50 > 2000 µM) although their phenyl side chains are substituted with electron-withdrawing groups, i.e. 4-chloro and 4-bromo groups, respectively, deducing that the oversized Br atom enfeebled the BChE inhibition of 2o. Meanwhile, compound 2m with its phenyl side chain substituted with electron-donating methoxy (OCH 3 ) group is inactive towards BChE inhibition up to a concentration of 2000 µM, which emphasizes the importance Figure 6. Crystal packing relationship in some studied compounds. Red boxes and yellow arrows indicate 3D and 2D similarities, respectively.

Compounds
Substituents Tacrine 0.24 ± 0.04 0.03 ± 0.01 Molecular docking studies of the potent imidazo[1,2-a]pyridine-based derivatives were performed to provide putative binding modes within the cholinesterase enzymes. The validation of the docking accuracy was carried out by docking the native co-crystallized ligand (tacrine) back into the parent enzymes and the conformation of the best-scored pose was compared with the bound ligand in the native crystal. The overall structure of human's BChE is similar to a human's AChE with phenylalanines (Phe) in the acyl binding pocket of AChE being replaced by leucine (Leu) and valine (Val) in BChE 24 .
In the protein crystal structure of 1ACJ, the inhibitor tacrine was bound to the catalytic active site of the AChE protein via π···π interactions, involving residues Trp84 and Phe330 (Fig. 8a). The docking models of active compounds 2i, 2g, 2c, 2h and 2f are illustrated in Fig. 8b-f, respectively. The molecule of 2i and 2g, which bear substituted imidazo[1,2-a]pyridine ring at the position R 2 were not penetrated to bound with any residue in the catalytic active site. Instead, the molecules of 2i and 2g only managed to bind only with the peripheral anionic site, which lies essentially on the surface of AChE 25 by π···π and C-H···O interactions (involving residues Trp279, Tyr 334, Tyr112 and Tyr70 for 2i, and residues Trp279, Tyr121 and Tyr70 for 2g) (Fig. 8b,c) and resulted in weaker inhibitory effect. Similar to tacrine, the imidazole ring of molecule 2c was bound to the catalytic active site through π···π interactions with residues Trp84 and Phe330. Meanwhile, its methyl group substituted at the position R 3 was favored to bind with the residues Trp437 and Ile439 of the hydrophobic pocket via C-H···π and C-H···alkyl interactions (Fig. 8d). However, molecule 2c failed to bind with any residues at the peripheral anionic site, which eventually leads to a moderate inhibition against AChE. As the most potent AChE inhibitors, the molecules of 2h and 2f were bound in the AChE enzyme with similar binding mode, involving both catalytic active site and peripheral anionic site. Their biphenyl side chains tend to form π···π interaction with residues Trp279 and Tyr334 at the peripheral anionic site 26 . The imidazo[1,2-a]pyridine rings of molecules 2h and 2f with unsubstituted position R 2 were bound to the catalytic active site and oxyanion hole of the AChE enzyme through C-H···O and amide···π interactions, involving residues Trp84, Gly117 and Gly118 (Fig. 8e,f) and results in a stronger AChE inhibitory effect as compared to aforementioned compounds 2i, 2g and 2c.
As a summary, the molecular docking results showed that imidazo[1,2-a]pyridine ring was able to bind in the catalytic active site via π···π and C-H···O interactions (Trp84 in AChE; Trp82 and His438 in BChE). The finding reveals that AChE inhibitory effects are attributed to the biphenyl ring, which enable the imidazo[1,2-a] pyridine-based derivatives to bind in both catalytic active site and peripheral anionic site through π···π and C-H···O interactions. On the other hand, BChE inhibitory effects are mainly contributed by halogen···O interaction between the dichloro-substituted phenyl side chain and the residue of acyl pocket of BChE.

Material and Methods
The reagents and solvents were obtained commercially from Sigma Aldrich and Merck Corporation, and were used without additional purification. Melting points were determined on Stuart (UK) SMP10 apparatus. 1 H and 13 C nuclear magnetic resonance (NMR) spectra were determined in CDCl 3 at 500 MHz and 125 MHz, respectively, using Bruker Avance III 500 spectrometer. Fourier transform infrared spectroscopy (FTIR) spectra were recorded on Perkin Elmer Frontier FTIR spectrometer equipped with attenuated total reflection (ATR). Mass spectra were recorded using Agilent 5975C TAD Gas Chromatograph/Mass Selective Detector System with HP-5MS GC column. The X-ray analysis of all target compounds excluding 2d were performed using Bruker APEX II DUO CCD diffractometer employing Mo Kα radiation (λ = 0.71073 Å) with ϕ and ω scans. Data reduction and absorption correction were performed using SAINT and SADABS programs 28 . All structures were solved by direct methods and refined by full-matrix least-squares techniques on F 2 using SHELXTL software package 29 . All non-hydrogen atoms were refined anisotropically (except for the minor disordered component of 2k with site-occupancy less than 0.2) and all C-bound H atoms were calculated geometrically with isotropic displacement parameters set to 1.2 (or 1.5 for methyl group) times the equivalent isotropic displacement parameters of the parent carbon atoms. Crystallographic data are collected in Supplementary Information.

Synthesis.
Target compounds were synthesized via a two-step reaction (Fig. 10). First, various substituted-ethanones were refluxed with N-bromosuccinimide (NBS) and petroleum ether in methanol at 333 K for two to five hours. Progress of the reactions were monitored by thin layer chromatography (TLC) plate using silica gel, with acetone:benzene (1:1) solvent system. The substituted-2-bromoethan-1-one (1) precipitate was filtered and recrystallized with ethanol. After that, 1 (0.002 mol) was reacted with 2-aminopyridines (0.003 mol) with the presence of potassium carbonate in DMF (8 ml) and stirred at room temperature for about three hours. The reaction progress was monitored by TLC method. After the reaction completed, the reaction mixture was poured into crushed ice. The precipitate formed, 2(a-o) were filtered out and recrystallized using acetone after dried 30 . All target compounds were synthesized in good yield and high purity, their chemical structures were characterized by using FTIR, NMR, mass spectroscopy and single-crystal X-ray diffraction analysis.                14  Cholinesterase inhibition assay. Ellman's colorimetric method with minor modification was employed in the determination of AChE and BChE inhibitory activities 16 . Electric eel acetylcholinesterase (AChE, EC 3.1.1.7), equine serum butyrylcholinesterase (BChE, EC 3.1.1.8), 5,5′-Dithiobis(2-nitrobenzoic acid) (DTNB) and sodium phosphate dibasic (Na 2 HPO 4 ) were purchased from Sigma-Aldrich (USA). S-butyrylthiocholine chloride (BTCC) from Sigma-Aldrich (Switzerland), acetylthiocholine iodide (ATCI) from Sigma-Aldrich (United Kingdom) and sodium phosphate monobasic (NaH 2 PO 4 ) from Sigma-Aldrich (Germany) were purchased. Tacrine hydrochloride was purchased from Cayman Chemical Company and absolute ethanol from Merck Emsure ® Germany was used as the solvent. Ethanol was used to replace the sample as a control in this assay 31 . Tacrine hydrochloride was used as the reference drug in the inhibition assays of AChE and BChE.

Statistical analysis.
The results of cholinesterase inhibitory assay were expressed as mean ± standard deviation (SD) and were labeled if p < 0.05 by using ANOVA of IBM SPSS Statistics for Windows, Version 23.0 (IBM, New York).

Docking protocol. Molecular docking was performed for the most potent inhibitors using Genetic
Optimization for Ligand Docking (GOLD) package 5.4.1 [32][33][34] . Crystal structures of AChE (PDB ID: 1ACJ 35 ) and BChE (PDB ID: 4BDS 36 ) isolated from Teronarce californica and Homo sapiens, respectively, with tacrine were obtained from Protein Data Bank. Ligand-protein binding space and the conformational flexibility of ligand inside the protein were explored by Genetic Algorithm (GA). A spherical binding site with a radius of 15 Å was used, it lays around the binding site of tacrine covering both catalytic anionic and peripheral anionic sites. 100 GA runs were carried out and the top 100 ranked docking poses were scored using the Piecewise Linear Potential (PLP) scoring function. Default values were used for all other parameters. The intermolecular interactions of the best scored pose of each ligand were analyzed and illustrated using Discovery Studio 4.5 software (Discovery studio v4.5.0.15071. Accelrys Inc 2015).

Conclusion
A series of imidazo[1,2-a]pyridine-based derivatives 2(a-o) were synthesized and characterized by FTIR, NMR and GCMS spectroscopy analyses. The 3D structures of these derivatives were further confirmed by single crystal X-ray diffraction studies and the result showed that they tend to adapt periplanar conformation with its substituted side chain at the position R 1 . With the presence of imidazo[1,2-a]pyridine ring and the low degree-of-freedom, the introduction of adamantyl or biphenyl does not affect their structural occupancies. Meanwhile, the present compounds tend to pack in a similar pattern with π···π and C-H···π interactions, involving the imidazo[1,2-a]pyridine ring, these had led to redundant 2D or 3D structural similarity in imidazo[1,2-a] pyridine-based derivatives. The AChE inhibitory effects of imidazo[1,2-a]pyridine-based derivatives increased in the order of phenyl <adamantyl < biphenyl side chain substituted at the position R 1 . Those derivatives with methyl substituent at the position R 2 of the imidazo[1,2-a]pyridine ring exhibited weaker AChE and BChE inhibitory effects. The molecular docking studies of the imidazo[1,2-a]pyridine-based derivatives revealed that the π···π or halogen···O interaction between the active compounds and their surrounding residues in the peripheral anionic site or acyl pocket of the target enzymes are dominant as compared to the hydrophobic interaction in the catalytic active site.

Data Availability
The crystallographic data for 2(a, b, c and e-o) were deposited with the Cambridge Crystallographic Data Centre as supplementary publications with CCDC no. 1055984-1055988 and 1541130-1541138, respectively. Copies of available materials can be obtained free of charge on application to CCDC, 12 Union Road, Cambridge CB2 1EZ, UK, (Fax: + 44-(0)1223-336033 or e-mail: deposit@ccdc.cam.ac.uk).