Cholinesterase Inhibition Activity, Alkaloid Profiling and Molecular Docking of Chilean Rhodophiala (Amaryllidaceae)

Amaryllidaceae plants are the commercial source of galanthamine, an alkaloid approved for the clinical treatment of Alzheimer’s disease. The chemistry and bioactivity of Chilean representatives of Rhodophiala genus from the family of Amaryllidaceae have not been widely studied so far. Ten collections of five different Chilean Rhodophiala were analyzed in vitro for activity against enzymes such as acetylcholinesterase (AChE) and butyrylcholinesterase (BuChE) as well as for their alkaloid composition by GC-MS. To obtain an insight into the potential AChE and BuChE inhibitory activity of the alkaloids identified in the most active samples, docking experiments were carried out. Although galanthamine was found neither in aerial parts nor in bulbs of R. splendens, these plant materials were the most active inhibitors of AChE (IC50: 5.78 and 3.62 μg/mL, respectively) and BuChE (IC50: 16.26 and 14.37 μg/mL, respectively). Some 37 known alkaloids and 40 still unidentified compounds were detected in the samples, suggesting high potential in the Chilean Amaryllidaceae plants as sources of both novel bioactive agents and new alkaloids.


Introduction
The vast structural and chemical diversity of natural products gives them a significant role in drug discovery [1]. Alkaloids are of particular interest in biomedicine and drug discovery research [2] due to their structural diversity and specific biological potential [3]. The Amaryllidaceae is a plant family that contains an exclusive, large and still expanding alkaloid group known as the Amaryllidaceae alkaloids, which are characterized by unique skeleton arrangements and a broad spectrum of biological activities [4,5]. Amaryllidaceae plants have been used in folk medicine for their therapeutic and Acetylcholinesterase (AChE) and butyrylcholinesterase (BuChE) are involved in the hydrolysis of the neurotransmitter acetylcholine (ACh) [14]. Acetylcholinesterase is highly selective for ACh hydrolysis, and BuChE can metabolize different substrates [15]. In the brain of AD patients, AChE activity tends to decrease while that of BuChE increases [16]. Consequently, cholinesterase inhibitors that suppress both AChE and BuChE may provide a better therapeutic response rather than AChE-selective agents [17].
Little is known on the chemistry and bioactivity of the South American endemic Amaryllidaceae genus Rhodophiala. The Rhodophiala species have high ornamental potential due to its attractive red, yellow, white or orange flowers [18][19][20][21][22]. The approximately 40 species described at present occur in Argentina, Bolivia, southern Brazil, Chile and Uruguay [22]. Rhodophiala plants have a tunicate bulb of 4-6 cm diameter, which is set 20-30 cm underground, a single umbel holding up to six flowers, each flower being 4-6 cm wide and a flower stem 35 to 50 cm long [18].
Rhodophiala species are well known ornamental plants. Propagation of the Chilean species was reported [19] as well as phylogenetic [20] and morphological studies [21,22]. The renewed interest on galanthamine sources including the search for additional alkaloids with inhibitory activity towards Acetylcholinesterase (AChE) and butyrylcholinesterase (BuChE) are involved in the hydrolysis of the neurotransmitter acetylcholine (ACh) [14]. Acetylcholinesterase is highly selective for ACh hydrolysis, and BuChE can metabolize different substrates [15]. In the brain of AD patients, AChE activity tends to decrease while that of BuChE increases [16]. Consequently, cholinesterase inhibitors that suppress both AChE and BuChE may provide a better therapeutic response rather than AChE-selective agents [17].
Little is known on the chemistry and bioactivity of the South American endemic Amaryllidaceae genus Rhodophiala. The Rhodophiala species have high ornamental potential due to its attractive red, yellow, white or orange flowers [18][19][20][21][22]. The approximately 40 species described at present occur in Argentina, Bolivia, southern Brazil, Chile and Uruguay [22]. Rhodophiala plants have a tunicate bulb of 4-6 cm diameter, which is set 20-30 cm underground, a single umbel holding up to six flowers, each flower being 4-6 cm wide and a flower stem 35 to 50 cm long [18].
Rhodophiala species are well known ornamental plants. Propagation of the Chilean species was reported [19] as well as phylogenetic [20] and morphological studies [21,22]. The renewed interest on galanthamine sources including the search for additional alkaloids with inhibitory activity towards the enzymes AChE and BuChE has prompted research work on the South American species of this family. Brazilian [23] and Argentinean wild Amaryllidaceae [24] showed high structural diversity in the alkaloid patterns and encouraged work on the Chilean species of this family. The knowledge of the chemical composition and anti-cholinesterase activity of the Chilean species belonging to genus Rhodophiala is limited. The comparative studies of their chemical composition are needed to identify the best potential sources of bioactive alkaloids in the South American species.
The aim of this work was to disclose the potential of Chilean Rhodophiala species as inhibitors of the enzymes AChE and BuChE as well as to analyze the alkaloid content and composition of native species by gas chromatography coupled to mass spectrometry (GC-MS) looking for galanthamine sources and searching for other compounds with effect of AChE and BuChE. Molecular docking studies were also carried out to investigate the affinity of the alkaloids identified in the most promising sample at the active sites of AChE and BuChE.

Results and Discussion
Twenty alkaloid extracts from five different Chilean Rhodophiala species (Figures 2 and 3) were assessed for inhibitory activity towards the enzymes AChE and BuChE. The alkaloid composition of the extracts was analyzed by GC-MS and the single alkaloid content was quantified and reported as mg GAL/g alkaloid extract (AE). The extracts have been obtained as described in Section 2.2. Briefly, the (fresh) plant material was lyophilized to obtain the dry tissues content (values included in the the enzymes AChE and BuChE has prompted research work on the South American species of this family. Brazilian [23] and Argentinean wild Amaryllidaceae [24] showed high structural diversity in the alkaloid patterns and encouraged work on the Chilean species of this family. The knowledge of the chemical composition and anti-cholinesterase activity of the Chilean species belonging to genus Rhodophiala is limited. The comparative studies of their chemical composition are needed to identify the best potential sources of bioactive alkaloids in the South American species. The aim of this work was to disclose the potential of Chilean Rhodophiala species as inhibitors of the enzymes AChE and BuChE as well as to analyze the alkaloid content and composition of native species by gas chromatography coupled to mass spectrometry (GC-MS) looking for galanthamine sources and searching for other compounds with effect of AChE and BuChE. Molecular docking studies were also carried out to investigate the affinity of the alkaloids identified in the most promising sample at the active sites of AChE and BuChE.

Results and Discussion
Twenty alkaloid extracts from five different Chilean Rhodophiala species (Figures 2 and 3) were assessed for inhibitory activity towards the enzymes AChE and BuChE. The alkaloid composition of the extracts was analyzed by GC-MS and the single alkaloid content was quantified and reported as mg GAL/g alkaloid extract (AE). The extracts have been obtained as described in Section 2.2. Briefly, the (fresh) plant material was lyophilized to obtain the dry tissues content (values included in the Table) of the different plant parts investigated. The extraction yields are calculated and reported from the lyophilized plant material.

AChE and BuChE Inhibitory Activities
The crude alkaloid samples from the different Chilean Rhodophiala were tested in vitro for AChE and BuChE-inhibitory activity. The percent dry weight of the samples, w/w extraction yields, and cholinesterase inhibition are summarized in Table 1. Galanthamine was used as a control and presented AChE and BuChE inhibition with IC50 values of 0.48 ± 0.07 and 3.70 ± 0.24 μg/mL, respectively. All the alkaloid extracts tested were active against AChE. The highest AChE inhibitory potential was found in bulbs of R. pratensis (sample Q) followed by R. splendens (sample S) with IC50 values of 3.32 ± 0.26 and 3.62 ± 0.02 μg/mL, respectively. Lowest activity was measured for the aerial part of R. pratensis (sample N) (IC50 value: 102.27 ± 6.61 μg/mL). Nearly 50% of the samples showed some activity against BuChE, with better effect for the bulbs of R. splendens (sample S) (IC50 14.37 ± 1.94 μg/mL).
The bulb (I) and leaf (J) extracts of R. montana presented moderate activity against both enzymes with better effect of (I) against AChE and (J) towards BuChE. The differences in the chemical profiles of I and J could explain these results. However, the high number of unknown alkaloids in the extracts precludes further discussion. The bulb (Q) and leaf (R) extracts of white flowering R. pratensis were active towards AChE, with IC50 of 3.32 and 8.39 μg/mL, respectively but with mild to low effect against BuChE (Table 1), reducing the pharmacological interest of this species. The high AChE and BuChE inhibitory activity of R. splendens bulb (S) and aerial parts/leaves (T) renders this plant as the most promising species in the search for active molecules for AD therapy (Table 1).

AChE and BuChE Inhibitory Activities
The crude alkaloid samples from the different Chilean Rhodophiala were tested in vitro for AChE and BuChE-inhibitory activity. The percent dry weight of the samples, w/w extraction yields, and cholinesterase inhibition are summarized in Table 1. Galanthamine was used as a control and presented AChE and BuChE inhibition with IC 50 values of 0.48 ± 0.07 and 3.70 ± 0.24 µg/mL, respectively. All the alkaloid extracts tested were active against AChE. The highest AChE inhibitory potential was found in bulbs of R. pratensis (sample Q) followed by R. splendens (sample S) with IC 50 values of 3.32 ± 0.26 and 3.62 ± 0.02 µg/mL, respectively. Lowest activity was measured for the aerial part of R. pratensis (sample N) (IC 50 value: 102.27 ± 6.61 µg/mL). Nearly 50% of the samples showed some activity against BuChE, with better effect for the bulbs of R. splendens (sample S) (IC 50 14.37 ± 1.94 µg/mL).
The bulb (I) and leaf (J) extracts of R. montana presented moderate activity against both enzymes with better effect of (I) against AChE and (J) towards BuChE. The differences in the chemical profiles of I and J could explain these results. However, the high number of unknown alkaloids in the extracts precludes further discussion. The bulb (Q) and leaf (R) extracts of white flowering R. pratensis were active towards AChE, with IC 50 of 3.32 and 8.39 µg/mL, respectively but with mild to low effect against BuChE (Table 1), reducing the pharmacological interest of this species. The high AChE and BuChE inhibitory activity of R. splendens bulb (S) and aerial parts/leaves (T) renders this plant as the most promising species in the search for active molecules for AD therapy (Table 1).

Alkaloid Identification by GC-MS
The activity of the Chilean Rhodophiala towards acetylcholinesterase is a consequence of the chemical composition of the extracts. Therefore, the alkaloid composition is a key point to understand the chemical diversity of these plants as a source of potential therapies for AD. The alkaloids occurring in the different extracts were identified by comparing their GC-MS spectra and Kovats Retention Index (RI) values with those of authentic samples. Thirty-seven known alkaloids were identified in these samples ( Figure 4). About 50% of them belong to three different alkaloid types, namely: lycorine, haemanthamine and crinine. The others belong to six different alkaloid types: tazettine, homolycorine, galanthamine, montanine, mesembrenone and narciclasine. Two unusual alkaloids known as ismine and galanthindole were also found. The occurrence and quantification of the alkaloids in the samples is summarized in Table 2 (5) 313 (14) 415        (14) 415  (16) 273            The highest alkaloid concentration was detected in the aerial parts of R. andicola (F) and in the aerial parts of R. pratensis (R) (311.1 and 274.1 mg GAL/g AE, respectively). Lowest content was found in the aerial parts of R. andicola (B) and in the aerial parts of R. pratensis (P) (133.2 and 138.1 mg GAL/g AE, respectively). In 70% of the samples, lycorine-, haemanthamine/crinine-and tazettine-type alkaloids were predominant. Lycorine-type alkaloids were present in all species with higher content in R. araucana (G) and R. montana bulbs (I) (79.0 and 78.6 mg GAL/g AE, respectively) and lowest values in the aerial parts of R. andicola (B) and (D) (8.6 and 7.5 mg GAL/g AE, respectively).
Haemanthamine/crinine-alkaloids occur in all samples except the aerial parts of R. andicola (B). However, the higher content was found in the bulbs of R. andicola from the same collection place (A) and in the aerial parts of the plant collected at Volcan Lonquimay (F). Compounds from the tazettine-type were not detected in the aerial parts and bulbs of R. montana (I and J). From the different Amaryllidaceae alkaloids groups, tazettine-type alkaloids were the main compounds in several samples, occurring in highest concentration in R. andicola (E, F) with values of 85.9 and 95.8 mg GAL/g AE of tazettine-type alkaloids in bulbs and aerial parts, respectively.
Galanthamine-type alkaloids were detected in low quantities in three species, namely R. andicola, R. araucana and R. montana (samples C, D, E, F, G, H and J) ranging between 5.1 to 19.0 mg GAL/g AE. Montanine-type alkaloids were present in all species, except R. andicola. The highest level of montanine-type compounds was detected in R. pratensis (K) (41.7 mg GAL/g AE), which presented three different alkaloids: pancratinine C, montanine and pancracine (5.3, 29.8 and 6.6 mg GAL/g AE, respectively). Mesembrenone-type was the least representative alkaloid-type. It was represented by demethylmesembrenol, detected in low quantities in three different samples of R. pratensis (K, M and Q) (7.2, 5.2 and 5.1 mg GAL/g AE, respectively). Narciclasine-type occurs in most samples in a range between 5.1 mg GAL/g AE in bulbs of R. splendens (S) to 44.5 and 32.2 mg GAL/g AE in aerial parts of R. pratensis with red flowers and leaves (N) and aerial parts of R. splendens (T), respectively. All species investigated presented ismine and/or galanthindole alkaloids, except R. montana. Forty structures occurring in the extracts could not be identified using the available databases. Three of the unidentified compounds were highly representative among the samples. The highest content of non-identified alkaloids was detected in the aerial parts of R. montana (J) and in the bulbs of R. pratensis (red flowers and without leaves) collected in the sand dunes at the sea shore (K) (126.9 and 92.6 mg GAL/g AE, respectively). The lowest content was detected in aerial parts of R. pratensis with red flowers and without leaves (P) and in aerial parts of R. splendens (T) (20.6 and 21.8 mg GAL/g AE, respectively).

Molecular Docking
In this study, R. splendens was the most active inhibitor of AChE and BuChE. Alkaloid analysis by GC-MS allowed the identification of 17 compounds in the leaf extract of R. splendens (T) including two unidentified constituents ( Table 2). The 15 alkaloids identified in the extract were evaluated for their theoretical AChE and BuChE inhibitory potential by molecular docking (Tables 3 and 4). As expected, no alkaloid identified in sample T presented better theoretical AChE inhibitory activity than galanthamine. Molecular simulation of six alkaloids identified in sample T on the 4BDS structure theoretically showed higher enzymatic inhibition against BuChE than galanthamine by 0.80 kcal/mol. 184.68 ± 11.58 >200 galanthamine (28) 0.48 ± 0.07 3.70 ± 0.24 To gain further insight into the molecular docking results, an experiment was carried out to check the AChE and BuChE inhibitory activities of 11-hydroxyvittatine (20), lycorine (9), 8-O-demethylmaritidine (16), hamayne (20b), deacetylcantabricine (17) and haemanthamine (18a) ( Table 3). The best AChE and BuChE inhibitory activities were obtained for lycorine (9) (IC 50 101.70 ± 23.79 µg/mL) and hamayne (20b) (IC 50 48.40 ± 1.13 µg/mL), respectively. However, their theoretical BuChE inhibition was not supported by the experimental assays. The difference in origin of the BuChE structure used in the molecular docking (human) and experimental assays (equine serum), as well as the inability of these compounds to arrive at the BuChE active site of the enzyme could help to explain the difference between theoretical and practical results.
Two important regions in the active sites of the hBuChE enzymes have been located: the first corresponding to the catalytic triad composed by the residues His438, Ser198, and Glu325 [25], while the second corresponds to a choline binding site (α-anionic site), composed principally by the residues Trp82 and Phe329 [25]. A graphical representation of molecular binding of 11-hydroxyvittatine (20a) and hamayne (20b) alkaloids with the hBuChE protein is presented in Figure 10. The alkaloid 11-hydroxyvittatine (20a) shows two strong interactions, hydrogen bonds, with the residues Trp82 and Trp430; however, this molecule does not present any interactions close to the catalytic triad His438, Ser198, and Glu325. On the other hand, hamayne (20b) shows one hydrogen bond interaction with the residue Gly115, an amino acid located close to the catalytic triad His438, Ser198, and Glu325. In the case of the interactions at the choline binding site (α-anionic site), both alkaloids show the same π-π stacking interaction with the residue Trp82. These molecular interactions suggest that the β-orientation of the hydroxyl group at C-3 in 11-hydroxyvittatine (20a) could theoretically increase the BuChE inhibition on the 4BDS structure by 0.49 kcal/mol, compared to the α-orientation of the hydroxyl group at the C-3 position in hamayne (20b). However, in the experimental assays, hamayne (20b) showed BuChE inhibitory activity (48.40 ± 1.13 µg/mL). It can be hypothesized that the β-orientation of the hydroxyl group at C-3 in 11-hydroxyvittatine (20a) probably makes it difficult for the compound to arrive at the catalytic triad in the active site of the BuChE. could theoretically increase the BuChE inhibition on the 4BDS structure by 0.49 kcal/mol, compared to the α-orientation of the hydroxyl group at the C-3 position in hamayne (20b). However, in the experimental assays, hamayne (20b) showed BuChE inhibitory activity (48.40 ± 1.13 μg/mL). It can be hypothesized that the β-orientation of the hydroxyl group at C-3 in 11-hydroxyvittatine (20a) probably makes it difficult for the compound to arrive at the catalytic triad in the active site of the BuChE. Studies on alkaloid composition associated with cholinesterase inhibition and binding-mode prediction have been reported [26,27]. A work on Argentinean Amaryllidaceae [24] reported the composition and acetylcholinesterase inhibition of four wild growing species, including Rhodophiala mendocina. Two R. mendocina samples collected in different locations presented similar activity towards AChE with IC50 values of 2.0 μg/mL but with relevant differences in the qualitative and quantitative alkaloid composition. The sample from the Provincia de San Juan showed high content of haemanthamine/crinamine (31.2%), tazettine (32.9%) and lycorine (13.3%) while the plant collected in the Provincia de Neuquen presented 6.8% haemanthamine/crinamine and 20.4% of lycorine, respectively. Galanthamine was found in both samples with 0.6 and 0.8% for the San Juan and Neuquen plants, respectively. In a report from acetylcholinesterase inhibitory alkaloids from Brazilian Amaryllidaceae [23] the bulbs of Rhodophiala bifida (Herb.) Traub were investigated. The activity on AChE was moderate with an IC50 value of 8.45 μg/mL, being lower than that from R. mendocina [24]. The alkaloid extract of R. bifida bulbs contained high amounts of montanine (91.94%). The alkaloid composition of the Chilean Rhodophiala ananuca (formerly: Hippeastrum ananuca) was described [28,29]. The bulbs contained phenantridine alkaloids, including hippeastidine and epi-homolycorine.
The alkaloid montanine isolated from R. bifida showed activity towards a panel of eight human cancer cell lines. According to [30], montanine at 2.5 μg/mL was more active than doxorubicine on the multi-drug resistant breast cell line NCLADR. Montanine also showed antimicrobial effect with MIC of 5 μg/mL against S. aureus ATCC 6538 and E. coli ATCC 24922 and 20 μg/mL against P. aeruginosa ATCC 27853, respectively [31]. In a screening towards Trichomonas vaginalis, dichloromethane and n-butanol extracts from Brazilian Hippeastrum species and Rhodophiala bifida showed activity against the protozoa [32]. The most active fractions contained the alkaloids lycorine and lycosinine.
In a study on the alkaloids of Zephyranthes robusta (Amaryllidaceae), the compounds isolated were evaluated as inhibitors of human cholinesterases [33]. The authors used human erythrocye AChE and serum BuChE. The compounds were tested in a range of 0.5-500 μg/mL and the inhibition was reported as IC50 values in μMolar concentration. While the activity of the reference compound galanthamine was similar in both studies, 11-hydroxyvitattine, lycorine and haemanthamine were not active on the human AChE and BuChE. 8-O-demethylmaritidine and Studies on alkaloid composition associated with cholinesterase inhibition and binding-mode prediction have been reported [26,27]. A work on Argentinean Amaryllidaceae [24] reported the composition and acetylcholinesterase inhibition of four wild growing species, including Rhodophiala mendocina. Two R. mendocina samples collected in different locations presented similar activity towards AChE with IC 50 values of 2.0 µg/mL but with relevant differences in the qualitative and quantitative alkaloid composition. The sample from the Provincia de San Juan showed high content of haemanthamine/crinamine (31.2%), tazettine (32.9%) and lycorine (13.3%) while the plant collected in the Provincia de Neuquen presented 6.8% haemanthamine/crinamine and 20.4% of lycorine, respectively. Galanthamine was found in both samples with 0.6 and 0.8% for the San Juan and Neuquen plants, respectively. In a report from acetylcholinesterase inhibitory alkaloids from Brazilian Amaryllidaceae [23] the bulbs of Rhodophiala bifida (Herb.) Traub were investigated. The activity on AChE was moderate with an IC 50 value of 8.45 µg/mL, being lower than that from R. mendocina [24]. The alkaloid extract of R. bifida bulbs contained high amounts of montanine (91.94%). The alkaloid composition of the Chilean Rhodophiala ananuca (formerly: Hippeastrum ananuca) was described [28,29]. The bulbs contained phenantridine alkaloids, including hippeastidine and epi-homolycorine.
The alkaloid montanine isolated from R. bifida showed activity towards a panel of eight human cancer cell lines. According to [30], montanine at 2.5 µg/mL was more active than doxorubicine on the multi-drug resistant breast cell line NCLADR. Montanine also showed antimicrobial effect with MIC of 5 µg/mL against S. aureus ATCC 6538 and E. coli ATCC 24922 and 20 µg/mL against P. aeruginosa ATCC 27853, respectively [31]. In a screening towards Trichomonas vaginalis, dichloromethane and n-butanol extracts from Brazilian Hippeastrum species and Rhodophiala bifida showed activity against the protozoa [32]. The most active fractions contained the alkaloids lycorine and lycosinine.
In a study on the alkaloids of Zephyranthes robusta (Amaryllidaceae), the compounds isolated were evaluated as inhibitors of human cholinesterases [33]. The authors used human erythrocye AChE and serum BuChE. The compounds were tested in a range of 0.5-500 µg/mL and the inhibition was reported as IC 50 values in µMolar concentration. While the activity of the reference compound galanthamine was similar in both studies, 11-hydroxyvitattine, lycorine and haemanthamine were not active on the human AChE and BuChE. 8-O-demethylmaritidine and hamayne were inactive on human BuChE but presented activity on erythrocyte AChE. The differences in the results can be explained by the biological source of the enzymes (human cholinesterases for [33] and electric eel AChE and horse (equine serum) BuChE in this work. For a better comparison of results, the use of enzymes from the same biological source should be recommended.
In summary, the AChE and BuChE inhibitory activity of the Chilean Rhodophiala species investigated led to an interesting source of inhibitors that do not contain the alkaloid galanthamine. Our results suggest that Chilean Rhodophiala could be a promising source of new alkaloids with effect towards cholinesterases. The difficulty in finding a species with high activity against AChE and BuChE, the similarity of the AChE and BuChE inhibitory values, the low complexity of the alkaloid profile of aerial parts of R. splendens, together with the absence of galanthamine-type alkaloids in this sample, prompted us to further explore the results.

Plant Material
The samples were collected in central-southern Chile and were identified following the reference [19,34].

Extraction
The freshly collected plant material was cleaned and separated into bulbs and aerial parts, frozen and lyophilized before extraction. The dry weight percentage was determined. The lyophilized plant material was extracted with MeOH under sonication for 10 min (3×) changing the solvent each time. The plant to solvent ratio ranged from 1:10 to 1:60 and was selected according to the volume of plant material for extraction. The combined MeOH solubles were taken to dryness under reduced pressure to afford the crude extracts. The crude extracts were then acidified to pH 3 with diluted H 2 SO 4 (2%, v/v) and the neutral material was removed with Et 2 O. The aqueous solutions were basified up to pH 9-10 with NH 4 OH (25%, v/v) and extracted with EtOAc to provide the alkaloid extracts which were used for all experiments (enzyme inhibition assays and chemical analysis by GC-MS).

Acetylcholinesterase (AChE) and Butyrylcholinesterase (BuChE) Inhibitory Activity
Cholinesterase inhibitory activities were determined according to [35] with some modifications [36]. Stock solutions with 518U of AchE from Electrophorus electricus (Merck, Darmstadt, Germany) and BuChE from equine serum (Merck, Darmstadt, Germany), respectively, were prepared and kept at −20 • C. Acetylthiocholine iodide (ATCI), S-butyrylthiocholine iodide (BTCI) and 5,5 -dithiobis (2-nitrobenzoic acid) (DTNB) were obtained from Merck (Darmstadt, Germany). Fifty microliters of AChE or BuChE (both enzymes used at 6.24 U) in phosphate buffer (8 mM K 2 HPO 4 , 2.3 mM NaH 2 PO 4 , 0.15 NaCl, pH 7.5) and 50 µL of the sample dissolved in the same buffer were added to the wells. The plates were incubated for 30 min at room temperature. Then, 100 µL of the substrate solution (0.1 M Na 2 HPO 4 , 0.5 M DTNB, and 0.6 mM ATCI or 0.24 mM BTCI in Millipore water, pH 7.5) was added. After 10 min, the absorbance was read at 405 nm in a Labsystem microplate reader (Thermo Fischer, Waltham, MA, USA). Enzyme activity was calculated as percent compared to a control using buffer without any inhibitor. Galanthamine served as positive control. In a first step, samples were assessed at 10, 100 and 200 µg/mL towards both enzymes. Samples with an IC 50 > 200 µg/mL were considered inactive. Samples with an IC 50 < 200 µg/mL were further analyzed to determine the IC 50 values. The cholinesterase inhibitory data were analyzed with the software Microsoft Office Excel 2010 (Microsoft, Redmond, WA, USA).

Equipment
The equipment used for the identification and quantification of the alkaloids was a GC-MS 6890N apparatus (Agilent Technologies, Santa Clara, CA, USA) coupled with MSD5975 inert XL operating in the electron ionization (EI) mode at 70 eV. A Sapiens-X5 MS column (30 m × 0.25 mm i.d., film thickness 0.25 µm) was used. The temperature gradient was as follows: 12 min at 100 • C, 100-180 • C at 15 • C/min, 180-300 • C at 5 • C/min and 10 min hold at 300 • C. The injector and detector temperatures were 250 and 280 • C, respectively, and the flow-rate of carrier gas (He) was 1 mL/min. Two mg of each alkaloid extract was dissolved in 1 mL of MeOH:CHCl 3 (1:1, v/v) and 1 µL was injected using the splitless mode. Codeine (0.05 mg/mL) was used as an internal standard in all the samples.

Alkaloids Identification
Amaryllidaceae alkaloids occurring in the samples were identified by comparison of the Rt, fragmentation patterns and data interpretation of the spectra. The database used was built using single alkaloids previously isolated and identified by spectroscopic and spectrometric methods (NMR, UV, CD, IR, MS) in the Natural Products Laboratory, Universidad de Barcelona, the NIST 05 Database and literature data [37][38][39][40][41][42].

Alkaloid Quantification
To quantify the single constituents, a calibration curve of galanthamine (10, 20, 40, 60, 80 and 100 µg/mL) was used. The same amount of codeine (0.05 mg/mL) was added to each solution sample as an internal standard. The peak areas were manually obtained considering selected ions for each compound (usually the base peak of their MS, i.e., m/z at 286 for galanthamine, at 299 for codeine). The ratio between the values obtained for galanthamine and codeine in each solution was plotted against the corresponding concentration of galanthamine to obtain the calibration curve and its equation (y = 0.0224x − 0.2037; R 2 = 0.9977). All data were standardized to the area of the internal standard (codeine) and the equation obtained for the calibration curve of galanthamine was used to calculate the amount of each alkaloid. Results are expressed as mg GAL, which was finally related to the alkaloid extract weight. As the peak area does not only depend on the corresponding alkaloid concentration but also on the intensity of the mass spectra fragmentation, the quantification is not absolute. However, the methodology is considered suitable to compare the specific alkaloid amount between samples [39,42].

Molecular Docking
The molecular docking simulations for the alkaloids identified in Rhodophiala splendens, the most promising species found in this study, were performed to investigate the binding mode into the active site of two different enzymes, namely Torpedo californica AChE (TcAChE) and hBuChE: proteins with PDB codes 1DX6 [43] and 4BDS [44], respectively. The 3D-structures of the alkaloids were drawn using the Chemcraft program [45], and then submitted to a geometrical optimization procedure at PBE0 [46]/6-311+g* [47] level of theory using the Gaussian 09 program [48]. All optimized alkaloids were confirmed as a minimum on the potential energy surface. The docking simulations for the set of optimized ligands were performed using the AutoDock v.4.2 program [49]. AutoDock combines a rapid energy evaluation through pre-calculated grids of affinity potentials with a variety of search algorithms to find suitable binding positions for a ligand on a given macromolecule. To compare the results from the docking simulations, the water molecules, cofactors, and ions were excluded from each X-ray crystallographic structure. Likewise, the polar hydrogen atoms of the enzymes were added, and the non-polar hydrogen atoms were merged. Finally, the enzyme was treated as a rigid body. The grid maps of interaction energy for various atom types with each macromolecule were calculated by the auxiliary program AutoGrid, choosing a grid box with dimensions of 70 × 70 × 70 Å around the active site, which was sufficiently large to include the most important residues of each enzyme. The docking searches for the best orientations of the ligands binding to the active site of each protein were performed using the Lamarckian Genetic Algorithm (LGA) [50]. The LGA protocol applied a population size of 2000 individuals, while 2,500,000 energy evaluations were used for the 50 LGA runs. The best conformations were chosen from the lowest docked energy solutions in the cluster populated by the highest number of conformations. The best docking complex solutions (poses) were analyzed according to the potential intermolecular interactions (ligand/enzyme), such as hydrogen bonding and the cation-π, π-π stacking.