Design, Synthesis, and In Vitro, In Silico and In Cellulo Evaluation of New Pyrimidine and Pyridine Amide and Carbamate Derivatives as Multi-Functional Cholinesterase Inhibitors

Alzheimer disease is an age-linked neurodegenerative disorder representing one of the greatest medical care challenges of our century. Several drugs are useful in ameliorating the symptoms, even if none could stop or reverse disease progression. The standard approach is represented by the cholinesterase inhibitors (ChEIs) that restore the levels of acetylcholine (ACh) by inhibiting the acetylcholinesterase (AChE). Still, their limited efficacy has prompted researchers to develop new ChEIs that could also reduce the oxidative stress by exhibiting antioxidant properties and by chelating the main metals involved in the disease. Recently, we developed some derivatives constituted by a 2-amino-pyrimidine or a 2-amino-pyridine moiety connected to various aromatic groups by a flexible amino-alkyl linker as new dual inhibitors of AChE and butyrylcholinesterase (BChE). Following our previous studies, in this work we explored the role of the flexible linker by replacing the amino group with an amide or a carbamic group. The most potent compounds showed higher selectivity against BChE in respect to AChE, proving also to possess a weak anti-aggregating activity toward Aβ42 and tau and to be able to chelate Cu2+ and Fe3+ ions. Molecular docking and molecular dynamic studies proposed possible binding modes with the enzymes. It is noteworthy that these compounds were predicted as BBB-permeable and showed low cytotoxicity on the human brain cell line.


Introduction
Alzheimer's disease (AD) is a neurodegenerative disorder with a chronic and progressive course, exhibited by memory loss and impaired cognitive abilities. As is well known, the major risk factor associated with the onset of dementia is age. The elderly population is constantly growing in the world and life expectancy is increasing at a steady rate [1]. The weight of this aging contributes to confirming the estimates of numerous international epidemiological studies that count 55 million cases of people with dementia in the world and predict, in 2030, a growth of more than 78 million, which could add, in the next twenty years, a figure exceeding 139 million people, with the vast majority concentrated in developing countries [2].
The pathogenesis of AD is complex and, since the currently available therapy brings limited benefits, finding new molecules that act on more specific targets is one of the most promising strategies in the search for effective treatments [3][4][5][6]. The disease is characterized by intracellular neurofibrillary tangles, due to the hyperphosphorylation and aggregation of tau protein, and by senile plaques, extracellular deposits of neurotoxic insoluble fibrils of β-amyloid peptide (Aβ) [7][8][9]. The latter is a metalloprotein with high binding affinity for biometals essential for cellular homeostasis, such as Cu 2+ , Zn 2+ and Fe 3+ , and high standard redox potential that involves the reduction of these ions and the trapping of molecular oxygen, giving rise to radical species causing oxidation of proteins, DNA and RNA and increased lipid peroxidation [10].
These lesions severely impair cholinergic transmission, so the pharmacological approach is primarily based on restoring acetylcholine (ACh) levels through the use of cholinesterase inhibitors (ChEIs). In the enzymatic pocket of cholinesterase (ChEs) there are two binding sites: the catalytic active site (CAS) and the peripheral anionic site (PAS) [11][12][13]. The involvement of acetylcholinesterase (AChE), through PAS, in the fibrillar aggregation of Aβ has given new meaning to the therapeutic use of AChE inhibitors with mixed or non-competitive mechanisms, as they act both by restoring ACh levels and interfering with Aβ aggregation [14][15][16]. Several studies also reported the role of butyrylcholinesterase (BChE) in the deposition of Aβ plaques and in the compensation of AChE function in ongoing AD patients [17][18][19][20].
Other important factors in the complex etiology of AD are metal ions (mainly copper, zinc and iron), probably considered responsible for the protein misfolding, and the oxidative stress [21][22][23].
During the last few years, scientific advances have inspired researchers to study new pharmacological therapies focused more on the pathophysiological events of the AD [24]. In particular, one of the central choices of medicinal chemistry research is the development of multitarget or multifunctional ChEIs [25][26][27]. Some of these studies suggest that a possible way forward could be the development of new therapeutic agents with anticholinesterasic activity that also mitigate oxidative stress both directly, by exerting antioxidant functions, and indirectly, through the chelation of the main metal ions involved in the generation of reactive oxygen species [28,29].
In our previous works, we synthetized and studied new molecules, characterized by two small aromatic moieties separated by various functionalized linkers, designed as new dual binding AChE and BChE inhibitors, able to bind both CAS and PAS of the enzymes, as well as being endowed with other properties, such as anti-aggregating activity, chelating and antioxidant abilities, which are potentially interesting in the light of the multifunctional nature of AD [30][31][32]. From these molecules, we recently developed some derivatives constituted by a 2-amino-pyrimidine or a 2-amino-pyridine moiety connected to various aromatic groups by a flexible amino-alkyl linker [32]. In this way, we obtained mixed ChEs inhibitors, in which the two aromatic moieties and the amine group inserted into the linear alkyl chain are able to interact with the amino acids located in the CAS, PAS and mid-gorge cavities by means of π−π stacking and π−cation interactions, acting as dual binding site inhibitors, as highlighted by in silico analyses. Moreover, the 2-amino-pyrimidinic or 2-amino-pyridinic moieties confer chelating activity to the compounds, due to the presence of the two adjacent nitrogen atoms. Noteworthy, these molecules are also characterized by low toxicity toward human brain cells and good predicted physicochemical and ADME parameters.
As a development of these works, we designed new pyrimidine and pyridine derivatives, replacing the amino group inserted into the alkyl chain with an amide or a carbamic group (Figure 1), and we studied these molecules in vitro, in silico and in cellulo to assay their ChEs inhibitory activity, chelating ability and anti-aggregating properties. The amide or carbamic groups were chosen to confer a major rigidity to the alkyl chain, as well as to form hydrogen bonds with the amino acids of the enzymes, with the aim to obtain a better and more stable orientation of the two aromatic moieties into the enzymatic gorge. Moreover, the carbamic group could establish covalent interaction with the serine of the catalytic site, forming a carbamoyl-cholinesterase complex like other ChEs carbamoylating inhibitors, contributing directly to the inhibitory action of the molecules. Based on the previous results, the length of the alkyl chain was chosen was six methylene units, to cover the distance between CAS and PAS. Among the aromatic groups introduced on one side of the aliphatic chain, besides the unsubstituted or substituted phenyl rings already studied for the amine derivatives of the previous work, bulkier groups, such as 4-phenoxybezyl or cinnamic one, were selected with the purpose to study the effects of an increase in the size of the molecules in the capacity to inhibit both ChEs or, selectively, BChE. due to the presence of the two adjacent nitrogen atoms. Noteworthy, these molecules are also characterized by low toxicity toward human brain cells and good predicted physicochemical and ADME parameters. As a development of these works, we designed new pyrimidine and pyridine derivatives, replacing the amino group inserted into the alkyl chain with an amide or a carbamic group (Figure 1), and we studied these molecules in vitro, in silico and in cellulo to assay their ChEs inhibitory activity, chelating ability and anti-aggregating properties. The amide or carbamic groups were chosen to confer a major rigidity to the alkyl chain, as well as to form hydrogen bonds with the amino acids of the enzymes, with the aim to obtain a better and more stable orientation of the two aromatic moieties into the enzymatic gorge. Moreover, the carbamic group could establish covalent interaction with the serine of the catalytic site, forming a carbamoyl-cholinesterase complex like other ChEs carbamoylating inhibitors, contributing directly to the inhibitory action of the molecules. Based on the previous results, the length of the alkyl chain was chosen was six methylene units, to cover the distance between CAS and PAS. Among the aromatic groups introduced on one side of the aliphatic chain, besides the unsubstituted or substituted phenyl rings already studied for the amine derivatives of the previous work, bulkier groups, such as 4-phenoxybezyl or cinnamic one, were selected with the purpose to study the effects of an increase in the size of the molecules in the capacity to inhibit both ChEs or, selectively, BChE. The synthesized compounds were screened against Electrophorus electricus AChE (EeAChE) and equine BChE (eqBChE) by means of Ellman's spectrophotometric method in order to evaluate their inhibition potency and selectivity. The studied compounds appear more selective towards eqBChE than EeAChE, and the most potent eqBChE inhibitors were further studied to identify inhibition mechanism and inhibition constants. Molecular docking and molecular dynamic studies were carried out with the aim of identifying the possible interactions among these molecules and the amino acids located in the CAS, PAS and mid-gorge region of both enzymes and to explain the selectivity towards eqBChE. Further studies were conducted on the most active compounds with the aim to evaluate their Cu 2+ and Fe 3+ chelation ability by means of UV-Vis spectroscopic analyses, while in cellulo studies on E. coli cells were carried out to reveal their anti-aggregating activity toward Aβ42 and tau. Finally, for these compounds, the ability to permeate the BBB was predicted by in silico simulation and the cytotoxicity on a human brain cell line was also evaluated. The synthesized compounds were screened against Electrophorus electricus AChE (EeAChE) and equine BChE (eqBChE) by means of Ellman's spectrophotometric method in order to evaluate their inhibition potency and selectivity. The studied compounds appear more selective towards eqBChE than EeAChE, and the most potent eqBChE inhibitors were further studied to identify inhibition mechanism and inhibition constants. Molecular docking and molecular dynamic studies were carried out with the aim of identifying the possible interactions among these molecules and the amino acids located in the CAS, PAS and mid-gorge region of both enzymes and to explain the selectivity towards eqBChE. Further studies were conducted on the most active compounds with the aim to evaluate their Cu 2+ and Fe 3+ chelation ability by means of UV-Vis spectroscopic analyses, while in cellulo studies on E. coli cells were carried out to reveal their anti-aggregating activity toward Aβ 42 and tau. Finally, for these compounds, the ability to permeate the BBB was predicted by in silico simulation and the cytotoxicity on a human brain cell line was also evaluated. Scheme 1. Synthesis of pyrimidine and pyridine intermediates 1 and 2. Reagents and conditions: a TEA (1 eq), CH 3 OH, reflux, 18h; b CuI (0.05 eq), Cs 2 CO 3 (2 eq), 2-isobutyrylcyclohexanone (0.2 eq), DMF, N 2 , 40 • C, 21 h; c TFA (20 eq), CH 2 Cl 2 , rt, 3.5 h.
The pyrimidine and pyridine amide compounds 3-15 were prepared as described in Scheme 2. Initially, the opportune carboxylic acid was activated using carbonyldiimidazole (CDI) in acetonitrile at reflux; the reaction was monitored through IR spectroscopy, evaluating the appearance of two carbonyl stretching bands around 1735 and 1700 cm −1 , concerning the CDI-carboxylic acid intermediate. At this point, pyrimidine intermediate 1 or pyridine intermediate 2 was added and the reaction was refluxed overnight. The compounds were purified by column chromatography on silica gel and/or by crystallisation. The structures of final compounds were confirmed by spectroscopic analysis; as an example, in the 1 H-NMR spectra of amide derivatives the shifts of the signals due to the -CH 2 -NH-groups were observed from 2.70 ppm (-CH 2 -NH 2 in intermediates 1 and 2) to 3.38-2.99 ppm (-CH 2 -NH-C=O in amide derivatives).
The pyrimidine and pyridine amide compounds 3-15 were prepared as described in Scheme 2. Initially, the opportune carboxylic acid was activated using carbonyldiimidazole (CDI) in acetonitrile at reflux; the reaction was monitored through IR spectroscopy, evaluating the appearance of two carbonyl stretching bands around 1735 and 1700 cm −1 , concerning the CDI-carboxylic acid intermediate. At this point, pyrimidine intermediate 1 or pyridine intermediate 2 was added and the reaction was refluxed overnight. The compounds were purified by column chromatography on silica gel and/or by crystallisation. The structures of final compounds were confirmed by spectroscopic analysis; as an example, in the 1 H-NMR spectra of amide derivatives the shifts of the signals due to the -CH2-NH-groups were observed from 2.70 ppm (-CH2-NH2 in intermediates 1 and 2) to 3.38-2.99 ppm (-CH2-NH-C=O in amide derivatives).
The detailed synthetic procedures, the analytical and spectroscopic data of the synthesized compounds, are reported in the experimental section and agree with the proposed structures.

Enzymatic Assays
Regarding the synthesized compounds enzymatic inhibition studies were carried out on Electrophorus electricus AChE (EeAChE) and equine BChE (eqBChE), according to Ellman's spectrophotometric method [33]. Initially, for each compound, except those insoluble under experimental conditions (see Table 1), the percentages of inhibition were determined at the inhibitor concentration equal to 9 µ M and, for the most potent compounds, also at 900 nM, in the presence of 0.0833 U/mL of enzyme and 100 µ M of acetylthiocholine as substrate. The detailed synthetic procedures, the analytical and spectroscopic data of the synthesized compounds, are reported in the experimental section and agree with the proposed structures.

Enzymatic Assays
Regarding the synthesized compounds enzymatic inhibition studies were carried out on Electrophorus electricus AChE (EeAChE) and equine BChE (eqBChE), according to Ellman's spectrophotometric method [33]. Initially, for each compound, except those insoluble under experimental conditions (see Table 1), the percentages of inhibition were determined at the inhibitor concentration equal to 9 µM and, for the most potent compounds, also at 900 nM, in the presence of 0.0833 U/mL of enzyme and 100 µM of acetylthiocholine as substrate.
In Table 1 the percentages of inhibition of amide derivatives 3-15 are reported. These data show that pyrimidine amide derivatives 3-11 are weak inhibitors of both enzymes, EeAChE and eqBChE. If compared with the corresponding pyrimidine amino compounds previously investigated [32], it is evident that the replacement of the amino group with the amide has led to a drastic reduction in the inhibitory potency towards both cholinesterases. As an example, pyrimidine amide 3, at a concentration equal to 9 µM, shows about 3% and 9% of inhibition, respectively, towards EeAChE and eqBChE, while its corresponding diamine derivative, at the same concentration, showed 93% of inhibition towards EeAChE and 26% towards eqBChE. Ellman's spectrophotometric method [33]. Initially, for each compound, except those insoluble under experimental conditions (see Table 1), the percentages of inhibition were determined at the inhibitor concentration equal to 9 µ M and, for the most potent compounds, also at 900 nM, in the presence of 0.0833 U/mL of enzyme and 100 µ M of acetylthiocholine as substrate. In Table 1 the percentages of inhibition of amide derivatives 3-15 are reported. These data show that pyrimidine amide derivatives 3-11 are weak inhibitors of both enzymes, EeAChE and eqBChE. If compared with the corresponding pyrimidine amino compounds previously investigated [32], it is evident that the replacement of the amino group with the amide has led to a drastic reduction in the inhibitory potency towards both cholinesterases. As an example, pyrimidine amide 3, at a concentration equal to 9 µ M, shows about 3% and 9% of inhibition, respectively, towards EeAChE and eqBChE, while its corresponding diamine derivative, at the same concentration, showed 93% of inhibition towards EeAChE and 26% towards eqBChE.
The pyridine amide derivatives 12-15 have percentages of inhibition significantly higher on both cholinesterases than the corresponding pyrimidine amide compounds 3- In Table 1 the percentages of inhibition of amide derivatives 3-15 are reported. These data show that pyrimidine amide derivatives 3-11 are weak inhibitors of both enzymes, EeAChE and eqBChE. If compared with the corresponding pyrimidine amino compounds previously investigated [32], it is evident that the replacement of the amino group with the amide has led to a drastic reduction in the inhibitory potency towards both cholinesterases. As an example, pyrimidine amide 3, at a concentration equal to 9 µ M, shows about 3% and 9% of inhibition, respectively, towards EeAChE and eqBChE, while its corresponding diamine derivative, at the same concentration, showed 93% of inhibition towards EeAChE and 26% towards eqBChE.
The pyridine amide derivatives 12-15 have percentages of inhibition significantly higher on both cholinesterases than the corresponding pyrimidine amide compounds 3-  In Table 1 the percentages of inhibition of amide derivatives 3-15 are reported. These data show that pyrimidine amide derivatives 3-11 are weak inhibitors of both enzymes, EeAChE and eqBChE. If compared with the corresponding pyrimidine amino compounds previously investigated [32], it is evident that the replacement of the amino group with the amide has led to a drastic reduction in the inhibitory potency towards both cholinesterases. As an example, pyrimidine amide 3, at a concentration equal to 9 µ M, shows about 3% and 9% of inhibition, respectively, towards EeAChE and eqBChE, while its corresponding diamine derivative, at the same concentration, showed 93% of inhibition towards EeAChE and 26% towards eqBChE.
The pyridine amide derivatives 12-15 have percentages of inhibition significantly higher on both cholinesterases than the corresponding pyrimidine amide compounds 3- In Table 1 the percentages of inhibition of amide derivatives 3-15 are reported. These data show that pyrimidine amide derivatives 3-11 are weak inhibitors of both enzymes, EeAChE and eqBChE. If compared with the corresponding pyrimidine amino compounds previously investigated [32], it is evident that the replacement of the amino group with the amide has led to a drastic reduction in the inhibitory potency towards both cholinesterases. As an example, pyrimidine amide 3, at a concentration equal to 9 µ M, shows about 3% and 9% of inhibition, respectively, towards EeAChE and eqBChE, while its corresponding diamine derivative, at the same concentration, showed 93% of inhibition towards EeAChE and 26% towards eqBChE.
The pyridine amide derivatives 12- 15  In Table 1 the percentages of inhibition of amide derivatives 3-15 are reported. These data show that pyrimidine amide derivatives 3-11 are weak inhibitors of both enzymes, EeAChE and eqBChE. If compared with the corresponding pyrimidine amino compounds previously investigated [32], it is evident that the replacement of the amino group with the amide has led to a drastic reduction in the inhibitory potency towards both cholinesterases. As an example, pyrimidine amide 3, at a concentration equal to 9 µ M, shows about 3% and 9% of inhibition, respectively, towards EeAChE and eqBChE, while its corresponding diamine derivative, at the same concentration, showed 93% of inhibition towards EeAChE and 26% towards eqBChE.
The pyridine amide derivatives 12-15 have percentages of inhibition significantly  In Table 1 the percentages of inhibition of amide derivatives 3-15 are reported. These data show that pyrimidine amide derivatives 3-11 are weak inhibitors of both enzymes, EeAChE and eqBChE. If compared with the corresponding pyrimidine amino compounds previously investigated [32], it is evident that the replacement of the amino group with the amide has led to a drastic reduction in the inhibitory potency towards both cholinesterases. As an example, pyrimidine amide 3, at a concentration equal to 9 µ M, shows about 3% and 9% of inhibition, respectively, towards EeAChE and eqBChE, while its corresponding diamine derivative, at the same concentration, showed 93% of inhibition towards EeAChE and 26% towards eqBChE.
The pyridine amide derivatives 12-15 have percentages of inhibition significantly  In Table 1 the percentages of inhibition of amide derivatives 3-15 are reported. These data show that pyrimidine amide derivatives 3-11 are weak inhibitors of both enzymes, EeAChE and eqBChE. If compared with the corresponding pyrimidine amino compounds previously investigated [32], it is evident that the replacement of the amino group with the amide has led to a drastic reduction in the inhibitory potency towards both cholinesterases. As an example, pyrimidine amide 3, at a concentration equal to 9 µ M, shows about 3% and 9% of inhibition, respectively, towards EeAChE and eqBChE, while its corresponding diamine derivative, at the same concentration, showed 93% of inhibition towards EeAChE and 26% towards eqBChE.
The pyridine amide derivatives 12-15 have percentages of inhibition significantly  In Table 1 the percentages of inhibition of amide derivatives 3-15 are reported. These data show that pyrimidine amide derivatives 3-11 are weak inhibitors of both enzymes, EeAChE and eqBChE. If compared with the corresponding pyrimidine amino compounds previously investigated [32], it is evident that the replacement of the amino group with the amide has led to a drastic reduction in the inhibitory potency towards both cholinesterases. As an example, pyrimidine amide 3, at a concentration equal to 9 µ M, shows about 3% and 9% of inhibition, respectively, towards EeAChE and eqBChE, while its corresponding diamine derivative, at the same concentration, showed 93% of inhibition towards EeAChE and 26% towards eqBChE.
The pyridine amide derivatives 12- 15  In Table 1 the percentages of inhibition of amide derivatives 3-15 are reported. These data show that pyrimidine amide derivatives 3-11 are weak inhibitors of both enzymes, EeAChE and eqBChE. If compared with the corresponding pyrimidine amino compounds previously investigated [32], it is evident that the replacement of the amino group with the amide has led to a drastic reduction in the inhibitory potency towards both cholinesterases. As an example, pyrimidine amide 3, at a concentration equal to 9 µ M, shows about 3% and 9% of inhibition, respectively, towards EeAChE and eqBChE, while its corresponding diamine derivative, at the same concentration, showed 93% of inhibition towards EeAChE and 26% towards eqBChE. In Table 1 the percentages of inhibition of amide derivatives 3-15 are reported. These data show that pyrimidine amide derivatives 3-11 are weak inhibitors of both enzymes, EeAChE and eqBChE. If compared with the corresponding pyrimidine amino compounds previously investigated [32], it is evident that the replacement of the amino group with the amide has led to a drastic reduction in the inhibitory potency towards both cholinesterases. As an example, pyrimidine amide 3, at a concentration equal to 9 µ M, shows about 3% and 9% of inhibition, respectively, towards EeAChE and eqBChE, while its corresponding diamine derivative, at the same concentration, showed 93% of inhibition towards EeAChE and 26% towards eqBChE. CH 34.0 ± 2.1 nd b 80.5 ± 1. In Table 1 the percentages of inhibition of amide derivatives 3-15 are reported. These data show that pyrimidine amide derivatives 3-11 are weak inhibitors of both enzymes, EeAChE and eqBChE. If compared with the corresponding pyrimidine amino compounds previously investigated [32], it is evident that the replacement of the amino group with the amide has led to a drastic reduction in the inhibitory potency towards both cholinesterases. As an example, pyrimidine amide 3, at a concentration equal to 9 µ M, shows about 3% and 9% of inhibition, respectively, towards EeAChE and eqBChE, while its corresponding diamine derivative, at the same concentration, showed 93% of inhibition towards EeAChE and 26% towards eqBChE.  In Table 1 the percentages of inhibition of amide derivatives 3-15 are reported. These data show that pyrimidine amide derivatives 3-11 are weak inhibitors of both enzymes, EeAChE and eqBChE. If compared with the corresponding pyrimidine amino compounds previously investigated [32], it is evident that the replacement of the amino group with the amide has led to a drastic reduction in the inhibitory potency towards both cholinesterases. As an example, pyrimidine amide 3, at a concentration equal to 9 µ M, shows about 3% and 9% of inhibition, respectively, towards EeAChE and eqBChE, while its corresponding diamine derivative, at the same concentration, showed 93% of inhibition towards EeAChE and 26% towards eqBChE.  In Table 1 the percentages of inhibition of amide derivatives 3-15 are reported. These data show that pyrimidine amide derivatives 3-11 are weak inhibitors of both enzymes, EeAChE and eqBChE. If compared with the corresponding pyrimidine amino compounds previously investigated [32], it is evident that the replacement of the amino group with the amide has led to a drastic reduction in the inhibitory potency towards both cholinesterases. As an example, pyrimidine amide 3, at a concentration equal to 9 µ M, shows about 3% and 9% of inhibition, respectively, towards EeAChE and eqBChE, while its corresponding diamine derivative, at the same concentration, showed 93% of inhibition towards EeAChE and 26% towards eqBChE. CH 45.4 ± 2.5 nd b 64.7 ± 2.9 31.5 ± 3.1 Tacrine  100 97.6 ± 0.1 100 99.7 ± 0.3 a Data are the average of three replicates; b nd stands for not determined; c ns stands for not soluble; d na stands for not active.
The pyridine amide derivatives 12-15 have percentages of inhibition significantly higher on both cholinesterases than the corresponding pyrimidine amide compounds 3-6. In particular, these molecules are more potent towards eqBChE than towards EeAChE. The comparison between pyridine amide 12 and the corresponding diamine compound [32]  shows that, in this case, the replacement of the amino group with the amide reduced the inhibitory potency towards EeAChE (73% for amine derivatives and 34% for 12, both at 9 µM), but preserved the activity towards eqBChE (72% for amine derivatives and 81% for 12, both at 9 µM), leading to a more selective BChE inhibitor. Among pyridine amides 12-15, the most interesting for their inhibitory activity towards eqBChE are 12 and 14, respectively, with a phenyl or 4-phenoxybenzyl ring adjacent to the carbonyl group. Compound 12 shows an 81% inhibition at 9 µM and 22% at 900 nM, while 14 has an 88% inhibition at 9 µM and 45% at 900 nM.
For carbamate derivatives 16-19, whose data are reported in Table 2, a time-dependent inhibition assay was carried out to evaluate if they could act as carbamoylating inhibitor, forming a carbamoyl-cholinesterase complex by covalent interaction with the serine in the catalytic site. For this purpose, the percentages of inhibition were measured both at a timepoint of zero and following an incubation of 1 h. An increase in the percentages of inhibition towards both cholinesterases was observed after the incubation period, suggesting that these compounds could act as carbamoylating agents. Moreover, these derivatives are more selective towards eqBChE than EeAChE. In eqBChE, a more significant increase in the percentage of inhibition was observed following the incubation of the phenylcarbamates 16 and 18, compared to the corresponding benzylcarbamates 17 and 19; this could be explained by the fact that the phenate ion, being a better leaving group in comparison to the benzyl alcohol ion, can more easily allow the carbamoylation of the serine in the catalytic site. Moreover, in this case, as for amide derivatives, pyridine carbamates 18 and 19 show greater inhibitory potency on both cholinesterases, both at zero time and after incubation, compared to the corresponding pyrimidine derivatives 16 and 17. the inhibitory potency towards EeAChE (73% for amine derivatives and 34% for 12, both at 9 µ M), but preserved the activity towards eqBChE (72% for amine derivatives and 81% for 12, both at 9 µ M), leading to a more selective BChE inhibitor. Among pyridine amides 12-15, the most interesting for their inhibitory activity towards eqBChE are 12 and 14, respectively, with a phenyl or 4-phenoxybenzyl ring adjacent to the carbonyl group. Compound 12 shows an 81% inhibition at 9 µ M and 22% at 900 nM, while 14 has an 88% inhibition at 9 µ M and 45% at 900 nM. For carbamate derivatives 16-19, whose data are reported in Table 2, a time-dependent inhibition assay was carried out to evaluate if they could act as carbamoylating inhibitor, forming a carbamoyl-cholinesterase complex by covalent interaction with the serine in the catalytic site. For this purpose, the percentages of inhibition were measured both at a time-point of zero and following an incubation of 1 h. An increase in the percentages of inhibition towards both cholinesterases was observed after the incubation period, suggesting that these compounds could act as carbamoylating agents. Moreover, these derivatives are more selective towards eqBChE than EeAChE. In eqBChE, a more significant increase in the percentage of inhibition was observed following the incubation of the phenylcarbamates 16 and 18, compared to the corresponding benzylcarbamates 17 and 19; this could be explained by the fact that the phenate ion, being a better leaving group in comparison to the benzyl alcohol ion, can more easily allow the carbamoylation of the serine in the catalytic site. Moreover, in this case, as for amide derivatives, pyridine carbamates 18 and 19 show greater inhibitory potency on both cholinesterases, both at zero time and after incubation, compared to the corresponding pyrimidine derivatives 16 and 17. In general, these carbamates are better inhibitors of both cholinesterases if compared with the corresponding amide compounds. On the other hand, if compared with the corresponding pyrimidine and pyridine amino derivatives previously investigated [32], the replacement of the amino group with the carbamic one decreases the inhibitory potency towards EeAChE but enhances it towards eqBChE, leading to more selective BChE inhibitors.
To monitor the progressive increase over time of the eqBChE inhibition exerted by compound 18, the most potent eqBChE inhibitor among these carbamate derivatives, the percentages of inhibition were determined every fifteen minutes for two hours, at the inhibitor concentration equal to 700 nM, in presence of 0.0833 U/mL of enzyme and 100 µ M the inhibitory potency towards EeAChE (73% for amine derivatives and 34% for 12, both at 9 µ M), but preserved the activity towards eqBChE (72% for amine derivatives and 81% for 12, both at 9 µ M), leading to a more selective BChE inhibitor. Among pyridine amides 12-15, the most interesting for their inhibitory activity towards eqBChE are 12 and 14, respectively, with a phenyl or 4-phenoxybenzyl ring adjacent to the carbonyl group. Compound 12 shows an 81% inhibition at 9 µ M and 22% at 900 nM, while 14 has an 88% inhibition at 9 µ M and 45% at 900 nM. For carbamate derivatives 16-19, whose data are reported in Table 2, a time-dependent inhibition assay was carried out to evaluate if they could act as carbamoylating inhibitor, forming a carbamoyl-cholinesterase complex by covalent interaction with the serine in the catalytic site. For this purpose, the percentages of inhibition were measured both at a time-point of zero and following an incubation of 1 h. An increase in the percentages of inhibition towards both cholinesterases was observed after the incubation period, suggesting that these compounds could act as carbamoylating agents. Moreover, these derivatives are more selective towards eqBChE than EeAChE. In eqBChE, a more significant increase in the percentage of inhibition was observed following the incubation of the phenylcarbamates 16 and 18, compared to the corresponding benzylcarbamates 17 and 19; this could be explained by the fact that the phenate ion, being a better leaving group in comparison to the benzyl alcohol ion, can more easily allow the carbamoylation of the serine in the catalytic site. Moreover, in this case, as for amide derivatives, pyridine carbamates 18 and 19 show greater inhibitory potency on both cholinesterases, both at zero time and after incubation, compared to the corresponding pyrimidine derivatives 16 and 17. In general, these carbamates are better inhibitors of both cholinesterases if compared with the corresponding amide compounds. On the other hand, if compared with the corresponding pyrimidine and pyridine amino derivatives previously investigated [32], the replacement of the amino group with the carbamic one decreases the inhibitory potency towards EeAChE but enhances it towards eqBChE, leading to more selective BChE inhibitors.
To monitor the progressive increase over time of the eqBChE inhibition exerted by compound 18, the most potent eqBChE inhibitor among these carbamate derivatives, the percentages of inhibition were determined every fifteen minutes for two hours, at the inhibitor concentration equal to 700 nM, in presence of 0.0833 U/mL of enzyme and 100 µ M the inhibitory potency towards EeAChE (73% for amine derivatives and 34% for 12, both at 9 µ M), but preserved the activity towards eqBChE (72% for amine derivatives and 81% for 12, both at 9 µ M), leading to a more selective BChE inhibitor. Among pyridine amides 12-15, the most interesting for their inhibitory activity towards eqBChE are 12 and 14, respectively, with a phenyl or 4-phenoxybenzyl ring adjacent to the carbonyl group. Compound 12 shows an 81% inhibition at 9 µ M and 22% at 900 nM, while 14 has an 88% inhibition at 9 µ M and 45% at 900 nM. For carbamate derivatives 16-19, whose data are reported in Table 2, a time-dependent inhibition assay was carried out to evaluate if they could act as carbamoylating inhibitor, forming a carbamoyl-cholinesterase complex by covalent interaction with the serine in the catalytic site. For this purpose, the percentages of inhibition were measured both at a time-point of zero and following an incubation of 1 h. An increase in the percentages of inhibition towards both cholinesterases was observed after the incubation period, suggesting that these compounds could act as carbamoylating agents. Moreover, these derivatives are more selective towards eqBChE than EeAChE. In eqBChE, a more significant increase in the percentage of inhibition was observed following the incubation of the phenylcarbamates 16 and 18, compared to the corresponding benzylcarbamates 17 and 19; this could be explained by the fact that the phenate ion, being a better leaving group in comparison to the benzyl alcohol ion, can more easily allow the carbamoylation of the serine in the catalytic site. Moreover, in this case, as for amide derivatives, pyridine carbamates 18 and 19 show greater inhibitory potency on both cholinesterases, both at zero time and after incubation, compared to the corresponding pyrimidine derivatives 16 and 17. In general, these carbamates are better inhibitors of both cholinesterases if compared with the corresponding amide compounds. On the other hand, if compared with the corresponding pyrimidine and pyridine amino derivatives previously investigated [32], the replacement of the amino group with the carbamic one decreases the inhibitory potency towards EeAChE but enhances it towards eqBChE, leading to more selective BChE inhibitors.
To monitor the progressive increase over time of the eqBChE inhibition exerted by compound 18, the most potent eqBChE inhibitor among these carbamate derivatives, the percentages of inhibition were determined every fifteen minutes for two hours, at the inhibitor concentration equal to 700 nM, in presence of 0.0833 U/mL of enzyme and 100 µ M the inhibitory potency towards EeAChE (73% for amine derivatives and 34% for 12, both at 9 µ M), but preserved the activity towards eqBChE (72% for amine derivatives and 81% for 12, both at 9 µ M), leading to a more selective BChE inhibitor. Among pyridine amides 12-15, the most interesting for their inhibitory activity towards eqBChE are 12 and 14, respectively, with a phenyl or 4-phenoxybenzyl ring adjacent to the carbonyl group. Compound 12 shows an 81% inhibition at 9 µ M and 22% at 900 nM, while 14 has an 88% inhibition at 9 µ M and 45% at 900 nM. For carbamate derivatives 16-19, whose data are reported in Table 2, a time-dependent inhibition assay was carried out to evaluate if they could act as carbamoylating inhibitor, forming a carbamoyl-cholinesterase complex by covalent interaction with the serine in the catalytic site. For this purpose, the percentages of inhibition were measured both at a time-point of zero and following an incubation of 1 h. An increase in the percentages of inhibition towards both cholinesterases was observed after the incubation period, suggesting that these compounds could act as carbamoylating agents. Moreover, these derivatives are more selective towards eqBChE than EeAChE. In eqBChE, a more significant increase in the percentage of inhibition was observed following the incubation of the phenylcarbamates 16 and 18, compared to the corresponding benzylcarbamates 17 and 19; this could be explained by the fact that the phenate ion, being a better leaving group in comparison to the benzyl alcohol ion, can more easily allow the carbamoylation of the serine in the catalytic site. Moreover, in this case, as for amide derivatives, pyridine carbamates 18 and 19 show greater inhibitory potency on both cholinesterases, both at zero time and after incubation, compared to the corresponding pyrimidine derivatives 16 and 17. In general, these carbamates are better inhibitors of both cholinesterases if compared with the corresponding amide compounds. On the other hand, if compared with the corresponding pyrimidine and pyridine amino derivatives previously investigated [32], the replacement of the amino group with the carbamic one decreases the inhibitory potency towards EeAChE but enhances it towards eqBChE, leading to more selective BChE inhibitors.
To monitor the progressive increase over time of the eqBChE inhibition exerted by compound 18, the most potent eqBChE inhibitor among these carbamate derivatives, the percentages of inhibition were determined every fifteen minutes for two hours, at the inhibitor concentration equal to 700 nM, in presence of 0.0833 U/mL of enzyme and 100 µ M  the inhibitory potency towards EeAChE (73% for amine derivatives and 34% for 12, both at 9 µ M), but preserved the activity towards eqBChE (72% for amine derivatives and 81% for 12, both at 9 µ M), leading to a more selective BChE inhibitor. Among pyridine amides 12-15, the most interesting for their inhibitory activity towards eqBChE are 12 and 14, respectively, with a phenyl or 4-phenoxybenzyl ring adjacent to the carbonyl group. Compound 12 shows an 81% inhibition at 9 µ M and 22% at 900 nM, while 14 has an 88% inhibition at 9 µ M and 45% at 900 nM. For carbamate derivatives 16-19, whose data are reported in Table 2, a time-dependent inhibition assay was carried out to evaluate if they could act as carbamoylating inhibitor, forming a carbamoyl-cholinesterase complex by covalent interaction with the serine in the catalytic site. For this purpose, the percentages of inhibition were measured both at a time-point of zero and following an incubation of 1 h. An increase in the percentages of inhibition towards both cholinesterases was observed after the incubation period, suggesting that these compounds could act as carbamoylating agents. Moreover, these derivatives are more selective towards eqBChE than EeAChE. In eqBChE, a more significant increase in the percentage of inhibition was observed following the incubation of the phenylcarbamates 16 and 18, compared to the corresponding benzylcarbamates 17 and 19; this could be explained by the fact that the phenate ion, being a better leaving group in comparison to the benzyl alcohol ion, can more easily allow the carbamoylation of the serine in the catalytic site. Moreover, in this case, as for amide derivatives, pyridine carbamates 18 and 19 show greater inhibitory potency on both cholinesterases, both at zero time and after incubation, compared to the corresponding pyrimidine derivatives 16 and 17. In general, these carbamates are better inhibitors of both cholinesterases if compared with the corresponding amide compounds. On the other hand, if compared with the corresponding pyrimidine and pyridine amino derivatives previously investigated [32], the replacement of the amino group with the carbamic one decreases the inhibitory potency towards EeAChE but enhances it towards eqBChE, leading to more selective BChE inhibitors.
To monitor the progressive increase over time of the eqBChE inhibition exerted by compound 18, the most potent eqBChE inhibitor among these carbamate derivatives, the percentages of inhibition were determined every fifteen minutes for two hours, at the inhibitor concentration equal to 700 nM, in presence of 0.0833 U/mL of enzyme and 100 µ M  In general, these carbamates are better inhibitors of both cholinesterases if compared with the corresponding amide compounds. On the other hand, if compared with the corresponding pyrimidine and pyridine amino derivatives previously investigated [32], the replacement of the amino group with the carbamic one decreases the inhibitory potency towards EeAChE but enhances it towards eqBChE, leading to more selective BChE inhibitors.
To monitor the progressive increase over time of the eqBChE inhibition exerted by compound 18, the most potent eqBChE inhibitor among these carbamate derivatives, the percentages of inhibition were determined every fifteen minutes for two hours, at the  For compounds 12 and 14, the most potent amide derivatives, such as the eqBChE inhibitors, the inhibition constant (Ki) and the corresponding inhibition mechanism, were determined according to Dixon's method [34], and the reciprocal of the hydrolysis rate versus the inhibitor concentrations at a fixed concentration of substrate were reported in graph. The recorded data were analysed with the enzyme kinetic module of SigmaPlot in order to find the best fitting model of inhibition, using the linear regression analysis. The reference kinetic models used in the regression analysis were: competitive, non-competitive, uncompetitive and mixed. Each determination was repeated five times and incorrect values were discarded to reduce the standard deviation to within the limit of 5%. In this way, the regression lines obtained have a linear regression coefficient (R 2 ) higher than 0.98. The amide derivatives 12 and 14 revealed competitive inhibition mechanism towards eqBChE, with Ki respectively equal to 2.988 ± 0.190 µ M and 0.621 ± 0.043 µ M ( Table 3). The Dixon's plots of tested compounds are reported in the Supplementary Materials ( Figures  S1 and S2).  For compounds 12 and 14, the most potent amide derivatives, such as the eqBChE inhibitors, the inhibition constant (K i ) and the corresponding inhibition mechanism, were determined according to Dixon's method [34], and the reciprocal of the hydrolysis rate versus the inhibitor concentrations at a fixed concentration of substrate were reported in graph. The recorded data were analysed with the enzyme kinetic module of SigmaPlot in order to find the best fitting model of inhibition, using the linear regression analysis. The reference kinetic models used in the regression analysis were: competitive, non-competitive, uncompetitive and mixed. Each determination was repeated five times and incorrect values were discarded to reduce the standard deviation to within the limit of 5%. In this way, the regression lines obtained have a linear regression coefficient (R 2 ) higher than 0.98. The amide derivatives 12 and 14 revealed competitive inhibition mechanism towards eqBChE, with K i respectively equal to 2.988 ± 0.190 µM and 0.621 ± 0.043 µM ( Table 3). The Dixon's plots of tested compounds are reported in the Supplementary Materials ( Figures S1 and S2). For pyridine phenylcarbamate 18 the IC 50 on eqBChE, following incubation for one hour, was determined by virtue of a time dependent inhibition. The IC 50 value was obtained by plotting the percentages of inhibition towards eqBChE versus the concentration of inhibitor expressed in logarithmic scale, at fixed substrate concentration (100 µM). The recorded data were analyzed with the enzyme kinetic module of SigmaPlot (Systat Software, Palo alto, CA, USA). Each measurement was replicated three times and the IC 50 value  (Table 4). The IC 50 plot of the tested compound is reported in the Supplementary Materials ( Figure S3). Moreover, it was evaluated whether the inhibition mechanism of carbamate 18 towards eqBChE is reversible, irreversible or pseudo-irreversible. For this purpose, the percentages of inhibition were determined both after one hour of the incubation of compound 18 (5.4 µM) with eqBChE (1.25 U), both following the removal of the inhibitor and the washing of the enzyme, carried out through ultrafiltration devices for the concentration of biological samples [35,36]. To validate the method, the same study was performed with carbaryl (54 µM), a known carbamoylating cholinesterase inhibitor.
After the incubation phase, the solution containing the inhibitor not bound to the enzyme was removed by centrifugation, retaining eqBChE on the membrane of the ultrafiltration devices. Then the enzyme was washed with phosphate buffer, to remove the bound inhibitor and to restore the enzyme activity. For carbamate 18, a partial restoration of the enzymatic activity after washing was recorded, with a reduction in the percentage of inhibition from 89.7 to 56.7%. For carbaryl, a greater restoration of the enzymatic activity was recorded, with a reduction in the percentage of inhibition from 97.4 to 29.9%. Therefore, a pseudo-irreversible inhibition mechanism was observed for both compounds, with the difference that the decarbamoylation phase for carbaryl is more rapid than for carbamate 18.

In Silico Studies
Computational approach was performed to gain structural insights into the binding mode of pyrimidine and pyridine amide and carbamate derivatives to human AChE (hAChE) and human BChE (hBChE). For the current analysis, the recent X-ray crystallographic structures of the hAChE (PDB code: 6ZWE) and the hBChE (PDB code: 6ZWI) in complex with a triphenylphosphonium conjugate were used [37].
Even if the apo-crystal structures of EeAChE are deposited in the protein data bank (PDB), they cannot be considered to be good quality models. By superposing the EeAChE (PDB code: 1C2O) and hAChE structures (PDB code: 4EY7) with the model (PDB code: 6ZWE), the root-mean-square deviation (RMSD) values of 0.79 and 0.55 Å, respectively, were obtained (Supplementary Materials, Figure S4a), thus providing a reliable hAChE receptor for the molecular recognition analysis. Instead, regarding the choice of the BChE receptor, no crystal structure of eqBChE was available in PDB. The sequence of eqBChE derived from the Uniprot Database shares 90% sequence identity with that of the adopted hBChE, above all comparing residues of the binding site ( Figure S5). By the overlapping of the X-ray structure of the previous work (PDB code: 5NN0) and the selected model (PDB code: 6ZWI) ( Figure S4b), it was found the RMSD value equal to 0.26 Å, so justifying the selection of hBChE for the computational study. According to the criteria for the identification of reliable complexes reported in the literature [38], and the availability of the resolved crystal structures of both enzymes in the same experimental conditions, these models were considered as a good starting point for a better comparison by in silico approaches.
The molecular recognition protocol was validated by docking cocrystallized ligands into the binding site. The RMSD values between the native pose of cocrystallized ligands against hAChE and hBChE and the related best redocked conformations were found to be 0.96 Å and 0.76 Å, respectively, thus revealing the assurance of docking protocol ( Figure S6). According with previously published studies [31,32,39] no linear correlation between docking score and the experimental data is expected (Table S1). Despite the structural differences of the two ChEs cavities, the alkyl linker of six methylene units between the pyrimidine or pyridine moiety and the amide or carbamate group of the 3-19 compounds ensured the appropriate flexibility to accommodate into the three domains lining the gorge of the ChEs. Indeed, analysing the best docking poses, all compounds well recognized both the pockets of the enzymes. The average molecular surface area of all docked ligands, calculated by using the "ink-blot" method, is equal to 342.56 Å 3 and 354.13 Å 3 , into hAChE (Figure 3a,c) and hBChE (Figure 3b,d) receptors, respectively (Table S2). models were considered as a good starting point for a better comparison by in silico approaches.
The molecular recognition protocol was validated by docking cocrystallized ligands into the binding site. The RMSD values between the native pose of cocrystallized ligands against hAChE and hBChE and the related best redocked conformations were found to be 0.96 Å and 0.76 Å , respectively, thus revealing the assurance of docking protocol ( Figure  S6). According with previously published studies [31,32,39] no linear correlation between docking score and the experimental data is expected (Table S1). Despite the structural differences of the two ChEs cavities, the alkyl linker of six methylene units between the pyrimidine or pyridine moiety and the amide or carbamate group of the 3-19 compounds ensured the appropriate flexibility to accommodate into the three domains lining the gorge of the ChEs. Indeed, analysing the best docking poses, all compounds well recognized both the pockets of the enzymes. The average molecular surface area of all docked ligands, calculated by using the "ink-blot" method, is equal to 342.56 Å 3 and 354.13 Å 3 , into hAChE (Figure 3a,c) and hBChE (Figure 3b,d) receptors, respectively (Table S2).  Taking into account the obtained experimental data, we performed the molecular recognition of compounds 12, 14 and 18 towards the two ChEs receptors in order to better understand their binding modes.
Regarding hAChE, the pyridine group of the compound 12 ( Figure 4a) formed two π-π stacking interactions with His447 and Trp86, while the amine reached Ser203 through a water bridge. Instead, the phenyl ring bound Trp286 by a π-π stacking and the amide portion contributed to the accommodation at the mid-gorge pocket by means of two H- bonds with Phe295 and Arg296 backbone atoms. The growth size of the compound 14 due to the presence of the 4-phenoxybenzyl ring (Figure 4b) led to an opposite orientation with respect to the others binding poses, characterized by three π-π stacking interactions with Trp86, Tyr337 and Tyr341. Conversely, the amine group engaged a H-bond with Glu292. Concerning compound 18 (Figure 4c), the carbamate group pointed towards the backbone of Arg296, the pyridine moiety was stabilized by a π-π stacking interaction with His447 and the amine portion interacted with Ser203 by means of a water bridge.
Taking into account the obtained experimental data, we performed the molecular recognition of compounds 12, 14 and 18 towards the two ChEs receptors in order to better understand their binding modes.
Regarding hAChE, the pyridine group of the compound 12 (Figure 4a) formed two π-π stacking interactions with His447 and Trp86, while the amine reached Ser203 through a water bridge. Instead, the phenyl ring bound Trp286 by a π-π stacking and the amide portion contributed to the accommodation at the mid-gorge pocket by means of two Hbonds with Phe295 and Arg296 backbone atoms. The growth size of the compound 14 due to the presence of the 4-phenoxybenzyl ring (Figure 4b) led to an opposite orientation with respect to the others binding poses, characterized by three π-π stacking interactions with Trp86, Tyr337 and Tyr341. Conversely, the amine group engaged a H-bond with Glu292. Concerning compound 18 (Figure 4c), the carbamate group pointed towards the backbone of Arg296, the pyridine moiety was stabilized by a π-π stacking interaction with His447 and the amine portion interacted with Ser203 by means of a water bridge.  Analysing the behaviour of these three compounds towards hBChE enzyme, an increased theoretical binding affinity for compounds 12, 14 and 18 to hBChE with respect to the hAChE pocket was observed. These findings were associated with the presence of a more hydrophilic and charged region in the hBChE binding site, which probably led to the selectivity in regard to hBChE ( Figure S7). Specifically, the pyridine moiety of compound 12 (Figure 5a) created a π-π interaction with Trp82, while the amine portion engaged a hydrogen bond with His438. Furthermore, two π-π interactions were formed between the phenyl group and Trp231 and Phe329 residues, and an additional H-bond between the N atom of the amide and Pro285 backbone was observed. Analysing the behaviour of these three compounds towards hBChE enzyme, an increased theoretical binding affinity for compounds 12, 14 and 18 to hBChE with respect to the hAChE pocket was observed. These findings were associated with the presence of a more hydrophilic and charged region in the hBChE binding site, which probably led to the selectivity in regard to hBChE ( Figure S7).
Specifically, the pyridine moiety of compound 12 (Figure 5a) created a π-π interaction with Trp82, while the amine portion engaged a hydrogen bond with His438. Furthermore, two π-π interactions were formed between the phenyl group and Trp231 and Phe329 residues, and an additional H-bond between the N atom of the amide and Pro285 backbone was observed.  Compound 14 (Figure 5b) interacted with Tyr332 by means of the 4-phenoxybenzyl ring, while the pyridine group pointed to Trp231 in CAS region. Its ammine portion was able to engage a H-bond with Ser198.
As regards the compound 18, the most interesting among the carbamate derivatives, the pyridine moiety established a π-π interaction with Trp82, the amino group engaged

Molecular Dynamics Studies
The best docked poses of compounds 12, 14 and 18 into the binding pockets of both hAChE and hBChE were submitted to 300 ns of molecular dynamics (MDs) simulations. The RMSD of both the proteins' Cα atoms and the ligands were calculated in order to evaluate the structural conformational changes of the protein with respect to the initial frame and the stability of the ligand, in regard to both the receptors and their binding pockets ( Figure 6). As shown, the RMSD plot indicated that the three compounds in complex to hAChE (Figure 6a) and hBChE (Figure 6b)  Compound 14 (Figure 5b) interacted with Tyr332 by means of the 4-phenoxybenzyl ring, while the pyridine group pointed to Trp231 in CAS region. Its ammine portion was able to engage a H-bond with Ser198.
As regards the compound 18, the most interesting among the carbamate derivatives, the pyridine moiety established a π-π interaction with Trp82, the amino group engaged a H-bond with His438, while the phenylcarbamate portion formed a closer network of interactions with Trp231 and Pro285 (Figure 5c).

Molecular Dynamics Studies
The best docked poses of compounds 12, 14 and 18 into the binding pockets of both hAChE and hBChE were submitted to 300 ns of molecular dynamics (MDs) simulations. The RMSD of both the proteins' Cα atoms and the ligands were calculated in order to evaluate the structural conformational changes of the protein with respect to the initial frame and the stability of the ligand, in regard to both the receptors and their binding pockets ( Figure 6). As shown, the RMSD plot indicated that the three compounds in complex to hAChE (Figure 6a) and hBChE (Figure 6b)   In order to obtain further insights about the selectivity of the compounds 12, 14 and 18 towards BChE with respect to AChE, the single contributions of hydrophobic, water In order to obtain further insights about the selectivity of the compounds 12, 14 and 18 towards BChE with respect to AChE, the single contributions of hydrophobic, water bridge and hydrogen-bonding interactions, pivotal interactions into ChEs active sites, were considered.
In detail, regarding the hAChE enzyme, it was observed that the compound 12 benzamido group interacted mainly with Phe295 for 32% of the MDs time by means a water bridge, while engaged hydrophobic contacts with Trp86, Trp286 and Tyr441 residues. The presence of a 4-phenoxybenzyl ring led to compound 14 being well accommodated into CAS and mid-gorge pockets thus ensuring several interactions in the enzyme. Indeed, during the MDs, the amide group provided a H-bond with Tyr124 (79%) at the CAS site and a water bridge with Phe295 (60%), located at the mid-gorge region. The amino portion was directly related to Ser293 (30%), while the 4-phenoxybenzyl ring engaged a π-π stacking interaction with Tyr337 for half of the MDs. Finally, the amino group and the phenylcarbamate of compound 18 participated in major interactions with His447 (36%) and Phe295 (75%), respectively.
Analysing the dynamic behaviour of the investigated compounds towards the hBChE enzyme, during 300 ns of MDs, it was observed that in compound 12, the benzamido group was anchored to Ser287 (43%), Trp231 (37%) and Phe329 (81%) in the CAS and mid-gorge pockets. However, the amino-pyridine portion engaged in a hydrophobic interaction with Phe73 (25%) and an H-bond with Thr284 (37%), respectively, in PAS and mid-gorge sites.
The steric hindrance brought about by the 4-phenoxybenzyl ring in compound 14 allowed it to orient the pyridine moiety versus the catalytic region, which established a water bridge with His438 (53%) and a π-π stacking interaction with Trp231 (79%), whilst the amino group was able to form a hydrogen bond with Pro285 (67%) in the CAS pocket. Lastly, the phenylcarbamate moiety of compound 18 engaged hydrophobic interactions with Trp231 (53%) and a H-bond with Ser287 (75%). Presumably, this pivotal interaction in the catalytic site could represent the main reason for the selectivity of compound 18 towards BChE, carrying it effortlessly to the carbamoylation of the serine.
Computational analysis revealed that the multiple interactions of compound 12, 14 and 18 with specific gorge hot spots may trap the enzyme in a peculiar conformation, decreasing the ability of the enzyme-inhibitor complex to access protein fluctuations, thus improving the compounds inhibitory potency. Taking into account that the conformational heterogeneity in ChEs is related to specific classical or non-classical functions, a promising inhibitor should be able to freeze the "classical" enzyme structure in a particular conformational state, thus preventing acetylcholine hydrolysis and non-classical protein functions. Our results may confirm that the pyridine moiety placed into the alkyl chain of six methylene units represents a crucial pharmacophore feature for ChEs inhibition, as detected in our previous work [32]. Still, the stronger inhibition capability of compound 18 may suggest that the presence of phenylcarbamate insertion is more suitable for BChE enzyme interaction, ensuring the selectivity toward its binding pocket.

Chelation Studies
For compounds 12, 14 and 18, the most potent as inhibitors of eqBChE, the chelating ability on the biometals Fe 3+ , Cu 2+ and Zn 2+ was evaluated, in order to define their possible multifunctional profile.
Initially, the UV-Vis spectrum of the ligand was recorded and compared with the spectra obtained by adding an excess of metal to the ligand solution, maintaining the same concentration of the ligand (ligand/metal ratio 1:3, 1:5 or 1:10). The variation of the UV-Vis spectrum of the ligand in presence of metal ions is indicative of the formation of the complex. Based on this, it has been observed that all the tested pyridine derivatives have the ability to chelate Fe 3+ and Cu 2+ ions but not Zn 2+ . As previously supposed, the chelating activity of these compounds is due to the presence of two adjacent nitrogen atoms in the 2-amino-pyridine moiety.
UV-Vis titrations of these compounds were carried out with the metal ions with which they form complexes, recording first the UV-Vis spectrum of the ligand and then the spectra obtained by mixing solutions of ligand and metal according to increasing metal/ligand molar ratios (Figure 7a). mole fractions of the ligand are in abscissa. The mole fraction X, which causes the maximum variation of absorbance, is extracted from the graph and then the value of the coefficient n, which corresponds to the number of ligand molecules per cation, is obtained (Figure 7b). The UV-Vis titrations spectra and the Job's plots of all tested compound are reported in the Supplementary Materials (Figures S8-S17), while the data extrapolated from Job's plots are summarized in Table 5.   (12,14,18) in presence of Fe 3+ and Cu 2+ . The table shows the wavelength values (λ) in which the absorbance variation was measured, the molar ratios (X) that cause the maximum variation of absorbance and the coordination values (n) which correspond to the number of ligand molecules per cation. The stoichiometries of metal-ligand complexes were determined through Job's method [40,41].
In the Job's plot, the values of ∆A, measured at the selected wavelengths, where evident absorbance variations were observed in the titration spectra, are in ordinate and the mole fractions of the ligand are in abscissa. The mole fraction X, which causes the maximum variation of absorbance, is extracted from the graph and then the value of the coefficient n, which corresponds to the number of ligand molecules per cation, is obtained (Figure 7b). The UV-Vis titrations spectra and the Job's plots of all tested compound are reported in the Supplementary Materials (Figures S8-S17), while the data extrapolated from Job's plots are summarized in Table 5. Table 5. Job's plot data for the tested compounds (12,14,18) in presence of Fe 3+ and Cu 2+ . The table shows the wavelength values (λ) in which the absorbance variation was measured, the molar ratios (X) that cause the maximum variation of absorbance and the coordination values (n) which correspond to the number of ligand molecules per cation.

Inhibition of Amyloid and Tau Aggregation
For compounds 12, 14 and 18, the inhibitory activity against Aβ 42 and tau aggregation was evaluated. The anti-aggregating effect of tested compounds was monitored by a cell-based assay in intact Escherichia coli cells that overexpress either Aβ 42 peptide or tau protein, which, upon expression, form insoluble inclusion bodies that were stained with thioflavin-S [42]. The percentages of inhibition towards Aβ 42 and tau aggregation were determined at an inhibitor concentration of 100 µM and the obtained results are reported in Table 6. In general, these compounds show a weak anti-aggregating activity; among these, the pyridine carbamate 18 is the best inhibitor of aggregation of both Aβ 42 and tau. More assays were performed to obtain a SEM < 5% with a maximum of 10 independent assays. b The known antiamyloid drug, DP-128 (100 µM), was used as a reference compound.

Computation of Physicochemical Descriptors and ADME Parameters
Physicochemical descriptors and ADME parameters of compounds 12, 14 and 18, the most interesting by virtue of their inhibitory activity towards BChE, as well as their chelating ability, were predicted by means of the SwissADME public server [43] and the obtained data are reported in Table 7. All studied compounds fit the Lipinski's rule of five (MW ≤ 500; MLogP ≤ 4.15; H bond acceptor ≤ 10; H bond donor ≤ 5) [44]. Table 7. Physicochemical properties and ADME parameters: Molecular Weight (MW); number H-bond acceptors (HBA); number H-bond donors (HBD); number of heavy atoms (heavy atoms); number rotatable bonds (RB); topological polar surface area in Å 2 (TPSA); octanol/water partition coefficient (MLogP); water solubility (LogS ESOL); water solubility class (Sol class); gastrointestinal absorption (GI); blood-brain barrier permeation (BBB); number of Lipinski's rule of five violations (Lipinski viol They should be soluble or moderately soluble in water. Moreover, they should have high gastrointestinal absorption after oral administration, and they should be BBB accessible.

Cell Viability Studies
Compounds 12, 14 and 18 were tested to evaluate the cytotoxic effects on U-87 MG Cell Line from human brain (glioblastoma astrocytoma) at concentration ranging from 1 to 50 µM. The obtained data, represented in the histogram of Figure 8, suggest that the tested compounds are characterized by low toxicity towards the studied cells; in particular, compounds 12 and 18 showed a similar dose related toxicity profile, providing a cell viability over 90% up to a concentration equal to 10 µM and above 35% at 50 µM; instead, compound 14 showed a cell viability above 65% at 10 µM and 30% at 50 µM. to 50 µ M. The obtained data, represented in the histogram of Figure 8, suggest that the tested compounds are characterized by low toxicity towards the studied cells; in particular, compounds 12 and 18 showed a similar dose related toxicity profile, providing a cell viability over 90% up to a concentration equal to 10 µ M and above 35% at 50 µ M; instead, compound 14 showed a cell viability above 65% at 10 µ M and 30% at 50 µ M.

Figure 8.
Cell viability assay of U-87 MG cells in the presence of an increasing concentration of the studied compounds was evaluated using the MTT method. Data represent the absorbance relative to untreated cells in the same experiment and standardized to 100%. All data points were performed in triplicate and at least in three independent experiments.

Chemistry
All reagents and solvents were of analytical grade and were purchased from Sigma-Aldrich (Milano, Italy) or from Fluorochem (Hadfield, UK). Triethylamine was freshly purified by distillation over potassium hydroxide. Column chromatographies were performed on silica gel (Merck; 63−200 μm particle size). 1 H-NMR and 13 C-NMR spectra were acquired at 25 °C on a Bruker AVANCE-400 spectrometer at 9.4 Tesla operating at 400 MHz ( 1 H-NMR) and 100 MHz ( 13 C-NMR) or on a Bruker AVANCE-200 spectrometer at 4.7 Tesla operating at 50 MHz ( 13 C-NMR) (Bruker, Billerica, MA, USA); chemical shift values (δ) are given in ppm, relative to TMS, using the solvent as the internal reference; coupling constants are given in Hz. The following abbreviation were used: s = singlet, d = doublet, t = triplet, q = quartet, dd = double doublet, ddd = double double doublet, dt = double triplet, bs = broad singlet, bt = broad triplet, m = multiplet. Mass spectra were recorded on a ThermoFinnigan LCQ Classic (San Jose, CA, USA) LC/MS/MS ion trap equipped with an ESI source and a syringe pump; samples (10 −4 -10 −5 M in MeOH/H2O 80:20) were infused in the electrospray system at a flow rate of 5-10 µL·min −1 ; when necessary, 50 µ L of 10 −2 M HCOOH or 10 −2 M NH3 were added to the sample solutions, in order to promote the analyte ionisation; the ESI-MS data are given as m/z, with mass expressed in amu. Melting points were determined on a FALC Mod. 360 D apparatus and are uncorrected. Infrared spectra were recorded on a Perkin Elmer Spectrum One FT-IR spectrometer equipped with an ATR system. The purity of the compounds was determined by elemental analyses obtained by a PE 2400 (Perkin-Elmer, Waltham, MA, USA) analyzer and the analytical results were within ±0.4% of the theoretical values for all compounds.

Chemistry
All reagents and solvents were of analytical grade and were purchased from Sigma-Aldrich (Milano, Italy) or from Fluorochem (Hadfield, UK). Triethylamine was freshly purified by distillation over potassium hydroxide. Column chromatographies were performed on silica gel (Merck; 63−200 µm particle size). 1 H-NMR and 13 C-NMR spectra were acquired at 25 • C on a Bruker AVANCE-400 spectrometer at 9.4 Tesla operating at 400 MHz ( 1 H-NMR) and 100 MHz ( 13 C-NMR) or on a Bruker AVANCE-200 spectrometer at 4.7 Tesla operating at 50 MHz ( 13 C-NMR) (Bruker, Billerica, MA, USA); chemical shift values (δ) are given in ppm, relative to TMS, using the solvent as the internal reference; coupling constants are given in Hz. The following abbreviation were used: s = singlet, d = doublet, t = triplet, q = quartet, dd = double doublet, ddd = double double doublet, dt = double triplet, bs = broad singlet, bt = broad triplet, m = multiplet. Mass spectra were recorded on a ThermoFinnigan LCQ Classic (San Jose, CA, USA) LC/MS/MS ion trap equipped with an ESI source and a syringe pump; samples (10 −4 -10 −5 M in MeOH/H 2 O 80:20) were infused in the electrospray system at a flow rate of 5-10 µL·min −1 ; when necessary, 50 µL of 10 −2 M HCOOH or 10 −2 M NH 3 were added to the sample solutions, in order to promote the analyte ionisation; the ESI-MS data are given as m/z, with mass expressed in amu. Melting points were determined on a FALC Mod. 360 D apparatus and are uncorrected. Infrared spectra were recorded on a Perkin Elmer Spectrum One FT-IR spectrometer equipped with an ATR system. The purity of the compounds was determined by elemental analyses obtained by a PE 2400 (Perkin-Elmer, Waltham, MA, USA) analyzer and the analytical results were within ±0.4% of the theoretical values for all compounds.

General Procedure for the Synthesis of Pyrimidine and Pyridine Amide Derivatives 3-15
To a solution of 1,1-carbonyldiimidazole (1 eq) in 15 mL of CH 3 CN the opportune carboxylic acid (1 eq) was added and the reaction mixture was stirred under reflux conditions for 5 h. Then intermediate 1 or 2 (1 eq), obtained as previously reported [31], was dissolved in 10 mL of warm CH 3 CN and added to the reaction mixture, which was then refluxed for 18 h. After this time, the solvent was removed under reduced pressure. The obtained residue was purified by column chromatography on silica gel and/or by crystallisation.

Percent Inhibition of EeAChE and eqBChE
Each experiment was repeated at least in triplicate. The potency of each compound to inhibit EeAChE or eqBChE activity was expressed as percent inhibition calculated using the following equation: where A i and A c represent the average of absorbance variation in presence of inhibitor and without the inhibitor, respectively.

Time Dependent Inhibition Assay for Carbamate Derivatives 16-19
For carbamate derivatives 16-19, the percentages of inhibition were measured both at zero time, according to the procedure described above, and after incubation of 1 h at 25 • C. In this case, for each compound, three cuvettes were prepared with 3.0 mL of a solution in 0.1 M phosphate buffer (pH = 7.4) containing DTNB (0.25 mM) and EeAChE (0.083 UmL −1 ) or eqBChE (0.083 UmL −1 ) and with 1 µL of inhibitor solution in DMSO, in order to obtain final compound concentrations equal to 9 µM for EeAChE and equal to 900 nM for eqBChE; in parallel, three control cuvettes were prepared with 3.0 mL of a solution in 0.1 M phosphate buffer containing DTNB (0.25 mM) and enzyme (0.0833 UmL −1 ), without the inhibitor, to determine the 100% of enzyme activity. The cuvettes were prepared consecutively at intervals of 3 min to allow equal incubation times.
Each cuvette was incubated for 1 h at 25 • C. Then the blank was made, and the reaction was activated by adding 30 µL of a solution in 0.1 M phosphate buffer (pH = 7.4) of acetylthiocholine (10 mM), in order to obtain a final concentration of substrate equal to 100 µM. The absorbance variation was measured at 412 nm at 25 • C, between 0.5 and 1.5 min, and the percentage of inhibition after incubation was calculated.
For carbamate 18, the percentages of inhibition toward eqBChE were determined at increasing incubation times. The same protocol described above was followed, but with different incubation times, every 15 min for 2 h. For each incubation time three cuvettes were prepared as previously described, with a final compound concentration equal to 700 nM. In order to determine the inhibition constants (K i ) and inhibition mechanism the rates of hydrolysis at three different concentrations of substrate in the presence of five different concentration of inhibitor were measured. The recorded data were analyzed by the enzyme kinetic module of SigmaPlot (version 8.02), plotting the reciprocal of rate of hydrolysis (1/v, min/µM) vs. the concentration of inhibitor (nM), according to Dixon's method [34] in order to find the best fitting model of inhibition, as indicated by calculated linear regression coefficient R 2 .

Determination of IC 50 vs. eqBChE for Compound 18
For compound 18 the IC 50 value vs. eqBChE was calculated. Eight different inhibitor solutions in DMSO were prepared, having suitable concentrations, so as to add 1 µL of these solutions to the assay mixture to obtain a final concentration between 50 and 1500 nM, with a final DMSO content <0.033%. For each inhibitor concentration the same protocol described above for the determination of the percentage of inhibition towards eqBChE (0.0833 UmL −1 ) after incubation for 1 h at 25 • C and then activating the reaction with acetylthiocholine (100 µM) was followed. The recorded data were analyzed by the enzyme kinetic module of SigmaPlot, plotting the percentages of inhibition as a function of the different concentrations of inhibitor, expressed in logarithmic scale, in order to obtain the sigmoid graph, from which it was possible to determinate the IC 50 value of the tested compound. The IC 50 value was confirmed by repeating the experiment twice. • three were prepared with 2.5 mL of eqBChE solution (0.50 UmL −1 ) and 5 µL of inhibitor solution; • three were prepared with 2.5 mL of eqBChE solution (0.50 UmL −1 ) and 5 µL of DMSO.
The Vivaspin ® thus prepared were incubated for 1 h at 25 • C. After this time the percentage of inhibition was determined following the Ellman method. The blank was recorded by taking 500 µL from the Vivaspin ® and diluting them in the cuvette in 2.5 mL of DTNB solution (0.3 mM). The reaction was activated by adding 30 µL of acetylthiocholine solution (10 mM), so as to obtain a final concentration of 100 µM in the cuvette. The absorbance variation was measured between 0.5 and 1.5 min at 412 nm, repeating the measurement for each of the six Vivaspin ® .
After this procedure the first centrifugation of the Vivaspin ® was carried out for 60 min at 4000 rpm and subsequently the solutions that crossed the membrane and settled in the filtrate tube were removed. In this way the solutions present in the concentrator body were concentrated up to 100 µL. To these solutions 1 mL of 0.1 M phosphate buffer (pH = 7.4) was added and the second centrifugation was carried out for 30 min at 4000 rpm. The solutions present in the concentrator body were concentrated up to 200 µL. The second washing was carried out, adding 1 mL of buffer and centrifuging for 30 min at 4000 rpm. The solutions present in the concentrator body were concentrated up to 250 µL. A third wash was carried out adding 1 mL of buffer and centrifuging for 50 min at 4000 rpm. The solutions present in the concentrator body were concentrated up to 150 µL. At this point, to recover the concentrated solution, the filtrate tube was removed, the concentrator body was turned upside down and inserted inside the concentrate recovery cap. The Vivaspin ® was centrifuged for 10 min at 4000 rpm. After this time the recovered solution was taken up with 3 mL of DTNB solution (0.3 mM) and transferred into the cuvette: blank was recorded. The reaction was activated when 30 µL of substrate solution (10 mM) were added so as to have a final concentration of 100 µM in the cuvette. The absorbance variation between 0.5 and 1.5 min was measured at 412 nm, repeating the measurement for each of the recovered solutions. With these data the percentage of inhibition was calculated using the following equation: where A ri and A rc represent the average of absorbance variation in the solutions recovered from the Vivaspin ® , initially prepared with the inhibitor, and in the solutions recovered from the Vivaspin ® , initially prepared without the inhibitor, respectively. As a positive control carbaryl was used, for which the same procedure was followed, using a 27 mM carbaryl stock solution in DMSO.

Molecular Docking Studies and ADME Prediction
All in silico simulations were carried out by means of Schrödinger Suite 2018-1 [45]. The newest X-ray crystallographic structures of both hAChE and hBChE in complex with a triphenylphosphonium conjugate with the 6ZWE and 6ZWI PDB codes, respectively, were adopted for the computational studies [37]. hChEs structures were prepared by using the Maestro Protein Preparation Wizard [46] tool and were refined to optimize hydrogenbonds and energy minimized by using OPLS_2005 as force field at pH 7.4 [47,48].
The homology between the human and non-human ChEs isoforms used for the enzymatic assays was assessed by the sequence alignment performed by Prime [49].
All synthesized ligands were submitted to 5000 steps of Monte Carlo conformational search applied to all rotatable bonds. The water solvent effect was considered using the implicit model GB/SA [50]. The global minimum energy structures were used to carry out molecular docking simulations. Then, target binding sites were defined by means of a regular grid box of about 27,000 Å 3 centred on the catalytic serine residues. All docking simulations were computed using Glide [51] ligand flexible algorithm at standard-precision (SP) level [52]. The best docked poses of the most interesting inhibitors resulted through the enzymatic assays of pyridine amides and carbamate derivatives were submitted to 300 ns of molecular dynamics (MDs) simulations by using Desmond ver. 4.2 [53]. The system was solvated in TIP3P explicit solvent model and counter ions were added to neutralize the system net charge. After the optimisation of the solvated model, we relaxed the system with Martyna-Tobias_Klein isobaric-isothermal ensemble (MTK_NPT). This preliminary stage included two energy minimisations of 2000 steps: in the first run, the system was restrained with a force constant of 50 kcal·mol −1 ·A −1 , while in the second one all the system was released without any restrains. The following conditions for MDs were used: NPT ensemble, a temperature of 300 K, a pressure of 1 bar with the Berendsen thermostat-barostat, a recording interval equal to 250 ps both for energy and for trajectory collecting 1000 frames for each simulation [54]. ADME descriptors of the most active compounds were predicted using the Swis-sADME public server [43]. The following ADME parameters were selected: Molecular weight (MW); number H-bond acceptors (HBA); number H-bond donors (HBD); number of heavy atoms (heavy atoms); number rotatable bonds (RB); topological polar surface area in Å 2 (TPSA); Octanol/water partition coefficient (MLogP); water solubility (LogS ESOL); water solubility class (Sol class); gastrointestinal absorption (GI); number of Lipinski's rule of five violations (Lipinski viol) [44]. For the first measure, 150 µL of ligand solution were diluted in 2.85 mL of MeOH into sample cuvette, while only MeOH was placed in reference cuvette and the UV-Vis spectrum was recorded. For subsequent measurements, 2.85 mL of the opportune metal solution were placed both in a sample cuvette and in a reference cuvette; then 150 µL of ligand solution were added to a sample cuvette, while 150 µL of MeOH were added to the reference cuvette and the UV-Vis spectra were recorded. The spectra thus obtained were superimposed to observe the variations present between the ligand alone and in presence of increasing amount of metal.

Job's Plot
To determine the stoichiometric coefficients of the complexes, Job's method [41] was used, which requires the mixing, in appropriate proportions, of the equimolar solution of metal ion and ligand, so that the final volume and the total moles present in the cuvette are equal for each measurement. The absorbance values were recorded at the wavelengths, extrapolated from the UV-Vis titration spectra, where the maximum variation of absorbance was observed. At the same time two solutions in MeOH of equal concentration were prepared, one of ligand and one of metal ion, For each determination 9-20 measurements were made, introducing the appropriate volumes of metal and ligand solutions in the sample cuvette to obtain mole fractions of the ligand in the range 0.1-0.9, while appropriate volumes of ligand solution and MeOH were placed in the reference cuvette, to obtain the same concentration of ligand in the sample cuvette. The absorbance values of the metal ion, calculated by the Lambert-Beer relation, known the extinction coefficients ε of the metal ion at used wavelengths and the nominal concentrations of the metal ion for each measure, were algebraically subtracted from the recorded absorbance values, in order to obtain the exclusive absorbance variations due to the complex formation.
The resulting ∆A were reported in graph as a function of the mole fraction of the ligand and through the intersections of the linear regression lines, the mole fraction X, which caused the maximum variation in absorbance, was determined and used to calculate the value of the coefficient n, which corresponds to the number of ligand molecules per cation, applying the following equation: n = X 1 − X 3.5. Inhibition of Amyloid and Tau Aggregation 3.5.1. Cloning and Overexpression of the Aβ 42 Peptide E. coli competent cells BL21 (DE3) were transformed with the pET28a vector (Novagen, Inc., Madison, WI, USA) carrying the DNA sequence of Aβ 42 . Because of the addition of the initiation codon ATG in front of both genes, the overexpressed peptide contains an additional methionine residue at its N terminus. For overnight culture preparation, an amount of 10 mL of the M9 minimal medium containing 50 µg·mL −1 of kanamycin was inoculated with a colony of BL21 (DE3) bearing the plasmid to be expressed at 37 • C. For expression of the Aβ 42 peptide, the required volume of overnight culture to obtain 1:500 dilution was added into the fresh M9 minimal medium containing 50 µg·mL −1 of kanamycin and 250 µM Th-S. The bacterial culture was grown at 37 • C and 250 rpm.
When the cell density reached OD 600 = 0.6, an amount of 980 µL of culture was transferred into Eppendorf tubes of 1.5 mL with 10 µL of each compound to be tested in DMSO and 10 µL of isopropyl 1-thio-β-D-galactopyranoside (IPTG) at 100 mM. The final concentration of the drug was fixed at 100 µM. The samples were grown overnight at 37 • C and 1400 rpm using a ThermoMixer (Eppendorf, Hamburg, Germany). As a negative control (maximal amyloid presence), the same amount of DMSO without the drug was added in the sample. In parallel, non-induced samples (in the absence of IPTG) were also prepared and used as positive controls (non-amyloid presence). In addition, these samples were used to assess the potential intrinsic toxicity of the compounds and to confirm the correct bacterial growth.

Cloning and Overexpression of the Tau Protein
E. coli BL21 (DE3) competent cells were transformed with pTARA containing the RNApolymerase gene of T7 phage (T7RP), under the control of the promoter PBAD. E. coli BL21 (DE3), with pTARA competent cells, were transformed with the pRKT42 vector encoding four repeats of the tau protein in two inserts. For overnight culture preparation, 10 mL of the M9 medium containing 0.5% of glucose, 50 µg·mL −1 of ampicillin, and 12.5 µg·mL −1 of chloramphenicol was inoculated with a colony of BL21 (DE3) bearing the plasmids to be expressed at 37 • C. For the expression of the tau protein, the required volume of overnight culture to obtain 1:500 dilution was added to the fresh M9 minimal medium containing 0.5% of glucose, 50 µg·mL −1 of ampicillin, 12.5 µg·mL −1 of chloramphenicol, and 250 µM Th-S. The bacterial culture was grown at 37 • C and 250 rpm. When the cell density reached OD 600 = 0.6, an amount of 980 µL of culture was transferred into Eppendorf tubes of 1.5 mL with 10 µL of each compound to be tested in DMSO and 10 µL of arabinose at 25%. The final concentration of the drug was fixed at 100 µM. The samples were grown overnight at 37 • C and 1400 rpm using a ThermoMixer (Eppendorf, Hamburg, Germany). As a negative control (maximal presence of tau), the same amount of DMSO without the drug was added in the sample. In parallel, non-induced samples (in the absence of arabinose) were also prepared and used as positive controls (absence of tau). In addition, these samples were used to assess the potential intrinsic toxicity of the compounds and to confirm the correct bacterial growth.

Thioflavin S (Th-S) Steady-State Fluorescence
Th-S (T1892) and other chemical reagents were purchased from Sigma (St. Louis, MO, USA). The Th-S stock solution (2500 mM) was prepared in double-distilled water purified through a Milli-Q system (Millipore, Burlington, MA, USA). Th-S fluorescence and absorbance were tracked using a DTX 800 plate reader Multimode Detector equipped with a Multimode Analysis Software (Beckman-Coulter, Indianapolis, IN, USA). Filters of 430/35 and 485/20 nm were used for the excitation and emission wavelengths, respectively. A filter of 535/25 nm was also used for the absorbance determination. To normalize the Th-S fluorescence as a function of the bacterial concentration, OD 600 was obtained using a Shimadzu UV-2401 PC UV-Vis spectrophotometer (Shimadzu, Japan). Note that the fluorescence normalisation was carried out considering 100% of the Th-S fluorescence of the bacterial cells expressing the peptide or protein in the absence of the drug and 0% of the Th-S fluorescence of the bacterial cells not expressing the peptide or protein. A minimum of five independent assays (with three replicates for assay) was performed for each tested compound. More assays were performed to obtain a SEM <5% with a maximum of 10 independent assays.

Cytotoxicity Assay
The compounds 12, 14 and 18 were evaluated for cell viability effects using the MTT assay. Briefly, U-87 MG Cell Line from human brain (glioblastoma astrocytoma) (ATCC, Manassas, VA, USA) were seeded into 96-well microtiter plates (NunclonTM, Nunc, Germany) at a density of 1.5 × 10 4 cells/well. After 24 h cells were exposed to increased concentrations of tested compounds (1, 5, 10 and 50 µM) in cell culture medium. After an incubation time of 24 h, the medium was replaced by a new, fresh medium containing 0.5 mg/mL MTT (Sigma, Deisenhofen, Germany). After 2 h at 37 • C in 5%, CO 2 unreacted dye was removed and the cells were dissolved by DMSO (Merck, Darmstadt, Germany). Absorbance was read at 570 nm by a microtiter spectrophotometer plate reader. The data were expressed as absorbance relative to untreated cells in the same experiment and standardized to 100%. All data points were performed in triplicate and at least three independent experiments.

Conclusions
Based on the interesting results previously reported for some pyrimidine and pyridine diamine derivatives [32], new compounds, containing a 2-amino-pyrimidine or a 2-amino-pyridine moiety connected to an aromatic group by a flexible alkyl chain and an amide or carbamic group, were designed and tested both in vitro and in silico to evaluate their ChEs inhibitory ability and to propose possible binding mode with these enzymes. Regarding the synthesized compounds, the determination of the percentages of inhibition toward EeAChE and eqBChE were carried out at concentrations equal to 9−0.9 µM. In general, among amide derivatives, the pyrimidines are weak inhibitor of both ChEs, while the pyridines appear more potent inhibitors than the pyrimidine analogues, in particular toward eqBChE. Additionally, among carbamates derivatives, pyridine carbamates show greater inhibitory potency on both ChEs compared to the corresponding pyrimidine derivatives. Moreover, the carbamate derivatives are more selective towards eqBChE than EeAChE. The increase in the percentages of inhibition towards both enzymes over time, in particular for the phenylcarbamates 16 and 18, suggests that these compounds could act as carbamoylating agents. If compared with the corresponding amine derivatives already studied, the replacement of the amino group with the amide (pyridine 12) or the carbamate (pyrimidine 16 and pyridine 18) reduce the inhibitory activity on EeAChE, but enhance it toward eqBChE, leading to selective eqBChE inhibitors. Among the studied compounds, the most potent eqBChE inhibitors are the pyridine amides 12 (K i = 2.988 ± 0.190 µM) and 14 (K i = 0.621 ± 0.043 µM), with an unsubstituted phenyl ring or a 4-phenoxybenzyl ring on one side of the aliphatic chain, respectively, and the pyridine phenylcarbamate 18 (IC 50 = 0.454 ± 0.082 µM).
In silico analyses allowed to identify and rationalize the crucial pharmacophore features for BChE inhibition, which are represented by the pyridine moiety and the alkyl chain of six methylene units. Additionally, the presence of phenylcarbamate in compound 18 resulted in enhancing the binding affinity, owing to the long-lasting interaction with Ser287, observed by MDs run. This interaction seems to maximize the selectivity toward BChE enzyme.
As already evidenced for the pyrimidine and pyridine amine derivatives previously studied, metal chelation analyses confirmed that the tested compounds 12, 14 and 18 have the ability to chelate Fe 3+ and Cu 2+ ions, due to the presence of two adjacent nitrogen atoms in the 2-amino-pyridine moiety.
In cellulo studies on E. coli cells showed that the tested compounds have a weak anti-aggregating activity toward both Aβ 42 and tau proteins, while cytotoxicity assays revealed a low toxicity on human brain cells.
In conclusion, among the molecules studied in this work, compounds 12, 14 and 18 are selective inhibitors of eqBChE and retain the chelating ability, the low toxicity and the good predicted ADME parameters already shown for amine analogues.