Assessments of Ceanothanes Triterpenes as Cholinesterase Inhibitors: An Investigation of Potential Agents with Novel Inspiration for Drug Treatment of Neurodegenerative Diseases

The purpose of this study was to determine the inhibitory capacity of ceanothanes triterpenes isolate from Chilean Rhamnaceae on acetylcholinesterase (AChE) and butyrylcholinesterase (BChE) enzymes. Seven ceanothanes triterpenes were isolated from aerial parts of plant material by classical phytochemical methods or prepared by the hemisynthetic method. Structures were determined by the spectroscopic method (1H-NMR and 13C NMR) and mass spectrometry (MS). AChE and BChE activity were determined by the Ellmann method for all compounds. All tested compounds exerted a greater affinity to AChE than to BChE, where compound 3 has an IC50 of 0.126 uM for AChE and of >500 uM to BChE. Kinetic studies indicated that its inhibition was competitive and reversible. According to the molecular coupling and displacement studies of the propidium iodide test, the inhibitory effect of compound 3 would be produced by interaction with the peripheral anionic site (PAS) of AChE. The compounds tested (1–7) showed an important inhibitory activity of AChE, binding to PAS. Therefore, inhibitors that bind to PAS would prevent the formation of the AChE-Aβ complex, constituting a new alternative in the treatment of Alzheimer’s disease (AD).


Introduction
Acetylcholinesterase (AChE) is the enzyme responsible for the degradation of the neurotransmitter acetylcholine (ACh) and is considered an important target in the development of molecules with insecticidal potential. On the other hand, this enzyme is the target of drugs for the treatment of symptoms that are produced by a decrease in ACh levels in patients affected by neurodegenerative diseases of the central nervous system (CNS), for example, AD [1].
While the physiological role of AChE in the neuronal transmission is well known, it remains the focus of pharmaceutical research, directed at treatments for Myasthenia Gravis, Glaucoma, and Alzheimer's disease (AD). Cholinergic deficiency has been clarified to be associated with AD [2]; therefore, one of the main therapeutic strategies is to inhibit the biological activity of AChE and thus increase the level of acetylcholine in the brain. Currently, most of the drugs used to treat AD are AChE inhibitors, including the synthetic Cholinesterase inhibitors (also called acetylcholinesterase inhibitors) are a group of molecules that block the normal breakdown of acetylcholine. Acetylcholine is the main neurotransmitter found in the body and has functions in both the peripheral nervous system and the central nervous system. For example, acetylcholine is released by motor neurons to activate muscles; acetylcholine also plays an important role in arousal, attention, learning, memory, and motivation [10][11][12].
The main use of cholinesterase inhibitors is for the treatment of dementia in patients with Alzheimer's disease. People with Alzheimer's disease have reduced levels of acetylcholine in the brain. Cholinesterase inhibitors have been shown to have a modest effect on dementia symptoms such as cognition [1, [13][14][15]. Cholinesterase inhibitors tend to cause side effects such as vasodilation; the constriction of the pupils in the eyes; increased secretion of sweat, saliva, and tears; slow heart rate; mucus secretion in the respiratory tract; and constriction of the airways.
The amyloid cascade hypothesis, which was first proposed in 1992, and which continues to be the leading model of AD pathogenesis, points to the deposition of amyloidbeta (Aβ) plaques in the brain as the initiating step of AD pathogenesis, which in turn leads to the accumulation of neurofibrillary tangles composed of hyperphosphorylated tau, synaptic and neuronal dysfunction and loss, and cognitive decline [15,16].
Since 1992, the discovery of a wide range of molecular and cellular processes that play a critical role in the development of AD has led experts to revise and expand the original hypothesis [17][18][19][20]. Genetic studies have also provided further insight into the complex mechanisms and biological pathways underlying AD, including those involving Cholinesterase inhibitors (also called acetylcholinesterase inhibitors) are a group of molecules that block the normal breakdown of acetylcholine. Acetylcholine is the main neurotransmitter found in the body and has functions in both the peripheral nervous system and the central nervous system. For example, acetylcholine is released by motor neurons to activate muscles; acetylcholine also plays an important role in arousal, attention, learning, memory, and motivation [10][11][12].
The main use of cholinesterase inhibitors is for the treatment of dementia in patients with Alzheimer's disease. People with Alzheimer's disease have reduced levels of acetylcholine in the brain. Cholinesterase inhibitors have been shown to have a modest effect on dementia symptoms such as cognition [1, [13][14][15]. Cholinesterase inhibitors tend to cause side effects such as vasodilation; the constriction of the pupils in the eyes; increased secretion of sweat, saliva, and tears; slow heart rate; mucus secretion in the respiratory tract; and constriction of the airways.
The amyloid cascade hypothesis, which was first proposed in 1992, and which continues to be the leading model of AD pathogenesis, points to the deposition of amyloid-beta (Aβ) plaques in the brain as the initiating step of AD pathogenesis, which in turn leads to the accumulation of neurofibrillary tangles composed of hyperphosphorylated tau, synaptic and neuronal dysfunction and loss, and cognitive decline [15,16].
Since 1992, the discovery of a wide range of molecular and cellular processes that play a critical role in the development of AD has led experts to revise and expand the original hypothesis [17][18][19][20]. Genetic studies have also provided further insight into the complex mechanisms and biological pathways underlying AD, including those involving amyloid precursor protein (APP), tau, immune response and inflammation, lipid transport and endocytosis, synaptic function, cytoskeletal function, and axonal transport [21].
Different studies show that pentacyclic triterpenes (ursolic acid, oleanolic acid, and taraxerol), as well as some steroids (leucisterol), are capable of inhibiting AChE [22,23]. Yoo and Park in 2012 [24] demonstrate that ursolic acid inhibits the AChE competitive/non-competitive way. Lupeol and calenduladiol isolated from Chuquiraga erinacea showed good AChE inhibition, but greater inhibitory activity was achieved with a derivative (calenduladiol disulfate), which showed much greater inhibitory activity than its precursor [25]. The literature reviewed does not report studies on ceanothanes with AChE inhibitory activity.
According to the above, the search or synthesis of molecules with inhibitor activity on AChE is an interesting alternative to developing a new therapeutic alternative by treating the degenerative illness of the central nervous system. Many phytochemical studies are bio-directed to find biopesticides of botanical origins, which have effects on the AChE of insects. The three-dimensional structure of AChE is highly conserved evolutionarily. The folding is similar when comparing the AChE structures of Homo sapiens and Drosophila melanogaster, but there are several active site and peripheral anion site residue mutations in the D. melanogaster AChE structure compared to that of H. sapiens AChE [26]. Some taxa of the Americas, such as the Rhamnaceae family, are toxic to insects, fungi, and several bacteria strains. These effects have been associated with the presence of alkaloids, phenolics, and terpenes. Our studies of Chilean flora to develop botanical insecticides, mainly plants of family Rhamnaceae, have allowed isolated different compounds with a pentacyclic triterpenes skeleton and inhibitory activity on AChE. Some of these compounds isolated from Chilean Rhamnaceae family plants have the capacity to inhibit AChE by interaction with PAS.

Results and Discussion
From the selected species of Rhamnaceae family growing in Chile, it was possible to isolate a set of pentacyclic triterpenes with the ceanothane skeleton ( Figure 2). The literature on this family of plants indicates that these types of compounds are considered as taxonomic markers since they have not been reported in other botanical families [27]. On the other hand, the studies on the biological activity associated with this type of secondary metabolites indicate that it possesses antibacterial, cytotoxic, antiprotozoal, and insecticidal activity [27][28][29][30][31]. The analysis of the literature reveals that investigations concerning the impact of ceanothane triterpenes on the inhibition of AChE have not been undertaken so far. amyloid precursor protein (APP), tau, immune response and inflammation, lipid transport and endocytosis, synaptic function, cytoskeletal function, and axonal transport [21].
Different studies show that pentacyclic triterpenes (ursolic acid, oleanolic acid, and taraxerol), as well as some steroids (leucisterol), are capable of inhibiting AChE [22,23]. Yoo and Park in 2012 [24] demonstrate that ursolic acid inhibits the AChE competitive/non-competitive way. Lupeol and calenduladiol isolated from Chuquiraga erinacea showed good AChE inhibition, but greater inhibitory activity was achieved with a derivative (calenduladiol disulfate), which showed much greater inhibitory activity than its precursor [25]. The literature reviewed does not report studies on ceanothanes with AChE inhibitory activity.
According to the above, the search or synthesis of molecules with inhibitor activity on AChE is an interesting alternative to developing a new therapeutic alternative by treating the degenerative illness of the central nervous system. Many phytochemical studies are bio-directed to find biopesticides of botanical origins, which have effects on the AChE of insects. The three-dimensional structure of AChE is highly conserved evolutionarily. The folding is similar when comparing the AChE structures of Homo sapiens and Drosophila melanogaster, but there are several active site and peripheral anion site residue mutations in the D. melanogaster AChE structure compared to that of H. sapiens AChE [26]. Some taxa of the Americas, such as the Rhamnaceae family, are toxic to insects, fungi, and several bacteria strains. These effects have been associated with the presence of alkaloids, phenolics, and terpenes. Our studies of Chilean flora to develop botanical insecticides, mainly plants of family Rhamnaceae, have allowed isolated different compounds with a pentacyclic triterpenes skeleton and inhibitory activity on AChE. Some of these compounds isolated from Chilean Rhamnaceae family plants have the capacity to inhibit AChE by interaction with PAS.

Results and Discussion
From the selected species of Rhamnaceae family growing in Chile, it was possible to isolate a set of pentacyclic triterpenes with the ceanothane skeleton ( Figure 2). The literature on this family of plants indicates that these types of compounds are considered as taxonomic markers since they have not been reported in other botanical families [27]. On the other hand, the studies on the biological activity associated with this type of secondary metabolites indicate that it possesses antibacterial, cytotoxic, antiprotozoal, and insecticidal activity [27][28][29][30][31]. The analysis of the literature reveals that investigations concerning the impact of ceanothane triterpenes on the inhibition of AChE have not been undertaken so far.

Structural Elucidation of Ceanothane Triterpenes
Compounds 1-7 are pentacyclic triterpenes with ceanothane skeleton. All of them contain a 19-isopropylidene group (IR ca 1642, 881 cm −1 ; 1 H NMR δ ca 4.71 br s; 4.89 br s; and 1.75 s Me), and a 17-carboxylic acid function (δC ca 179.0). The 1 H and 13 C NMR spectral data revealed their structural differences in ring A. We described below the details of their structural elucidation.
Compound 2 was obtained as a white powder. The spectral data indicated it to be closely related to compound 1. The 1 H NMR spectra showed signals for an aldehydic group, a vinylic methyl, five angular methyls, and two olefinic protons. The main difference was the signal at δ 9.97 ppm, attributable to the aldehyde proton, which shifted 0.31 ppm downfield with respect to the same signal in 1, and another doublet at δ 4.16 ppm, which was assumed to correspond to the signal for the proton germinal to the 3β-hydroxyl group. The 13 C NMR showed one oxygenated methine (δ 81.08 ppm). Consequently, several 2D-NMR spectra (COSY, HMQC, and HMBC) were obtained. The results served to confirm the cyclopentane nature of ring A in 2 and assigned all the signals of both carbon and proton spectra. Thus, the doublet at 9.73 ppm of the aldehyde proton showed long-range correlations with the carbon signal assigned to its vicinal methane carbon (72.3 ppm) and to the methane supporting the hydroxyl function (79.7 ppm). The signal at 4.16, corresponding to the proton germinal to the hydroxyl group, was long-range coupled to that of the aldehyde carbonyl. Additionally, the two methyl signals at 0.80 and 0.87 ppm, assigned to a gem-dimethyl moiety due to its mutual hydroxylated carbon signal at 79.7 ppm, as well as the signal for the quaternary carbon (40.4 ppm) supported both methyl groups. Finally, the splitting (dd, J = 8.5 and 4.7 Hz) of the proton germinal to the aldehyde group (2.02 ppm), along with the long-range heteronuclear correlation between the methane at 72.3 ppm and the proton signal for the angular methyl group at 0.93 ppm, confirmed the rearrangement of ring A in 2. Accordingly, compound 2 is zyziberanalic acid, also called colubrinic acid.
Compound 3 was obtained by successive crystallization as a white powder. The IR spectrum showed absorption at 2500-3100 cm −1 and 1698 cm −1 for the carboxylic acid function. 1 H NMR showed C-2βH δ 3.18 (s) and C-3αH, δ 4.99 (bs). Ceanothic acid revealed two carboxyls: one at δ178.7 and another at δ177.9, assignable to C-28 and C-1, respectively. The signals to δ151.3 and δ110.8 are assigned to C-20 and 29 vinyl carbons. The signal at δ 65.69 is due to C-2, and the signal at δ85.0 is due to C-3. The 1 H NMR spectrum of zizyberanalic acid 2 showed that C-2αH is at 2.53 (dd, J = 4.7 Hz and J = 8.7 Hz), and C-3αH is at δ 4.32 (d, J = 8.7 Hz). The 1 H NMR of ceanothic 3 showed C-2βH δ 3.18 (s) and C-3αH δ 4.99. In the literature, the 1 H NMR of isoceanothic acid showed C-2αH δ 2.60 (d, J = 8.7 Hz) and C-3βH at δ 3.5 (m). The chemical shift and coupling constant of C-2δH of isoceanothic acid is similar to that of zizyberanalic acid and differs from that of ceanothic acid. The chemical shift of isoceanothic acid, C-3βH, is resonated at upfield (δ 3.50) compared to that of zizyberanalic acid (3βH, δ 4.32). This is due to the shielding effect of 2β-COOH in isoceanothic acid on 3βH. In zizyberanalic acid, the shielding by 2β-CHO effect on C-2αH is not prominent and hence appeared downfield (3βH, 4.32).
Compound 4 was obtained as a white amorphous powder. The IR spectrum showed absorption at 3067 cm −1 (CH=C), 1721 (C=O), and 1685 (C=O). 13 C NMR indicated the presence of the 29 carbons. Carbon multiplicity, deduced from HMQC and DEPT experiments, indicated the presence of five methyl groups, nine methylene, seven methine, and seven quaternary carbons. Its spectral data were found to be similar to those of 1-norceanotha-1(3),20(30)-diene-28-oic acid, a norceanothane derivative prepared by us from ceanothic acid [32]. A cross comparison with this reference compound indicated that the signal of C-27 methyl was absent. Except for the methyl signal (δ 1.66 ppm) of the isopropenylidene substituent at C-19, this compound contained four methyl singlets, one less than the compound 3. The 1 H NMR displayed two coupled olefinic proton signals at δ 5.99 and δ 5.45, JAX = 5.7 Hz, assignable to H-1 and H-3, which were confirmed by an NOE experiment, which enhanced H-3 (δ 5.45, d) upon irradiation at the frequency of H-23 (δ 0.97,s) or H-24 (δ 0.89,s). Another irradiation at the signal of H-25 (δ 0.88 ppm) enhancing H-1 (δ 5.99,d) and H-24 also confirmed the assignment of H-1 and H-24.
Compound 5 is a white amorphous powder recrystallized from n-hexane. The proton decoupled 13 C NMR spectrum shows the presence of 32 carbon atoms. A DEPT-90 sub-spectrum indicates five methine (CH) carbons, while DEPT-135 suggests ten methylene (CH 2 ) and eight methyl groups. Nine quaternary carbons were identified in the signals that appear additionally in the proton broadband decoupled 13 C NMR. The IR spectra showed an absorption at 3562 cm −1 that corresponds to an (O-H) for hydrogen bonded hydroxyl group. The signal at 1767 cm −1 is due to a carbonyl carbon (C=O, ester), whereas that at 1729 cm −1 is due to (C=O, acid).
Compound 6 was isolated as a white powder amorphous. The IR spectrum showed absorption at 3400-3250 (COOH), 2915, 2849, 1708 (C=O), 1588, 1381, and 1026 cm −1 . 1 H NMR is very similar to ceanothic acid 3, being able to observe the disappearance of the signal δ 4.99 ppm; it corresponds to the H-3 product of the oxidation of the hydroxyl group in the same position. This correlates with the appearance of a signal in the 13 C NMR at δ216.6 ppm corresponding to C=O in C-3. Compound 7 IR spectrum showed absorption at 3073 (CH=C), 1731 (ester), and 1681 (carbonyl) cm −1 . Carbon multiplicity, deduced from HMQC (heteronuclear multiple quantum coherence) and DEPT (distortionless enhancement by polarization transfer) experiment, indicated the presence of seven methyl groups, nine methylene groups, seven methine groups, and nine quaternary carbons. The 1 H NMR spectrum of 7 displayed seven three-proton singlets at δ 0.88, 0.98, 0.99, 1.06, 1.16, and 2.03 ppm, consistent with the methyl groups attached to the quaternary carbons. The presence of the acetyl groups was confirmed by the HMBC (heteronuclear multiple bond coherence) correlation of the methyl at δ 2.03 ppm with carbonyl, whereas an isopropenyl group was assigned from NMR signals corresponding to a methyl to a methyl (δH 1.68 ppm attached to a sp2-carbon (δC 152.0 ppm) showing HMBC correlation with two vinylic protons at δ 4.58 ppm and δ4.70 ppm.

Enzyme Inhibition and Kinetics Assays
The inhibitory capacity of these triterpenes on AChE and BuChE was studied using a colorimetric method as described in the methodology [33] and galantamine hydrobromide as the reference compound. The results of these assays are summarized in Figure 3 and Table 1. The compounds showed high affinity for AChE IC 50 = 0.125 µM for the commercial enzyme and 0.146 µM for the enzyme extracted from human blood and very low affinity for BuChE (IC 50 > 500 µM). In addition, as shown in Table 1, the results obtained when evaluating the compounds with commercial AChE or with the enzyme extracted from human blood do not differ significantly.    Graphical analysis of the Lineweaver-Burk plot gives information about the binding mode. As shown in Figure 4, the lines cross the first quadrant at the same point, and Vmax decreases as the concentration of compound 3 and 7 increases. The Lineweaver-Burk plot reveals that 3 and 7 are partially competitive AChE inhibitors. Graphical analysis of the Lineweaver-Burk plot gives information about the bindin mode. As shown in Figure 4, the lines cross the first quadrant at the same point, and Vma decreases as the concentration of compound 3 and 7 increases. The Lineweaver-Burk plo reveals that 3 and 7 are partially competitive AChE inhibitors.  With the Dixon graphs ( Figure 5), it is clear that the type of inhibition for ceanothane 3 and 7 is competitive since they intersect in the second quadrant, and an increase in the concentration of substrate (inhibitor) generates a line parallel to the X-axis, which indicates that the rate does not change, even when the inhibitor concentration changes. The other ceanothanes (1, 2, 4, 5, and 6) depict profiles that were fitted to a mixed competitive mode inhibition ( Figure 3).
The above was inferred from the calculation of alpha, whose value is greater than one, this means that these compounds can act directly on the enzyme or the substrate enzyme complex.

Docking Studies
The docking binding energy for the compounds (1-7) had a good correlation with the experimental IC 50 values (0.8831 Pearson products at the 95% confidence level) ( Table 1). The compounds that show most negative downlink energy (2, 3, and 5) was bound to the active site by two or three conventional hydrogen bonds with residues near to aromatic cavity such as Ser 286, Phe 288, Arg 289, and Phe 330, and π-sigma, π-alkyl, and alkyl interactions with other residues of PAS (e.g., Tyr 70, Tyr 121, Trp 279, Phe 331, and Tyr 334). In contrast, compounds with higher binding energy and IC 50 value present none or just one hydrogen bond; such is the case of compound 4. In the case of compound 7, an unfavorable negative-negative interaction with the residue Aps 285 was observed. It is a plausible explanation for its lower affinity.

EER REVIEW 7 of 18
With the Dixon graphs ( Figure 5), it is clear that the type of inhibition for ceanothane 3 and 7 is competitive since they intersect in the second quadrant, and an increase in the concentration of substrate (inhibitor) generates a line parallel to the X-axis, which indicates that the rate does not change, even when the inhibitor concentration changes. The other ceanothanes (1, 2, 4, 5, and 6) depict profiles that were fitted to a mixed competitive mode inhibition (Figure 3). The above was inferred from the calculation of alpha, whose value is greater than one, this means that these compounds can act directly on the enzyme or the substrate enzyme complex.

Docking Studies
The docking binding energy for the compounds (1-7) had a good correlation with the experimental IC50 values (0.8831 Pearson products at the 95% confidence level) ( Table  1). The compounds that show most negative downlink energy (2, 3, and 5) was bound to the active site by two or three conventional hydrogen bonds with residues near to aromatic cavity such as Ser 286, Phe 288, Arg 289, and Phe 330, and π-sigma, π-alkyl, and alkyl interactions with other residues of PAS (e.g., Tyr 70, Tyr 121, Trp 279, Phe 331, and Tyr 334). In contrast, compounds with higher binding energy and IC50 value present none or just one hydrogen bond; such is the case of compound 4. In the case of compound 7, an unfavorable negative-negative interaction with the residue Aps 285 was observed. It is a plausible explanation for its lower affinity.
The complex obtained for compound 3 was evaluated by molecular dynamics (MD) simulation. In the last 3 ns of simulation, it was observed that this compound interacts with residues of aromatic cavity Phe 288 and Phe 331 by a stronger conventional hydrogen bond of 2.10 ± 0.28 Å and 2.13 ± 0.28 Å, respectively, and a less strong unconventional hydrogen bond with Ile 287 2.74 ± 0.26 Å. These bonds play a crucial role in the inhibitory process. Another additional interaction is between Trp 279 and ring B, D and a methyl group by a triple perpendicular pi-alkyl interaction ( Figure 6). The role of water molecules that stabilize the complex by forming hydrogen bonds into the gorge cavity and in the external zone of the active site, mostly with carboxylates anions, is also important. Fur- The complex obtained for compound 3 was evaluated by molecular dynamics (MD) simulation. In the last 3 ns of simulation, it was observed that this compound interacts with residues of aromatic cavity Phe 288 and Phe 331 by a stronger conventional hydrogen bond of 2.10 ± 0.28 Å and 2.13 ± 0.28 Å, respectively, and a less strong unconventional hydrogen bond with Ile 287 2.74 ± 0.26 Å. These bonds play a crucial role in the inhibitory process. Another additional interaction is between Trp 279 and ring B, D and a methyl group by a triple perpendicular pi-alkyl interaction ( Figure 6). The role of water molecules that stabilize the complex by forming hydrogen bonds into the gorge cavity and in the external zone of the active site, mostly with carboxylates anions, is also important. Furthermore, this functional group works like an intermediator for long electrostatic contacts with other residues such as Trp 279:O39 in the external zone, and Phe 330:O24 in the opposite site, Figure 7. Another factor that influences the inhibitor action is the steric effects of the voluminous structure of compound 3 in the gorge, staying in the PAS by the effect of constriction over Phe330-Tyr 121, see Figure 4. In the relaxation time, the root mean square deviation (RMSD) was smaller than 1.5 Å.

Propidium Iodide Displacement Assay
The PAS-AChE binding capability of ceanothanes 1-7 was evaluated by propidium iodide displacement assay at the concentration 0.75 and 1.5 mM. Propidium iodide is known for ability to specifically binding to the PAS region of AChE. This assay is based on the change in fluorescence that undergoes a test solution as a result of competition between propidium iodide and the compounds to be tested.

Propidium Iodide Displacement Assay
The PAS-AChE binding capability of ceanothanes 1-7 was evaluated by propidium iodide displacement assay at the concentration 0.75 and 1.5 mM. Propidium iodide is known for ability to specifically binding to the PAS region of AChE. This assay is based on the change in fluorescence that undergoes a test solution as a result of competition between propidium iodide and the compounds to be tested.
All ceanothanes were examined for their ability to bind to the PAS of EeAChE and competitively displace propidium iodide. The results are presented in Table 2. The compounds decreased in fluorescence intensity by 24-34% at 0.75 µM and by 40-44% at 1.5 µM. The results of the assay indicate that the ceanothanes assayed are capable of interacting specifically with the PAS region of AChE, corroborating the information provided by the docking approach. Based on the molecular docking study of the most active compound 3, it was possibly inferred that this active compound displayed a significant binding interaction with the PAS. This biological profile highlights the importance of these molecules as a prototype for the development of new protective and regenerative drugs for the potential treatment of neurodegenerative diseases. The above mentioned in the text is based on the information available from the studies about the mechanisms of inhibition of AChE that indicated that there are two important sites where the inhibitors of this enzyme, such as the PAS and the catalytic site (CAS), are joined; there is even the possibility that some do so in both [11,33]. Studies in vitro have suggested that AChE may interact with beta-amyloid to promote the deposition of amyloid plaques in the brain of patients with AD [29]. This action of AChE is primarily mediated by the PAS, through which it co-localizes with the Aβ peptide and promotes Aβ fibrillogenesis by forming a stable AChE-Aβ complex [7,34]. The binding of ligands to PAS could limit the catalytic  (C). Human AChEs RMSF from MDS trajectory. These RMSF plots contain the low bar with the secondary structure f AChE, coil (green), sheet (yellow), and helix (magenta); the red line shows the more populated residues from the cluster analysis for each MDS and the compounds' interaction.
All ceanothanes were examined for their ability to bind to the PAS of EeAChE and competitively displace propidium iodide. The results are presented in Table 2. The compounds decreased in fluorescence intensity by 24-34% at 0.75 µM and by 40-44% at 1.5 µM. The results of the assay indicate that the ceanothanes assayed are capable of interacting specifically with the PAS region of AChE, corroborating the information provided by the docking approach. Based on the molecular docking study of the most active compound 3, it was possibly inferred that this active compound displayed a significant binding interaction with the PAS. This biological profile highlights the importance of these molecules as a prototype for the development of new protective and regenerative drugs for the potential treatment of neurodegenerative diseases. The above mentioned in the text is based on the information available from the studies about the mechanisms of inhibition of AChE that indicated that there are two important sites where the inhibitors of this enzyme, such as the PAS and the catalytic site (CAS), are joined; there is even the possibility that some do so in both [11,33]. Studies in vitro have suggested that AChE may interact with beta-amyloid to promote the deposition of amyloid plaques in the brain of patients with AD [29]. This action of AChE is primarily mediated by the PAS, through which it co-localizes with the Aβ peptide and promotes Aβ fibrillogenesis by forming a stable AChE-Aβ complex [7,34]. The binding of ligands to PAS could limit the catalytic efficiency of AChE via steric and electrostatic blockage of the inhibitors' trafficking, generating conformational changes in the active site [35]. Data are mean ± SEM (n ≥ 3 experiments).
The evidence suggests that the PAS, besides its role in the allosteric regulation of AChE-catalysed hydrolysis, also mediates heterologous protein associations that contribute to cell recognition and adhesion processes during synaptogenesis, and the nucleation of amyloid peptides during the onset of AD in humans and mammalian model systems.
The ligands or bonding by PAS and the subsequent penetration of the AChE-gorge are essential, implying the role of both the peripheral anionic site and the formation of cation-π interactions in the ligand entrance. In particular, the simulation with our molecules shows the important role of this residue in anchoring the ligand at the PAS of the enzyme and in its positioning before the gorge entrance. Once the ligand is properly oriented, the formation of specific and synchronized cation-π interactions with our molecules enables the gorge penetration. Eventually, the ligand is stabilized in a free energy basin through cation-π interactions with ceanothanes.
The inhibitory effect on AChE of terpene-related compounds has been previously reported. For instance, from Buxus baleraica Wild Sauvaitre et al., the tetracyclic triterpene N-3-isobutyrylcycloxobuxidine-F was isolated [36]. This compound was able to inhibit both CAS and PAS of AChE. In early works such as Eubanks et al. [37], using computational modeling of the THC-AChE interaction, it was discovered that THC from Cannabis also enters the binding pocket of AChE PAS. Several terpenes have affinity by the hydrophobic pocket of AChE [38]. However, to date, there are no reports regarding the proposed mechanism of AChE inhibition by ceanothanes triterpenes. It is known that the Rhamnaceae family is a rich source of pentacyclic triterpenes, particularly with the ceanothane skeleton [28,39]. These inhibitors had low micromolar IC 50 (0.126-0.188 µM) values for AChE (Table 1). Alkaloids that were isolated from the active extracts of Esenbeckia leiocarpa (Rutaceae), leptomerine and kokusaginine, with IC 50 values of 2.5 and 46 µM, respectively, were observed to elicit AChE inhibitory activity [40] or galangin with an IC 50 of 120 µM [41], showing that our values are in these ranges. The kinetic analysis demonstrated that the compounds tested exhibited a competitive-type of inhibition on AChE. It is noteworthy that the Ki of our compounds is in the same range as galantamine.
Pentacyclic triterpenoids generally exhibit low oral bioavailability; in particular, compounds with the oleanane, ursane, or ceanothane skeleton fall into class IV, according to the Biopharmaceutical Classification System, due to low aqueous solubility and poor intestinal permeability. Our studies in silico about the bioavailability of the compounds assayed using SwissASDME software show good gastrointestinal absorption but a low capacity for the cross of barrier hematoencephalic (BBB), a significant characteristic or requirement of compounds is that they act in SNC. Penetrating BBB may be achieved by modifying compounds' liposolubility by transforming the functional groups present in the molecules assayed [42][43][44].

Equipment and General Experimental Procedures
NMR spectra were recorded on a Bruker spectrometer 400 MHz (Palo Alto, CA, USA) using CDCl 3 as the solvent and TMS as the internal standard. Chemical shifts are reported in δ units (ppm) and coupling constants (J) in Hz. Optical rotations were carried out on an ATAGO POLAX-2 L semiautomatic polarimeter. IR spectra were recorded on a Shimadzu FTIR-8400 infrared spectrophotometer. Silica gel (Kieselgel-mesh 0.15/0.30, Merck, Darmstadt, Germany) was used for all liquid chromatography procedures (LC). For thin layer chromatography (TLC), silica gel GF 254 was used as the stationary phase with a plate dimension of 20 cm × 20 cm × 0.20 mm for analytical TLC (Merck, Darmstadt, Germany) and 20 cm × 20 cm × 0.25 mm for semi-preparative TLC (SPTLC) (Merck, Darmstadt, Germany). Spots on chromatogram were visualized under UV light and by spraying with 5% H 2 SO 4 in methanol, and then heating at 110 • C for 5 min. Melting points were measured with a Kofler hot-stage apparatus and are reported uncorrected. Acetylcholinesterase (from Electrophorus electricus), 5,50-dithiobis-(2-nitrobenzoic acid) (DTNB), acetylthiocholine iodide, and butyrylcholinesterase (from equine serum) were purchased from Sigma-Aldrich.

Plant Materials
The aerial parts and roots of Talguenea quinquenervia (Gill.et Hook) Johnston were collected on the roadside at a pass 4.7 km NW of Portezuelo on the road to Ninhue

Extraction and Isolation
The leaves and stems (1.5 kg) of each plant under study separality were milled and fourfold extracted for 48 h with MeOH, at room temperature, and the combined macerate was filtered and evaporated under reduced pressure. The crude extract (180 g) was partitioned by being dissolved in a mixture of MeOH/H 2 O (1:2), transferred to a separating funnel, and extracted with n-hexane (15 × 300 mL) and with ethyl acetate (EtOAc) (10 × 250 mL). Both n-hexane and EtOAc fractions were concentrated under reduced pressure, and the aqueous phase was concentrated by lyophilization.

Preparation of 3-Acetoxyceanothic Acid 7
Ceanothic acid 3 (200 mg, 0.41 mmol) was warmed with 1 mL of Ac 2 O in 0.25 mL of anhydrous-pyridine for 1 h, diluted with water, allowed to stand for 5 h, and extracted with Et 2 O. The removal of Et 2 O produced a residue that crystallized from MeOH as colorless fine needles of 7 (206 mg, yield 95%).

In Vitro AChE/BChE Inhibitory Activity Assay
The Ellman assay was used to test acetylcholinesterase (AChE from Electrophorus electricus, and human blood) and butyrylcholinesterase (BChE from equine serum) inhibition activity [33]. Human blood was drawn using sterilized syringes and stored in BD Vacutainer tubes with heparin as an anticoagulant at 4 • C for 1 h. To lyse the erythrocytes, 1 mL of a 1:50 dilution of the blood sample was prepared with the non-ionic detergent triton X-100. A mixture of the DTNB (125 µL), enzyme solution (25 µL), and compound solution (25 µL) was prepared and incubated at room temperature for 20 min. All the assays were under 0.1 M KH 2 PO 4 /K 2 HPO 4 buffer, pH 8.0. The substrate was added to start the enzymatic reaction. The absorbance (λ = 405 nm) was recorded at a controlled temperature of 30 • C for 5 min. All measurements were performed fivefold as triplicate. The compounds were assayed in a dilution interval of 15 to 500 µg/mL. Galantamine was used as positive control. The percentage of inhibition was determined as follows:

Kinetic Characterization of AChE Inhibition
To investigate the inhibition mechanism of the tested compounds on AChE, a kinetic analysis was performed. The experiments were carried out using a combination of four substrate concentrations and three inhibitor concentrations with the view to obtain a double reciprocal plot (Lineweaver-Burk), in which each point is the mean of three different experiments. A parallel control with no inhibitor in the mixture allowed adjusting activities to be measured at various times.

Propidium Displacement Assay
This trial was carried out in order to evaluate the interaction of the ceanothanestriterpenes understudy with the PAS of AChE. A solution of the test compound or standard donezepil was incubated with five units of EeAChE at 25 • C for 15 min. After incubation, 50 µL of 1 µM propidium iodide solution was added to make the final assay volume 200 µL. After 15 min, the fluorescence intensity was observed at an excitation wavelength (λex = 535 nm) and an emission wavelength (λem = 595 nm) using a fluorescence plate reader Perkin Elmer VictorX2 (Perkin Elmer, Singapur). The percentage inhibition was calculated by the following expression: 100-(IFi/IF0 × 100), where IF1 and IF0 correspond to fluorescence intensities with and without the test compound, respectively. Each assay was carried out in triplicate.

Ligand Construction
All structures of the ligands were constructed using Spartan'10 1.1.0 2011, and their geometries were optimized using Density Functional B3LYP 6-31**G ab initio methods in vacuum.

Molecular Docking
The crystallized structure of Homo sapiens AChE (hAChE) PDB ID: 4M0E was obtained from the protein data bank [44]. The molecular docking of hAChE and ligands 1-7 was carried out using Smina, a fork of Autodock Vina [45]. Dockings with donepezil, rivastigmine, and tacrine were used as controls. The geometries of ligands were optimized and minimized with Avogadro using the MMFF94 force field and the conjugate gradient algorithm. The grid was centered in the middle of the gorge of the catalytic site through the ligand bound to the crystallographic structure as a reference (12 Å in every direction). With regard to the macromolecule, all residues within 6 Å from the reference ligand were counted as flexible. Exhaustiveness and seed parameters had the settings 64 and 0, respectively. The RMSD threshold for multiple clustering was set to 1 Å. The result was analyzed by a ranked cluster and binding energy (∆G), where the lower-energy and more populated cluster was selected as the best protein-ligand complex for further analyses. To test the docking accuracy, the co-crystallized ligand was re-docked under the same conditions, and an RMSD of 0.62147 Å was obtained. All experiments were made with a physiological pH.

MD Protocol
To obtain the parameters and topology files of all compounds, the most stable docking pose of the complex was employed by using the generalized amber force field (GAFF) through the Antechamber software [46]. For the protein, it was necessary to add the missing residues with Modeller [47] before establishing the protonation state of the ionizables residues using PDB2PQR 3.3.1 software from the APBS web server [48]. The parameters derived from the ff19SB force field were used for protein [49]. The system was solvated using a water box (88.5 × 92.1 × 102.7 Å, 837,647.4 Åˆ3) of an optimum point charge (OPC) water model [50] using tleap from Ambertools21 [51]. Counter-ions (Na + ) were used to maintain electroneutrality. Energy minimizations were carried out following four successive stages of minimization, wherein each minimization cycle consisted of 5000 steps using the steepest descent algorithm, followed by 5000 steps using a conjugate gradient method. In the first two-stage, a harmonic positional restraint of 500 kcal/mol xÅ 2 was applied first over the all-protein atoms to accommodate the solvent and ions and then was applied only to heavy protein atoms to minimize the hydrogens. In the following two stages, minimizations were carried out successively, reducing the restraint from 10 to 0 kcal/mol x Å 2 , except for a restriction of 10 kcal/mol x Å 2 over ligands. The minimized systems were equilibrated under NVT conditions, heating the system from 0 K to 310.15 K using the weak Langevin thermostat [52], within 300 ps. Then, 20 ns-long equilibrations in NPT conditions were carried out for each system, keeping the restriction of 10 kcal/mol x Å to the ligands at a constant temperature of 310.15 K using the Langevin thermostat and constant pressure of 1 atm using the Berendsen barostat [51,52]. Lastly, productions of 200 ns were realized for each system under the same conditions of the NPT equilibration without restraints. All the simulations were performed using pmemd.cuda of the Amber20 software [48], using periodic boundary conditions with a time step of 2 fs. The SHAKE algorithm for bond length constraints involving hydrogen atoms was used. Non-bonded interactions were calculated using a cut-off of 8 Å, and the Particle Mesh Ewald method [48] was used for treating long-range electrostatic interactions. An analysis was carried out with cpptraj and pytraj software [48].

Conclusions
The present work demonstrated potent AChE inhibitory and very low BuChE inhibitory activity of seven ceanothane-triterpenes isolated from Rhamnaceae plants growing in Chile. Furthermore, these ceanothanes not only bind to the AChE catalytic site but also to the PAS and in this way could inhibit the formation of senile plaques typical of Alzheimer's disease. Currently, we are working to increase the interaction between the active site and the peripheral anionic site in the search for the enhancement of the potency of natural cholinesterase inhibitors.