Computer-assisted evaluation of plant-derived β-secretase inhibitors in Alzheimer’s disease

Alzheimer’s disease (AD) is a progressive neurodegenerative age-related dementia that results in memory loss of elderly people. Many hypotheses have been formally articulated till now to decipher the pathogenesis of this disease. According to the compelling amyloidogenic hypothesis, β-secretase is a key regulatory enzyme in AD development and is therefore considered as one of the major targets for the development of drugs to treat AD. In this study, 40 plant-derived phytocompounds, proven to have β-secretase inhibitory activity in different laboratory experiments, were evaluated using computational approaches in order to identify the best possible β-secretase inhibitor(s). Amentoflavone (IFD score: − 7.842 Kcal/mol), Bilobetin (IFD score: − 7.417 Kcal/mol), and Ellagic acid (IFD score: − 6.923 Kcal/mol) showed highest β-secretase inhibitory activities with high binding affinity among all the selected phytocompounds and interacted with key amino acids, i.e., Asp32, Tyr71, and Asp228 in the catalytic site of β-secretase. Moreover, these three molecules exhibited promising results in different drug potential assessment experiments and displayed signs of correlation with significant pharmacological and biological activities. Amentoflavone, Biolbetin, and Ellagic acid could be investigated further in developing β-secretase-dependent drug for the effective treatment of AD. However, additional in vivo and in vitro experiments might be required to strengthen the findings of this experiment.


Background
Alzheimer's disease (AD) is a perilous and progressive neurodegenerative disorder that worsens over time and ultimately culminates in dementia [1,2]. AD manifests different types of characteristic symptoms, such as cognitive dysfunction (memory loss and language difficulties), behavioral disturbances (depression and hallucinations, also known as non-cognitive symptoms), and impairment in daily activities in affected individuals [3][4][5]. AD predominantly results from genetic heritability, age, and abnormal deposition of proteins inside as well as outside of brain cells. These proteins include amyloid β (Aβ, extracellularly) which is derived from the proteolysis of amyloid precursor protein (APP). APP is not only an integral but also a type 1 transmembrane protein that is found in various types of tissues but predominantly in neurons. It has been suggested that APP might be responsible for the formation of synapses and repairment of neurons [6]. Extracellular deposition of amyloid β proteins called amyloid plaques constitute the fundamental pathogenic driver of AD. Generated by the proteolytic activity of the βsecretase enzyme, it is the prime drug target for the inhibition of Aβ production in AD [7]. Presently, it is a matter of great concern as AD is the sixth prevalent cause of death in the USA, killing more people than other diseases such as breast and prostate cancers. Astonishingly, deaths from AD have increased by 89% since 2000 such that it is now estimated that someone develops AD in every 33 s in the USA. Much to our dismay, the number will continue to rise as high as 16 million by 2050 [8].
Amyloidogenic and non-amyloidogenic pathways and role of β-secretase in AD β-secretase, also known as beta-site amyloid precursor protein cleavage enzyme 1 (BACE 1), is an aspartyl protease enzyme. Both APP and β-secretase mature in the endoplasmic reticulum (ER) and are later transported to the plasma membrane through the participation of the Golgi body ( Fig. 1) [9,10]. As endosomes internalize more than half of the APP, only a small amount is cleaved by α-secretase enzyme (located on the cell surface). Some of the APP are transported back to the plasma membrane, while some are degraded by lysosomes. The enzyme α-secretase competes with the βsecretase in the trans-Golgi network (TGN) for cleaving and processing APP. β-secretase is located in the TGN and endosomes [11]. It causes impairment in the trafficking of APP in the cell surface upon finding an endosomal dysfunction. The internalization of APP increases with the mediation of the β-secretase cleavage of APP [12,13]. The cleaved portion of APP by α or β-secretase is further cleaved by γ-secretase (found in the ER), TGN, and endosomes [11].
Several hypotheses propose that β-secretase enzyme generates Aβ protein in neurons. The active site of the enzyme βsecretase is located in the extracellular domain and contains two key aspartate residues: Asp332 and Asp28 [11,12].
In normal health condition (in the absence of βsecretase) and non-amyloidogenic pathway, amyloid precursor protein (APP) is cleaved sequentially. Firstly, the amino-terminal portion of the protein is cleaved by αsecretase. This secretes APP and produces carboxyterminal fragment (CTF83), which precludes the generation of Aβ protein. CTF83 is further cleaved by γ-secretase and produces P3, a soluble and amino-terminal APP intracellular domain (AICD). On the contrary, the βsecretase enzyme processes APP through another pathway named the amyloid pathway (for AD). In this pathway, β-secretase enzyme trims the APP generating secreted APP (sAPPβ) and carboxyterminal fragment (CTF99). They are further cleaved by γ-secretase and generates highly insoluble Aβ protein and AICD. These insoluble Aβ proteins lead to the formation of amyloid plaques and impair the activity of neurons. These plaques destroy synapses, thus contributing to the development of AD [13][14][15][16]. Fig. 1 Role of β-secretase in Alzheimer's disease and amyloidogenic and non-amyloidogenic pathways. In non-amyloidogenic pathway, cleavage of APP sequentially by αand γ-secretase produces sAPPα and soluble P3, respectively. On the other hand, in amyloidogenic pathway, β-and γsecretase produces sAPPβ and insoluble Aβ protein, respectively, leading to the formation of β amyloid plaques. Matured APP and β-secretase are transported to the plasma membrane by Golgi body. A portion of APP is degraded by lysosome

Treatment of Alzheimer's disease
Acetylcholinesterase inhibitor drugs are the most prescribed drugs for the treatment of Alzheimer's disease owing to their beneficial effects in cognitive, non-cognitive, and behavioral symptoms. These drugs include donepezil (reversible and selective inhibitor), rivastigmine (selective and pseudoreversible inhibitor), galantamine, tacrine (irreversible and nonselective inhibitor), metrifonate (pseudoreversible inhibitor), eptastigmine (reversible and selective inhibitor) as well as a newly approved drug, memantine, an NMDA (N-methyl-D-aspartate receptor) channel blocker [11,17,18]. Although these drugs are widely used to treat AD, several reports have shown that some of these drugs are not useful in mild cognitive impairment. Some of them are not effective in severe cognitive impairment either and have various side effects including nausea, vomiting, decreased appetite, and gastric acid production. In some cases, patients even show symptoms of hallucination, fatigue, dizziness, etc. Therefore, there are a few detrimental effects that result from the prolonged use of these drugs [19,20]. Complementary medicines, such as extracts from medicinal plants, also have beneficial effects on cognitive and non-cognitive symptoms of AD. These medicines are thought to have less toxicity or side effects [19].
In this study, we evaluated the potentiality of 40 plant-derived phytocompounds to inhibit βsecretase which have already shown β-secretase inhibition activity in different laboratory experiments. We selected 3 best inhibitors from those 40 compounds and subsequently analyzed them for their drug-like potentials (Fig. 2).

Selection of ligand molecules
An extensive literature survey was carried out to identify phytochemicals from plants that showed β-secretase inhibition in laboratory experiments. Forty ligands were selected from literature review that displayed a βsecretase inhibitory effect to varying extents in different in vitro assays. Table 1 enlists the selected compounds from various plant sources that were used in the next phases of this study.

Molecular docking study Protein preparation
The three-dimensional crystallographic structure of human β-secretase (PDB ID:2OHM) was downloaded in PDB format from Protein Data Bank (www.rcsb.org) [50]. The structure was then prepared and processed using the Protein Preparation Wizard in Maestro Schrödinger Suite (v11.4). Bond orders were assigned to the structures and hydrogens were added to heavy atoms. All of the water molecules were erased from the atoms, missing side chains were added to the protein backbone using Prime, and the states were generated with Epik at pH 7 ± 2 [51]. At last, the structures were refined and then minimized utilizing the Optimized Potentials for Liquid Simulations force field (OPLS_ 2005). Minimization was performed by setting the greatest substantial particle RMSD (root-meansquare deviation) to 30 Å, and any extraordinary water under 3H bonds to non-water was again eradicated during the minimization step.

Ligand preparation
A total of 40 selected ligand molecules were downloaded in SDF format from the PubChem database (https:// pubchem.ncbi.nlm.nih.gov/). These structures were then processed and prepared using the LigPrep wizard of Maestro Schrödinger suite [52]. Minimized 3D structures of ligands were generated using Epik2.2 within pH 7.0 ± 2.0 in the suite. Minimization was again carried out using the OPLS_2005 force field which generated a maximum of 32 possible stereoisomers depending on available chiral centers on each molecule.

Receptor grid generation
Grid usually restricts the active site to a particular area of the receptor protein for the ligand to dock specifically within that area. A grid was generated using default Van der Waals radius scaling factor 1.0 and charge cutoff 0.25, which was then subjected to OPLS_2005 force field for the minimized structure in Glide [53]. A cubic box was generated around the active site (co-crystallized reference ligand) of the target molecule. Then the grid box dimension was adjusted to 14 Å × 14 Å ×14 Å for docking to be carried out.

Glide standard precision (SP) and extra precision (XP) ligand docking
Usually, extra precision (XP) ligand docking is more accurate in predicting the binding affinity of small ligand libraries than standard precision (SP) ligand docking; the latter is generally recommended for large ligand collections and sometimes may come with imperfection [54]. However, both of the docking methods were applied for the selected ligand molecules to make a comparison between two different docking scores. The Van der Waals radius scaling factor and charge cutoff was set to 0.80 and 0.15 respectively for all the ligand molecules under study. The final score was assigned according to the pose of docked ligand within the active site of the receptor molecules. The best possible poses and types of ligand-receptor interactions were then analyzed using the Discovery Studio Visualizer (v4.5) [55].

Prime MM-GBSA rescoring
After SP and XP ligand docking, the ligands were then again subjected to MM-GBSA (Molecular mechanicsgeneralized born and surface area) rescoring using the Prime module of Maestro Schrödinger suite for further evaluation. This technique utilizes the docked complex and uses an implicit solvent, which then assigns a more accurate scoring function and improves the overall free binding affinity score upon the reprocessing of the complex [54,56]. It combines OPLS molecular mechanics energies (E MM ), a surface generalized born solvation model for polar solvation (G SGB ), and a nonpolar salvation term (G NP ) for total free energy (ΔG bind ) calculation. The total free energy of binding was calculated by the following equation:

Induced fit docking
Ten compounds were selected based on the MM-GBSA score for further evaluation as it is a more rigorous scoring method. After that, these ligands were subjected to induced fit docking (IFD), which is an even more accurate docking method to generate the native poses of the ligands from different sources [57]. Again, OPLS_2005 force field was applied after generating a grid around the co-crystallized ligand of the receptor. This time the ten best ligands docked rigidly. Receptor and Ligand Van Der Waals screening was set at 0.70 and 0.50, respectively. Residues within 2 Å were refined to generate the 2 best possible poses with extra precision. Three best ligands were selected based on the IFD score.

Pharmacological and biological activity prediction
Pharmacological activities of the three best ligand molecules were predicted using PASS online server (http:// www.pharmaexpert.ru/passonline/) and the biological activities with GPCR (G protein-coupled receptor) ligand, ion channel, enzymes, etc., were predicted again using Molinspiration Cheminformatics server (https:// www.molinspiration.com/cgi-bin/properties) [61,62]. Prediction of Activity Spectra for Substances (PASS) estimates the pharmacological activities based on the compound's native chemical structure. PASS predicts the function of a compound based on Structure-Activity Relationship Base (SAR Base), which assumes that the activity of a compound is related to its structure. It functions by comparing the 2D structure of an unknown compound to other well-known compounds possessing specific biological activities existing in the database with almost 95% accuracy [63].

Density functional theory (DFT) calculation
Minimized ligand structures obtained from LigPrep were used for DFT calculation using the Jaguar panel of Maestro Schrödinger Suite. Becke's three-parameter exchange potential and Lee-Yang-Parr correlation functional (B3LYP) theory with 6-31G* basis set were used for DFT calculation [64][65][66][67]. Quantum chemical properties such as surface properties (MO, density, potential) and multipole moments were calculated along with HOMO (highest occupied molecular orbital) and LUMO (lowest unoccupied molecular orbital) energy. Afterwards, the global frontier orbital was analyzed, and hardness (η) and softness (S) of selected molecules were calculated using the following equation as per Parr and Pearson interpretation and Koopmans theorem [68,69].

Molecular docking study
In this experiment, 40 selected ligand molecules that have been shown to have in vitro β-secretase inhibitory effect were docked against an intended target protein.
Initially, 10 ligand molecules were selected based on free binding energy (Table 2). A significant correlation among the IC 50 value, SP docking score, XP docking score, and free binding energy was observed for the molecules having molecular weight between 330 and 580 Da with minimal exception. However, a slight variation between SP and XP docking scores was observed among the different ligand molecules. A total of 10 ligand molecules were selected based on MM-GBSA docking scores (Table 2). Afterwards, the 10 best-selected ligands were subjected to IFD which is a powerful and accurate docking method for generating poses and assigning binding scores [70]. Finally, 3 best-performed ligands were selected from IFD study which were then analyzed further in the subsequent phases of different experiments.
Binding mode of amentoflavone, bilobetin, and ellagic acid with β-secretase Amentoflavone, bilobetin, and ellagic acid showed superior IFD score and XP Gscore among the 10 selected ligand molecules ( Table 2). Both the hydrogen bonds and hydrophobic interactions play vital roles in drug-receptor interaction by strengthening and specifying the interaction between ligand and target molecules [71]. All of the ligand molecules formed a significant number of hydrogen and hydrophobic interactions with respective amino acids inside the binding pocket of the target (Table 3 and Fig. 3). It has been observed that amentoflavone docked with β-secretase with an IFD score of − 823.501 kcal/mol and XP Gscore of − 7.842 kcal/mol interacted with 5 amino acids within the binding pocket and formed a total of 8 interactions (Table 5). It formed 4 conventional hydrogen bonds with Asp228, Phe108, Lys107, and Gly34 amino acid residues at 3.04-, 1.88-, 1.97-, and 2.34-Å distance apart respectively within the binding site of βsecretase. Moreover, amentoflavone also formed additional hydrophobic interactions, i.e., Pi-Pi stacked and Pi-Pi T shaped with Phe108 and Tyr71 amino acid residues (Fig. 3).
The bilobetin docked with β-secretase with an IFD score of − 821.327 kcal/mol and XP Gscore of − 7.417 kcal/mol and interacted with 6 amino acids within the binding pocket and formed a total of 9 interactions within the binding pocket (Table 5). It formed 3 conventional hydrogen bonds with Phe108, Lys107, and Tyr198 amino acid residues at 1.86-, 1.94-, and 2.03-Å distance apart, respectively. Moreover, amentoflavone also formed 2 non-conventional hydrogen bonds with Asp228, and Gly34 amino acid residues at 2.81-and 2.69-Å distance apart respectively, within the binding cleft. It also formed few hydrophobic interactions (i.e., Pi-Pi stacked and Pi-Pi T-shaped interactions) with Tyr71 and Phe108 amino acid residues (Fig. 3). Ellagic acid docked with β-secretase with an IFD score of − 819.500 kcal/mol and XP Gscore of − 6.923 kcal/ mol, interacted with 3 amino acids within the binding pocket and formed a total of 7 interactions within the binding pocket (Table 5). It formed 2 conventional hydrogen bonds with Thr72 amino acid residue at 1.99and 2.01-Å distance apart respectively. Moreover, amentoflavone also formed another conventional hydrogen bond with Asp32 amino acid residue at 1.97-Å distance apart respectively within the binding cleft. It was also reported to form other hydrophobic interactions (i.e., Pi-Pi stacked and Pi-Pi T-shaped interactions with only Tyr71 amino acid residue (Fig. 3).

ADME/toxicity profiling
The three best-ligand molecules (Table 4) were subjected to analyze their absorption, distribution, metabolism, excretion, and toxicity profiles. The result of ADME/toxicity analysis is represented in Table 5. All of the ligand molecules were reported to be highly absorbed in the intestine and have low Caco-2 permeability. Only ellagic acid exhibited high oral bioavailability.
None of them reported having the ability to cross the blood-brain barrier, and no ligand showed the sign to be P-glycoprotein substrate. Only bilobetin was reported to be a P-glycoprotein inhibitor among the three selected ligand molecules. Ellagic acid was shown to be neither a substrate nor an inhibitor of any of the selected enzymes of the cytochrome P450 family. Both amentoflavone and bilobetin were reported to be the substrate of only CYP3A4. None of the selected ligands was predicted to be OCT2 (organic cation transporter 2) substrate. All of the selected ligands were predicted to induce hepatotoxicity, and only bilobetin was reported to inhibit the hERG (human Ether-a-go-go-related gene) channel. Moreover, ellagic acid was reported to cause eye irritation. Bilobetin showed type III acute oral toxicity, whereas the other two ligand molecules were reported to have type II oral toxicity.

Pharmacological and biological activity prediction
The pharmacological connection of the ligand molecules with other neurological diseases, proteins, and enzymes involved in AD was analyzed using PASS online (Table 6). PASS predicts the result in the context of probability of activity (Pa) and inactivity (Pi) of a compound. This result varies between 0.000 and 1.000. The activity is considered possible for a compound only when Pa > Pi [72]. When Pa > 0.7, the compound is very likely to exhibit the activity, but its possibility of being an analogue to a known pharmaceutical is also high. When 0.5 < Pa < 0.7, the compound is likely to exhibit the activity, but the probability along with the chance of being a known pharmaceutical agent is also lower. When Pa < 0.5, the compound is less likely to exhibit the activity [73]. Ellagic acid was reported to have better pharmacological activity with the highest probability scores. Amentoflavone and bilobetin were predicted to have similar pharmacological activity. The three best ligands were analyzed to predict their involvement in biological activities with GPCR ligand, ion channels, enzyme, etc. The result of this experiment is summarized in Table 7. GPCRs (G protein-coupled receptors), ion channels, enzymes, nuclear receptors, etc., are the most potent drug targets in the human body. Among them, only GPCRs are the targets of almost 50% of currently available drugs [74,75]. Amentoflavone showed better biological activities followed by bilobetin and ellagic acid. Ellagic acid had the highest positive score. Amentoflavone and bilobetin also reported showing positive scores as enzyme inhibitors. The findings of biological activity prediction indicate the important therapeutic significance of these molecules, but as useless may raise the concern of adverse effects.

Analysis of frontier orbitals
Density functional theory calculation aids to understand the energetics, structure, and properties of molecules at a lower cost. Analysis of the frontier orbitals of molecules allows us to understand the pharmacological properties of the molecule under investigation [76,77]. The detailed HOMO (highest occupied molecular orbital) energy, LUMO (lowest unoccupied molecular orbital) energy, energy gap (HOMO-LUMO gap), hardness, and softness of the selected four compounds are summarized in Table 8. The occupation of HOMO and LUMO is illustrated in Fig. 4 for each compound. HOMO and LUMO orbitals help to understand the kinetic stability and chemical reactivity of a compound. HOMO usually occupies a small region in a molecule with a higher ability to donate one or more electrons during bond formation. LUMO represents the region capable of accepting electron(s) from another macromolecule (Fig. 4) [78]. The compound with higher HOMO-LUMO gap is considered energetically unfavorable to undergo a chemical reaction [79]. The highest gap was observed for Amentoflavone, whereas Ellagic acid was reported to have the lowest energy gap. According to the energy gap, the order of the compounds is amentoflavone > bilobetin > ellagic acid. The dipole moment of each compound was also calculated along with HOMO and LUMO energy, the scores of which are also reported in Table 8. According to the dipole moment score, the stability order of compounds is ellagic acid > amentoflavone > bilobetin.

Discussion
Molecular docking is one of the most commonly used computer-aided techniques in structure-based drug designing. It defines the best possible orientation of a small ligand molecule when bound to the binding site of a second target macromolecular target [80,81]. This technique utilizes a specific scoring algorithm and assigns binding energy to the ligand molecules based on the   poses that fit with the target. It reflects the binding affinity of a ligand molecule for the target. The low binding energy of a ligand-receptor complex indicates high stability of the ligand-receptor complex. It confers that the complex remains more time in contact [82]. Besides being a more rigorous scoring that guides the post-processing of docking experiment, MM-GBSA scoring predicts the most accurate free binding energy of ligand-receptor complex [83][84][85]. In this experiment, the three ligands-amentoflavone, bilobetin, and ellagic acid-were found to be the best inhibitors of β-secretase upon sequential molecular docking experiment (Tables 2 and 3) (Fig. 3).
In silico ADME/toxicity analysis allows rigorous pharmacokinetic property and toxicity testing. This analysis is required to confirm whether a drug should sustain the Phase I clinical trial or not and, in turn, assists in the in vitro assays of the candidate drug [86,87]. Bloodbrain barrier permeability is required to be confirmed for the drugs that primarily target the cells of the central nervous system (CNS). Since the oral delivery system is the most frequently used route of drug delivery and the delivered drug travels through the digestive tract into the intestine, it is expected that the drug is highly absorbed in the human intestinal tissue. P-glycoproteins are the embedded glycoproteins on the cell membrane that are responsible for facilitating the transport of many drugs through the cell membrane. As a result, their inhibition may affect the normal drug transport inside the human body. Caco-2 cell line is usually used for the assessment of drug permeability of a new candidate which, in turn, reflects the human intestinal tissue permeability [88][89][90][91][92]. Cytochrome P450 family of enzymes is responsible to control drug interaction, metabolism, and excretion inside the body. Inhibition of these enzymes may lead to acute drug toxicity, slow clearance, and eventually malfunction of the drug compound inside the human body [93][94][95]. The AMES toxicity examines the toxicity of chemicals [96,97]. hERG (Human ether-a-gogo-related gene) channels are the voltage-gated potassium ion channels that play key roles in potassium ion transport along the cell membrane. Different structurally and functionally unrelated drugs have been reported to block the hERG potassium channel, raising the concern of off-target drug interaction. Therefore, it is imperative to screen compounds for activity on hERG channels early in the lead optimization process of a drug discovery approach to reduce the risk of a drug candidate failing in preclinical safety studies due to the blockade of hERG channels [98]. Renal OCT2 (organic cation transporter 2) is important for drug and xenobiotic excretion through the kidney. The substrates of this transporter protein are thought to be excreted easily with urine [99]. All of the selected ligand molecules were reported to have similar ADME/T properties ( Table 5).
The best-selected molecules were also analyzed for respective pharmacological activity and biological activity. Ellagic acid and amentoflavone reported having better pharmacological and biological activity, respectively (Tables 6 and 7). The analysis of frontier orbitals revealed that the best-selected molecules were also stable to undergo chemical reactions (Table 8) (Fig. 4).
Finally, medicinal plants are potential sources of numerous phytocompounds of great therapeutic values, with many being considered in alleviating Alzheimer's disease conditions [100]. In this experiment, a total of 40 plant-derived phytochemicals were analyzed in stepwise computational molecular docking approaches to identify potential inhibitors of β-secretase. Eventually, 3 ligands (i.e., amentoflavone, bilobetin, and ellagic acid) were selected as the best inhibitors. A significant correlation among the IC 50 values, SP docking scores, XP docking scores, and free binding energies for the molecules having a molecular weight in a specific range was observed with minimal exception (Table 2).
Asp32 and Asp228 amino acid residues from the catalytic dyad inside the active site of β-secretase. Portray crucial contributions for the cleavage of the amyloid precursor protein. Again, Tyr71 residue within the active site responds to inhibitor binding by changing its conformation [101,102]. In this experiment, amentoflavone formed 1 conventional hydrogen bond with Asp228 and  3 hydrophobic interactions with Tyr71 amino acid residues (Table 3 and Fig. 4). Again, bilobetin formed 1 non-conventional hydrogen bond with Asp228 and 3 hydrophobic interactions with Tyr71 amino acid residues. Moreover, ellagic acid also formed 1 conventional hydrogen bond with Asp32 and 4 hydrophobic interactions with Tyr71 amino acid residues. Hence, these compounds are expected to interfere with the normal function of β-secretase. Later, they were predicted to perform almost similar when analyzed in different postscreening studies. Impermeability of the selected compounds to the blood-brain barrier might require further modification since this characteristic is a major concern for AD drug development (Table 5) [103]. Finally, the de novo drug discovery process involves prolonged time, multiple steps and elevated cost. Computer-aided drug designing nowadays offers great fidelity of prediction and has become increasingly popular in the last few decades. It is used in conjunction with in vitro drug discovery and greatly helps in reducing time and cost of any de novo drug discovery initiatives. In this study, we screened 40 known inhibitors (phytocompounds) of β-secretase and eventually found amentoflavone, bilobetin, and ellagic acid as the best inhibitors. The compounds were found to interact with the key amino acids in the active site of β-secretase. Moreover, the best molecules proved to be adequately harmless when they were analyzed in the different steps of drug-likeness property analysis experiment. Therefore, these molecules can be investigated further to develop potent anti-AD drugs. The other molecules of this study could also be investigated since they also performed very similarly in docking experiments. However, computational exploration is largely based on the modeling of the molecules and sometimes may generate faulty outcomes [104,105]. Therefore, further in vivo and in vitro studies might be required to strengthen the findings of this study.

Conclusion
The underlying mechanism of AD development still remains unclear. On top of that, the presence of the blood-brain barrier makes it difficult to design a drug to treat AD. As a result, not a single drug has yet been proven to treat this progressive neurological disease to the extent of any satisfactory margin. In this study, we analyzed 40 phytocompounds that showed a β-secretase inhibitory effect in laboratory studies via different computational experiments. Upon continuous exploration, we found amentoflavone, bilobetin, and ellagic acid as the most potent inhibitors of β-secretase, which could be the best possible drugs for β-secretase dependent AD treatment. These compounds also performed well in different post-screening studies unveiling the potential druggable properties. Nonetheless, amentoflavone, bilobetin, and ellagic acid could be investigated further for potential AD drug discovery. Additionally, the other selected compounds could also be analyzed further since they also performed well in different experiments of this study. This study should contribute to the development of an effective drug for AD treatment in the very near future. However, further supportive laboratory experiments and interventions might be required to support the findings of the study.