Amaryllidaceae Alkaloids of Norbelladine-Type as Inspiration for Development of Highly Selective Butyrylcholinesterase Inhibitors: Synthesis, Biological Activity Evaluation, and Docking Studies

Alzheimer’s disease (AD) is a multifactorial neurodegenerative condition of the central nervous system (CNS) that is currently treated by cholinesterase inhibitors and the N-methyl-d-aspartate receptor antagonist, memantine. Emerging evidence strongly supports the relevance of targeting butyrylcholinesterase (BuChE) in the more advanced stages of AD. Within this study, we have generated a pilot series of compounds (1–20) structurally inspired from belladine-type Amaryllidaceae alkaloids, namely carltonine A and B, and evaluated their acetylcholinesterase (AChE) and BuChE inhibition properties. Some of the compounds exhibited intriguing inhibition activity for human BuChE (hBuChE), with a preference for BuChE over AChE. Seven compounds were found to possess a hBuChE inhibition profile, with IC50 values below 1 µM. The most potent one, compound 6, showed nanomolar range activity with an IC50 value of 72 nM and an excellent selectivity pattern over AChE, reaching a selectivity index of almost 1400. Compound 6 was further studied by enzyme kinetics, along with in-silico techniques, to reveal the mode of inhibition. The prediction of CNS availability estimates that all the compounds in this survey can pass through the blood-brain barrier (BBB), as disclosed by the BBB score.


Introduction
Alzheimer's disease (AD) is a devastating neurodegenerative age-dependent disorder characterized by memory loss and progressive cognitive impairment [1]. AD currently affects over 50 million people worldwide and is estimated to nearly triple by 2050 [2]. AD influences not only the quality of life of AD patients but also their families, with an enormous economic burden as well.
The etiology of AD has historically been ascribed to disrupted cholinergic neurotransmission [3], with more up-to-date explanations comprising excessive production of amyloid-β (Aβ) toxic species [4], dysfunction and accumulation of post-translationally modified tau protein [4], metal ions imbalance [5], and oxidative stress [6]. Although several clinical candidates have been generated upon these hypotheses over the past three decades of intensive research, current therapy of AD's mild-to-moderate stage still relies almost exclusively on the administration of acetylcholinesterase (AChE; EC 3.1.1.7) inhibitors [7], represented by galantamine, donepezil, and rivastigmine [8]. Congruently, they halt the degradation of acetylcholine (ACh), thus elevating its levels at the synapses and facilitating the cholinergic receptors' stimulation. It is noteworthy that in the later stages of the disease, the hydrolysis of ACh is more preferentially under the control of another cholinesterase enzyme, namely butyrylcholinesterase (also known as pseudocholinesterase; BuChE; EC 3.1.1.8) [9]. Notably, the crosstalk between the cholinergic system and other sub-pathologies of AD has been proven, pointing out the privileged role of cholinesterase inhibitors. For example, the cholinergic depletion promoted Aβ burden and tauopathy, both significantly contributing to cognitive decline [10]. Some studies have disclosed that the treatment with cholinesterase inhibitors can shift amyloid precursor protein (APP) processing to the non-amyloidogenic pathway [11], thus generating non-toxic Aβ species. Professor Inestrosa's group first reported the association between AChE and assembly of Aβ [12]. This study constitutes the basis for intensive research in the field of so-called dual binding site AChE inhibitors that can slow down the assembly of Aβ via the interaction with the peripheral anionic site (PAS) of the enzyme [12]. BuChE has also been suggested to participate in the formation and maturation of Aβ plaques [13,14]. Moreover, some studies have shown that BuChE may contribute to the pathogenesis of type 2 diabetes mellitus (T2DM) by causing insulin resistance; the latter is considered as one of the major risk factors contributing to AD onset [15,16]. Recently, several new structural scaffolds have been designed, synthesized, and studied as new cholinesterase inhibitors [17][18][19].
Natural products (NPs) represent one of the major sources for the development of novel drugs for a variety of disorders, including neurodegenerative ones [20,21]. Drug development starting from NPs can be considered a highly challenging task given the NPs structural complexity, a number of chiral centers, laborious multi-step syntheses, and unfavorable ADMET properties [22]. Only some of the NPs are suitable bioactive pharmacophores to be directly used as clinical drugs. Indeed, after deciphering structureactivity relationships (SAR), NPs can be considered as lead compounds requiring further optimization like potency/selectivity/water solubility improvement, metabolic stability enhancement, and reducing side effects [23]. From the large family of NPs, a significant number of them have been identified as AChE/BuChE inhibitors [24][25][26][27]. With respect to neurodegenerative disorders, only galantamine, the Amaryllidaceae alkaloid isolated from the bulbs of snowdrops, e.g., Galanthus woronowii, is currently used in the management of AD symptoms.
In the present study, we report the synthesis of a pilot series (compounds 1-20; Figure 1) of highly selective hBuChE inhibitors. Their structure has been inspired by Amaryllidaceae alkaloids (AAs) of the belladine-type, namely carltonine A and B (Figure 1, [28]). We explored the SAR by establishing their AChE/BuChE inhibition properties. We have also determined the prolyloligopeptidase (POP; EC 3.4.21.2) inhibition ability of this class of compounds since its blockade is associated with cognitive improvement. The inhibition of POP mediates the signaling cascade by releasing Ca 2+ at the endoplasmatic reticulum membrane, which may induce memory enhancement [29]. Some highly selective POP inhibitors entered clinical trial testing for neurodegenerative disorders, but these were terminated, mostly due to their shortcomings related to the compound's pharmacokinetic profile. Besides, we have also predicted the CNS availability of the selected compounds. Molecular dynamic (MD) simulation was performed for the most active BuChE inhibitor (6) to reveal the structural requirements responsible for the ligand binding. Complex 6-BuChE was compared with BuChE inactive inhibitor and topological isomer 5. inhibition of POP mediates the signaling cascade by releasing Ca 2+ at the endoplasmatic reticulum membrane, which may induce memory enhancement [29]. Some highly selective POP inhibitors entered clinical trial testing for neurodegenerative disorders, but these were terminated, mostly due to their shortcomings related to the compound's pharmacokinetic profile. Besides, we have also predicted the CNS availability of the selected compounds. Molecular dynamic (MD) simulation was performed for the most active BuChE inhibitor (6) to reveal the structural requirements responsible for the ligand binding. Complex 6-BuChE was compared with BuChE inactive inhibitor and topological isomer 5. Figure 1. Design of novel cholinesterase inhibitors derived from hit compound 1 and norbelladine-type Amaryllidaceae alkaloids isolated from Narcissus pseudonarcissus cv. Carlton, namely carltonine A and B.

Design
AAs are an extensively studied class of NPs. As mentioned previously, one such drug from this family, namely galantamine, was approved for AD treatment in 2001 [8]. Lycorine, haemanthamine, pancratistatine, and narciclasine are known to possess antiproliferative properties [30][31][32]. Montanine governs antidepressive, anticonvulsive, anti-rheumatic, antimicrobial, and antitumor activities [33]. The emerging data and potential applicability of AAs spur the search for new alkaloids from this family.
Indeed, within our previous study, we isolated new alkaloids from Narcissus pseudonarcissus cv. Carlton, namely carltonine A and B (Figure 1), demonstrating highly selective in vitro hBuChE inhibition potency in the nanomolar range [28] (Figure 1). Unfortunately, these alkaloids are present in plant material only in trace amounts, and they cannot be isolated for either more detailed biological investigation or commercial use. Thus, we decided to preserve some of the crucial structural requirements from carltonine A/B that are plausibly responsible for high BuChE inhibition activity, i.e., the 4-[2-(benzylamino)ethyl]phenol moiety, and modified other molecular regions to elucidate SAR. Specifically, we were interested in (i) the role of the secondary or tertiary amino group (the presence of allyl group), (ii) etherification of the phenolic hydroxyl group at aromatic ring B, and (iii) position of alkoxy or aryloxy substituents at benzene ring A (Figure 1), all concerning cholinesterase inhibitory activity. Some structural aspects related to the design of the novel family of compounds investigated herein can be found in the previously reported studies on disease-modifying multifunctional anti-AD agents. Indeed, one such family of compounds carries a benzylamine moiety, the critical fragment for cholinesterase binding [34]. It has also been observed that various aromatic substitutions in the benzylamine region, along with the introduction of a basic moiety into alicyclic amines, can yield cholinesterase inhibitors endowed with additional activity to beta-secretase 1 (BACE-1) [35,36]. Taking the precedents mentioned above, in this study, we have designed and synthesized 20 novel potential cholinesterase inhibitors derived either from commercially available Obenzylvanillin, or O-benzylisovanillin tethered to tyramine. Further details regarding the structural modifications are outlined in Figure 1.
In general, a portfolio of 20 compounds with a norbelladine framework was constructed. All prepared derivatives were characterized by LC-MS, HRMS, and NMR analysis. The purity was determined higher than 96% for each synthesized compound (Supplementary Materials).

In Vitro Cholinesterase Inhibitory Activities of New Norbelladine Derivatives (1-20)
All compounds (1)(2)(3)(4)(5)(6)(7)(8)(9)(10)(11)(12)(13)(14)(15)(16)(17)(18)(19)(20) were initially screened for hAChE/hBuChE inhibition potency at a concentration of 100 µM according to a slightly modified Ellman's method [28,38]. Compounds displaying inhibition ability >60% against one or both cholinesterases at the screening concentration were selected for the determination of their IC 50 values (Table 1). Galantamine and eserine were used as positive controls in the hAChE/hBuChE assay. The results are summarized in Table 1. All the compounds displayed only weak to moderate hAChE inhibition potency (IC 50 > 20 µM) and were selective towards hBuChE. Indeed, all the novel compounds revealed hBuChE inhibition potency in the micromolar to low nanomolar range (Table 1). Compound 5 was the only one displaying no inhibition potency against either enzyme.  (Table 1). On the contrary, except for compound 5, all derivatives showed intriguing inhibition potency against hBuChE at the screening concentration; thus, we were able to determine their IC 50 values to draw the structureactivity relationship. Compound 6, originating from O-benzylisovanillin and tyramine with two allyl substitutions, was the most pronounced hBuChE inhibitor (IC 50 = 72 ± 5 nM). Most strikingly, 6 emerged as the most selective hBuChE with a SI value of almost 1400. Surprisingly, its close topological derivative 5 was completely inactive in the hBuChE assay (IC 50 > 100 µM). In this specific case, the complete loss of inhibition ability is associated with the opposite positioning of the methoxy-and benzyloxy-groups on the aromatic ring (further inspected by in silico techniques). The etherification/amination by allyl of 1 grad-ually decreases the inhibition ability of compounds 3 and 5, with compound 5, endowed with two allyl appendages, being completely inactive. When template compound 2 was used, the opposite effect was observed, showing that the double-allyl derivative 6 was one order of magnitude more active than compound 4 (Table 1). Structure modifications in position C-4' of ring B by OH, H, and OMe in compounds 1, 7, and 11, respectively, were associated with gradual loss of hBuChE inhibition potency. This is the case when either the hydroxy group in 1 is removed, like in 7, or replaced by a methoxy group, like in 11 (IC 50 = 0.36 ± 0.03 µM, 1.28 ± 0.05 µM, and 2.39 ± 0.27 µM, for 1, 7, and 11, respectively). A similar effect was observed in the series derived from O-benzylisovanillin (compounds 2, 8, and 12 with IC 50 values of IC 50 = 0.29 ± 0.02 µM, 1.10 ± 0.05 µM, and 1.12 ± 0.11 µM, respectively). To elucidate further the necessity of benzyl attachment, a small series stemmed from 3-ethoxy-4-methoxybenzaldehyde (15-20; Scheme 3) was built. Within this subset, 20 (IC 50 = 0.69 ± 0.03 µM) was classified as the top-ranked hBuChE inhibitor, possessing N-benzyl substitution at the central nitrogen atom. This pattern of behavior might be ascribed to better occupancy of the ligand in the bulky gorge of hBuChE [39]. Recently, close derivatives with a norbelladine framework (N-benzyl-2-phenylethan-1-amine congeners) have been developed by Carmona-Viglianco et al. [40] and screened for their AChE/BuChE potency. The compounds reported herein include the same fundamental unit as compound 1, but they possess different substitution patterns within the A-ring. Generally, these derivatives yielded moderate to weak inhibition potency against both cholinesterases with IC 50 values above 10 µM. Most importantly, only one compound displayed selectivity to BuChE higher than 10.

POP Inhibition Activity of Selected Norbelladine Derivatives
Recently, some POP inhibitors have been evaluated in clinical trials to reverse and/or attenuate memory loss. Inhibition of POP can represent an additional supporting approach in AD treatment; thus, the POP inhibition potential of the most active norbelladine derivatives (1-4, 6 and 14) was tested. Unfortunately, the low solubility of tested compounds in buffer allowed the determination of IC 50 only for compound 1. This compound showed slightly lower POP inhibition potency (IC 50 = 186 ± 14 µM) compared to used standard POP inhibitor berberine (IC 50 = 142 ± 21 µM) [41].
a Tested at 100 µM compound concentration; b Compound concentration required to decrease enzyme activity by 50%; the values are the mean ± SEM of three independent measurements, each performed in triplicate; c Selectivity index for hBuChE is determined as ratio hAChE IC 50 /hBuChE IC 50 ; d calculated using BBB score [42]; e Reference compound; f Due to low solubility of compounds in buffer, the presented values correspond to the highest tested concentration; n.s. stands for not studied.

Prediction of CNS Availability
Significant brain exposure is a critical factor in the design of novel drugs to confront AD. Pursuant to this prerequisite, we have applied in-silico calculation of the so-called blood-brain barrier (BBB) score [42]. It is a well-validated predictive tool recently developed to distinguish between CNS and non-CNS drugs. It builds upon several physicochemical properties, including molecular weight, topological polar surface area, pKa, number of aromatic rings, heavy atoms, number of hydrogen bond donors, and acceptors. Most importantly, the BBB score exceeded the reliability of other predictive algorithms with a high value of area under the receiver operating characteristic curves (AUC), reaching a value of 0.86. All the compounds with BBB-score values higher than 4.0 are postulated to enter the CNS area. Indeed, all the novel compounds in this study (1)(2)(3)(4)(5)(6)(7)(8)(9)(10)(11)(12)(13)(14)(15)(16)(17)(18)(19)(20) reached a BBB-score ranging between 4.53 to 5.60 (Table 1), showing their ability to permeate the BBB. It is noteworthy that the template drugs carltonine A (BBB score 4.52) and carltonine B (4.54) displayed comparable values to compound 1, whereas the rest of the family revealed higher BBB score values.

Enzyme Kinetic Analysis of Compound 6
Compound 6, the derivative with the most robust inhibitory activity, was subjected to enzyme kinetics analysis to determine the mode of hBuChE inhibition. Inhibition kinetics were elucidated from velocity curves that were measured at several concentrations of substrate butyrylthiocholine and test compound. The type of enzyme inhibition and corresponding kinetic parameters (K i and K i ) were determined by nonlinear regression analysis. Results for each type of model of inhibition (competitive, non-competitive, uncompetitive, and mixed) were compared by the sum-of-squares F-test. Statistical analysis showed a mixed type of inhibition (p < 0.05), which is in line with the Lineweaver-Burk plot, used for visualization of the obtained data ( Figure 2). The intersection of lines is located above the x-axis, which means a reversible binding mode to both free enzyme and enzyme-substrate complex, with a stronger affinity to the free enzyme (K i < K i ). K m was increased whereas V max was slightly reduced at a higher concentration of 6. A K i value of 60.3 ± 11.6 nM and K i of 347 ± 79 nM were measured against hBuChE.

Molecular Dynamic Simulation for Compounds 5 and 6
Molecular dynamic simulation was used to determine the binding differences between the topological isomers 5 (hBuChE inactive at 100 µM) and 6 (hBuChE IC 50 = 72 nM) in the hBuChE active site (PDB ID: 6QAA) [43]. As portrayed in Figure 3A,B, compound 5 is accommodated in a "wrapped" manner. The critical interaction can be classified as parallel π-π stacking between the A phenyl ring and Tyr332 (3.5 Å), and the salt bridge formed between the carboxyl of Asp70 and the protonated tertiary amino group (4.8 Å). Two hydrogen bonds with water molecules are also apparent. One can be observed with oxygen from the ligand's methoxy group (2.4 Å); the second is mediated to the protonated tertiary amino group (1.9 Å). The benzyloxy group is implicated in the formation of T-shaped π-π interactions with Trp82 (4.7 Å) and His438 (4.8 Å). The latter residue is part of the enzyme's catalytic machinery. The other two catalytic triad residues, namely Ser198 and Glu325, stand aside from the ligand anchoring. The allyl group seems to protrude outside the cavity gorge providing no specific interaction with the enzyme at all. Ring B of ligand 5 occupies the oxyanion hole of the enzyme flanked by Gly116, Gly117, and Ala199. The allyloxy appendage contacts the acyl binding pocket (Val288, Leu286) of the enzyme via hydrophobic interaction, and Trp231 by aliphatic-π interaction. For the 6-hBuChE complex ( Figure 3C,D), the ligand is accommodated in a more "straight" manner compared to the lodging of 5 with the enzyme cavity. The protonated amino center of 6 is anchored by cation-π to Tyr332 (5.6 Å) and one water molecule (1.9 Å). The benzyloxy moiety orchestrates several hydrophobic interactions, e.g., π-π stackings to Trp231 (T-shaped; 3.6 Å), Phe329 (distorted; 4.2 Å), and Phe398 (distorted; 3.7 Å). Besides, the benzyloxy moiety is also surrounded by residues from the acyl pocket (Leu286, Val288) and oxyanion hole residues (Gly116, Gly117, Ala199). Moreover, π-π contact to His438 (4.4 Å) and van der Waals forces with Ser198 from the catalytic triad can also be observed. Phenyl ring A faces Trp82 (4.1 Å) via π-π interaction. The methoxy group attached to phenyl ring A is implicated in a hydrogen bridge with one water molecule. Phenyl ring B is oriented distally being exposed to Ala277, Ile69, and Asp70 residues.
From the MD simulation, it can be concluded that the higher inhibition ability of 6 can be ascribed to the accommodation of its benzyloxy substituent, revealing several crucial interactions with the enzyme. The inability of 5 to lodge in the hBuChE active site, like 6, led to less energetically favored overall topology, i.e., the "wrapped" conformation. For the sake of clarity, we are providing the superimposed spatial view of 5 and 6 ( Figure 4), demonstrating their different orientation in the hBuChE cavity.

Cytotoxicity of Compound 6
The cytotoxicity profile was assessed for compound 6 using the human neuroblastoma SH-SY5Y cell line and MTT assay. At a concentration of 10 µM, compound 6 was completely devoid of cytotoxicity, with viability reaching 95.1 ± 2.6%.

MAO-A and MAO-B Inhibition Potency of Compound 6
Since monoamine oxidase (MAO; EC 1.4.3.4) inhibitors are considered as neuroprotective agents, compound 6 was screened at two concentrations (1 and 10 µM) for its inhibition potential against both isoforms of human MAOs. Compound 6 displayed only weak MAO-A and MAO-B inhibition activity (residual activity of enzymes after treatment by 6 was higher than 68% in both concentrations; IC 50 >10 µM).

Conclusions
The development of highly selective BuChE inhibitors represents an emerging field of research [44]. There is strong evidence that inhibiting BuChE could benefit AD patients in more progressive stages of the disease compared to AChE-selective inhibitors [45]. Moreover, selective BuChE inhibitors do not exert any adverse effects mediated by the enzyme inhibition. Indeed, it has been reported that selective inhibition of BuChE with cymserine analogs raises ACh levels in the brain and improves cognitive performance without any adverse parasympathetic side effects [45,46]. To date, only a few studies have reported the discovery of selective, reversible BuChE inhibitors [43,44,47,48]. Within this work, we have constructed 20 novel compounds based on the norbelladine framework using carltonine A and B as template scaffolds. The compounds were generated by applying straightforward synthetic routes with moderate to excellent yields. The majority of the developed compounds displayed a highly selective inhibition pattern for hBuChE with IC 50 values below 1 µM. The most selective of them, compound 6, exerted nanomolar scale inhibition with an IC 50 of 72 nM. The mode of hBuChE inhibition was minutely inspected using enzyme kinetic analysis in tandem with MD simulation. The latter pointed out the differences in the binding modes between hBuChE inactive inhibitor 5 and its isomer, topranked compound 6. We have also predicted the CNS availability by applying BBB score calculation, presuming that compounds can pass through the BBB. The observations made herein should pave the way for the structure-based optimization of further norbelladine analogues, potentially applicable in the treatment of neurodegenerative diseases.

General Experimental Procedures
All solvents were treated by using standard techniques before use. All reagents and catalysts were purchased from commercial sources (Sigma Aldrich, Prague, Czech Republic) and used without purification. The NMR spectra were obtained in CDCl 3 at ambient temperature on a VNMR S500 (Varian) spectrometer operating at 500 MHz for 1 H and 125.7 MHz for 13 C and on a JNM-ECZ600R (Jeol) instrument operating at 600 MHz for 1 H and 151 MHz for 13 C. Chemical shifts were recorded as δ values in parts per million (ppm) and were indirectly referenced to tetramethylsilane (TMS) via the solvent signal (CDCl 3 -7.26 ppm for 1 H and 77.0 ppm for 13 C). Coupling constants (J) are given in Hz. ESI-HRMS were obtained with a Waters Synapt G2-Si hybrid mass analyzer of a quadrupole-time-of-flight (Q-TOF) type, coupled to a Waters Acquity I-Class UHPLC system. TLC was carried out on Merck precoated silica gel 60 F254 plates. Compounds on the plate were observed under UV light (254 and 366 nm) and visualized by spraying with Dragendorff's reagent. 1, 2,  7, 8, 11, 12, 15, and 18 One hundred mg of the required aromatic aldehyde (0.41 mmol of O-benzylvanillin/Obenzylisovanillin; 0.55 mmol of 3-ethoxy-4-methoxybenzaldehyde) was dissolved in 5 mL of methanol, the corresponding amine (0.28-0.38 mmol) was added, and the reaction was stirred for one day. Sodium borohydride (0.28-0.38 mmol) was added to the reaction mixture at 0 • C and the mixture was stirred for a further 3 h at room temperature. The reaction mixture was evaporated to dryness and separated by preparative thin-layer chromatography (To:cHx:Et 2 NH 5:10:2). The NMR spectra for known compounds 1 and 2 [37,49] matched those in the literature. The characterization of all new compounds is given below. (1) spectra in accordance with literature [49]. Yield: 143 mg (95%), white amorphous solid.

hAChE and hBuChE Inhibition Assay
The hAChE/hBuChE activities were determined using a modified method of Ellman [38], as recently described by our group [50]. A detailed description of the assay can be found in Supplementary Materials.

Kinetic Study of hBuChE Inhibition
Kinetic studies employed the same procedure that was previously described [50]. Supplementary Materials provides details of this assay.

POP Inhibition Assay
The previously applied procedure was used [50]. Detailed description of the assay can be found in Supplementary Materials.

MAOs Inhibition Assay
MAO inhibition assay was performed as previously described [51]. Supplementary Materials provides details of this assay.
The MTT (3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyl-tetraziolium bromide (Sigma-Aldrich) reduction assay was used for measurement of the compounds' cytotoxicity, according to Mosmann [51]. MTT is a water-soluble tetrazolium salt that is converted to purple formazan by succinate dehydrogenase in mitochondria of viable cells [52,53]. Cell viability was detected after 24-h incubation with the tested substances. For the assay, SH-SY5Y cells were seeded into 96-well plates in 100 µL volume and at a density of 20 × 10 3 . Cells were allowed to attach overnight before the treatment. The stock solutions of tested compounds were prepared in dimethylsulfoxide (DMSO, Sigma-Aldrich), which were serially diluted in DMEM and added to the cells in a 96-well culture plate. The final concentration of DMSO was less than 0.25% per well. After 24 h incubation, the medium containing serially diluted substances was aspirated from each well and replaced by 100 µL of fresh medium containing MTT (0.5 mg/mL). Plates were subsequently incubated at 37 • C in a CO 2 incubator for 1 h. Medium containing MTT was then aspirated, and formazan dissolved in 100 µL of DMSO. The optical density of each well was measured using a Spark ® multimode microplate reader (Tecan Group Ltd., Männedorf, Switzerland) at 570 nm. The cell viability was expressed as a percentage of the untreated control. Each experiment was performed in triplicate and repeated three independent times.
The docking poses of 5 and 6 were improved by MD simulation. The receptor structure was prepared using the software Chimera. The best-scored docking pose was taken as the initial for MD. The force-field parameters for ligands were assessed by Antechamber [59], v. 20.0 using General Amber force-field 2 [60]. MD simulations were carried out by Gromacs, v. 2018.1 [61]. The complex receptor-ligand was solvated in the periodic water box using the TIP3P model [62]. The system was neutralized by adding Na + and Cl − ions to the concentration of 10 nM. The system energy was minimalized and equilibrated in a 100-ps isothermal-isochoric NVT and then a 100-ps isothermal-isobaric NPT phase. Then, a 10-ns MD simulation was run at the temperature of 300 K. The molecular docking and MD results were 3D visualized by the PyMOL Molecular Graphics System, Version 2.4.1, Schrödinger, LLC.