In Vitro and In Silico Anti-Acetylcholinesterase Activity from Macaranga tanarius and Syzygium jambos

Macaranga tanarius (MT) and Syzygium jambos (SJ) are pharmacologically reported to have anti-oxidant, anti-inflammatory, and anti-diabetic effects, and can be neuroprotective agents. Our previous work revealed that MT and SJ exhibited 76.32% and 93.81% inhibition against acetylcholinesterase (AChE) at 50 μg/mL final concentration in their ethyl acetate and hexane fractions, respectively. This study was aimed to investigate the bioactive constituents of MT and SJ and their molecular mechanism toward AChE inhibition. Bioassay-guided isolation afforded prenylflavonoids 1–3 from MT and anacardic acid derivatives 4 and 5 from SJ that were confirmed by NMR and MS data. Compound 5 exerted the strongest anti-AChE potential (IC50: 0.54 μM), followed by 1, 4, 3, and 2 (IC50: 1.0, 2.4, 6.8, and 33 μM, respectively). In silico molecular docking revealed 5 formed stronger molecular interactions including three H-bonds than its derivative 4 based on the saturation of their alkyl chains. The addition of a five carbon-prenyl chain in 1 increased the number of binding interactions, justifying its greater activity than derivatives 2 and 3. This research reflects the first report of AChE inhibitors from these species, thereby adding pharmacological values to MT and SJ as potential remedies in neuroprotection.


Introduction
Macaranga tanarius (L.) Mull.Arg. (MT), from the Euphorbiaceae family is a fastgrowing evergreen tree that can reach a height up to 20 m [1]. It is native to Malaysia as well as other tropical regions such as Africa, Madagascar, Southeast Asia, Australia, and the Pacific. MT traditionally has been used in folk medicine as anti-pyretic and anti-tussive agents to relieve fever and cough [2]. The anti-inflammatory properties of its leaves are used to treat wounds and swellings [2,3]. Macaranga species are rich in flavonoids and stilbenes, which are thought to be responsible for a variety of pharmacological activities such as anticancer, anti-oxidant, anti-microbial, anti-inflammatory, and anti-diabetic properties [4][5][6][7].
Syzygium jambos (L.) Alston (SJ), also known as rose apple is widely distributed in Central America, sub-Saharan Africa, and Asia, including Malaysia [8]. It is a member of the Myrtaceae family and has long been used in folk medicine for its anti-pyretic and anti-inflammatory properties. The leaves are decocted for diuretic, rheumatism, and sore eyes relief, while the seeds are used to treat diarrhea, dysentery, diabetes, and catarrh [8,9]. The bark extract is used to soothe asthma and bronchitis [8,9]. SJ was reported to comprise primarily gallic acid, vitamin C, cyanidins, tannins, and anthocyanins that contribute to its anti-diabetic, anti-oxidant, and anti-inflammatory properties [10][11][12].
Our ongoing search for acetylcholinesterase (AChE) inhibitors has been focused on Malaysian plants exhibiting anti-inflammatory, anti-oxidant, and anti-diabetic activities that are considered as neuroprotective agents [13]. AChE is crucial in both central and peripheral nervous systems, as it degrades the neurotransmitter acetylcholine homeostasis [14]. However, in patients inflicted with Alzheimer's disease, the progression of this illness is associated with low levels of acetylcholine, possibly due to increased hydrolysis of acetylcholine molecules by AChE [14]. Therefore, to achieve a homeostatic neurotransmitter equilibrium, AChE inhibitors play a critical role in preventing the AChE activity from breaking down more neurotransmitters. Physostigmine and tacrine are among the early discovered drugs used to inhibit cholinesterase, while rivastigmine, galanthamine, and donepezil are the currently approved AChE inhibitors [15].
In our previous report, MT and SJ leaves revealed strong AChE inhibition in the ethyl acetate (76.32%) and hexane (93.81%) extracts, respectively, at 50 µg/mL final concentration (Supplementary Materials Figure S1) [13]. As both species were reported to promote antiinflammatory, anti-oxidant, and anti-diabetic properties, their anti-AChE potential has yet to be addressed; this study was initiated to isolate the active principles from these extracts by anti-AChE assay-guided chromatography techniques. The structures were identified on the basis of NMR and MS spectral data, and the enzyme-ligand interactions were predicted based on Autodock 4.3 molecular docking.

Isolation and Identification of Active Constituents
Our previous screening study evaluated 177 Malaysian plant extracts for their anti-AChE potential [13]. Among the 18 plant extracts reported to show strong anti-cholinesterase activity at 50 µg/mL final concentration, Artocarpus altilis (AA) and MT exhibited more pronounced inhibitory activity in the ethyl acetate fraction, while SJ exhibited stronger activity in the hexane fraction ( Figure S1). Identification of AA's bioactive principles was hampered by the low amount of sample collected and, therefore, was not pursued. Moreover, the bioactive principles from butanol and water fractions constitute high molecular weight polyphenols [13]. Our interest to identify low molecular weight active constituents, which have not been reported from MT and SJ before, led to this bioassay-guided isolation.
The structures of isolated compounds ( Figure 1) were identified based on 1D and 2D NMR and MS spectral data ( Figures S4-S34). Their physicochemical properties including 1 H and 13 C NMR chemical shifts were compared and identical to the reported data (Tables S1-S6). Compounds 1, 2, and 3 from MT were deduced to have a flavanone backbone structure based on their UPLC-MS analysis revealing λ max around 233, 290, and 334 nm ( Figure S35). Compounds 4 and 5 from SJ showed UV characteristics of λ max at 243 and 311 nm ( Figure S36). Compound 1 was obtained as a light yellow amorphous solid. The 1 H NMR spectrum of 1 showed three singlet methyl signals (δ H = 1.54, 1.60, and 1.75), a doublet signal integrating two methyl protons (δ H = 1.65), one methylene α signal to the carbonyl (δ H = 2.59, dd and 3.09, dd), as well as seven methine signals above 5.00 ppm (Table S2). The observation of four methylene signals suggested the presence of side chains. Its 13 C NMR indicated 30 carbon signals including a carbonyl signal (δ C = 198.4) and one oxymethine signal (δ C = 77.9). These data were in concordant with the MS spectrum of 1. Its ESI-MS revealed quasi-molecular ions observed at m/z 493.8 [M+H] + and 491.9 [M-H] -( Figure S30), which suggested a formula of C 30 H 36 O 6 . The 2D NMR correlations suggested 1 had a flavanone skeleton with a B-ring substitution pattern. The presence of a geranyl chain and a prenyl chain were observed, where the methylene (δ H = 3.21, d) was correlated to A-ring carbons (δ C = 166.1, 109.8, and 162.6) and the methylene (δ H = 3.47, d) was correlated to B-ring carbons (δ C = 128.3 and 144.5). These characteristics proposed 1 was either nymphaeol C or solophenol A. It is evident that proton H-1 " was correlated to carbon C-7, C-8, and C-9, and proton H-6 was correlated to carbon C-5, C-7, C-8, and C-10, confirming the position of the geranyl chain that was attached to C-8. Compound 1 hence was identified as solophenol A. It was first reported in propolis collected from the Solomon Islands [16].
Compound 2 showed a similar characteristic of geranylation at C-2 on ring B to 1, but lacked a prenyl chain on ring A. Its ESI-MS spectrum showed quasi-molecular ions at m/z 425.  Figure S31). These characteristics indicated 2 was nymphaeol B, a known compound identified in a few Macaranga species [4]. The absence of ortho-coupling in 3 suggested it had no vicinal protons in ring B. Based on HSQC correlations of 3, the singlet proton (δ H = 6.79) integrating a total of two hydrogen atoms was assigned to C-2 and C-4 , while the singlet proton (δ H = 6.91) was assigned to C-6 ( Table S3). The doublet proton H-1" expressing HMBC correlations with C-5, C-6, and C-7 indicated that geranylation occurred at C-6 on ring A. Compound 3 showed similar MS characteristics as 2, in having a geranyl chain ( Figure S32 and Table S4), and therefore was identified as schizolaenone C. It was first isolated from Schizolaena hystrix [17]. Compounds 1 and 3 were newly identified from Macaranga species.

In Vitro Anti-AChE Activity
Compound 1 manifested the lowest IC 50 value against AChE at 1.0 µM among the prenylflavonoids, followed by 3 and 2, which showed moderate to low AChE inhibition ( Table 1). The different position of geranyl group appeared to influence the activity of 2 and 3, where geranylation on ring A in 3 might contribute to a better inhibition than on ring B. However, 1 showed greater activity by seven times than 3, which could be postulated due to the presence of a prenyl group in 1. On the other hand, 5 demonstrated the strongest anti-AChE activity among all isolated compounds with an IC 50 value of 0.5 µM (Table 1). Anacardic acid derivatives are the major constituents of cashew nutshell liquid (CNSL), whereby different degrees of saturation were found on the side chains [20]. It was suggested that the saturation level found on their aliphatic chain could affect biological activities. Compound 4, which contained an extra double bond on the aliphatic chain, exhibited a weaker AChE inhibition by four times than 5. Higher saturated anacardic acids such as 5 seemed to have a flexible structural conformation during enzyme inhibition, but further molecular docking is necessary to confirm.

In Silico Molecular Docking
Compounds 1-3 showed binding energies around −10.5 to −12.6 kcal/mol, while the binding energies of 4 and 5 fell within the range −7.0 to −8.0 kcal/mol, implying ideal values for the formation of a stable complex with the target enzyme were observed ( Table 2). Their varying molecular weights may explain the differences in their binding energy values.
Compound 1 established three polar interactions with TYR 70, TRP 84, and HIS 440, and seven non-polar interactions with TRP 84, TYR 121, TRP 279, PHE 330, and PHE 331 (Figure 2a; Table 2). In addition, all compounds 1-3 interacted with four similar key TcAChE residues: TRP 84, TYR 121, PHE 330, and PHE 331 (Figure 2a; Table 2). Compound 2 revealed seven interactions including two conventional H-bonds. On the other hand, 3 showed eight interactions including three H-bonds, and hence justified its superior inhibitory activity than 2. It is worth mentioning that 1 had shorter distances of H-bonds than 3. In addition, the prenyl group of 1 formed a few contacts with TcAChE residues, and therefore may help in creating stronger molecular interactions with the enzyme.
Compounds 4 and 5 formed three similar interactions with TRP 84 (non-polar), GLY 118, and GLY 119 (Figure 2b; Table 2). However, 5 made one extra polar H-bond with ALA 201 and a π-alkyl non-polar interaction with TYR 121. Its saturated alkyl chain was observed to create more flexible conformation than 4, allowing more interactions especially H-bonds with the enzyme residues. Anacardic acid derivatives formed interactions mostly at the oxyanion hole (OH), which is located near the base cavity. Conversely, prenylflavonoids favored binding interactions at the peripheral active site (PAS) located near the gorge mouth, where the entrance was blocked. Further molecular dynamics study is necessary, but these results rationalized their in vitro activity.
Positive control physostigmine was obtained from TCI Tokyo Chemical Industry (Tokyo, Japan). All other reagents acquired from commercial sources were of analytical grade.

In Vitro Anti-Acetylcholinesterase Assay
The plant fractions and isolated compounds were evaluated for their anti-AChE potential in triplicate based on the Ellman's method described previously [13]. Physostigmine was used as a positive control.

In Silico Molecular Docking
The crystal structure of Torpedo californica acetylcholinesterase (TcAChE) in complex with galanthamine [21] (PDB: 1W6R, 2.05 Å) was retrieved from the Research Collaboratory for Structural Bioinformatics (RCSB) Protein Data Bank (www.rcsb.org, accessed on 24 November 2021). TcAChE was selected in this study based on literature reviews [22,23], where it depicts closely the in vitro model. The water molecules and heteroatoms were removed using UCSF Chimera version 1.15. To add hydrogens, reconstruct missing atoms, and assign atomic charges, the PDB2PQR (https://server.poissonboltzmann.org/pdb2pqr, accessed on 24 November 2021) and MolProbity (http://molprobity.biochem.duke.edu, accessed on 24 November 2021) web services were used. The protonation state for the ionizable groups of the protein was set at 7.40 by using the most used empirical pKa predictor (PROPKA3). The protein was then added with hydrogen atoms and Kollman charges using AutoDock Tools 1.5.6. The 3D ligand structures were acquired from Pub-Chem (https://pubchem.ncbi.nlm.nih.gov, accessed on 24 November 2021), structurally minimized, and added with Gaisteger charges using UCSF Chimera. The ligand structure then was prepared with torsion using AutoDock Tools 1.5.6. A control docking using galanthamine (heteroatom GNT) was conducted, whereby the docking grid was set to 50, 50, and 50 of X-, Y-, and Z-dimensions with 0.375 Å grid point spacing, respectively. As a result, the coordinates of central grid points were set at x = 3.518, y = 65.122, z = 64.481. The protein-ligand docking simulation was performed using AutoDock 4.2, allowing 150 confirmations of genetic algorithm (GA) run with Lamarckian GA output. The docking simulation generating the lowest free energy of binding (FEB) with the highest cluster was selected and visualized using Discovery Studio Visualizer. The molecular interactions were observed and discussed based on our comparisons with the literature reviews [22,24].

Conclusions
This study highlights the pharmacological relevance of MT and SJ as AChE inhibitors in neuroprotection. Their first report of anti-AChE activity was represented here. Compound 5 exerted the strongest anti-AChE activity at 0.54 µM of IC 50 . The in silico molecular docking supported the in vitro data, where prenylation in prenylflavonoids and saturation of the alkyl chain in anacardic acid derivatives may modulate the molecular interactions with the enzyme.