Antioxidant, LC-MS Analysis, and Cholinesterase Inhibitory Potentials of Phoenix dactylifera Cultivar Khudari: An In Vitro Enzyme Kinetics and In Silico Study

We evaluated the therapeutic potentials of Khudari fruit pulp, a functional food and cultivar of Phoenix dactylifera, against neurological disorders. Our results demonstrate a good amount of phytochemicals (total phenolic content: 17.77 ± 8.21 µg GA/mg extract) with a high antioxidant potential of aqueous extract (DPPH assay IC50 = 235.84 ± 11.65 µg/mL) and FRAP value: 331.81 ± 4.56 µmol. Furthermore, the aqueous extract showed the marked inhibition of cell-free acetylcholinesterase (electric eel) with an IC50 value of 48.25 ± 2.04 µg/mL, and an enzyme inhibition kinetics study revealed that it exhibits mixed inhibition. Thereafter, we listed the 18 best-matched phytochemical compounds present in aqueous extract through LC/MS analysis. The computational study revealed that five out of eighteen predicted compounds can cross the BBB and exert considerable aqueous solubility. where 2-{5-[(1E)-3-methylbuta-1,3-dien-1-yl]-1H-indol-3-yl}ethanol (MDIE) indicates an acceptable LD50. value. A molecular docking study exhibited that the compounds occupied the key residues of acetylcholinesterase with ΔG range between −6.91 and −9.49 kcal/mol, where MDIE has ∆G: −8.67 kcal/mol, which was better than that of tacrine, ∆G: −8.25 kcal/mol. Molecular dynamics analyses of 100 ns supported the stability of the protein–ligand complexes analyzed through RMSD, RMSF, Rg, and SASA parameters. TRP_84 and GLY_442 are the most critical hydrophobic contacts for the complex, although GLU_199 is important for H-bonds. Prime/MM-GBSA showed that the protein–ligand complex formed a stable confirmation. These findings suggest that the aqueous extract of Khudari fruit pulp has significant antioxidant and acetylcholinesterase inhibition potentials, and its compound, MDIE, forms stably with confirmation with the target protein, though this fruit of Khudari dates can be a better functional food for the treatment of Alzheimer’s disease. Further investigations are needed to fully understand the therapeutic role of this plant-based compound via in vivo study.


Introduction
Alzheimer's disease (AD) stands as one of the top three contributors of neurological disorders mainly associated with old age, with around 130,000 cases reported only in Saudi Arabia [1].The symptoms of AD patients generally include memory loss, agitation, apathy, cognitive impairments, and aberrant motor behavior [2].The pathophysiology of AD is quite complex, and studies have reported that amyloid-beta (AB) deposition forms plaques, and the hyperphosphorylation of tau proteins forms neurofibrillary tangles during AD, which further disturbs redox balance and the cholinergic activity of neurons [3][4][5].Normally, cholinergic activity greatly depends on the release of acetylcholine (ACh), a neurotransmitter that further binds to ACh receptors and transmits the signals, or it can be hydrolyzed via acetylcholinesterase (AChE) [6].It was reported that, during AD, due to neuronal damage or the degradation of cholinergic neuron-rich region, there is a sudden decline in ACh levels, and cholinesterase enzymes further decrease its levels, which hinders the neuronal signal transmission [7].The elevated level of acetylcholine in synapses stimulates nicotinic and muscarinic receptors, which give beneficial relief for cognitive problems in AD [8,9].Therefore, cholinesterase inhibition can manage AD by reducing the catalysis of Ach, which ultimately supports the signal's transmission and improves the life quality of AD patients.
However, the most common FDA-approved medication for AD, i.e., cholinesterase inhibitors (tacrine, donepezil), exhibits several side effects, such as nausea, headache, decrease appetite, and vomiting [10,11].As a result, there is a great deal of interest in developing more effective treatments to alleviate illness symptoms and stop the loss of neurons, which causes the disease.Natural ingredients and their byproducts are known to have higher therapeutic potentials for a variety of metabolic and infectious problems, including diabetes, oxidative damage, high cholesterol levels, ulcer, neurological disorders, tumors, and infections caused by microbes [12][13][14][15][16][17][18][19][20][21][22].Because nature provides a plentiful supply of bioactive metabolites, it may be safe and cost-effective to develop new inhibitors of important regulatory enzymes for the treatment of metabolic diseases, including neurological illnesses [15,23].
The date palm (Phoenix dactylifera L.) represents the most grown fruit-bearing cultivated crop in the world's arid and semiarid areas.Throughout the millennia, Islamic physicians have carried out plenty of research related to Phoenix dactylifera in the treatment and management of various disease, as they believed that the Prophet Mohammed (Peace Be Upon Him) said: "There is no disease that Allah has created, except that Allah also has created its cure" [24].Phoenix dactylifera leaves, fruit, pits, and pollen have been utilized in the management of various health problems and disorders because they are a significant source of bioactive chemicals that are primarily responsible for their therapeutic activity [25].Fruits of Phoenix dactylifera are high in carbs, amino acids, vitamins, and minerals.Alkaloids, tannins, glycosides, flavonoids, steroids, terpenoids, and carotenoids are active phytoconstituents in fruits.The fruit's principal active ingredients, including flavonoids, exhibit various medicinal properties.The plant extract also contains antioxidants and has neuroprotective properties, and the fruits have the potential to help control Alzheimer's disease [26].Khudari, a cultivar of Phoenix dactylifera that is mostly cultivated in the Al-Qassim region, possesses a higher protein content (3.42 gm/100 gm dry weight), higher energy values (287.5 kcal/100 gm), and a lower amount of lipids (0.18 gm/100 gm DW) than most Phoenix dactylifera varieties [27,28].The total phenol content (TPC) and radical scavenging (DPPH) potentials of extracts of Khudari fruits were found to be nearly identical to Awja date varieties (the most researched and therapeutic variety of dates).This date (Khudari) was also shown to have a significantly reduced level of sugar (74.5 grammes/100 gm DW) [29,30].
Various natural bioactive compounds have been identified in fruit extract of Phoenix dactylifera, i.e., luteolin, protocatechuic acid, quercetin, gallic acid, resorcinol, chelidonic acid, and chlorogenic acid [31].The antioxidant potential of Phoenix dactylifera is possible due to the occurrence of phenolics, melatonin, carotenoids, and vitamins contents, e.g., p-coumaric acid [32].Hepatoprotective activity and protection against elevated alkaline phosphatase was also observed, and this is due to the presence of bioactive compounds i.e., ferulic acid, caffeic acid, quercetin, anthocyanins, and proanthocyanidins [33,34].There have been plenty of in silico, in vitro, and in vivo reports published, which exhibit the potential effect of fruit extract and their bioactive compounds in the treatment and management of AD via modulating various risk factors, i.e., antioxidant and hepatoprotective activity, memory and behavior improvement, and cerebral anti-ischemic, neuroprotective, anxiolytic, nootropic, and antipsychotic activities [35][36][37].There is no study available that demonstrates the inhibitory activity of AChE from Phoenix dactylifera fruits-derived bioactive compounds.Therefore, we selected Phoenix dactylifera fruit pulp extract and its bioactive compounds as a potential AChE inhibitor for investigation.
It is generally known that computational biology and in vitro inhibition of enzymes are the primary techniques for elucidating binding patterns and optimizing bioactive chemicals for target-based drug development [38][39][40][41][42].However, despite having beneficial phytoconstituents and antioxidant activity, this Khudari variety of dates is still the least explored for its therapeutic potentials against several pathologies.Considering these facts, we planned to explore the therapeutic potential of the fruit pulp of Khudari cultivars of Phoenix dactylifera and their secondary metabolites for the treatment of neurological disorders through cholinesterase inhibition via a variety of methods, such as cell-free enzyme inhibition kinetics, molecular docking, and MDS (molecular dynamics simulation) studies.

Collection and Preparation of Dates Extract
Khudari cultivars of Phoenix dactylifera dates were procured from a local market of the Qassim region in Saudi Arabia.The fruit pulp of dates (500 g) was separated, dried in a shed, and pulverized into powder form.The powdered form of the date weighed around 25 gm and was extracted using the optimum amount (1:10 ratio) of methanol and water solvent in the Soxhlet apparatus (Borosil), respectively, until the solvent became colorless.Once the solvent attained room temperature, it was filtered and air-dried, and the residue was scraped out and subsequently stored at −20 • C for further experimental use.
The method for calculating percent yield of different fractions is as follows: %yield = Dry weight of plant extract Dry weight of plant material × 100

Total Phenolic Content (TPC) and Phytochemical Analysis
Khudari date extracts were tested qualitatively for phytochemical composition using a standard approach [43].To measure the total phenolic content of the extracts, the Folin-Ciocalteu method was used [44].

Scavenging Activity of DPPH Radical
The methanolic and aqueous fraction of Khudari extracts were scrutinized for their potential to scavenge the DPPH radical using Brand-William's method.Briefly, 3 mL of DPPH solution was added to 100 µL of plant extract and the reference drug (ascorbic acid) and incubated at 37 • C for 30 min; thereafter, absorbance was measured at 517 nm with a Biospectrum kinetics spectrophotometer (Eppendorf) [45].To measure the percentage (%) of DPPH free radicals that were scavenged, the following equation was applied: %DPPH = ∆Absorbance of control − ∆absorbance of test sample ∆Absorbance of control × 100

Analysis of Ferric-Reducing Antioxidant Power
An adaptation of Benzie and Strain's method was used for assessment of the ferricreducing antioxidant potential (FRAP) of Khudari extracts [46].Concisely, A freshly prepared FRAP reagent can be attained by combining sodium acetate buffer (300 mM, pH 3.6), 10 mM TPTZ solution, and 20 mM FeCl 3 solution in a volume ratio of 10:1:1.To 3 mL of FRAP reagent, 100 microliters of sample were added at varying concentrations.The absorbance was determined at 593 nm after incubating it for 30 min at room temperature.Results were expressed as µmol Fe (II)/g dry weight of date material, based on plots of the standard curve using the FeSO 4 solution.

Acetylcholinesterase Inhibition Assay
We performed the acetylcholinesterase test according to Ellman's protocol with slight modifications [47].In a 2 mL cuvette, 33 µL of 10 mM-DTNB, 100 µL of 1 mM-AChI, and 767 µL of 50 mM-Tris HCl (calibrated at pH = 8.0), along with 100 µL of extracts at different concentrations, were added.The above-mentioned cuvette was blank, whereas in another cuvette, 300 µL of buffer was substituted with an equal volume of 0.28 U ml −1 of acetylcholinesterase enzyme.To carry out the competitive analysis, standard drug tacrine was used.The following reaction was monitored continuously in a Biospectrum kinetics spectrophotometer (Eppendorf) for 20 min at 405 nm wavelength.The experiment was performed in triplicate.The following equation was used to estimate the percentage (%) of acetylcholinesterase enzyme activity: %inhibition = ∆Absorbance of control − ∆Absorbance of sample ∆Absorbance of control × 100

Spectrometric Study of the Enzyme Inhibition Kinetics
Kinetic analysis was carried out at 30 • C using four distinct concentrations (0.5, 1.0, 1.5, and 2.0 mM) of substrate (AChI) to determine the AChE inhibition mechanism via absence/presence of various inhibitory concentrations of standard (tacrine) and aqueous fractions of Khudari extract."AChE catalyzed-hydrolysis" of AChI was analyzed kinetically using spectrophotometry (Eppendorf) at a wavelength of 405 nm for a total of 20 min, and the OD values were measured at 60 s intervals.Additionally, a Lineweaver-Burk, Dixon, and Km and V max plot was devised to define the inhibitory mechanism [48].

LC-MS (Liquid Chromatography-Mass Spectrometry) Analysis
The fraction that exhibited the greatest inhibition of AChE activity was further explored for bioactive compounds through LC-MS analysis, performed at the Sophisticated Analytical Instrument Facility, CSIR-Central Drug Research Institute, Lucknow.The LC-MS and MS/MS simulations were carried out using a Waters Alliance E2695/HPLC-TQD mass spectrometer with a mass range of 100-2000 Da.Separation was achieved on an Agilent Poroshell 120 EC C-18 column (150 mm × 4.6 mm, 2.7 µm).Two mobile phases were used: A-0.1% formic acid in water and B-90% of acetonitrile in water, at a flow rate of 1.500 mL/min.The gradient started with A-5% and B-95%, a linear increase in A% and decrease in B% after 6 min at A-30% and B-70%; during 12-16 min, A-60% B-40%; 20 min, A-80% and B-20%; and 26-30, min A-5% and B-95%.The chromatograms were compared to the available literature to identify the peaks, and then, the output file was subjected to analyze the peaks via Mnova software version 14.2.1 offered by Mestrelab Research, chemistry software solutions.

Computational Study 2.5.1. Platform of In Silico Computational Study
The computer-added drug design studied was performed on a Lenovo Notebook equipped with Intel(R) Core (TM) i5-8265U @1.80 GHz powered machine with 2 GB graphics card (GeForce MX250 NVIDIA) and 8 GB of RAM.The software for molecular modelling studies, such as AutoDock 4.2 32-bit (http://autodock.scripps.edu/(accessed on 16 April 2023)), Cygwin, SPDBV, Discovery studio (DS) visualizer, ChemDraw, Notepad++ version 7.9.2., and Open Babel GUI, was installed.The Desmond (Schrodinger-2020, LLC, New York, NY, USA) tool was used for the molecular dynamics simulation study, carried out on workstation equipped with NVIDIA Quadro P5000 graphic card, 28 GB RAM, and 3.50 GHz processor.

Retrieval and Preparation of 3D Structure of Ligands
The LC-MS-predicted compound's 3D structures were retrieved through the PubChem (http://pubchem.ncbi.nlm.nih.gov(accessed on 16 May 2023)) database as an SDF file, which was later converted to a PDB file format using Discovery studio visualizer, and the energy was minimized by applying CHARMM force field using the steepest descent method for 500 steps and an RMS gradient of 0.01.Thereafter, the ligands were converted to a ".pdbqt" file format [42,49].

Retrieval and Preparation of Target Protein
The target protein 3D structure of AChE (ID: 1ACJ) was downloaded from the protein data bank (PDB) (https://www.rcsb.org/search(accessed on 16 May 2023)) and visualized by Discovery Studio Visualizer 2020 (https://discover.3ds.com/discovery-studiovisualizer-download(accessed on 3 October 2020)) [50].The resolution of the selected AChE crystal structure was 2.80 Å.The active site of the protein was obtained through an online tool, the Play-Molecule (https://www.playmolecule.com/(accessed on 16 May 2023)) DEEPSITE [51].The structure of the target protein was cleaned by removal of heteroatoms and polar hydrogen and addition of Kollmann charges using AutoDock tools [52].After that, the final structure for docking was prepared by converting it to a PDBQT file format through AutoDock 4.2 tool [42,[52][53][54].

Molecular Docking Studies of Selected Ligands against AChE
Molecular docking was performed by using the AutoDock 4.2 version as described earlier [55,56].Briefly, Hetatms were removed, and Kollman charges were assigned to protein (1ACJ), and the processed protein was saved in a ".pdbqt" file for docking study.The grid box size was maintained at 60 × 60 × 60 points (x, y, and z, respectively) with 0.608 Å grid spacing, and the grid center was selected at dimensions x, y, and z: −7.36, 64.04, and 134.87, respectively.A total of forty independent docking runs were executed for each ligand.The conformational clustering of the poses of the target-ligand complexes was determined by considering a root mean square deviation tolerance <2.0 Å.The complexes (target-ligand) with the least free energy of binding (∆G) were considered the most stable and indicated the inhibitory effects of ligands [53,57].Further, the protein-ligand complex interaction was studied using Discovery studio visualizer version 2020.

Pharmacokinetics and Drug-Likeness Profiling of Ligands
All the predicted compounds (ligands) were evaluated for their pharmacokinetics and physicochemical parameters via SWISS ADME tool [58].They were also evaluated for their drug-like characteristics based on their agreement with Lipinski's rule of five (http://www.swissadme.ch(accessed on 18 May 2023)) [59].

Predicted Toxicity Assessment of the Selected Ligands
The ProTox-II web-server (https://tox-new.charite.de/protox_II/index.php?site= compound_input (accessed on 18 May 2023)), an online tool that predicts the toxicity of small organic molecules, was used for the toxicity assessment.It provides extensive facts concerning the toxicity of the chemicals, such as LD 50 , carcinogenicity, immunotoxicity, mutagenicity, cytotoxicity, and hepatotoxicity [60].

Molecular Dynamics Simulation Analysis
The complex of the best-scoring ligand (best hit) and AChE protein (PDB id: 1ACJ) was selected for the molecular dynamics simulation study using Schrodinger-2020, LLC, New York, NY, USA tools [42].The complex was imported to Maestro interface for pre-processing (optimization and minimization).
The protein-ligand complex developed with AutoDock Vina was imported into the Maestro interface in Schrodinger's software 2020.The optimization and minimization of the protein-ligand complex was performed using Maestro's Protein Processing Wizard.All the systems were built using the System Builder tool.As a solvent model with an orthorhombic box (10 Å), TIP3P was selected.The simulation used the OPLS_2005 force field [61].The models were neutralized by using counterions.The addition of sodium chloride (NaCl) 0.15 M simulated physiological conditions.The normal pressure test (NPT) ensemble with a pressure of one atmosphere and 300 K temperature was employed for the whole simulation.The trajectories were saved for analysis every 100 ps.Measures including the root mean square deviation (RMSD), the percentage of secondary structural components, root mean square fluctuation (RMSF), and protein-ligand interactions were analyzed to assess the stability of the protein-ligand complexes.The findings of the three independent research studies are shown as mean ± standard deviation.

Calculations of Free Energy (Prime-MM/GBSA)
Using the Prime module (Schrodinger, LLC, New York, NY, USA) and the MM-GBSA approach, the interaction's free energy value for the protein-ligand complex was calculated as previously described [42].Once equilibration was achieved, free energy was estimated on every 10 ns MD simulation trajectory up to 100 ns.In conclusion, the docked complexes were first optimized locally in Prime using molecular mechanics (MM), and then, their energies were reduced using the OPLS-AA (2005) force field and the generalized Born surface area (GBSA) continuum solvent model.

Phytochemical Estimation of Khudari Fruit Pulp
Our results showed the qualitative presence/absence of different phytochemicals, such as alkaloids, carbohydrate, flavonoids, glycosides, tannins, terpenoids, phenols, cardiac glycosides, and steroids, in both the MeOH and aqueous fraction of Khudari extract (Table 1).The percentage yield was found to be MeOH 48%, and aqueous 36%.Phytochemicals are easily accessible, bioavailable, and far less toxic than synthetic medications and show significant promise for preventing and treating AD.They also have several synergistic benefits, such as improving cognitive and cholinergic functioning.A great amount of scientific research has demonstrated the anti-Alzheimer's disease benefits of polyphenols, including cell damage protection, neuroprotective effects, AChE inhibitory action, and amyloid aggregation [62,63].

Scavenging Property of DPPH Radical
The comparatively stable DPPH radical is commonly employed to assess the free radical scavenging ability of natural antioxidants such as fruit and plant extracts.The percent scavenging capacity of DPPH radical by the MeOH and aqueous fraction of Khudari extract is shown in Figure 1.The strongest antioxidant activity was found in the aqueous fraction (IC 50 value: 235.84 ± 11.24 µg/mL) in a dose-dependent style, where the ascorbic acid used as reference standard showed an IC 50 value of 15.58 ± 0.76 µg/mL.Our results are in agreement with previous reports of free radical scavenging potentials of Khudari fruit extract [28,29].

Scavenging Property of DPPH Radical
The comparatively stable DPPH radical is commonly employed to assess the free radical scavenging ability of natural antioxidants such as fruit and plant extracts.The percent scavenging capacity of DPPH radical by the MeOH and aqueous fraction of Khudari extract is shown in Figure 1.The strongest antioxidant activity was found in the aqueous fraction (IC50 value: 235.84 ± 11.24 µg/mL) in a dose-dependent style, where the ascorbic acid used as reference standard showed an IC50 value of 15.58 ± 0.76 µg/mL.Our results are in agreement with previous reports of free radical scavenging potentials of Khudari fruit extract [28,29].

Ferric-Ions-Reducing Potential and Total Phenolic Content
The FRAP (ferric-reducing antioxidant property) assay was performed to evaluate the total antioxidant capacity of the Khudari extract fractions.The assay is based on the extract's ability to reduce ferric ions to ferrous ones.Our findings elucidated that the aqueous fraction has significantly higher values of FRAP than the MeOH, which are 331.81± 4.56 µmol and 51.57± 1.183 µmol Fe (II)/g, respectively.The total content of phenols of both extracts was measured by the standard-procedure Folin-Ciocalteu method.The result shows that the highest content of phenols was found in the aqueous (17.77 ± 8.21 µg GA/mg extract) rather than the MeOH (7.70 ± 0.45 µg GA/mg extract) fraction.A linear relationship test was also used to examine the relationship between TPC and FRAP.Our findings indicated that as the phenolic content rose, the ferric-reducing potential also increased.Based on our data, we can conclude that TPC and FRAP have a positive association (Figure 2).

Ferric-Ions-Reducing Potential and Total Phenolic Content
The FRAP (ferric-reducing antioxidant property) assay was performed to evaluate the total antioxidant capacity of the Khudari extract fractions.The assay is based on the extract's ability to reduce ferric ions to ferrous ones.Our findings elucidated that the aqueous fraction has significantly higher values of FRAP than the MeOH, which are 331.81± 4.56 µmol and 51.57± 1.183 µmol Fe (II)/g, respectively.The total content of phenols of both extracts was measured by the standard-procedure Folin-Ciocalteu method.The result shows that the highest content of phenols was found in the aqueous (17.77 ± 8.21 µg GA/mg extract) rather than the MeOH (7.70 ± 0.45 µg GA/mg extract) fraction.A linear relationship test was also used to examine the relationship between TPC and FRAP.Our findings indicated that as the phenolic content rose, the ferric-reducing potential also increased.Based on our data, we can conclude that TPC and FRAP have a positive association (Figure 2).

AChE Inhibition Property and Mode of Inhibition
Based on the antioxidant property of the aqueous fraction, we further examined its AChE inhibitory activity.The aqueous fraction of Khudari extract showed the marked inhibition of AChE in a dose-dependent manner, with an IC50 value of 48.25 ± 2.04 µg/mL.We observed 3.59% inhibition at 12.5 µg/mL, and the percentage inhibition increased to 34.72% at 25 µg/mL and 51.28% at 50 µg/mL, whereas percentage inhibition for the highest concentration used, which was 100 µg/mL, was almost 66.26% inhibition of AChE activity.The reference standard showed an IC50 value of 17.31 nM.We further determined the mode of inhibition (enzyme inhibition kinetic) of the aqueous fraction via the AChIhydrolysis inhibition at different concentrations (Figure 3A-D).The standard drug tacrine binds to the allosteric site, a site that differs from the active site of the enzyme, resulting in decreased enzyme affinity without competing with its substrate binding site.The fundamental method for distinguishing non-competitive inhibition from competitive inhibition is by the decrease in Vmax and the unchanged Km (Figure 3A, B).On the other hand, enzyme kinetic studies revealed that the aqueous fraction of Khudari extracts at different concentrations do not intersect on any axis, nor is their slope same; thus, they exhibit a mixed type of inhibition that might be a ributed to the presence of several secondary metabolites in the plant extract (Figure 3C).Considering the decline in both the values of Km and Vmax by increasing the dose concentration, it is evident that the plant extract exhibits non-competitive inhibition.The plant extract is composed of several secondary metabolites, and therefore, some of them are involved in a non-competitive mode of inhibition and some in a mixed type of inhibition (competitive and noncompetitive) (Figure 3D).Therefore, from the data, we can clearly state that there was a mixed type of AChE inhibition demonstrated by the aqueous fraction.Our findings are consistent with earlier reports in which tacrine was described as a reversible noncompetitive inhibitor towards the AChE enzyme [64,65] as well as a mixed-type inhibitor with a high non-competitive component [66].

AChE Inhibition Property and Mode of Inhibition
Based on the antioxidant property of the aqueous fraction, we further examined its AChE inhibitory activity.The aqueous fraction of Khudari extract showed the marked inhibition of AChE in a dose-dependent manner, with an IC 50 value of 48.25 ± 2.04 µg/mL.We observed 3.59% inhibition at 12.5 µg/mL, and the percentage inhibition increased to 34.72% at 25 µg/mL and 51.28% at 50 µg/mL, whereas percentage inhibition for the highest concentration used, which was 100 µg/mL, was almost 66.26% inhibition of AChE activity.The reference standard showed an IC 50 value of 17.31 nM.We further determined the mode of inhibition (enzyme inhibition kinetic) of the aqueous fraction via the AChI-hydrolysis inhibition at different concentrations (Figure 3A-D).The standard drug tacrine binds to the allosteric site, a site that differs from the active site of the enzyme, resulting in decreased enzyme affinity without competing with its substrate binding site.The fundamental method for distinguishing non-competitive inhibition from competitive inhibition is by the decrease in Vmax and the unchanged Km (Figure 3A, B).On the other hand, enzyme kinetic studies revealed that the aqueous fraction of Khudari extracts at different concentrations do not intersect on any axis, nor is their slope same; thus, they exhibit a mixed type of inhibition that might be attributed to the presence of several secondary metabolites in the plant extract (Figure 3C).Considering the decline in both the values of K m and V max by increasing the dose concentration, it is evident that the plant extract exhibits non-competitive inhibition.The plant extract is composed of several secondary metabolites, and therefore, some of them are involved in a non-competitive mode of inhibition and some in a mixed type of inhibition (competitive and non-competitive) (Figure 3D).Therefore, from the data, we can clearly state that there was a mixed type of AChE inhibition demonstrated by the aqueous fraction.Our findings are consistent with earlier reports in which tacrine was described as a reversible non-competitive inhibitor towards the AChE enzyme [64,65] as well as a mixed-type inhibitor with a high non-competitive component [66].

LC/MS Analysis for Identification of Phytochemical Compounds in Aqueous Fraction of Phoenix dactylifera (Khudari Cultivar)
In the present work, bioactive composition from the aqueous fraction of fruit pulp of Phoenix dactylifera (Khudari cultivar) was carried out on a Waters Alliance E2695/HPLC-TQD mass spectrometer in positive-and negative-ionization mode.The output file was analyzed in the Mnova software offered by Mestrelab Research, chemistry software solutions, and the peaks were matched by using the inbuilt molecular match plugin tool.The reference compounds for the molecular match analysis were downloaded from the PubChem classification browser> spectral information> Mass Spectrometry> LC-MS, subcategory, in which 18,804 compounds were enlisted.According to the Mnova analysis, a total of the 18 best-matched phytochemical compounds are listed in Tables 2 and 3.
The compounds identified in this study were tentatively characterized using MS data, and the understanding of the observed MS/MS spectra as compared with reference compounds was obtained from PubChem using Mnova software.

LC/MS Analysis for Identification of Phytochemical Compounds in Aqueous Fraction of Phoenix dactylifera (Khudari Cultivar)
In the present work, bioactive composition from the aqueous fraction of fruit pulp of Phoenix dactylifera (Khudari cultivar) was carried out on a Waters Alliance E2695/HPLC-TQD mass spectrometer in positive-and negative-ionization mode.The output file was analyzed in the Mnova software offered by Mestrelab Research, chemistry software solutions, and the peaks were matched by using the inbuilt molecular match plugin tool.The reference compounds for the molecular match analysis were downloaded from the PubChem classification browser> spectral information> Mass Spectrometry> LC-MS, sub-category, in which 18,804 compounds were enlisted.According to the Mnova analysis, a total of the 18 best-matched phytochemical compounds are listed in Tables 2 and 3.
The compounds identified in this study were tentatively characterized using MS data, and the understanding of the observed MS/MS spectra as compared with reference compounds was obtained from PubChem using Mnova software.
The positive-and negative-ionization modes were used to identify the compounds in the aqueous fraction, with five parameters: molecular weight, retention duration, m/z [M − H] − , m/z [M − H] + , molecular formula, and fragmentation pattern [67].The representations of phytochemicals in the form of chromatograms from LC-MS data are shown in  The positive-and negative-ionization modes were used to identify the compounds in the aqueous fraction, with five parameters: molecular weight, retention duration, m/z [M − H] − , m/z [M − H] + , molecular formula, and fragmentation pa ern [67].The representations of phytochemicals in the form of chromatograms from LC-MS data are shown in Figure 4.

ADME Pharmacokinetic Profiling of the Detected Compounds
Our target is AChE, which primarily maintains the neurotransmi er homeostasis via hydrolyzing the acetylcholine (Ach) into acetic acid and choline.Neuronal activity is directly connected with the brain, so the oral bioactive compounds must penetrate the blood-brain barrier (BBB) to a ain their activity.The pharmacological measurements, as well as the BBB-penetration power of chemicals via using computational methods, are the most effective and plausible steps in the field of drug discovery [68].The BBB-penetration power was predicted by the ADME tool in the form of a boiled egg, where white region denoted the probability of gastrointestinal absorption, and yellow (yolk) showed the BBBpassing capability [69].According to the SWISS ADME analysis, there were five predicted compounds (2-Aminooctadecane, Isoxsuprine, 8-Hydroxyloxapine, 2-{5-[(1E)-3methylbuta-1,3-dien-1-yl]-1H-indol-3-yl} ethanol (MDIE), and Methyl hydroxysterpurate ethylidene acetal (MHEA)) within the yolk region, which means they can cross the BBB.

ADME Pharmacokinetic Profiling of the Detected Compounds
Our target is AChE, which primarily maintains the neurotransmitter homeostasis via hydrolyzing the acetylcholine (Ach) into acetic acid and choline.Neuronal activity is directly connected with the brain, so the oral bioactive compounds must penetrate the blood-brain barrier (BBB) to attain their activity.The pharmacological measurements, as well as the BBB-penetration power of chemicals via using computational methods, are the most effective and plausible steps in the field of drug discovery [68].The BBBpenetration power was predicted by the ADME tool in the form of a boiled egg, where white region denoted the probability of gastrointestinal absorption, and yellow (yolk) showed the BBB-passing capability [69].According to the SWISS ADME analysis, there were five predicted compounds (2-Aminooctadecane, Isoxsuprine, 8-Hydroxyloxapine, 2-{5-[(1E)-3methylbuta-1,3-dien-1-yl]-1H-indol-3-yl} ethanol (MDIE), and Methyl hydroxysterpurate ethylidene acetal (MHEA)) within the yolk region, which means they can cross the BBB.The absorption capacity across the BBB was high in 2-Aminooctadecane, Isoxsuprine, and MDIE but moderate in 8-Hydroxyloxapine and MHEA [70,71] (Figure 6).The absorption capacity across the BBB was high in 2-Aminooctadecane, Isoxsuprine, and MDIE but moderate in 8-Hydroxyloxapine and MHEA [70,71] (Figure 6).All five compounds exert considerable aqueous solubility, ranging from 2.32 to 629.45 mg/L.The permeability of Caco-2 and MDCK cells has been identified as a critical parameter in the drug development process [71][72][73].The pre-ADME analysis evaluated the permeability score of Caco-2 and MDCK cells, ere ranging from 15.32-57.68nm/s and 0.23 to 231.82 nm/s, respectively.The appropriate range for the permeability score of Caco-2 and MDCK cells is classified as low (<4 nm/s), moderate (>4-70 nm/s), and high (>70 nm/s).[74].Additionally, the skin permeability (Log Kp) of these compounds was also assessed and found appropriate for the drug candidate, ranging from −0.77 to −3.75 cm/h [75].Plasma protein binding (PPB) was found to be the most important parameter for distribution of the drug from the plasma to the target tissue, which ultimately effects the drug efficacy [76,77].Our compounds showed different efficiencies of PPB, in which 2-Aminooctadecane and Methyl hydroxysterpurate ethylidene acetal (MHEA) were found to have the highest PPB efficiency, that is, 100% (Table 4).It is commonly known that drug diffusion is inversely proportional to the efficiency of PPB.We can clearly see that the higher the PPB efficiency, the lower the drug diffusion, and the lower the PPB efficiency, the higher the drug diffusion throughout the body.All five compounds exert considerable aqueous solubility, ranging from 2.32 to 629.45 mg/L.The permeability of Caco-2 and MDCK cells has been identified as a critical parameter in the drug development process [71][72][73].The pre-ADME analysis evaluated the permeability score of Caco-2 and MDCK cells, ere ranging from 15.32-57.68nm/s and 0.23 to 231.82 nm/s, respectively.The appropriate range for the permeability score of Caco-2 and MDCK cells is classified as low (<4 nm/s), moderate (>4-70 nm/s), and high (>70 nm/s).[74].Additionally, the skin permeability (Log Kp) of these compounds was also assessed and found appropriate for the drug candidate, ranging from −0.77 to −3.75 cm/h [75].Plasma protein binding (PPB) was found to be the most important parameter for distribution of the drug from the plasma to the target tissue, which ultimately effects the drug efficacy [76,77].Our compounds showed different efficiencies of PPB, in which 2-Aminooctadecane and Methyl hydroxysterpurate ethylidene acetal (MHEA) were found to have the highest PPB efficiency, that is, 100% (Table 4).It is commonly known that drug diffusion is inversely proportional to the efficiency of PPB.We can clearly see that the higher the PPB efficiency, the lower the drug diffusion, and the lower the PPB efficiency, the higher the drug diffusion throughout the body.
Similarly, human intestinal absorption (HIA) and cytochrome P450-2D6 (CYP 2D6) are responsible for the drug absorption and excretion in the body [78][79][80].Our results demonstrated that two compounds, namely 2-Aminooctadecane and Isoxsuprine, showed inhibitory activity against CYT P450 2D6, while the rest were non-inhibitors.The inhibition of the CYTP450 2D6 via these compounds results in increased bioavailability, and clearance from circulation was also delayed.This may eventually lead to increased therapeutic efficacy of these drugs.Contrarily, CYT P450 2D6 was shown as a substrate by three compounds: 2-Aminooctadecane, Isoxsuprine, and 8-Hydroxyloxapine.This reflects their rapid metabolism, immediate therapeutic activities, and clearance from the body.

Drug-Likeness Properties of Predicted Compounds
A set of distinct factors is involved in the assessment of the drug-likeness property of a chemical entity.To date, there are several in silico approaches used to access the drug-likeness property to prevent wasting time, money, and labor; hence, we used a computational tool to analyze the ADMET parameters [68,81].Our selected compounds were tested for the same to assess Lipinski's rule of five (RO5) and the physiochemical and drug-like properties.The RO5 is a guideline for determining whether a chemical molecule with a specific pharmacological or biological activity possesses qualities that would make it a likely orally active medication in humans [82].The chemical entity to be utilized as a medicine for oral delivery, if a biologically active molecule, must meet five requirements.Poor absorption or permeation are most likely if the molecular weight >500 Da, hydrogen bond donors >5, hydrogen bond acceptor > 10, and logP (octanol water partition coefficient) >5 [82].An orally active drug-like compound should have no more than one violation of the following criteria as stated by RO5.Based on the drug-likeness analysis, our selected compounds obeyed all the parameters of Lipinski's rule of five except one compound, i.e., 2-Aminooctadecane, which violated one Lipinski's rule (LogP > 5).However, this violation is not enough to exclude this compound from the study (Table 5).These characteristics indicate that the chosen compounds passed all the drug-likeness tests and could therefore be evaluated for additional pharmacological actions.

Assessment of Toxicity of Selected Compounds
Toxicology evaluation using AI-based methods is a vital step in the drug discovery process that helps to cut down on time, resources, labor, and animal experiments.The toxicity of the selected ligands was assessed by an online web server, i.e., ProToxII.The server predicted the hepatotoxicity, carcinogenicity, immunogenicity, mutagenicity, as well as cytotoxicity along with LD 50 .The output results are tabulated in Table 6.
Our results indicated that 2-Aminooctadecane and MHEA successfully qualified all the parameters of toxicity, while MDIE indicated suitable LD 50 but is active against immunotoxicity.On the other hand, Isoxsuprine and 8-Hydroxyloxapine (toxic against immune system) showed low LD 50 , falling under class III and class II categories, respectively.Note: NT, non-toxic; T, toxic.
Biomolecules 2023, 13, x FOR PEER REVIEW 18 of 26 (Figure 8).All the other selected ligands' binding energies and amino acid residues interactions with the active pocket of AChE are given in tabular form in Table 7.The in-silico findings suggest that the presence of the natural compounds in the aqueous fraction of Phoenix dactylifera Khudari cultivar have anti-acetylcholinesterase activity, as they showed strong binding affinity towards the active pocket of AChE when compared to the tacrine.

Molecular Dynamics Simulation (MDS) Analysis
Schrödinger LLC's Desmond application was utilized to model molecular dynamics for 100 ns [83].The initial phase of protein and ligand complexes for molecular dynamics modelling was docking studies.Molecular docking studies can forecast the ligand binding state under static conditions.As a result of the fact that docking provides a static representation of a molecule's binding posture at a protein's active site [84], MD simulations generally reproduce the motion of atoms over time by integrating Newton's classical theory of motion, where the physiological environment's ligand binding state is predicted using simulations [85,86].

RMSD and RMSF Data Analysis
As a function of simulation time, the RMSD was calculated by deviating from the structure of the protein-ligand starting posture [67].In this experiment, we assessed the RMSD of the target proteins (AChE) in the presence of the substance "2-5-[(1E)-3methylbuta-1,3-dien-1-yl]-1H-indol-3-ylethanol" (PubChem CID: 90675402) during a simulation duration of 100 ns (Figure 9).During the initial few seconds and throughout the duration of the simulation, the RMSD of AChE grew.However, in the presence of 90675402, the RMSD of AChE varied from 0 to 50 ns before reaching equilibrium for 50 to 100 ns.The AChE's average RMSD values were 2.52 0.06 and 3.22 0.25, respectively, during 50 and 100 ns in the presence and absence of 90675402 (Figure 9).Over time, ligand-bound proteins' C-alpha atoms' RMSD values varied.The proteins in the complex a ained stability at 20 ns, as shown in the RMSD plot.After that, the simulation's RMSD value fluctuations remained within 1.0 angstrom, which is acceptable.Until 25 ns, the RMSD values of ligands fi ed to proteins fluctuated within 1.0 angstrom, after which The in-silico findings suggest that the presence of the natural compounds in the aqueous fraction of Phoenix dactylifera Khudari cultivar have anti-acetylcholinesterase activity, as they showed strong binding affinity towards the active pocket of AChE when compared to the tacrine.

Molecular Dynamics Simulation (MDS) Analysis
Schrödinger LLC's Desmond application was utilized to model molecular dynamics for 100 ns [83].The initial phase of protein and ligand complexes for molecular dynamics modelling was docking studies.Molecular docking studies can forecast the ligand binding state under static conditions.As a result of the fact that docking provides a static representation of a molecule's binding posture at a protein's active site [84], MD simulations generally reproduce the motion of atoms over time by integrating Newton's classical theory of motion, where the physiological environment's ligand binding state is predicted using simulations [85,86].

RMSD and RMSF Data Analysis
As a function of simulation time, the RMSD was calculated by deviating from the structure of the protein-ligand starting posture [67].In this experiment, we assessed the RMSD of the target proteins (AChE) in the presence of the substance "2-5-[(1E)-3-methylbuta-1,3dien-1-yl]-1H-indol-3-ylethanol" (PubChem CID: 90675402) during a simulation duration of 100 ns (Figure 9).During the initial few seconds and throughout the duration of the simulation, the RMSD of AChE grew.However, in the presence of 90675402, the RMSD of AChE varied from 0 to 50 ns before reaching equilibrium for 50 to 100 ns.The AChE's average RMSD values were 2.52 0.06 and 3.22 0.25, respectively, during 50 and 100 ns in the presence and absence of 90675402 (Figure 9).Over time, ligand-bound proteins' C-alpha atoms' RMSD values varied.The proteins in the complex attained stability at 20 ns, as shown in the RMSD plot.After that, the simulation's RMSD value fluctuations remained within 1.0 angstrom, which is acceptable.Until 25 ns, the RMSD values of ligands fitted to proteins fluctuated within 1.0 angstrom, after which there was a flip in the ligand mode, with regained stability at 50 ns that remained consistent throughout the simulation.To detect the local conformational changes in a protein's side chains brought on by the binding of a ligand, RMSF was frequently monitored [87].Here, in the presence and absence of the chemical ligand (90675402), the RMSF of AChE was calculated (Figure 10).Proteins' greater degrees of flexibility at their N-and C-termini cause fluctuations there.In the presence of 90675402, the average RMSF values for AChE were 1.196 and 0.67.The residues with higher maxima belong to loop regions or the N and C termini, according to MD trajectories.Low RMSF values of putative binding residues show the stability of the ligand binding to the protein.The compactness of the protein-ligand complex was quantified by the radius of gyration, and its exposure to solvent molecules was determined by the solvent-accessible surface area.Each of these traits provides insight into the stability of the protein-ligand To detect the local conformational changes in a protein's side chains brought on by the binding of a ligand, RMSF was frequently monitored [87].Here, in the presence and absence of the chemical ligand (90675402), the RMSF of AChE was calculated (Figure 10).Proteins' greater degrees of flexibility at their N-and C-termini cause fluctuations there.In the presence of 90675402, the average RMSF values for AChE were 1.196 and 0.67.The residues with higher maxima belong to loop regions or the N and C termini, according to MD trajectories.Low RMSF values of putative binding residues show the stability of the ligand binding to the protein.To detect the local conformational changes in a protein's side chains brought on by the binding of a ligand, RMSF was frequently monitored [87].Here, in the presence and absence of the chemical ligand (90675402), the RMSF of AChE was calculated (Figure 10).Proteins' greater degrees of flexibility at their N-and C-termini cause fluctuations there.In the presence of 90675402, the average RMSF values for AChE were 1.196 and 0.67.The residues with higher maxima belong to loop regions or the N and C termini, according to MD trajectories.Low RMSF values of putative binding residues show the stability of the ligand binding to the protein.The compactness of the protein-ligand complex was quantified by the radius of gyration, and its exposure to solvent molecules was determined by the solvent-accessible  The compactness of the protein-ligand complex was quantified by the radius of gyration, and its exposure to solvent molecules was determined by the solvent-accessible surface area.Each of these traits provides insight into the stability of the protein-ligand combination during the simulation [87].In this study, the Rg of AChE was determined in the presence of 90675402 during 0-100 ns, as shown in Figure 11A.The Rg of the AChE-90675402 complex varied in the range of 22.84-23.32Å, with an average value of 23.15 ± 0.06 Å.Similarly, the SASA of AChE was determined in the presence of 90675402 during 0-100 ns (Figure 11B).The SASA of AChE-90675402 complex fluctuated in the range of 20,216.44-22,252.01Å 2 , with an average value of 21,186.77± 384.72 Å 2 .These results demonstrate that changes in the Rg and SASA of the target protein in the presence of 90675402 did not significantly differ, supporting the continued stability of the proteinligand complexes.For the purpose of determining if a protein has undergone any structural changes as a result of ligand binding, it is important to track changes in the secondary structure elements (SSE) of the protein in a protein-ligand complex [88].The total SSE (α-helix + βsheet) of AChE in the presence of 90675402 was assessed here as a function of simulation

SSEs, i.e., Secondary Structural Elements Study
For the purpose of determining if a protein has undergone any structural changes as a result of ligand binding, it important to track changes in the secondary structure elements (SSE) of the protein in a protein-ligand complex [88].The total SSE (α-helix + β-sheet) of AChE the presence of 90675402 was assessed here as a function of simulation duration (Figure 12).AChE's SSE was 42.01 2.84% on average when it was in a complex with 90675402 (-helix = 28.63 2.96% and -sheet = 13.38 0.24%).The results reveal a stable protein structure in the complex of protein and ligand since the total SSE for all the targeted proteins did not significantly change when 90675402 was present.
Biomolecules 2023, 13, x FOR PEER REVIEW 21 of 26 duration (Figure 12).AChE's SSE was 42.01 2.84% on average when it was in a complex with 90675402 (-helix = 28.63 2.96% and -sheet = 13.38 0.24%).The results reveal a stable protein structure in the complex of protein and ligand since the total SSE for all the targeted proteins did not significantly change when 90675402 was present.

Interactions Analysis
The majority of significant ligand-protein interactions discovered by MD are hydrogen bonds and hydrophobic interactions, as seen in Figure 13.TRP_84 and GLY_442 are the most critical hydrophobic contacts for the complex, although GLU_199 is important for H bonds.The stacked bar graphs were standardized along the trajectory; for instance, a value of 1.0 denotes that the precise interaction was maintained for 100% of the simulation time.Values above 1.0 are achievable simply owing to the possibility of numerous interactions of the same subtype involving a few protein residues and the ligand (Figure 13A,B).

Prime/MM-GBSA Free Energy Analysis
A flexible and reliable way to assess the stabilization of a protein-ligand complex is free energy (Prime/MM-GBSA) computation [89].After every 10 ns, free energy and its components were calculated using the Prime/MM-GBSA method, and the results are

Interactions Analysis
The majority of significant ligand-protein interactions discovered by MD are hydrogen bonds and hydrophobic interactions, as seen in Figure 13.TRP_84 and GLY_442 are the most critical hydrophobic contacts for the complex, although GLU_199 is important for H bonds.The stacked bar graphs were standardized along the trajectory; for instance, a value of 1.0 denotes that the precise interaction was maintained for 100% of the simulation time.Values above 1.0 are achievable simply owing to the possibility of numerous interactions of the same subtype involving a few protein residues and the ligand (Figure 13A,B).duration (Figure 12).AChE's SSE was 42.01 2.84% on average when it was in a complex with 90675402 (-helix = 28.63 2.96% and -sheet = 13.38 0.24%).The results reveal a stable protein structure in the complex of protein and ligand since the total SSE for all the targeted proteins did not significantly change when 90675402 was present.

Interactions Analysis
The majority of significant ligand-protein interactions discovered by MD are hydrogen bonds and hydrophobic interactions, as seen in Figure 13.TRP_84 and GLY_442 are the most critical hydrophobic contacts for the complex, although GLU_199 is important for H bonds.The stacked bar graphs were standardized along the trajectory; for instance, a value of 1.0 denotes that the precise interaction was maintained for 100% of the simulation time.Values above 1.0 are achievable simply owing to the possibility of numerous interactions of the same subtype involving a few protein residues and the ligand (Figure 13A,B).

Prime/MM-GBSA Free Energy Analysis
A flexible and reliable way to assess the stabilization of a protein-ligand complex is free energy (Prime/MM-GBSA) computation [89].After every 10 ns, free energy and its components were calculated using the Prime/MM-GBSA method, and the results are

Conclusions
Cholinesterase inhibition can manage Alzheimer's disease (AD) by reducing the catalysis of acetylcholine, which ultimately supports the signals transmission and improves the life quality of AD patients.Khudari, a cultivar of Phoenix dactylifera that is mostly cultivated in the Al-Qassim region, was found to possess a higher total phenol content and good radical scavenging (DPPH) potential.The Khudari aqueous extract showed the significant inhibition of acetylcholinesterase in a dose-dependent manner, with an enzyme inhibition kinetic study revealing that it exhibits a mixed type of inhibition.Further, for phytochemicals investigation, we used LC/MS analysis, which revealed the presence of eighteen major compounds predicted through the PubChem database.Molecular docking and ADMET analysis confirmed that compound, 2-{5-[(1E)-3-methylbuta-1,3-dien-1-yl]-1H-indol-3-yl} ethanol (MDIE) has better binding affinity than tacrine against acetylcholinesterase and was found to be safer, exert considerable aqueous solubility and acceptable LD 50 , and to easily penetrate the BBB.Molecular dynamics analysis of 100 ns and Prime/MM-GBSA supports the stability of the protein-ligand complexes.These findings suggest that aqueous extract of Khudari fruit pulp has significant antioxidant and acetylcholinesterase inhibition potentials and its compound, MDIE, forms a stable confirmation with the target protein, though the fruit of Khudari dates can be a better functional food for the management of Alzheimer's disease.Further in vivo investigation is required to fully characterize the pharmacokinetic properties, optimization of dose administration, and efficacy of this plant-based natural compound.

Figure 1 .
Figure 1.Radical (DPPH) scavenging potentials of aqueous and methanolic extract of Khudari fruit pulp.The data are represented as mean ± S.D. in triplicate value.

Figure 1 .
Figure 1.Radical (DPPH) scavenging potentials of aqueous and methanolic extract of Khudari fruit pulp.The data are represented as mean ± S.D. in triplicate value.

Figure 2 .
Figure 2. Correlation graph of TPC and FRAP values of aqueous extract of Khudari fruit pulp.

Figure 2 .
Figure 2. Correlation graph of TPC and FRAP values of aqueous extract of Khudari fruit pulp.

Figure 3 .
Figure 3. Lineweaver-Burk double reciprocal and K m and V max plot of 1/v versus 1/[S] of (A,B) tacrine and (C,D) Khudari aqueous extract against AChE.

Figure 6 .
Figure 6.The boiled-egg representation of the selected compounds.Five compounds were in the yellow region, i.e., PubChem IDs 3126, 3783, 43655, 90675402, and 139583859, indicating that they can cross the BBB.

Figure 6 .
Figure 6.The boiled-egg representation of the selected compounds.Five compounds were in the yellow region, i.e., PubChem IDs 3126, 3783, 43655, 90675402, and 139583859, indicating that they can cross the BBB.

Figure 8 .
Figure 8.In silico binding pa erns of standard drug tacrine (PubChem ID: 1935) against the active pocket of AChE (PDB ID: 1ACJ).(A1, A2) Two-dimensional interactions of tacrine within the active pocket of AChE and A2 panel showed the interaction of tacrine with AChE surrounded by α-helix and β-sheet conformations of the active domain of AChE.

Figure 8 .
Figure 8.In silico binding patterns of standard drug tacrine (PubChem ID: 1935) against the active pocket of AChE (PDB ID: 1ACJ).(A 1 ,A 2 ) Two-dimensional interactions of tacrine within the active pocket of AChE and A 2 panel showed the interaction of tacrine with AChE surrounded by α-helix and β-sheet conformations of the active domain of AChE.
Biomolecules 2023, 13, x FOR PEER REVIEW 19 of 26there was a flip in the ligand mode, with regained stability at 50 ns that remained consistent throughout the simulation.

Biomolecules 2023 ,
13, x FOR PEER REVIEW 19 of 26there was a flip in the ligand mode, with regained stability at 50 ns that remained consistent throughout the simulation.

Figure 12 .
Figure 12.Protein secondary distribution of structural elements by residue index across the complexed protein structures.Alpha helices are denoted by red columns, whereas beta strands are denoted by blue columns.

Figure 13 .
Figure 13.(A) Protein-ligand interaction histogram of the complex; (B) Residues that interacted with the ligand.

Figure 12 .
Figure 12.Protein secondary distribution of structural elements by residue index across the complexed protein structures.Alpha helices are denoted by red columns, whereas beta strands are denoted by blue columns.

Figure 12 .
Figure 12.Protein secondary distribution of structural elements by residue index across the complexed protein structures.Alpha helices are denoted by red columns, whereas beta strands are denoted by blue columns.

Figure 13 .
Figure 13.(A) Protein-ligand interaction histogram of the complex; (B) Residues that interacted with the ligand.

Figure 13 .
Figure 13.(A) Protein-ligand interaction histogram of the complex; (B) Residues that interacted with the ligand.

Table 4 .
Pharmacokinetic profiling of the selected compounds predicted in the Khudari aqueous extract.

Table 4 .
Pharmacokinetic profiling of the selected compounds predicted in the Khudari aqueous extract.

Table 5 .
Physicochemical properties of the predicted compounds.

Wt. TPSA Hydrogen Bond Acceptor Hydrogen Bond Donor No. of Violations No. of Rotatable Bond
* Standard drug.

Table 6 .
Toxicity assessment of the selected five compounds.

Table 7 .
The binding pattern of potential Phoenix-dactylifera-derived natural selected compounds against AChE.
[87]olecules 2023, 13, x FOR PEER REVIEW 20 of 26 combination during the simulation[87].In this study, the Rg of AChE was determined in the presence of 90675402 during 0-100 ns, as shown in Figure11A.The Rg of the AChE-90675402 complex varied in the range of 22.84-23.32Å,withanaveragevalue of 23.15 ± 0.06 Å.Similarly, the SASA of AChE was determined in the presence of 90675402 during 0-100 ns (Figure11B).The SASA of AChE-90675402 complex fluctuated in the range of 20,216.252.01Å 2 , with an average value of 21,186.77± 384.72 Å 2 .These results demonstrate that changes in the Rg and SASA of the target protein in the presence of 90675402 did not significantly differ, supporting the continued stability of the proteincomplexes.

Table 8 .
Before and after simulations, MM-GBSA calculations were performed.