Design, Synthesis and Biological Evaluation of Biscarbamates as Potential Selective Butyrylcholinesterase Inhibitors for the Treatment of Alzheimer’s Disease

As butyrylcholinesterase (BChE) plays a role in the progression of symptoms and pathophysiology of Alzheimer’s disease (AD), selective inhibition of BChE over acetylcholinesterase (AChE) can represent a promising pathway in treating AD. The carbamate group was chosen as a pharmacophore because the carbamates currently or previously in use for the treatment of AD displayed significant positive effects on cognitive symptoms. Eighteen biscarbamates with different substituents at the carbamoyl and hydroxyaminoethyl chain were synthesized, and their inhibitory potential toward both cholinesterases and inhibition selectivity were determined. The ability of carbamates to cross the blood–brain barrier (BBB) by passive transport, their cytotoxic profile and their ability to chelate biometals were also evaluated. All biscarbamates displayed a time-dependent inhibition with inhibition rate constants within 10−3–10−6 M−1 min−1 range for both cholinesterases, with generally higher preference to BChE. For two biscarbamates, it was determined that they should be able to pass the BBB by passive transport, while for five biscarbamates, this ability was slightly limited. Fourteen biscarbamates did not exhibit a cytotoxic effect toward liver, kidney and neuronal cells. In conclusion, considering their high BChE selectivity, non-toxicity, ability to chelate biometals and pass the BBB, compounds 2 and 16 were pointed out as the most promising compounds for the treatment of middle and late stages of AD.


Introduction
Alzheimer's disease (AD) is the most common and widespread neurodegenerative disease characterized by memory and judgment ability loss and personality changes [1]. Today, it affects more than 50 million people worldwide and has a tendency to continuously grow as a result of the aging of the world's population [2]. AD is a multifactorial disease with various pathological features: deficiency of the neurotransmitter acetylcholine (ACh), accumulation of beta-amyloid (Aβ) plaques, hyperphosphorylation of tau protein, overstimulation of N-methyl-D-aspartate receptor, changes in the homeostasis of biometals, increased MAO-B enzyme activity and oxidative stress [3,4].
Although the multifactorial nature of disease points at the existence of a number of possible targets, the existing treatment for AD is based mainly on increasing the concentration of ACh by inhibiting the action of enzymes responsible for its hydrolysis, acetylcholinesterase (AChE) and butyrylcholinesterase (BChE), or using N-methyl-D-aspartate (NMDA) receptor antagonists; both directions are directed at the restoration of cognitive functions of patients and alleviating the symptoms of the disease. The new entry in AD treatment is the use of a monoclonal IgG1 antibody [5], aducanumab, which is the only drug aimed to cure or stop the progression of the disease. It exerts its mechanism of action by selectively targeting and binding aggregated soluble oligomers and insoluble fibrils conformations of Aβ plaques to reduce their levels in the brain [5]. Three of the five currently approved drugs for AD treatment are cholinesterases inhibitors. Due to AChE's crucial role in neurotransmission and consequently the development and progression of AD, drugs aimed at increasing the levels of ACh in the brain were developed primarily as AChE inhibitors. The currently approved [6] cholinesterase inhibitors donepezil and galantamine are selective AChE inhibitors, while carbamate rivastigmine is a non-selective inhibitor of AChE and BChE [7,8].
An important role of the BChE in the regulation of brain ACh levels in late AD [9,10] was pointed out based on findings that, during disease progression, AChE activity decreases to about 33-45% of its normal activity, while BChE activity progressively increases by about 40-90% of its normal activity [8]. Animal studies on rodents have shown that selective inhibition of BChE, with respect to AChE, has a beneficial effect on cognitive abilities of rodents with AD and reduces accumulations of amyloid plaques in their brains [11,12]. Moreover, recent studies have shown that selective inhibition of BChE reduces the occurrence of side effects like tremor, anxiety, hypersalivation and sweating seen with the AChE or nonselective cholinesterase inhibitors currently in use [13]. Consequently, selective inhibition of BChE has evolved into a promising new approach in the treatment of middle and advanced AD.
To date, a number of selective BChE inhibitors with different structural elements have been synthesized and tested. Many of them had tacrine as a structural scaffold (thienothiazines, thiazoles, quinuclidines, benzofuranes, quinolines, etc.), but many other structural scaffolds were also used, as well as hybrids containing galantamine, donepezil and rivastigmine as moieties [14][15][16][17][18][19][20]. However, so far, no drug that is a selective BChE inhibitor has been approved for the treatment of AD.
Very promising results were obtained for compounds with a carbamate group as a functional scaffold due to the similarity of mechanism of cholinesterases hydrolysis of carbamates with the mechanism of AChE hydrolysis of its physiological substrate ACh [21,22]; the difference lies in the stability of the carbamylated cholinesterases which decarbamylate more slowly than acetylated AChE [22]. Several drugs which are currently in use or had been used for the treatment of AD were carbamates ( Figure 1) [23] of which rivastigmine is a non-selective cholinesterase inhibitor currently in use [24,25]. Physostigmine was the first carbamate clinically used as a cholinesterase inhibitor, but its use today is reduced due to its poor bioavailability and adverse side effects [8]. Many physostigmine derivatives entered clinical trials: development of eptastigmine as a drug was stopped due the development of neutropenia in patients [26]; phenserine tartrate entered phase III of clinical trials, but its further development failed due to certain concerns regarding the documentation of its clinical trials [27][28][29]; cymserine inhibits BChE 15 times faster than AChE, but its development has been discontinued due to unacceptable side effects caused by its toxic metabolite eseroline [30][31][32]. After cymserine, several cymserine derivatives have been synthesized with greater selectivity for BChE than cymserine, and several of them were tested in animals demonstrating the ability to increase brain ACh levels and produce nootropic effects, as well as reducing levels of amyloid precursor protein and Aβ [30], but so far only bisnorcymserine has entered phase I clinical trials for the treatment of AD [33].
Bambuterol, a biscarbamate ester of terbutaline, a drug for the treatment of asthma, has been shown to be highly selective to BChE by inhibiting it 20,000 times faster than AChE [22,34]. Considering the fact that it is already in use as a drug without side effects related to the cholinergic system, bambuterol stood out as a promising candidate for repurposing into an agent for AD treatment. This was additionally supported by the fact that the monocarbamate derivative of bambuterol, released after the decarbamylation of cholinesterase inhibited by bambuterol, also inhibits cholinesterases [22] which enables the prolonged action of bambuterol as an AD drug targeting BChE. However, repurposing bambuterol as an AD drug failed due to its poor ability to cross the blood-brain barrier Pharmaceuticals 2022, 15, 1220 3 of 23 (BBB) [35]. Wu and colleagues synthesized a series of bambuterol analogues with different amino parts of the molecule with the aim to achieve a potent and BChE selective carbamatebased inhibitor where derivatives with 2-methylbutan as a substituent on the amine were pointed out as the most promising for further evaluation as AD drugs [36,37]. Bambuterol, a biscarbamate ester of terbutaline, a drug for the treatment of asthma, has been shown to be highly selective to BChE by inhibiting it 20,000 times faster than AChE [22,34]. Considering the fact that it is already in use as a drug without side effects related to the cholinergic system, bambuterol stood out as a promising candidate for repurposing into an agent for AD treatment. This was additionally supported by the fact that the monocarbamate derivative of bambuterol, released after the decarbamylation of cholinesterase inhibited by bambuterol, also inhibits cholinesterases [22] which enables the prolonged action of bambuterol as an AD drug targeting BChE. However, repurposing bambuterol as an AD drug failed due to its poor ability to cross the bloodbrain barrier (BBB) [35]. Wu and colleagues synthesized a series of bambuterol analogues with different amino parts of the molecule with the aim to achieve a potent and BChE selective carbamate-based inhibitor where derivatives with 2-methylbutan as a substituent on the amine were pointed out as the most promising for further evaluation as AD drugs [36,37].
In this paper, we synthesized eighteen new biscarbamates with different substituents at the carbamoyl and dihydroxyaminoethyl part of the molecule with carbamate moieties in the meta position on the benzene ring. Their inhibition potency and selectivity toward human BChE or AChE were determined and analyzed using molecular modelling. As newly synthesized biscarbamates were synthesized with the aim of using them as central nervous system (CNS) active drugs, we evaluated their physicochemical properties and estimated their ability to cross the blood-brain barrier (BBB) by passive transport and tested their cytotoxicity on cells that represent models of individual organs. Additionally, the ability of tested biscarbamates to chelate biometals (Zn, Cu and Fe) was evaluated as several studies have shown that dyshomeostasis of biometals may contribute to AD pathology.

Design and Synthesis of Compounds
Biscarbamates were designed using bambuterol as a structural scaffold retaining the meta-position of carbamate groups on the benzene ring. The choice of meta disposition of carbamate groups on benzene was rationalized based on our previous study which determined that the inhibition potency and selectivity of biscarbamates towards BChE over AChE is dictated by the disposition of carbamate groups on the benzene ring, where a meta-position was preferred over an ortho-position [37,38]. In this study, we attempted to explore the impact of the carbamate group size and the size of the terminal amino group on inhibition potency and selectivity of biscarbamates to BChE over AChE. Eighteen dimethyl, ethylmethyl, diethyl, 1-pyrrolidine, 1-piperidine and methylphenyl biscarbamates with various alkyl, cyclic or aromatic terminal amines ( Figure 2) were synthesized as racemates starting from commercially available 3,5dihydoxyacetophenone, applying a slightly modified protocol described previously (Scheme 1) [36]. Briefly, the five-step synthesis included: carbamylation of 3,5dihydroxyacetophenone using carbamoyl chloride; the introduction of bromine into the In this paper, we synthesized eighteen new biscarbamates with different substituents at the carbamoyl and dihydroxyaminoethyl part of the molecule with carbamate moieties in the meta position on the benzene ring. Their inhibition potency and selectivity toward human BChE or AChE were determined and analyzed using molecular modelling. As newly synthesized biscarbamates were synthesized with the aim of using them as central nervous system (CNS) active drugs, we evaluated their physicochemical properties and estimated their ability to cross the blood-brain barrier (BBB) by passive transport and tested their cytotoxicity on cells that represent models of individual organs. Additionally, the ability of tested biscarbamates to chelate biometals (Zn, Cu and Fe) was evaluated as several studies have shown that dyshomeostasis of biometals may contribute to AD pathology.

Design and Synthesis of Compounds
Biscarbamates were designed using bambuterol as a structural scaffold retaining the meta-position of carbamate groups on the benzene ring. The choice of meta disposition of carbamate groups on benzene was rationalized based on our previous study which determined that the inhibition potency and selectivity of biscarbamates towards BChE over AChE is dictated by the disposition of carbamate groups on the benzene ring, where a meta-position was preferred over an ortho-position [37,38]. In this study, we attempted to explore the impact of the carbamate group size and the size of the terminal amino group on inhibition potency and selectivity of biscarbamates to BChE over AChE. Eighteen dimethyl, ethylmethyl, diethyl, 1-pyrrolidine, 1-piperidine and methylphenyl biscarbamates with various alkyl, cyclic or aromatic terminal amines ( Figure 2) were synthesized as racemates starting from commercially available 3,5-dihydoxyacetophenone, applying a slightly modified protocol described previously (Scheme 1) [36]. Briefly, the five-step synthesis included: carbamylation of 3,5-dihydroxyacetophenone using carbamoyl chloride; the introduction of bromine into the α-position in relation to the biscarbamate keto group; reduction of αbromoketone to β-bromoalcohol; addition of the corresponding amine on β-bromoalcohol; and the hydrochlorination of biscarbamate β-aminoalcohol. A more detailed description of the synthetic route is available in Supporting Information I. α-position in relation to the biscarbamate keto group; reduction of α-bromoketone to βbromoalcohol; addition of the corresponding amine on β-bromoalcohol; and the hydrochlorination of biscarbamate β-aminoalcohol. A more detailed description of the synthetic route is available in Supporting Information I. α-position in relation to the biscarbamate keto group; reduction of α-bromoketone to βbromoalcohol; addition of the corresponding amine on β-bromoalcohol; and the hydrochlorination of biscarbamate β-aminoalcohol. A more detailed description of the synthetic route is available in Supporting Information I.   All of the synthesized biscarbamates were obtained at a 10-70% yield as oils and had a chemical purity ≥95%. Fourteen out of eighteen compounds were new compounds. Structures of the prepared compounds were confirmed with NMR and HRMS spectra.

Kinetic Studies
The ability of biscarbamates to reduce the activity of human BChE and AChE was tested for all of the synthesized compounds. All of the compounds displayed a time-Pharmaceuticals 2022, 15, 1220 5 of 23 dependent inhibition of both, BChE and AChE, which indicated that, during the binding of the inhibitor to the active site of both enzymes, a covalent bond between the catalytic serine and the carbamate group of biscarbamate was formed. The experiments were designed so that the maximum preincubation time of the enzyme with carbamates was 30 min to assure that during that period no spontaneous decarbamylation would occur [39][40][41]. Inhibition of both cholinesterases by all of the tested carbamates followed first-order kinetics at any given inhibitor concentration ( Figure 3A,C). The inhibition potency of the compounds was expressed by the overall inhibition rate constant (k i ) that represents the first step in biscarbamate hydrolysis (enzyme inhibition scheme in Section 4) determined from the relation between the first-order rate constant (k obs ) and inhibitor concentration. For the majority of compounds, relation of k obs vs.
[CC ] deviated from linearity ( Figure 3B), allowing determination of intrinsic kinetic constants: the maximal first-order rate constant of carbamylation, k max , and dissociation constant, K i ( Table 1). For the rest of the compounds, the k obs constant was a linear function of biscarbamate concentration ( Figure 3D). a chemical purity ≥95%. Fourteen out of eighteen compounds were new compounds. Structures of the prepared compounds were confirmed with NMR and HRMS spectra.

Kinetic Studies
The ability of biscarbamates to reduce the activity of human BChE and AChE was tested for all of the synthesized compounds. All of the compounds displayed a timedependent inhibition of both, BChE and AChE, which indicated that, during the binding of the inhibitor to the active site of both enzymes, a covalent bond between the catalytic serine and the carbamate group of biscarbamate was formed. The experiments were designed so that the maximum preincubation time of the enzyme with carbamates was 30 min to assure that during that period no spontaneous decarbamylation would occur [39][40][41]. Inhibition of both cholinesterases by all of the tested carbamates followed first-order kinetics at any given inhibitor concentration ( Figure 3A,C). The inhibition potency of the compounds was expressed by the overall inhibition rate constant (ki) that represents the first step in biscarbamate hydrolysis (enzyme inhibition scheme in Section 4) determined from the relation between the first-order rate constant (kobs) and inhibitor concentration. For the majority of compounds, relation of kobs vs.
[CC′] deviated from linearity ( Figure  3B), allowing determination of intrinsic kinetic constants: the maximal first-order rate constant of carbamylation, kmax, and dissociation constant, Ki (Table 1). For the rest of the compounds, the kobs constant was a linear function of biscarbamate concentration ( Figure  3D).  Points in panels (A,C) indicate the logarithm of residual activity, while the slope of the lines determines the observed first-order rate constant k obs . On panels (B,D), constants k obs were plotted as a function of biscarbamate concentration wherefrom the second-order rate constant, k i , was calculated. Compound 1 displayed a nonlinear dependence of the k obs on biscarbamate concentration (panel (B)) allowing determination of the maximal first-order rate constant of carbamylation, k max , and dissociation constant of Michaelis type of complex between enzyme and inhibitor, K i , while k i is defined as their ratio. For biscarbamate 8, k obs was a linear function of biscarbamate concentration (panel (D)), and the slope defined the k i .

Inhibition of Butyrylcholinesterase
The tested carbamates were divided into six groups according to the substituents on carbamoyl nitrogen (Table 1). Within a particular group, the compounds differed in the size of the hydroxyaminoethyl chain attached to the benzene ring.
All of the compounds inhibited BChE with k i constants in the range 0.0144-38.0·10 6 M −1 min −1 , which makes these compounds fast or very fast BChE inhibitors. The fastest inhibition was obtained by compound 16 with piperidine in the carbamoyl and hydroxyaminoethyl chain of the compound, being almost a 10 times faster inhibitor than bambuterol. Compounds 6 and 13 were equally potent inhibitors and ten compounds (2, 4, 5, 8, 9, 12, 14, 15, 17 and 18) had an about two times lower inhibition potency of BChE than bambuterol. The tested carbamates were divided into six groups according to the substituents on the carbamoyl nitrogen. Within a particular group, the compounds differed in the size of the hydroxyaminoethyl chain attached to the benzene ring. Generally, no clear correlation between the length and/or the size of groups on the carbamoyl nitrogen and inhibition potency of compounds was detected It can be said that a lower inhibition potency can be expected for compounds with an aromatic amine on the hydroxyaminoethyl chain since the lowest inhibition potency, about 2400 times lower than that of compound 16, was determined for compounds 10 and 11 with phenyl as a substituent on the hydroxyaminoethyl chain combined with alkyl groups on carbamoyl nitrogen. But, when more detailed structure-activity analysis was performed, analyzing the influence of different substituents on carbamoyl nitrogen with the aniline on the hydroxyaminoethyl chain, it could be seen that replacement of aliphatic substituents does not affect the inhibition potency of the compounds; the increase in inhibition potency occurred when piperidine was used. That observation is in accordance with the structure-activity analysis of the same class of compounds reported earlier by Wu et al., where it was demonstrated that replacing dimethyl groups with ethylmethyl on the carbamoyl part of molecules did not affect the inhibition potency of the biscarbamates [36].
The BChE active site did not discriminate binding of structural isomers 10 and 11, which differed in the substituent on the hydroxyaminoethyl chain.
Detailed inhibition studies revealed that fifteen out of eighteen biscarbamates displayed a nonlinear dependence of the first-order rate constant (k obs ) on biscarbamate concentration allowing for the determination of the maximal first-order rate constant of carbamylation, k max , and dissociation constant of the enzyme-carbamate Michaelis type of complex, K i . For three biscarbamates 8, 10 and 11, k obs was a linear function of biscarbamate concentration. As the reciprocal value of K i represents the affinity of the enzyme to the Pharmaceuticals 2022, 15, 1220 7 of 23 compounds, it seems that the low inhibition potency of the compound 3 can be a result of BChE's low affinity to that compound. In line with that, the high inhibition potency of compounds 2, 5 and 16 can be due to BChE's high affinity towards them. On the other hand, though generally k max had a lower impact on the value of k i compared to K i , higher impact can be considered for compounds 4, 9, 13 and 14 with higher k max values than the rest of the compounds.
The inhibition potency of the tested compounds was compared to the inhibition potency of rivastigmine as an FDA-approved drug for the treatment of ADs currently in use. The majority of tested carbamates were equally or more potent (even up to 690 times) BChE inhibitors than rivastigmine. Compounds 3, 10 and 11 were the only whose inhibition potency was about three times lower than the inhibition potency of rivastigmine. The k i constant of rivastigmine determined in the study was about two-three times lower than that determined earlier [24]. The tested compounds inhibited human BChE at generally lower concentrations than dimethyl, and methylethyl biscarbamates, synthesized by Wu, inhibited horse BChE [36]. The inhibition rate constant determined for compound 5 corresponded to the value determine earlier by Wu [36].

Inhibition of Acetylcholinesterase
The tested carbamates inhibited AChE with k i constants in the range of (0.00330-4.56) 10 6 M −1 min −1 . The fastest inhibition was obtained by compound 8 with diethyl groups in the carbamoyl and tert-dimethylethyl substituent in the hydroxyaminoethyl chain of the compound. The tested carbamates were 15 to 20.727 times more potent AChE inhibitors than bambuterol [34,42]. Besides compound 8, compounds 4, 6 and 9 inhibited AChE with inhibition rate constants within the 10 6 M −1 min −1 range, compound 17 in the 10 5 M −1 min −1 range, eight compounds in 10 4 and five compounds below the 10 4 M −1 min −1 range. No relationship between inhibition potency and structure of tested compounds was observed. The AChE active site did not discriminate binding of structural isomers 10 and 11.
Detailed inhibition studies revealed that twelve compounds displayed a nonlinear dependence of the first-order rate constant (k obs ) on biscarbamate concentration, while for the rest of the carbamates, k obs was a linear function of biscarbamate concentration. According to the values of K i constants, it seems that the inhibition potency of the tested carbamates can be related to the AChE's affinity to the compounds; high inhibition potency of the compounds 8, 4, 6 and 9 is due to AChE's high affinity to these compounds. The inhibitory potential of the rest of the tested carbamates followed changes in the affinity of AChE towards the carbamates.
The majority of the tested carbamates were equally or more potent (up to 2054 times) AChE inhibitors than rivastigmine. The k i constant determined for rivastigmine in the study corresponded to that determined earlier [24]. Compounds 4, 6, 8 and 9 were about equally potent AChE inhibitors as physostigmine, a carbamate compound that was used as an FDA-approved drug for the treatment of AD, and neostigmine, a carbamate in use for the treatment of myasthenia gravis [43][44][45].

Selectivity of Inhibition
The inhibition selectivity of the tested biscarbamates was evaluated as the ratio of overall inhibition rate constants determined for BChE and AChE (Table 1). Generally, the tested carbamates were selective BChE inhibitors. The most selective were compounds 13 and 16, which inhibited BChE 1288 and 1087 times faster than AChE, respectively. Seven more compounds were more than 50 times selective to BChE than to AChE. The BChE/AChE selectivity ratio of the compound 17 was 12, and that of the rest of the compounds was lower than 10, whereas compounds 4 and 8 were up to 2.5 times more selective to AChE than to BChE or non-selective as compound 9. Regarding the groups on carbamoyl nitrogen, groups with cyclic amines as the amine part of the carbamate were more BChE selective than groups having alkyl amines. The lowest selectivity was that of compounds with a di-Pharmaceuticals 2022, 15, 1220 8 of 23 ethyl group on carbamoyl nitrogen regardless of the substituent on the hydroxyaminoethyl chain. Generally, none of the compounds reached the selectivity of bambuterol. Selectivity of compounds 13, 14, 15, 16 and 18 was up to 51 times higher than selectivity of rivastigmine [24].

Docking Analysis
Docking studies were conducted to gain insight into the structural features governing the observed differences in inhibition potency among the tested compounds. As rate of carbamylation of BChE depends on the entry of carbamate into the BChE binding site and the formation of the reversible BChE-carbamate complex close to the catalytic serine, we analyzed non-bonding interactions contributing to the stabilization of BChE-carbamate complexes and positioning of carbamates in a way that makes it susceptible to the nucleophilic attack of catalytic serine oxygen (O Ser198 ) on carbonyl group (C=O) and facilitates the release of the carbamate leaving group. A flexible docking protocol using the crystal structure of free BChE superimposed onto a network of conserved water molecules as the receptor structure was utilized to predict the binding mode of the tested compounds prior to carbamylation, i.e., the BChE-carbamate Michaelis-like complex. The protocol was applied as reported previously [45]. The predicted model complexes between the selected compounds and free BChE were analyzed for non-bonding interactions. Figure 4A illustrates the active site of the model complex between the most active carbamate from the tested series, 16, and free BChE. The predicted binding mode of carbamate 16 in the BChE active site showed that it occupied active site gorge with an orientation that fulfils prerequisites for the carbamylation reaction to occur once the carbamate slides deeper into the active site. Namely, both of its carbamate groups are oriented towards the bottom of the active site approximately occupying the choline-binding pocket and acyl-binding pocket, respectively, while the leaving group is oriented towards the gorge entrance. Carbamate 16 is predicted to engage in multiple non-bonding interactions with neighboring residues outlining the BChE active site gorge and conserved water molecules. These include water hydrogen bonds and hydrophobic-alkyl and hydrophobic-π-alkyl interactions with residues from the choline-binding pocket (Ala328, Phe329) and acyl-binding pocket (Trp231, Leu286) (the full list of predicted non-bonding interactions for all compounds can be found in Supplementary Materials). The predicted binding poses are a satisfactory illustration of one of the many transient complexes existing in the dynamic process of carbamate advancement towards the bottom of the active site gorge.
On the other hand, the predicted binding mode of one of the least active carbamates from the series, 10 ( Figure 4B), shows it is almost equally involved in non-bonding interactions with neighboring residues and conserved water molecules. These include water hydrogen bonds and hydrophobic-alkyl and hydrophobic-π-alkyl interactions predominantly with residues from the acyl-binding pocket (Trp231, Leu286, Val288). However, its predicted orientation reveals one significant difference. Namely, only one of its two carbamate groups is directed towards the bottom of the active site (the other is directed to the gorge entrance) reducing the probability of a successful carbamylation reaction, which could be the reason for its low activity. Moreover, such a spatial orientation of benzene substituents in 10 inside the BChE active site leads to the formation of an intramolecular carbon hydrogen bond between carbamate carbonyl oxygen and hydrogen attached to the carbon adjacent to cyclohexane moiety. This interaction locks-up carbamate in an unfavorable position for a carbamylation reaction, adding to its low activity.  On the other hand, the predicted binding mode of one of the least active carbamates from the series, 10 ( Figure 4B), shows it is almost equally involved in non-bonding interactions with neighboring residues and conserved water molecules. These include water hydrogen bonds and hydrophobic-alkyl and hydrophobic-π-alkyl interactions predominantly with residues from the acyl-binding pocket (Trp231, Leu286, Val288). However, its predicted orientation reveals one significant difference. Namely, only one of its two carbamate groups is directed towards the bottom of the active site (the other is directed to the gorge entrance) reducing the probability of a successful carbamylation reaction, which could be the reason for its low activity. Moreover, such a spatial orientation of benzene substituents in 10 inside the BChE active site leads to the formation of an intramolecular carbon hydrogen bond between carbamate carbonyl oxygen and hydrogen attached to the carbon adjacent to cyclohexane moiety. This interaction locksup carbamate in an unfavorable position for a carbamylation reaction, adding to its low activity.

Decarbamylation Process
Inhibition by biscarbamates occurs by the formation of a covalent bond between the carbamate group of the compound and the catalytic serine of the enzyme during which a carbamylated enzyme is formed. The decomposition of carbamylated enzyme is spontaneous and occurs by the action of water, where a free enzyme and monocarbamate were released. The recovery of enzyme activity, i.e., a spontaneous reactivation named rate of decarbamylation, is characterized by the first-order rate constant (decarbamylation rate constant; kdecarb) determined by following the spontaneous recovery of enzyme activity in time. The rate of decarbamylation depends on the substituents on carbamoyl nitrogen: rate of spontaneous decarbamylation of an enzyme inhibited by Nunsubstituted carbamates is faster than by N-substituted carbamates, and it also depends

Decarbamylation Process
Inhibition by biscarbamates occurs by the formation of a covalent bond between the carbamate group of the compound and the catalytic serine of the enzyme during which a carbamylated enzyme is formed. The decomposition of carbamylated enzyme is spontaneous and occurs by the action of water, where a free enzyme and monocarbamate were released. The recovery of enzyme activity, i.e., a spontaneous reactivation named rate of decarbamylation, is characterized by the first-order rate constant (decarbamylation rate constant; k decarb ) determined by following the spontaneous recovery of enzyme activity in time. The rate of decarbamylation depends on the substituents on carbamoyl nitrogen: rate of spontaneous decarbamylation of an enzyme inhibited by N-unsubstituted carbamates is faster than by N-substituted carbamates, and it also depends on the length and branching of substituents on the nitrogen atom [46]. In terms of drug development, the rate of spontaneous decarbamylation is important because it can be used for the estimation of the time of a drug's action.
Decarbamylation rate constants were determined for BChE and AChE, for all of the tested compounds (Table S1 in Supplementary Materials). For BChE, k decarb ranged from 0.0720 to 0.192 h −1 (k decarb = 0.136 ± 0.013 h −1 ), and 0.138 to 0.222 h −1 (0.173 ± 0.014 h −1 ) for AChE. In general, the rate of decarbamylation of BChE or AChE was very similar for all of the tested compounds, regardless of the substituents on the carbamoyl nitrogen or hydroxyaminoethyl chain. The values of k decarb constants were consistent with the literature data for N-disubstituted carbamates [40,[46][47][48][49]. Moreover, the rates of decarbamylation of BChE and AChE were similar, except for reactivation of enzymes carbamylated with cyclic carbamates, where decarbamylation of BChE was slightly faster than that of AChE.

Metal Chelating Ability
All of the synthesized biscarbamates were tested for the ability to chelate biometal ions Fe 2+ , Zn 2+ and Cu 2+ . The absorbance spectra of biscarbamates and metal mixture were recorded following 1, 30 and 60 min after mixing them and all showed changes in the spectra compared to the spectra of biscarbamates, which has been previously reported as evidence of formation of metal-compound complex: changes in absorbance intensity, bathochromic shifts of peaks or appearance of additional peaks [50,51]. Spectra recorded after a 30 min incubation period were chosen for differential spectra analysis, as after 30 min the spectra of biscarbamate-metal mixture did not change with time. Differential spectra ( Figure 5) confirmed the existence of a biscarbamate-metal complex, i.e., the biscarbamate chelating ability. than that of AChE.

Metal Chelating Ability
All of the synthesized biscarbamates were tested for the ability to chelate biometal ions Fe 2+ , Zn 2+ and Cu 2+ . The absorbance spectra of biscarbamates and metal mixture were recorded following 1, 30 and 60 min after mixing them and all showed changes in the spectra compared to the spectra of biscarbamates, which has been previously reported as evidence of formation of metal-compound complex: changes in absorbance intensity, bathochromic shifts of peaks or appearance of additional peaks [50,51]. Spectra recorded after a 30 min incubation period were chosen for differential spectra analysis, as after 30 min the spectra of biscarbamate-metal mixture did not change with time. Differential spectra ( Figure 5) confirmed the existence of a biscarbamate-metal complex, i.e., the biscarbamate chelating ability.  Furthermore, for the most selective BChE inhibitors 1, 2, 4, 5, 13, 14, 15, 16, and for most selective AChE inhibitor 8, the stoichiometry of the biscarbamate-metal complexes was evaluated using the mole fraction method. The changes of absorbance were recorded at 245 nm and plotted against the metal mole fraction (x M , Figure 6). Formation of biscarbamate-metal complexes in a ratio 1:1 with all three biometals was determined for biscarbamates 1 and 15, while biscarbamates 13 and 16 formed complexes with two biometals (compound 13 with Cu 2+ and Fe 2+ , and compound 16 with Zn 2+ and Cu 2+ ). Biscarbamates 2, 4 and 5 chelated one biometal in ratio 1:1 (complexes 2−Fe 2+ , 4−Cu 2+ and 5−Fe 2+ ). The observed stoichiometry of complexes for biscarbamates 2, 12 and 14 with two biometals (2 with Zn 2+ and Cu 2+ , 12 with Cu 2+ and Fe 2+ , and 14 with Zn 2+ and Fe 2+ ) differed from 1:1, indicating that chelation probably involved more than one molecule of biscarbamate and one biometal ion, as is the case for complex 12−Cu 2+ in which the metal-biscarbamate ratio was 3:2 or complex 14−Cu 2+ with a ratio 7:3. A similar situation was determined for biscarbamate 8, which formed complexes with Cu 2+ . For three biscarbmates, the stoichiometry of biscarbamate-metal complexes with two out of three biometals (4 with Zn 2+ and Fe 2+ , 5 with Zn 2+ and Cu 2+ , and compound 8 with Zn 2+ and Fe 2+ ) could not be determined using mole fraction method, indicating that those biscarbamates had a lower capacity to form complex metals and formed complexes only in excess concentrations of metals. The same was the case with three biscarbamates in combination with one of the biometals (12 with Zn 2+ , 13 with Zn 2+ , and 16 with Fe 2+ ). As an example of stoichiometry determination, the stoichiometry of the most selective BChE inhibitor, compound 13, is presented ( Figure 6). For the complex of 13-Cu 2+ , a sharp increase in absorbance followed by a sharp decline resulted in two straight lines for both absorption changes intersected at the metal mole fraction of 0.5 (Figure 6), pointing to the fact that Cu 2+ and compound 13 bonded stoichiometrically in a 1:1 ratio. The same trend was observed when the 13-Fe 2+ complex was analyzed. However, for the 13−Zn 2+ complex, although spectral changes and differential spectrum Figure 5B (green line) indicated that 13 could chelate Zn 2+ , the stoichiometry could not be determined.
concentrations of metals. The same was the case with three biscarbamates in combination with one of the biometals (12 with Zn 2+ , 13 with Zn 2+ , and 16 with Fe 2+ ). As an example of stoichiometry determination, the stoichiometry of the most selective BChE inhibitor, compound 13, is presented ( Figure 6). For the complex of 13-Cu 2+ , a sharp increase in absorbance followed by a sharp decline resulted in two straight lines for both absorption changes intersected at the metal mole fraction of 0.5 (Figure 6), pointing to the fact that Cu 2+ and compound 13 bonded stoichiometrically in a 1:1 ratio. The same trend was observed when the 13-Fe 2+ complex was analyzed. However, for the 13−Zn 2+ complex, although spectral changes and differential spectrum Figure 5B (green line) indicated that 13 could chelate Zn 2+ , the stoichiometry could not be determined.

The BBB Penetration Ability of Biscarbamates
The ability of biscarbamates to cross the BBB was estimated by comparing the calculated values of six physicochemical descriptors of compounds with the recommended values obtained for known CNS-active drugs [52][53][54][55]. CNS-active drugs generally have a molecular weight lower than 500 g mol −1 , moderate hydrophobicity (logP < 5), less than five hydrogen bonds donors (HBD), less than ten hydrogen bond acceptors (HBA), and less than ten rotatable bonds (RB) and are less polar (polar surface area (PSA) < 90 Å 2 ) than drugs that are not active in the CNS. The calculated values of molecular descriptors for the biscarbamates and recommended values [54][55][56] for CNS active drugs are shown in Table S2 in Supplementary Materials.
For the two compounds with the same amine substituent in the hydroxyaminoethyl chain, piperidine, i.e., compounds 2 (bisdimethyl carbamate) and 23 (bispiridinyl carbamate), in silico determined physicochemical descriptors were in the range of the upper recommended values, based on which it was estimated that they should be able to pass the BBB . For five biscarbamates 1, 3, 4, 12, and 13, the PSA value showed a minimal deviation from the upper recommended value (∆ = 0.134 Å 2 ), and it can be expected that these compounds could be able to pass the BBB since, according to the literature, some compounds with an even higher PSA can penetrate the BBB [52,56]. The physicochemical properties of the rest of the compounds showed a deviation of two or three physicochemical descriptors and are unlikely to be able to cross the BBB.
The values of the physicochemical properties for the tested biscarbamates were compared to those calculated for rivastigmine, carbamate currently in use for the treatment of AD and bambuterol. It seems that all of the tested carbamates are more likely to penetrate through the BBB than bambuterol, whose low BBB penetration can be accounted to its low logP value, lower than recommended for CNS active drugs [57]. All of the biscarbamates had higher values of molecular descriptors than rivastigmine.

Cytotoxicity
The cytotoxic effect of carbamates was evaluated on HepG2, HEK293 and SH-SY5Y in a 24 h exposure period. The concentration range was selected to correspond to the concentration range of carbamates used in in vitro kinetic experiments. The obtained results are given in Table 2. Fourteen out of eighteen biscarbamates did not exhibit hepatotoxicity, nephrotoxicity or neurotoxicity in the concentration range in which they displayed BChE or AChE inhibition activity. Four biscarbamates (10, 14, 17 and 18) were toxic to all of the tested cell lines at the concentrations in which they showed inhibition activity.

General Discussion
This study has shown that biscarbamates with a meta disposition of carbamate groups on the benzene ring are a promising structural base for the design of novel AD drugs aimed to alleviate the symptoms of the disease. This is particularly true for agents aimed to be used in the middle and late stages of the disease as the majority of tested carbamates were more than 10 times more selective to BChE than AChE (two of them more than 1000 times more selective), a feature that is considered beneficial due to the fact that ACh is mainly hydrolyzed by AChE in the early stage and by BChE in the late stage. Beside selectivity, the tested biscarbamates are very fast BChE inhibitors whose inhibition potency is equally or up to 690 times more potent than rivastigmine's, an ethylmethyl carbamate currently in use for treatment of AD and approved by the FDA. Generally, structure-activity analysis did not detect clear correlation between the length and/or the size of groups on carbamoyl nitrogen and inhibition potency of compounds. The size of aliphatic substituents did not affect the inhibition potency of the compounds, but the increase in inhibition potency can be expected when cyclic substituents were introduced. That observation is in line with that of Wu et al., who synthesized biscarbamates with aliphatic substituents on carbamoyl nitrogen [36].
The additional beneficial feature of the here presented compounds is their ability to chelate Fe 2 +, Cu 2+ and Zn 2+ , biometals whose dishomeostasis is related to the pathophysiology of AD through the production of free radicals or formation of toxic metal-Aβ plaques. In this respect, this study has identified the tested compounds as potential multi-target-directed ligands in AD treatment able to inhibit cholinesterases and chelate at least one of the biometals. The ability to chelate biometals is especially interesting in terms of ferroptosis, an iron-dependent mechanism of regulated cell death associated with an increase in oxidative stress generated by free radicals formed via the Fenton reaction. Due to its correlation to the etiopathology of AD, ferroptosis is proposed as a promising a new target for the treatment of AD [58]. The non-toxicity of most of them towards neural, liver and kidney cells in the concentration range in which they displayed BChE or AChE inhibition activity speaks in favor of the possibility of using the investigated biscarbamates as drugs. For two compounds, it was estimated that they should be able to pass the BBB by passive transport, and for five more, this ability is expected to be slightly reduced.
Development of biscarbamates as AD drugs provides the possibility of developing long-acting drugs in two ways. First, a slow decarbamylation process assures a prolonged action per se, and second, the product of the first decarbamylation process is a monocarbamate that can also inhibit one or both cholinesterases, as was earlier determined for the monocarbamate of bambuterol [22,37].
Considering all of the beneficial features, this study has singled out compounds 2 and 16 as the most promising compounds for the treatment of AD. They strongly and preferentially inhibit BChE, are non-toxic, have the potential to cross the BBB and, compared to rivastigmine, possess the ability to chelate biometals.

Synthesis of Compounds
All chemicals, reagents and solvents for the synthesis of biscarbamates were purchased from commercial sources and were used without further purification. Biscarbamates were synthesized starting from 3, 5-dihroxyacethophenone using a slightly modified protocol as previously described in Wu [37]. For the synthesis, two methods were applied: Method 1 without hydrochloride formation and Method 2 with hydrochloride formation, described in detail in Supplementary Materials. The reactions were monitored using thin-layer chromatography on Silica Gel aluminum sheets and visualized with a UV lamp. NMR spectra were recorded on Bruker AV 600 MHz and 300 MHz spectrometers, operating at 150.92 or 75.47 MHz for 13 C and 600.13 or 300.13 MHz for 1 H nuclei. Chemical shifts are quoted in ppm and are referenced to SiMe 4 as internal standard unless stated otherwise. Multiplets are abbreviated as follows: br-broad; s-singlet; d-doublet; t-triplet; qquartet; m-multiplet. Mass spectra were recorded on a 4800 Plus MALDI-TOF/TOF mass spectrometer (Applied Biosystems Inc., Foster City, CA, USA) equipped with a 200 Hz, 355 nm Nd: YAG laser. Compounds were purified using column chromatography with silica gel as the stationary phase and dichloromethane/methanol mixtures as the eluent system.
Sources of AChE and BChE were native human erythrocytes and native human plasma, respectively.
Final concentrations of carbamates were in the range of 0.01-200 µM, while substrates were 1.0 mM and 4.0 mM for ATCh and PTCh, respectively. The final dilution of AChE and BChE was 500 and 300 times, respectively.

Inhibition by Biscarbamates
In inhibition experiments, enzyme samples were incubated for up to 30 min with biscarbamates in the absence of substrate using 3-5 inhibitor concentrations ranging from 0.5-200 µM. An inhibitor was added to the reaction mixture containing DTNB, buffer and enzyme. After designated period of time, the inhibition reaction was stopped by the addition of substrate (4.0 mM PTCh for BChE, or 1.0 mM ATCh for AChE). The extent of inhibition was determined by measuring the residual activity of the enzyme. To measure the enzyme activity at zero time of inhibition, the enzyme was added to a reaction mixture containing DTNB, buffer, inhibitor and substrate immediately before the start of the measurement. With the inhibited probes, the activities of the control probes, which did not contain an inhibitor, were also measured. For each carbamate and enzyme, at least three independent experiments were performed.
Enzyme inhibition proceeds according to the following scheme:

Inhibition by Biscarbamates
In inhibition experiments, enzyme samples were incubated for up to 30 min with biscarbamates in the absence of substrate using 3-5 inhibitor concentrations ranging from 0.5-200 µM. An inhibitor was added to the reaction mixture containing DTNB, buffer and enzyme. After designated period of time, the inhibition reaction was stopped by the addition of substrate (4.0 mM PTCh for BChE, or 1.0 mM ATCh for AChE). The extent of inhibition was determined by measuring the residual activity of the enzyme. To measure the enzyme activity at zero time of inhibition, the enzyme was added to a reaction mixture containing DTNB, buffer, inhibitor and substrate immediately before the start of the measurement. With the inhibited probes, the activities of the control probes, which did not contain an inhibitor, were also measured. For each carbamate and enzyme, at least three independent experiments were performed.
Enzyme inhibition proceeds according to the following scheme: where E, CC′, ECC′, EC and C′ stand for free enzyme, inhibitor, Michaelis type complex between enzyme and inhibitor, carbamylated enzyme and monocarbamate, respectively. k+1, k−1, and kmax are rate constants of the respective reactions, while ki is the overall secondorder inhibition rate constant. The first-order rate constants (kobs) were calculated by linear regression analysis at any given inhibitor concentration [CC′]: where v0 and vi stand for the enzyme activity in the absence and in the presence of inhibitor at time t.
When the kobs was a linear function of [CC′], the slope represented the second-order inhibition rate constant (ki) When dependence of kobs vs.
[CC′] was not linear, indicating the presence of reversible enzyme-inhibitor complex, the maximum first-order inhibition rate constant (kmax) and the dissociation constants of enzyme-inhibitor complex (Ki) were determined from: where E, CC , ECC , EC and C stand for free enzyme, inhibitor, Michaelis type complex between enzyme and inhibitor, carbamylated enzyme and monocarbamate, respectively. k +1 , k −1 , and k max are rate constants of the respective reactions, while k i is the overall second-order inhibition rate constant. The first-order rate constants (k obs ) were calculated by linear regression analysis at any given inhibitor concentration [CC ]: where v 0 and v i stand for the enzyme activity in the absence and in the presence of inhibitor at time t. When the k obs was a linear function of [CC ], the slope represented the second-order inhibition rate constant (k i ) When dependence of k obs vs. [CC ] was not linear, indicating the presence of reversible enzyme-inhibitor complex, the maximum first-order inhibition rate constant (k max ) and the dissociation constants of enzyme-inhibitor complex (K i ) were determined from: Then the k i constant was determined as the ratio: All kinetic parameters were calculated using the statistical package GraphPadPrism (Graph Pad Inc. San Diego, CA, USA).

Spontaneous Decarbamylation
The rate of spontaneous recovery of activity of cholinesterase inhibited by biscarbamates, i.e., the rate of spontaneous decarbamylation, was determined by monitoring the time course of return of cholinesterase activity.
The spontaneous decarbamylation proceeds according to scheme:

Spontaneous Decarbamylation
The rate of spontaneous recovery of activity of cholinesterase inhibited by biscarbamates, i.e., the rate of spontaneous decarbamylation, was determined by monitoring the time course of return of cholinesterase activity.
The spontaneous decarbamylation proceeds according to scheme: where the carbamylated enzyme reacts with water to form the free enzyme (EH) and corresponding alcohol (COH). Before spontaneous decarbamylation measurements, BChE or AChE were incubated with 100 µM biscarbamates (incubation mixture) for 60 min to obtain the 90-100% inhibition of enzyme activity [48,49]. To avoid the possible reinhibition of the free enzyme released during the decarbamylation process, excess of biscarbamate was removed by filtering the reaction mixture through a filtration column (Strata ® C18-E, Phenomenex, Torrance, CA, USA). The time immediately after filtration was denoted as "zero" time of the onset of spontaneous decarbamylation. After a designated period of time (up to 8 h), aliquots of incubation mixture were added to a reaction mixture (buffer and 0.3 mM DTNB), and measurement of activity was started by the addition of substrate (ATCh for AChE and PTCh for BChE). An identical incubation mixture was prepared with enzyme and buffer instead of inhibitor to measure the control values of enzyme activity.
The decarbamylation rate constant (kdecarb) was calculated according to the equation: where v0 denotes the activity of the enzyme in the absence of inhibitor and vt the activity of the enzyme incubated by the inhibitor in time t. The -kdecarb is represented by the slope of the line [60][61][62][63]. All kinetic parameters were calculated using the statistical package GraphPadPrism (Graph Pad Inc. San Diego, CA, USA).

Docking Studies
For docking ligands into the enzyme receptors, the Flexible Docking protocol was used [64]. Ligands to be docked in the enzyme structures were created with ChemBio3D Ultra 13.0 (PerkinElmer, Inc., Waltham, MA, USA) and minimized using the CHARMm force field and Smart Minimizer minimization method of Minimize Ligands protocol implemented in Biovia Discovery Studio Client v18.1. (Dassault Systèmes, Vélizy-Villacoublay, France). Before starting molecular docking protocol, Prepare Ligands protocol was used to prepare ligands with regards to possible different protonation states, isomers and tautomers at pH 7.4.
The enzyme structures were prepared starting from the crystal structures of free BChE (PDB ID: 1P0I) [65]. The binding site within BChE was defined by the sphere surrounding the residues that outline the active sites gorge including those that were selected as flexible: Asn68, Asp70, Trp82, Gln119, Trp128, Glu197, Ser198, Trp231, Pro285, Leu286, Ser287, Glu325, Phe329, Phe398 and His438. The representative pose of each of where the carbamylated enzyme reacts with water to form the free enzyme (EH) and corresponding alcohol (COH).
Before spontaneous decarbamylation measurements, BChE or AChE were incubated with 100 µM biscarbamates (incubation mixture) for 60 min to obtain the 90-100% inhibition of enzyme activity [48,49]. To avoid the possible reinhibition of the free enzyme released during the decarbamylation process, excess of biscarbamate was removed by filtering the reaction mixture through a filtration column (Strata ® C18-E, Phenomenex, Torrance, CA, USA). The time immediately after filtration was denoted as "zero" time of the onset of spontaneous decarbamylation. After a designated period of time (up to 8 h), aliquots of incubation mixture were added to a reaction mixture (buffer and 0.3 mM DTNB), and measurement of activity was started by the addition of substrate (ATCh for AChE and PTCh for BChE). An identical incubation mixture was prepared with enzyme and buffer instead of inhibitor to measure the control values of enzyme activity.
The decarbamylation rate constant (k decarb ) was calculated according to the equation: where v 0 denotes the activity of the enzyme in the absence of inhibitor and v t the activity of the enzyme incubated by the inhibitor in time t. The -k decarb is represented by the slope of the line [60][61][62][63]. All kinetic parameters were calculated using the statistical package GraphPadPrism (Graph Pad Inc. San Diego, CA, USA).

Docking Studies
For docking ligands into the enzyme receptors, the Flexible Docking protocol was used [64]. Ligands to be docked in the enzyme structures were created with ChemBio3D Ultra 13.0 (PerkinElmer, Inc., Waltham, MA, USA) and minimized using the CHARMm force field and Smart Minimizer minimization method of Minimize Ligands protocol implemented in Biovia Discovery Studio Client v18.1. (Dassault Systèmes, Vélizy-Villacoublay, France). Before starting molecular docking protocol, Prepare Ligands protocol was used to prepare ligands with regards to possible different protonation states, isomers and tautomers at pH 7.4.
The enzyme structures were prepared starting from the crystal structures of free BChE (PDB ID: 1P0I) [65]. The binding site within BChE was defined by the sphere surrounding the residues that outline the active sites gorge including those that were selected as flexible: Asn68, Asp70, Trp82, Gln119, Trp128, Glu197, Ser198, Trp231, Pro285, Leu286, Ser287, Glu325, Phe329, Phe398 and His438. The representative pose of each of the docked ligands was selected based on the highest Consensus score calculated from the scoring functions used to estimate binding affinity, as implemented in the Biovia Discovery Studio Client v18.1. Score Ligand Poses protocol. Among those with highest Consensus score, the ones closest to fulfilling prerequisites for successful carbamylation reaction in terms of O Ser198 -C=O distance and alignment were finally chosen as representative poses. The parameters used in in silico protocols were configured as reported in a previously published work [45].

Metal Chelation Studies
The ability of carbamates to chelate biometals was tested using metal salts: zinc (II) chloride (ZnCl 2 ), copper chloride dihydrate (CuCl 2 ·2H 2 O) and iron dichloride tetrahydrate (FeCl 2 ·4H 2 O). All metal salts and biscarbamates were dissolved in methanol. Briefly, a fixed amount of biscarbamate (30 µM) was mixed with a fixed amount of metal salt (60 µM). The absorbance spectrum (200 to 600 nm) of such a mixture was measured in the first minute and 30 and 60 min after mixing. The spectra of metal salt and of biscarbamate were also measured [50,66].
Interaction between biscarbamates and metal was indicated by the changes of spectra of biscarbamates compared to that of biscarbamate-metal mixture and confirmed by the differential UV-VIS spectra. Differential UV-VIS spectra were obtained by numerical subtraction of the spectra of the metal and the biscarbamate from the spectra of the mixture of the metal and biscarbamate. By visual inspection of the differential spectra, 245 nm was denoted as the wavelength of maximum absorption corresponding to the formation of the biscarbamate-metal complexes.
The stoichiometry of the biscarbamate-metal complexes was determined using the mole fraction method [50,67]. The changes of absorbance (∆A) were recorded at 245 nm and plotted against the metal mole fraction (x M ). The intersection point in the plot corresponds to the mole fraction of metal in the biscarbamate-metal complex [50,51,66,67].
The study of metal chelation was performed in 1 cm quartz cuvette (final volume 1 mL) at 25 • C using UV-VIS spectrophotometer (Cary 300 spectrophotometer Varian, Inc., Melbourne, Australia). All presentations of the spectrum and stoichiometry of binding were carried out in GraphPad Prism, GraphPad Software.

In Silico Prediction of Blood-Brain Barrier (BBB) Penetration
The ability of synthesized biscarbamates to penetrate the BBB was estimated by calculating molecular descriptors important for passive transport [53,54] (optimal values of lipophilicity (log P), molecular weight (MW), a polar surface area (PSA), optimal numbers of H-bond donors (HBD) and H-bond acceptors (HBA)) and molecular flexibility characterized by the number of rotatable bonds (RB) using the Chemicalize protocol [68].The obtained results were compared to the upper and lower recommended values obtained for known CNS-active drugs [52][53][54][55][56][57].
The cytotoxic profiles of tested carbamates were determined by measuring the succinate dehydrogenase mitochondrial activity of cells exposed to them [69]. The commercially available MTS detection reagent assay was used (CellTiter 96 ® AQueous One Solution Cell Proliferation Assay, Promega, Madison, WI, USA). The procedure followed slightly modified manufacturer protocol described previously [72]. Data were presented as percentage of the inhibited cells to control untreated cells, i.e., percentage of cytotoxicity. Staurosporine was used as positive control.
Supplementary Materials: The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/ph15101220/s1, Synthesis scheme, Detailed description of synthesis method I and II, Figure S1: Structures of biscarbamates, Table S1: Decarbamylation rate constants, Table S2: Physiochemical properties of biscarbamates, 1H and 13C spectra of the compounds, HRMS spectra of the compounds, Tables S3-S20: The predicted non-bonding interactions between biscarbamtes and hBChE.