Next Article in Journal
Mechanisms of Resistance to Rituximab Used for the Treatment of Autoimmune Blistering Diseases
Previous Article in Journal
Palmitoylethanolamide (PEA) for Prevention of Gastroesophageal Inflammation: Insights from In Vitro Models
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Network Pharmacology, Molecular Dynamics and In Vitro Assessments of Indigenous Herbal Formulations for Alzheimer’s Therapy

by
Oluwafemi Adeleke Ojo
1,2,*,
Omolola Adenike Ajayi-Odoko
3,
Gideon Ampoma Gyebi
4,5,
Damilare IyinKristi Ayokunle
6,
Akingbolabo Daniel Ogunlakin
1,2,
Emmanuel Henry Ezenabor
1,2,
Adesoji Alani Olanrewaju
7,
Oluwatobi Deborah Agbeye
1,2,
Emmanuel Tope Ogunwale
8,
Damilare Emmanuel Rotimi
9,10,
Dalia Fouad
11,
Gaber El-Saber Batiha
12 and
Oluyomi Stephen Adeyemi
1,2,13
1
Good Health and Wellbeing Research Clusters (SDG 03), Bowen University, Iwo 232102, Nigeria
2
Phytomedicine, Molecular Toxicology, and Computational Biochemistry Research Group, Biochemistry Programme, Bowen University, Iwo 232102, Nigeria
3
Microbiology Programme, Bowen University, Iwo 232102, Nigeria
4
Natural Products and Structural (Bio-Chem)-Informatics Research Laboratory (NpsBC-RI), Department of Biochemistry, Bingham University, Karu 961105, Nigeria
5
Department of Biotechnology and Food Science, Faculty of Applied Sciences, Durban University of Technology, P.O. Box 1334, Durban 4000, South Africa
6
Pure and Applied Biology Programme, Bowen University, Iwo 232102, Nigeria
7
Chemistry and Industrial Chemistry Programme, Bowen University, Iwo 232102, Nigeria
8
Animal Science Programme, Bowen University, Iwo 232102, Nigeria
9
Department of Pharmacology and Pharmaceutical Sciences, Alfred E. Mann School of Pharmacy and Pharmaceutical Sciences, University of Southern California, Los Angeles, CA 90089, USA
10
Department of Biochemistry, Landmark University, Omu-Aran 251101, Nigeria
11
Department of Zoology, College of Science, King Saud University, P.O. Box 22452, Riyadh 11495, Saudi Arabia
12
Department of Pharmacology and Therapeutic, Faculty of Veterinary Medicine, Damanhour University, Damanhour 22511, AlBeheira, Egypt
13
Laboratory of Sustainable Animal Environment Systems, Graduate School of Agricultural Sciences, Tohoku University, Sendai 980-8579, Japan
*
Author to whom correspondence should be addressed.
Life 2024, 14(10), 1222; https://doi.org/10.3390/life14101222
Submission received: 14 August 2024 / Revised: 11 September 2024 / Accepted: 18 September 2024 / Published: 25 September 2024

Abstract

:
Alzheimer’s disease (AD) is an age-associated neurodegenerative condition marked by amyloid plaques, synaptic dysfunction, and neuronal loss. Besides conventional medical care, herbal therapies, both raw and refined, have attracted researchers for their potential therapeutic effects. As a proof-of-concept, our study combined HPLC-DAD analysis of bioactive constituents, network pharmacology, molecular dynamics (MD), molecular docking, post-MD analysis, and experimental verification to investigate the mechanisms of crude drug formulations as a therapeutic strategy for AD. We identified nine bioactive compounds targeting 188 proteins and 1171 AD-associated genes. Using a Venn diagram, we found 47 overlapping targets, forming “herb-compound-target (HCT)” interaction networks and a protein‒protein interaction (PPI) network. Simulations analyzed binding interactions among the three core targets and their compounds. MD assessed the stability of the best-ranked poses and beneficial compounds for each protein. Among the top 22 hub genes, AChE, BChE, and MAO, ranked 10, 14, and 34, respectively, were selected for further analysis. Two tetraherbal formulations, Form A and Form B, showed notable activity against AChE. Form A exhibited significant (p < 0.0001) inhibitory activity (IC50 = 114.842 ± 2.084 µg/mL) compared to Form B (IC50 = 142.829 ± 4.258 µg/mL), though weaker than galantamine (IC50 = 27.950 ± 0.122 µg/mL). Form B had significant inhibitory effects on BChE (IC50 = 655.860 ± 32.812 µg/mL) compared to Form A (IC50 = 679.718 ± 20.656 µg/mL), but lower than galantamine (IC50 = 23.126 ± 0.683 µg/mL). Both forms protected against Fe2+-mediated brain injury by inhibiting MAO. Docking identified quercetin (−10.2 kcal/mol) and myricetin (−10.1 kcal/mol) for AChE; rutin (−10.6 kcal/mol) and quercetin (−9.7 kcal/mol) for BChE; and kaempferol (−9.1 kcal/mol) and quercetin (−8.9 kcal/mol) for MAO. These compounds were thermodynamically stable based on MD analysis. Collectively, the results offer a scientific rationale for the use of these specifically selected medicinal herbs as AD medications.

1. Introduction

Alzheimer’s disease (AD) is an age-associated, progressive neurodegenerative condition characterized by amyloid plaques, neuronal cell death, synaptic disorders, and neurofibrillary tangles. One of the main factors contributing to AD pathogenesis is the accumulation of β-amyloid, which detrimentally affects synaptic activity, leading to neurodegeneration [1]. This condition predominantly affects elderly people, representing a significant global health concern, with dementia impacting 47 million people in 2015, a number projected to increase to 131 million by 2050. The risk of developing Alzheimer’s disease is significantly influenced by dietary practices, as various studies have highlighted the potential preventive effects of bioactive compounds from diverse food sources to potentiate drug development [2].
Current AD treatments mainly include N-methyl-D-aspartate (NMDA) receptor inhibitors and acetylcholinesterase inhibitors (AChE). However, these medicines often target a single pathway, exhibit low resistance to drug effects, and can cause adverse side effects, resulting in suboptimal clinical outcomes [3]. In light of the challenges associated with conventional treatments, there has been a growing interest in traditional herbal remedies as alternatives that address the limitations of synthetic drugs [4].
Herbal medicine offers several promising strategies for slowing AD progression and managing symptoms, primarily through mechanisms such as reducing oxidative stress, modulating inflammatory pathways, inhibiting acetylcholinesterase activity, and promoting neuroprotection. For instance, the antioxidative properties of certain herbs can mitigate the oxidative stress implicated in AD pathogenesis, while others may exert anti-inflammatory effects by downregulating pro-inflammatory cytokines. The production and sale of plant-derived medications are increasing, underscoring their increasing scientific and financial importance in healthcare [5].
This research investigated three beneficial herbal plants, Beta vulgaris (root and leaves), Persea americana (avocado) seeds, and Syzygium aromaticum (cloves), to identify novel treatment strategies for AD. Beta vulgaris is a member of the Chenopodiaceae family and is valued for its medicinal qualities and juice value [6]. It contains biologically active components such as flavonoids, ascorbic acid, saponins, nitrate, and polyphenols that contribute to its health benefits, including antioxidant, anti-inflammatory, and antihypertensive effects [7,8].
Persea americana, or avocado, belongs to the Lauraceae family and is well known worldwide because of its nutritionally rich fruit [9]. Seeds contain phytocompounds, including flavonoids, glycosides, alkaloids, phenols, and tannins, which have various medicinal uses, including for the treatment of diarrhea, dysentery, and hypertension [10].
Syzygium aromaticum, also referred to as clove, is a member of the Myrtaceae family and is indigenous to parts of India and Indonesia [11]. Cloves are recognized for their anti-inflammatory, antimicrobial, and antioxidant qualities, attributed to their bioactive constituents, including vitamin A and beta-carotene [12,13].
To advance drug discovery, the concept of network pharmacology has been introduced. This approach has been increasingly used over the past decade to understand the mechanisms of herbal formulas by examining the complex interactions between ligands, targets, and biological systems [14]. Molecular docking, a computational approach, is also used to determine the optimal binding configuration of a ligand and a molecular target. Network pharmacology aims to methodically and systematically examine the regulatory impacts of drugs on biomolecular networks, thereby advancing experimental studies conducted in vivo, ex vivo, and in vitro [15]. This study investigated the effects of distinct plant formulations combined with active bioactive compounds in combating AD through the use of a combination of network pharmacology and molecular dynamics modeling.

2. Methods

2.1. Plant Materials

The selected medicinal plants, namely, Beta vulgaris root and leaves, Persea americana seeds, and Syzygium aromaticum, were collected from different areas in Nigeria. Specifically, the plants were collected from the following locations: Beta vulgaris root and leaves from Jos, Plateau State (latitude 80°24′ N and longitude 80°32′), Persea americana seeds from Ado-Ekiti, Ekiti State (latitude 7.6124° N and longitude 5.2371° E), and Syzygium aromaticum from Omu-Aran, Kwara State (latitude 8.1402° N and longitude 5.0963° E).
The herbaria were prepared and identified with the help of a taxonomist from the Forestry Research Institute of Nigeria (FRIN), Ibadan. The identification sample of each collected medicinal plant was deposited in the herbarium department of FRIN (FHI 114105, FHI 113162, and FHI 114106). The plant samples were immediately rinsed with water and cut into small pieces. After that, the samples were completely air-dried. To determine whether the samples still contained moisture, a weight variation test was performed at various intervals. Following drying to a constant weight, the samples were ground into a powder using a grinder and sieved.

2.2. Flavonoid-Rich Extraction and Crude Drug Formulations

Fifty (50) grams of each powdered sample of B. vulgaris leaf, Persea americana seeds, Beta vulgaris root, and Syzygium aromaticum were steeped in methanol at a concentration of 80% for 72 h to obtain a crude methanolic extract. Flavonoid-rich extraction was performed following a method described elsewhere [16]. In addition, two formulations were created by combining ground crude extracts high in flavonoids at the same ratios, as indicated in Table 1. The obtained crude formulations A and B were used for further experiments.

2.2.1. Analysis of Crude Formulations A and B Using High-Performance Liquid Chromatography (HPLC–DAD)

HPLC-DAD asnalysis of crude formulations A and B was performed following methods described elsewhere [16,17].

2.2.2. Determination of the Cholinesterase Activity of Crude Formulations A and B

The inhibitory effects of crude formulations A and B, which are rich in flavonoids, on acetylcholinesterase (AChE) and butyrylcholinesterase (BChE) activities were assessed using Ellman’s method [18].

2.3. Ex Vivo Studies

2.3.1. Animals and Preparation of the Brain

Male healthy Wistar rats weighing between 150 and 200 g each were procured from the Department of Biochemistry, Bowen University, Nigeria. Brain preparation was performed following the method described by [16]. The study adhered to the approved protocols of the Bowen University Research Ethics Committee (approval number: BUI/BCH/2024/0005), which were reported using the ARRIVE guidelines.

2.3.2. Ex Vivo Induction of Brain Injury

Using Fe2+, brain injury was induced ex vivo using the techniques outlined by Erukainure et al. [19] and Ojo et al. [16].

2.3.3. Evaluation of the Activity of MAO, Also Called Monoamine Oxidase

The monoamine oxidase (MAO) activity was assessed using the method described by Green and Haughton [20].

3. Network Pharmacology and Molecular Dynamics Studies

3.1. Exploring Flavonoid-Rich Extracts of Crude Formulations A and B for Possible Target Genes for Biologically Active Compounds

The components in the flavonoid-rich extract of crude formulations A and B, which were identified as potential targets by HPLC, were predicted using SwissTarget Prediction (http://www.swisstargetprediction.ch/, retrieved on 12 July 2024) [21] while using a probability filter above zero in the “Homo sapiens” mode [22] and PharmMapper (http://www.lilab-ecust.cn/pharmmapper/, retrieved on 12 July 2024). The structures of these molecular bioactive compounds were obtained in SMILES format from the PubChem database (https://pubchem.ncbi.nlm.nih.gov/), which was accessed on 12 July 2024 [23], and uploaded to the servers. We used the species Homo sapiens to filter the findings in the SwissTarget Prediction platform, and only targets with a probability of ≥0.1 were considered. Hit target pharmacophore models, Human Protein Targets Only, and Maximum Generated Conformations 300 were the parameters utilized in PharmMapper. The models were sorted by normalized fit score, with scores below 0.5 being eliminated. Finally, the UniProt database (https://www.uniprot.org/, accessed on 12 July 2024) was used to acquire the unique matching five gene names and UniProt IDs [24].

3.2. Building the AD Target Database

The genes linked to AD were acquired by means of the DisGeNet database (https://www.disgenet.org/, retrieved on 12 July 2024) [25], the MalaCards database (https://www.malacards.org/, retrieved on 12 July 2024) [26], and the Online Mendelian Inheritance in Man (OMIM, https://www.omim.org/) database retrieved on 12 July 2024) [27]. “Alzheimer’s Disease” was used as the search term. One of the largest public databases of genes and variations linked to human diseases may be found on the discovery site DisGeNET. The inclusion of targets with a score of ≥0.1 was considered. Using the UniProt database, we succeeded in discovering the final gene list and establishing standard names after merging the three acquired targets and eliminating duplicate targets. (https://www.uniprot.org/, accessed on 12 July 2024) [24].

3.3. Identification of Bioactive Compound Targets for AD

Once all of the targets associated with AD and the active ingredients in the bioactive compounds were integrated, the overlapping targets for additional analysis were extracted using the VENNY 2.1 platform of the R language package (https://bioinfogp.cnb.csic.es/tools/venny/ accessed on 26 March 2024).

3.4. Construction of Bioactive Compounds from the Flavonoid-Rich Extracts of Crude Formulations A and B and the AD Target Network

Using the STRING database (https://string-db.org/, retrieved on 12 July 2024) [28], we constructed a protein–protein interaction network, limiting its scope within the genus Homo sapiens with a confidence level above 0.9 [28], a graphical user interface that is freely accessible and can be used to import, explore, and analyze biomolecular interaction connections visually. The nodes in the network represented the active substances and their target genes, while the edges displayed the interactions between the active substances and their target genes. We analyzed the network using the CytoHubba plug-in [29], which provides an easy-to-use interface for identifying important nodes in biological networks.

3.5. Pathway and Functional Enrichment Analysis

Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were used to evaluate the intersecting genes using the Shiny GO 0.77 tool (http://bioinformatics.sdstate.edu/go/, retrieved on 12 July 2024) [30]. The top ten results are displayed with a cutoff value of p < 0.05 and an FDR of <0.05.

3.6. Molecular Docking Analysis

Based on the results of the PIP-network interaction, important targets were selected for molecular docking analysis. Although MOA, which had a far lower score, was considered because it formed a strong connection with both AChE and BChE and has a great deal of potential as a viable AD treatment target, AChE and BChE were chosen because they possessed significantly higher values in the PPI network.

3.7. Molecular Docking Studies of Bioactive Compounds Against Target Substances

3.7.1. Preparation of Protein Structure

The three target protein 3D structures that were obtained from the Protein Data Bank (http://www.rcsb.org, accessed on 26 March 2024) are human monoamine oxide B complexed with safinamide (PDBID: 2V5Z), human acetylcholinesterase (hAChE) complexed with donepezil (PDBID: 4EY7), and human butyrylcholinesterase (hAChE) complexed with decamethonium (PDBID: 6EP4). Using MGL-AutoDockTools (ADT, v1.5.6), missing hydrogen atoms were added to all crystal structures, while the existing ligands and water molecules were eliminated [31].

3.7.2. Ligand Preparation

The HPLC-detected phytocompounds from the flavonoid-rich extract of crude formulations A and B as well as the structure data format (SDF) of the reference inhibitors (donepezil, decamethonium, and safinamide) were retrieved from the PubChem database (www.pubchem.ncbi.nlm.nih.gov, accessed on 26 March 2024). Open Babel was utilized to further transform the chemicals into the pdb chemical format [32]. The carbon atoms were combined with nonpolar hydrogen molecules, and Gasteiger-type polar hydrogen charges were assigned to the atoms. Additionally, AutoDock Tools was utilized to convert the ligand molecules into the dockable PDBQT format.

3.7.3. Molecular Docking Protocol Validation

The extracted crystallized ligand from the two proteins was superimposed with the docked poses of the cocrystallized ligands (donepezil) that had the lowest binding affinity according to the docking evaluation to verify the procedure employed for the docking evaluation. This procedure was used to validate the docking methodology used for the docking of the bioactive substances identified by HPLC. Using Discovery Studio Visualizer (BIOVIA discovery studio, 2020), the RMSD was calculated following the superimposition.

3.7.4. Molecular Docking of Phytochemicals with Targeted Active Sites

Using AutoDock Vina integrated into PyRx 0.8, active site target molecular docking of the reference inhibitors and the HPLC-DAD-discovered compounds to the binding site of the three target proteins was carried out [33]. The bioactive chemicals were imported using Open Babel and integrated into PyRx 0.8 prior to docking analysis [32]. Moreover, Open Babel further reduced the number of bioactive components. The universal force field (UFF) and the optimization procedure were utilized as the energy minimization parameter and conjugate gradient descent, respectively. Mapping the amino acid residues surrounding the natural ligand binding site allowed the identification of the binding site coordinates of the target proteins. Table 2 displays the dimensions of the grid boxes that were created around the targets’ actives. Discovery Studio Visualizer version 16 was used to perform an interactive analysis on the chosen conformer from the docking analysis.

3.7.5. Molecular Dynamics

For a 100 ns molecular dynamics simulation, the complexes of the top two phytochemicals with 4ey7 and 1b2y were also determined. The GROMACS 2019.2 and GROMOS96 43a1 force fields were used in the investigation. Charmm GUI was used to create topology files for the proteins and ligands [16,34,35]. The solvation system, periodic boundary conditions, physiological conditions, system minimization, and equilibration under constant number of atoms, pressure, and temperature (NPT) in the simulation were consistent with those employed in our previous study [16,36,37,38]. The temperature was maintained at 310 K, and the pressure was maintained at 1 atm using a Parrinello–Rahman barostat and velocity rescales. A leap-frog integrator was utilized with a 2-femtosecond time step. Every system had a 100 ns simulation, during which 1000 frames—a snapshot every 0.1 ns—were captured for every system. H-bonds, ROG, SASA, RMSD, and RMSF were obtained from the MD trajectories.

3.7.6. Calculation of the Binding Free Energy with MM-GBSA

The Molecular Mechanics Generalized Born Surface Area (MM-GBSA) method, along with decomposition analysis using the gmx MMPBSA package, was employed to obtain the binding energies of amino acids within 0.5 nm of the ligand. This analysis aimed to determine the binding free energy of the two top-docked phytochemicals from the initial docking study [39,40]. The techniques employed were the same as those described in our earlier reports [16,37,38].

3.7.7. Data Analysis

The data were analyzed using one-way ANOVA, and the mean ± standard deviation (n = 3) represents the output. We subsequently used Tukey’s multiple range post hoc test to determine statistically significant differences (p < 0.05). For graph charting, GraphPad Prism version 9.0 was used.

4. Results

4.1. HPLC-DAD Analyses of Flavonoid-Rich Extracts of Crude Formulations A and B

HPLC-DAD analysis of crude formulation A revealed the presence of six (6) bioactive flavonoids, namely, caffeic acid, syringic acid, gallic acid, kaempferol, quercetin, and rutin, while Form B contained seven bioactive flavonoids, namely, myricetin, rutin, gallic acid, caffeic acid, quercetin, p-coumaric acid, and methyl gallate (Table 3, Figures S1 and S2).

4.2. Inhibitory Action of Butyrylcholinesterase and Acetylcholinesterase

We assessed the potential of the formulations to inhibit AChE and BChE, pivotal enzymes in modulating acetylcholine levels crucial for nerve transmission. The percentage inhibition at varying concentrations is depicted in Figure 1 and Figure 2, illustrating the dose-dependent inhibitory effects on both enzymes. The flavonoid-rich extracts of crude formulations A (IC50 = 112.842 ± 2.084 µg/mL) and B (IC50 = 142.829 ± 4.258 µg/mL) exhibited notable inhibitory activity (p < 0.0001) against AChE, despite being less potent than the typical control, galantamine (IC50 = 27.950 ± 0.122 µg/mL), as shown in Figure 1. Similarly, the flavonoid-rich extracts of crude drug formulations A (IC50 = 679.718 ± 20.656 µg/mL) and B (IC50 = 655.860 ± 32.812 µg/mL) had notable (p < 0.0001) inhibitory effects on BChE, although the effects were weaker than those of galantamine (IC50 = 23.126 ± 0.683 µg/mL), as indicated in Figure 2.

4.3. Activity of Monoamine Oxidase

The inhibitory effects of flavonoid-rich crude drug formulations A and B on monoamine oxidase activity (MAO) in the oxidized brain were evaluated, and the results are shown in Figure 3. The experiment showed that the untreated rats had increased MAO activity (p < 0.0001). However, we observed a notable decrease in rat brains treated with various doses of the formulations (p < 0.0001). The inhibitory effect of the extracts rich in flavonoids from crude drug formulations A and B on enzyme activity was concentration-dependent; the effect increased as the concentration of the crude formulations decreased, with the greatest inhibitory effect recorded at 1000 µg/mL.

4.4. Screening of Bioactive Compounds and Databases for AD Targets

The screening of gallic acid, caffeic acid, syringic acid, rutin, quercetin, kaempferol, methyl gallate, p-coumaric acid, and myricetin in the SwissTarget and PharmMapper databases revealed 188 targets, while after the duplicates were removed, we obtained 1171 genes associated with AD from the DisGeNet, MalaCards, and Online Mendelian Inheritance in Man Databases. Using a Venn diagram, common AD genes and targets associated with bioactive compounds were identified. Forty-seven putative anti-AD genes were chosen and regarded as important targets. This schematic is displayed in Figure 4.

4.5. Analysis of the Target Protein‒Protein Interaction Network

Figure 5 displays the diagrammatic form of the PPI network. There were 326 edges and 47 nodes, with an average node degree of 13.9. The network was then loaded into Cytoscape 3.8.2 to identify important subnetworks. The network clustering process was carried out using proximity and degree metrics based on the scores determined by CutyoHubba and CytoNca. According to the values of genes that were greater than the median across all the results, the top 22 targets were filtered out. The top five targets are the PIK3CA, APP—amyloid-beta precursor, AKT1 GSK3B, and TLR4 proteins. Also, among the several targets that were identified, AChE, BChE, and MAO, which were ranked as 10, 14, and 34 based on the closeness parameter, respectively. These proteins have been identified to be inhibited by the formulation in the in vitro assays that were used to validate the activities of the formulations; hence, they were selected for further analysis.

4.6. Enrichment Analysis of Overlapping Targets

Gene Ontology (GO) enrichment revealed the top 22 out of the 47 overlapping targets for each enrichment. According to our biological process (BP) analysis (Figure 6A), the functions of the bioactive compounds were primarily focused on the response to compounds containing oxygen and the response to peptides containing oxygen. The majority of the genes encoded proteins in the endoplasmic reticulum lumen, cell membrane ratio and microdomains, mitochondria, and, notably, synapses (Figure 6B). The molecular function (MF) elements (Figure 6C) primarily included identical proteins, binding to enzymes, cholinesterase and nuclear estrogen receptor activity, and, notably, amyloid-beta binding. The bioactive compounds were primarily linked to endocrine resistance, lipids, cancers, and most notably, AD, according to Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis (Figure 6D).

4.7. Molecular Docking of Bioactive Compounds with Important Targets

We used molecular docking based on the PPI network to assess whether AD-related core targets and bioactive compounds could bind. The basis for the selection of target proteins stems from their enrichment in the PPI and related to the prior in vitro assay that was used to validate the activity of the formulation. Prior to the bioactive compounds being docked to the target proteins’ active sites, the process was verified, and the superimposed docked conformations of donepezil and safinamide were superimposed against the crystalized conformation. The computed RMSD for both compounds were 0.4642 and 1.0496 Å, respectively (Figure 7). The binding energies of the bioactive compounds with respect to the three target proteins are shown in Figure 7. Donepezil, decamethonium, and safinamide, the native ligands of the 4EY7, 6EP4, and 2V5Z targets, presented binding energies of −12.2, −5.4, and −9.9 kcal/mol, respectively. The two compounds with the lowest binding energies to the targets were chosen. Quercetin (−10.2 kcal/mol) and myricetin (−10.1 1 kcal/mol) are the top two compounds for 4EY7; rutin (−10.6 kcal/mol) and quercetin (−9.7 kcal/mol) are for 6EP4; and kaempferol (−9.1 kcal/mol) and quercetin (−8.9 kcal/mol) are for 2V5Z. The top two compounds of 6EP4 had lower binding energies than did the reference compound (decamethonium, −5.4 kcal/mol), but the binding energies for the 4EY7 and 2V5Z targets were rather close to those of the reference compounds. Moreover, at least two of the protein targets showed repeated high binding inclinations to quercetin, and rutin 16 (Figure 8).

4.8. Amino Acid Interactions of the Top Two Compounds from the Docking Analysis and Reference Molecules with the Five Protein Targets

The connections between the reference molecule and the top two compounds from the docking study with residues from the target protein binding site are shown in Figure 9, Figure 10 and Figure 11. A validation study revealed that donepezil was stretched in the long, narrow, and hydrophobic gorge of 4ey7 in a binding configuration similar to that of the natural ligand. There was a single hydrogen bond between Phe295 and the carbonyl oxygen of the indenone ring. Two Pi-alkyl interactions were generated by Tyr337 and Tyr341 of 4ey7 with the donepezil piperidine ring. With Donepezil’s 1-benzyl unit, Trp86 and His447 participate in an aromatic pi-pi stacking interaction. It was discovered that the 5-methoxy unit of inden-1-one and the piperidine rings Trp286 and Phe338 have two pi-sigma bonds (Figure 9).
As shown in Table 4, myricetin and quercetin, the best bioactive compounds docked to the 4EY7 target, were docked similarly to each other in the gorge’s active site, generating many hydrogen bonds and hydrophobic contacts. The interaction of myricetin in the active site was stabilized by 12 hydrogen bonds, 3 carbon hydrogen bonds, 2 pi-pi T shapes, and pi-pi stacking, while quercetin was stabilized by 12 hydrogen bonds, 3 carbon hydrogen bonds, 2 pi-pi T shapes, and pi-pi stacking. Decamethonium was stabilized in the active site of 6EP4 with major hydrophobic interactions such as pi-alkyl interactions with His438 and Trp82 and an attractive charge with Asp70 without hydrogen bonds. Quercetin and rutin, on the other hand, were stabilized by several hydrogen and hydrophobic contacts, including pi-sigma (Ala328), pi-pi stacking (Phe329), and pi-pi T-shaped contacts with Trp430 of 6EP4 (Figure 10). The reference compound for 2V5Z, safinamide, was stabilized with catalytic residues through several hydrophobic contacts, including pi-sigma (Leu171), pi-sulfur (Cys172), and pi-pi T-shaped interactions and one hydrogen bond with Gln206 of 2V5Z. In addition, the two top docked bioactive compounds produced multiple hydrophobic interactions and additional hydrogen bonds with the catalytic residues. Kaempferol forms eight hydrogen bonds with 2V5Z and hydrophobic contacts, including pi-sigma (Leu171 and Ile199), pi-pi T-shaped (Tyr326), amide-pi stacking (Ile198), and pi-alkyl (Ile316, Leu171 and Ile199) interactions. Quercetin was stabilized with four hydrogen bonds in the active site of 2V5Z and several hydrophobic contacts, including pi-sigma (Ile199 and Leu171), pi-pi T-shaped (Tyr398), and pi-pi stacking (Tyr326) (Figure 11).

4.9. Molecular Dynamics Simulations

To compare the thermodynamic stability of the top two bioactive compounds complexed with the 4EY7 protein with that of the reference chemical-bound system, the number of intramolecular hydrogen bonds as well as the RMSD, ROG, and SASA were assessed during the simulation. Table 5 displays the averages and standard deviations of all thermodynamic parameter values, while Figure 12 displays the parameter’s time-dependent spectra. The RMSD trajectory illustrates how a complex’s structure deviates over time from its unbound structure. Generally, an RMSD value of less than 3 Å is acceptable [41]. The 4EY7 complex RMSD graphs indicated equilibrium before 10 ns, and for the duration of the run, the system fluctuated very little (Figure 12a). By measuring the RMSF, the thermodynamic flexibility of amino acid residues following the binding of myricetin and quercetin to the 4EY7 active site was examined (Figure 12b).
During an MD simulation, the ROG, another thermodynamic stability metric, measures a complex’s time-dependent compactness; the lower the value is, the more stable and compact the complex is. The complex systems of the top-docked compounds presented similar mean ROG values to those of donepezil (23.20 Å, 23.19 Å, and 23.22 Å) (Figure 12c). It is commonly assumed that the binding of the top-docked compounds significantly impacts the structural integrity of the protein, which could lead to the unfolding of the protein structure. A strong positive correlation between the RMSD and ROG usually indicated a distortion of the protein structure. In this study, even though a higher RMSD value was recorded, a close mean ROG value was reported for the top docked compounds compared to the reference compound (donepezil).
The SASA is a measure of protein folding and surface area changes during a simulation; greater SASA values imply an increase in protein volume ([42]). Along with the mean SASA values of donepezil, those of the top-docked complexes were similar. The findings from the SASA analysis corroborate those of the mean ROG, indicating that the binding of the top-docked compound did not cause unfolding of the protein (Figure 12d). The average number of H-bonds created in the molecules as a whole did not vary significantly during the experiment. A close number of hydrogen bonds were observed in the ligand-bound complexes (Figure 12e).

4.10. Molecular Mechanics Generalized Born Surface Area (MM-GBSA) Analysis

The binding free energy of the two topmost phytochemicals docked to the protein 4ey7 was determined using the MMGBSA technique. Among the top-docked compounds, myricetin had the highest negative binding free energy of −25.78 ± 4.04 kcal/mol to 4EY7; this binding energy was greater than that of donepezil (−19.92 ± 3.62). Notably, the binding free energies of both phytocompounds are greater than those of donepezil. Table 6 lists the different parts that add up to the overall binding free energy. Using decomposition analysis, the contributing amino acids that make up the overall binding energy were examined and are shown in Figure 13. The majority of the binding free energy was attributed to the interacting residues during static docking.

5. Discussion

Many Alzheimer’s patients are currently interested in complementary or alternative therapies that use readily available herbal products. In the interim, a variety of elements have a direct impact on the selection criteria for herbal formulations used in illness treatment. The stage at which the disease progresses, the kinds of comorbidities that are present, the accessibility and cost of the herbs, and their safety profile are all taken into consideration when choosing herbal treatments for Alzheimer’s disease [43]. Three medicinal plants were chosen for the study based on these facts as well as the use of herbal medicines for their neuroprotective properties by different indigenous and ethnic groups in different parts of Nigeria.
The medicinal value of the majority of plants is attributable to their respective secondary metabolites. Secondary metabolites are organic substances found in plants that are produced from primary metabolites and have chemically distinct structures. These secondary metabolites primarily belong to the following classes: phenolics, terpenoids, alkaloids containing nitrogen, and substances containing sulfur [44]. Secondary metabolites found in plants are important radical scavengers. According to a study by Yu-Jie et al. [45], the two primary MD simulation categories of antioxidant phytochemicals are carotenoids and polyphenols. The five types of dietary and herbal polyphenols are as follows: tannins, coumarins, stilbenes, phenolic acid, and flavonoids [45]. Additional classifications for flavonoids include isoflavonoids, anthocyanidins, flavonols, flavones, flavanols, flavanones, and flavonols, all of which have been linked to antioxidant properties.
Acetylcholinesterase (AChE) and butyrylcholinesterase (BChE) are enzymes that degrade acetylcholine, leading to reduced neurotransmission and progressive cognitive decline. Inhibiting these enzymes can enhance cholinergic transmission, thereby alleviating Alzheimer’s disease (AD) symptoms, including memory loss, and reducing mortality risk [46]. Consequently, cholinesterase inhibitors are currently the only approved therapy for AD and neurodegenerative dementia [47]. Cholinergic dysfunction is marked by increased AChE activity, making it a viable target for therapeutic interventions [48]. Additionally, AChE inhibition not only improves cholinergic signaling but also reduces amyloid-beta peptide formation and aggregation in AD [49]. In our study, the flavonoid-rich crude formulations A and B exhibited notable inhibition of AChE and BChE, suggesting their potential neuroprotective effects on the management of AD. This inhibition implies the ability of the crude drug formulations to mitigate AD-related enzyme activity, thereby showing promise as therapeutic agents. Our findings are consistent with previous reports [16].
The inhibition of monoamine oxidase (MAO) is a crucial biomarker in the management of AD. The ability of the flavonoid-rich crude formulations A and B to inhibit this enzyme suggests potential therapeutic benefits for AD. This reduction in MAO activity, attributed to the crude drug formulations, could increase the levels of amine neurotransmitters such as dopamine and serotonin, while also preventing the degradation of amines by reactive oxygen species (ROS) [50].
Our study provides insights from the combination of molecular docking, network pharmacology, and MD to assess the possible active ingredients of flavonoid-rich crude drug formulations A and B against AD in intricate processes. The screening results revealed 1171 genes linked to AD and 188 AD targets from the DisGeNet, MalaCards, and Online Mendelian Inheritance in Man databases. Through 47 potential anti-AD gene targets, the “active ingredient-target” interaction network and Venn diagram allowed for an intuitive understanding of the precise relationship between crude drug formulations and AD. Eleven related active bioactive compounds were identified by our network analysis, all of which may play a relatively significant role in the anti-Alzheimic effect of crude drug formulations on AD.
A protein–protein interaction (PPI) network was constructed and then put into Cytoscape using the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) database. In this network, the 38 target proteins are nodes, and the connections between proteins are edges. The complex bioactive-target network reflected that the bioactive compounds in the flavonoid-rich crude drug formulations A and B corresponded to multiple bioactive constituents, thus revealing multicomponent properties. The target proteins with high confidence, defined as those with a score higher than 0.9, were chosen. The highest correlation suggested a strong correlation between the targeted genes, suggesting that each of these genes may be a significant target [51].
The intersection targets were subjected to GO and KEGG pathway enrichment analysis using the Shiny GO 0.77 tool. The 47 overlapping targets of the top 22 enriched targets were identified. Among the GO-enriched categories, peptides and their reactions to oxygen-containing chemicals were the primary subjects of BP enrichment. The majority of genes encode proteins in the rat cell membrane, microdomains, endoplasmic reticulum lumen, mitochondria, and, notably, synapses. Additionally, among other factors, the enhanced items in the molecular function (MF) category were primarily linked to the binding of cholinesterase, amyloid-beta, and protein kinase.
The intersecting proteins involved in the extracellular environment, such as membrane rafts, membrane microdomains, microtubule cytoskeletons, and synapses, were classified as enriched in cellular components (CCs).
Overall, our GO enrichment results showed that targets in the cytoplasm, cell membrane, or extracellular space that are closely related to the regulation of amyloid beta, the cellular response, and metabolism during the onset and progression of AD may bind to anti-AD targets of crude drug formulations. KEGG enrichment analysis was then carried out to identify possible signaling pathways involved in the effects of crude medication formulations on AD by merging genomes with cellular and species data. According to our research, endocrine resistance pathways in cancer and AD signaling pathways are the primary targets shared by AD and crude drug formulations.
We used molecular docking based on the PPI network to assess whether AD-related core targets and bioactive drugs could bind. Prior to the bioactive chemicals being docked to the target proteins’ active sites, the process was verified, and the docked conformations of donepezil and safinamide were superimposed on the crystallized conformation. For AChE (4EY7), quercetin and myricetin were the top two compounds; for BChE (6EP4), rutin and quercetin; and for monoamine oxidase A (2V5Z), kaempferol and quercetin were the top two compounds. The top two compounds of 6EP4 had higher binding energies than did the reference compound (decamethonium, −5.4 kcal/mol), but their binding energies for the 4EY7 and 2V5Z targets were rather close to those of the reference compounds.
A validation study revealed that donepezil was stretched in the long, narrow, hydrophobic gorge of 4ey7 in a binding configuration akin to that of the natural ligand. There was one hydrogen bond between Phe295 and the carbonyl oxygen of the indenone ring. Two Pi-alkyl interactions were generated by Tyr337 and Tyr341 of 4ey7 with the donepezil piperidine ring. With donepezil’s 1-benzyl unit, Trp86 and His447 participate in an aromatic pi-pi stacking interaction.
The mean RMSD values of the complexes of quercetin and myricetin fluctuated greatly, although they were still quite close to those of the reference compound donepezil. It can be deduced from the RMSD analysis of the trajectories that the binding of myricetin and quercetin did not generally affect the thermostructural stability of the protein [52]. The ability of atoms and residues in a protein structure to form stable intra- and intermolecular bonds is taken into account by the RMSF value, which measures the average fluctuation of these elements over the course of a simulation. The stronger the bonds are and the greater the affinity of the ligand for the protein is, the less fluctuation there is, especially in the active site where ligand binding and catalysis occur [53]. Compared to that of the reference inhibitor, the binding of the top-docked ligands in the present study did not result in a change in the protein residues. The greatest variation was observed for amino acid residues 243 and 375. Terminal fluctuations are the cause of the fluctuations at the end of the spectrum. In this study, the top-docked substances showed similar mean ROG values compared to the reference compound (donepezil), despite the greater RMSD value being recorded. During a simulation, the SASA measures changes in surface area and protein folding; higher SASA values indicate higher protein volumes [42]. Compared to the mean SASA values of donepezil, those of the top-docked complexes were similar. The results of the SASA analysis support the mean ROG results, suggesting that protein unfolding was not triggered by the binding of the top-docked molecule. The stability of a protein structure is mostly dependent on intramolecular hydrogen bonds and distance, which can be assessed to ascertain how ligand binding impacts a protein’s stability during simulation [41].
The binding free energy measures the energy differential between the ligand and receptor that are bound and unbound in a complex; the larger the negative value is, the greater the ligand’s affinity for the protein [41]. The binding free energy estimates offer comprehensive data on the binding mechanisms of the top-docked compounds that are critical to the preliminary stages of drug discovery and development [54].

6. Conclusions

This study provides compelling evidence that flavonoid-rich crude drug formulations have significant potential as therapeutic agents in the treatment of Alzheimer’s disease (AD). By integrating molecular docking, network pharmacology, and molecular dynamics (MD) simulations, we identified key bioactive compounds from the formulations that target crucial AD-related proteins. The inhibitory effects of these formulations on acetylcholinesterase (AChE), butyrylcholinesterase (BChE), and monoamine oxidase (MAO) highlight their potential in modulating cholinergic transmission and reducing oxidative stress, which are central to AD pathogenesis. The molecular interactions observed in docking studies and the stability of the compounds during MD simulations further reinforce the therapeutic promise of these formulations. These findings offer a scientific rationale for the use of selected medicinal herbs in managing AD and underscore the importance of exploring natural compounds as viable alternatives to conventional AD treatments. The top enriched targets that were not considered in the presented studies are suggested for further study.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/life14101222/s1, Figure S1: Chromatogram of crude formulation A; Figure S2: Chromatogram of crude formulation B.

Author Contributions

The study was conceptualized and designed by the O.A.O. The first draft was written by O.A.O., E.H.E., G.A.G. and O.D.A., O.A.O., O.A.A.-O., D.I.A., A.A.O., A.D.O., E.T.O., D.F., G.E.-S.B. and O.S.A. conducted the experiments. O.A.O. and G.A.G. looked at and analyzed the information. O.S.A. oversaw the experiment. O.A.O., E.H.E., G.A.G., O.D.A., O.A.A.-O., D.I.A., A.A.O., A.D.O., E.T.O., D.E.R., D.F., G.E.-S.B. and O.S.A. reviewed the final draft of the manuscript. The manuscript’s final form was approved by all the authors. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Bowen University via the Bowen University Research Grant (BURG) (Grant No: BURG/2024/006) for 2023–2024. In addition, the research was also funded by King Saud University (Riyadh, Saudi Arabia) through Researchers supporting project number RSPD2024R965.

Institutional Review Board Statement

The study adhered to the approved protocols of the Bowen University Research Ethics Committee (approval number: BUI/BCH/2024/0005), which were reported using the ARRIVE guidelines.

Informed Consent Statement

Not applicable.

Data Availability Statement

The datasets generated during and/or analyzed during the current study are available from the corresponding author upon reasonable request.

Acknowledgments

The authors would like to thank Bowen University for funding this research via the Bowen University Research Grant (BURG) (Grant No: BURG/2024/006) for 2023–2024. In addition, the authors would like to extend their gratitude to King Saud University (Riyadh, Saudi Arabia) for funding this research through Researchers supporting project number RSPD2024R965.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Raafat, K. Identification of phytochemicals from North African plants for treating Alzheimer’s diseases and of their molecular targets by in silico network pharmacology approach. J. Tradit. Complement. Med. 2020, 11, 268–278. [Google Scholar] [CrossRef] [PubMed]
  2. Bordoloi, S.; Pathak, K.; Devi, M.; Saikia, R.; Das, J.; Kashyap, V.; Das, D.; Ahmad, M.; Abdel-Wahab, B. Some promising medicinal plants used in Alzheimer’s disease: An ethno pharmacological perspective. Discov. Appl. Sci. 2024, 6, 215. [Google Scholar] [CrossRef]
  3. Tao, P.; Xu, W.; Gu, S.; Shi, H.; Wang, Q.; Xu, Y. Traditional Chinese medicine promotes the control and treatment of dementia. Front. Pharmacol. 2022, 13, 1015966. [Google Scholar] [CrossRef] [PubMed]
  4. Park, J.; Kim, H.; Yang, M. Network pharmacology based study on the efficacy and mechanism of Lonicera japonica Thunberg. Appl. Sci. 2022, 12, 9122. [Google Scholar] [CrossRef]
  5. Rao, V.; Descamps, O.; John, V.; Bredesen, E. Ayurvedic medicinal plants for Alzheimer disease: A review. Alzheimer’s Res. Ther. 2012, 4, 22. [Google Scholar] [CrossRef]
  6. Ndour, P.; Fall, J.; Loum, A.; Sagne, M.; Diagne, N.; Samb, S.; Ndong, D.; Diouf, M. Effects of Kigelia africana, Beta vulgaris and Ricinodendron heudelotti as feed additive on growth performance, survival and whole-body composition of Nile tilapia (Oreochromis niloticus) Linnaeus, 1758) fry. Int. J. Fish. Aquat. Stud. 2021, 9, 398–402. [Google Scholar] [CrossRef]
  7. Patel, M.; Patel, P.; Pithawala, E. Antioxidant and Antibacterial Properties of biopigment from Beta vulgaris. Biosci. Biotechnol. Res. Asia 2023, 20, 1327–1339. [Google Scholar] [CrossRef]
  8. Afiomah, S.; Iwuozor, O. Nutritional and phytochemical properties of Beta vulgaris Linnaeus (Chenopodiaceae): A review. Niger. J. Pharm. Appl. Sci. Res. 2020, 9, 38–44. [Google Scholar]
  9. Olasunkanmi, M.; Ogunyemi, O. Phytochemical constituents and antioxidant activity of Persea americana leave. Int. J. Chem. Res. 2023, 7, 1–4. [Google Scholar] [CrossRef]
  10. Egbuonu, C.; Opara, C.; Onyeabo, C.; Uchenna, O. Proximate, functional, anti-nutrient and antimicrobial properties of avocado pear (Persea americana) seeds. J. Nutr. Health Food Eng. 2018, 1, 00260. [Google Scholar]
  11. Musbau, S.; Yushau, M.; Dalha, T. Antibacterial activities and quantitative phytochemical screening of the hexane and ethanolic oil extracts of Syzygium aromaticum against ESBL producing bacteria isolates. Eur. J. Med. Res. 2024, 2, 123–127. [Google Scholar] [CrossRef] [PubMed]
  12. Khaled, S. Comparison the effect of plant extracts and antibiotics on gram-negative bacteriaisolated from patients with urinary tract infection in Tikrit city. GSC Biol. Pharm. Sci. 2024, 27, 60–67. [Google Scholar] [CrossRef]
  13. Bhowmik, D.; Kumar, S.; Yadav, A.; Srivastava, S.; Paswan, S.; Dutta, S. Recent trends in Indian traditional herbs Syzygium aromaticum and its health benefits. J. Pharmacogn. Phytochem. 2012, 1, 13–21. [Google Scholar]
  14. Sachdeo, R.; Khanwelkar, C.; Shete, A. In silico exploration of berberine as a potential wound healing agent via network pharmacology, molecular docking, and molecular dynamics simulation. Int. J. Appl. Pharm. 2024, 16, 188–194. [Google Scholar] [CrossRef]
  15. Shi, L.; Deng, Y.; Luo, D.; Li, L.; Kuang, X.; Qi, A.; Fu, B. Exploration of the possible mechanisms of Ling Gui Zhu Gan decoction innephrotic syndrome based on network pharmacology, molecular docking and molecular dynamics simulation. Medicine 2023, 102, e34446. [Google Scholar] [CrossRef]
  16. Ojo, O.A.; Gyebi, G.A.; Ezenabor, E.H.; Iyobhebhe, M.; Emmanuel, D.A.; Adelowo, O.A.; Olujinmi, F.E.; Ogunwale, T.E.; Babatunde, D.E.; Ogunlakin, A.D.; et al. Exploring beetroot (Beta vulgaris L.) for diabetes mellitus and Alzheimer’s disease dual therapy: In vitro and computational studies. RSC Adv. 2024, 14, 19362–19380. [Google Scholar] [CrossRef]
  17. Araujo-León, J.A.; Cantillo-Ciau, Z.; Ruiz-Ciau, D.V.; Coral-Martínez, T.I. HPLC profile simultaneous quantitative analysis of tingenone pristimerin in four Celastraceae species using, HPLC-UV-DAD-MS. Rev. Bras. De Farmacogn. 2019, 29, 171–176. [Google Scholar] [CrossRef]
  18. Perry, N.S.; Houghton, P.J.; Theobald, A.; Jenner, P.; Perry, E.K. In-Vitro Inhibition of Human Erythrocyte Acetylcholinesterase by Salvia lavandula Efolia Essential Oil and Constituent Terpenes. J. Pharm. 2000, 52, 895–902. [Google Scholar] [CrossRef]
  19. Erukainure, O.L.; Chukwuma, C.I.; Matsabisa, M.G.; Salau, V.F.; Koorbanally, N.A.; Islam, M.S. Buddleja saligna Willd (Loganiaceae) inhibits angiotensin-converting enzyme activity in oxidative cardiopathy with concomitant modulation of nucleotide hydrolyzing enzymatic activities and dysregulated lipid metabolic pathways. J. Ethnopharmacol. 2020, 248, 112358. [Google Scholar] [CrossRef]
  20. Green, A.L.; Haughton TM, A. Colorimetric Method for the Estimation of Monoamine Oxidase. Biochem. J. 1961, 78, 172–176. [Google Scholar] [CrossRef]
  21. Daina, A.; Michielin, O.; Zoete, V. SwissTargetPrediction: Updated data and new features for efficient prediction of protein targets of small molecules. Nucleic Acids Res. 2019, 47, W357–W364. [Google Scholar] [CrossRef] [PubMed]
  22. Gfeller, D.; Grosdidier, A.; Wirth, M.; Daina, A.; Michielin, O.; Zoete, V. SwissTargetPrediction: A web server for targe prediction of bioactive small molecules. Nucleic Acids Res. 2014, 42, W32–W38. [Google Scholar] [CrossRef]
  23. Kim, S.; Chen, J.; Cheng, T.; Gindulyte, A.; He, J.; He, S.; Li, Q.; Shoemaker, B.A.; Thiessen, P.A.; Yu, B. PubChem 2023 update. Nucleic Acids Res. 2023, 51, D1373–D1380. [Google Scholar] [CrossRef] [PubMed]
  24. The UniProt Consortium. UniProt UniProt: The universal protein knowledgebase in 2023. Nucleic Acids Res. 2023, 51, D523–D531. [Google Scholar] [CrossRef] [PubMed]
  25. Piñero, J.; Ramírez-Anguita, J.M.; Saüch-Pitarch, J.; Ronzano, F.; Centeno, E.; Sanz, F.; Furlong, L.I. The DisGeNET knowledge platform for disease genomics: 2019 update. Nucleic Acids Res. 2020, 48, D845–D855. [Google Scholar] [CrossRef]
  26. Rappaport, N.; Nativ, N.; Stelzer, G.; Twik, M.; Guan-Golan, Y.; Iny Stein, T.; Bahir, I.; Belinky, F.; Morrey, C.P.; Safran, M. MalaCards: An integrated compendium for diseases and their annotation. Database 2013, 2013, bat018. [Google Scholar] [CrossRef]
  27. Amberger, J.S.; Bocchini, C.A.; Schiettecatte, F.; Scott, A.F.; Hamosh, A. OMIM.org: Online Mendelian Inheritance in Man (OMIM®), an online catalog of human genes and genetic disorders. Nucleic Acids Res. 2015, 43, D789–D798. [Google Scholar] [CrossRef]
  28. Szklarczyk, D.; Kirsch, R.; Koutrouli, M.; Nastou, K.; Mehryary, F.; Hachilif, R.; Gable, A.L.; Fang, T.; Doncheva, N.T.; Pyysalo, S. The STRING database in 2023: Protein–protein association networks and functional enrichment analyses for any sequenced genome of interest. Nucleic Acids Res. 2023, 51, D638–D646. [Google Scholar] [CrossRef]
  29. Shannon, P.; Markiel, A.; Ozier, O.; Baliga, N.S.; Wang, J.T.; Ramage, D.; Amin, N.; Schwikowski, B.; Ideker, T. Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genome Res. 2003, 13, 2498–2504. [Google Scholar] [CrossRef]
  30. Ge, S.X.; Jung, D.; Yao, R. ShinyGO: A graphical gene-set enrichment tool for animals and plants. Bioinformatics 2020, 36, 2628–2629. [Google Scholar] [CrossRef]
  31. Morris, G.M.; Huey, R.; Lindstrom, W.; Sanner, M.F.; Belew, R.K.; Goodsell, D.S.; Olson, A.J. AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibility. J. Comput. Chem. 2009, 30, 2785–2791. [Google Scholar] [CrossRef] [PubMed]
  32. O’Boyle, N.; Banck, M.; James, C.A.; Morley, C. Vandermeersch T. Hutchison GR Open Babel: An open chemical toolbox. J. Cheminf. 2011, 3, 33. [Google Scholar] [CrossRef] [PubMed]
  33. Trott, O.; Olson, A.J. AutoDock Vina: Improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. J. Comput. Chem. 2010, 31, 455–461. [Google Scholar] [CrossRef]
  34. Lee, J.; Cheng, X.; Swails, J.M.; Yeom, M.S.; Eastman, P.K.; Lemkul, J.A.; Wei, S.; Buckner, J.; Jeong, J.C.; Qi, Y. CHARMM-GUI input generator for NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM simulations using the CHARMM36 additive force field. J. Chem. Theory Comput. 2016, 12, 405–413. [Google Scholar] [CrossRef]
  35. Lee, J.; Hitzenberger, M.; Rieger, M.; Kern, N.R.; Zacharias, M.; Im, W. CHARMM-GUI supports the Amber force fields. J. Chem. Phys. 2020, 153, 035103. [Google Scholar] [CrossRef] [PubMed]
  36. Ogunyemi, O.M.; Gyebi, G.A.; Ibrahim, I.M.; Esan, A.M.; Olaiya, C.O.; Soliman, M.M.; Batiha, G.E.-S. Identification of promising multitargeting inhibitors of obesity from Vernonia amygdalina through computational analysis. Mol. Divers. 2023, 27, 1–25. [Google Scholar] [CrossRef]
  37. Ogunyemi, O.M.; Gyebi, G.A.; Ibrahim, I.M.; Olaiya, C.O.; Ocheje, J.O.; Fabusiwa, M.M.; Adebayo, J.O. Dietary stigmastane-type saponins as promising dual-target directed inhibitors of SARS-CoV-2 proteases: A structure-based screening. As RSC Adv. 2021, 11, 33380–33398. [Google Scholar] [CrossRef]
  38. Gyebi, G.A.; Ogunyemi, O.M.; Ibrahim, I.M.; Afolabi, S.O.; Adebayo, J.O. Dual targeting of cytokine storm and viral replication in COVID-19 by plant-derived steroidal pregnanes: An in silico perspective. Comput. Biol. Med. 2021, 134, 104406. [Google Scholar] [CrossRef]
  39. Miller, B.R., III; McGee, T.D., Jr.; Swails, J.M.; Homeyer, N.; Gohlke, H.; Roitberg, A.E. MMPBSA. py: An efficient program for end-state free energy calculations. J. Chem. Theory Comput. 2012, 8, 3314–3321. [Google Scholar] [CrossRef]
  40. Valdés-Tresanco, M.S.; Valdés-Tresanco, M.E.; Valiente, P.A.; Moreno, E. gmx_MMPBSA: A new tool to perform end-state free energy calculations with GROMACS. J. Chem. Theory Comput. 2021, 17, 6281–6291. [Google Scholar] [CrossRef]
  41. Du, X.; Li, Y.; Xia, Y.-L.; Ai, S.-M.; Liang, J.; Sang, P.; Ji, X.-L.; Liu, S.-Q. Insights into protein–ligand interactions: Mechnisms, models, and methods. Int. J. Mol. Sci. 2016, 17, 144. [Google Scholar] [CrossRef] [PubMed]
  42. Leimkuhler, B.; Matthews, C. Molecular Dynamics; Interdisciplinary application mathematics series; Springer: Cham, Switzerland, 2015; Volume 39. [Google Scholar]
  43. Alzheimer’s Association. Revised Criteria for Diagnosis and Staging of Alzheimer’s Disease: Alzheimer’s Association Workgroup. J. Alzheimer’s Assoc. 2024, 20, 5143–5169. [Google Scholar] [CrossRef]
  44. Choudhury, H.; Pandey, M.; Hua, C.K.; Mun, C.S.; Jing, J.K.; Kong, L.; Ern, L.Y.; Ashraf, N.A.; Kit, S.W.; Yee, T.S.; et al. An update on natural compounds in the remedy of diabetes mellitus: A systematic review. J. Tradit. Complement. Med. 2018, 8, 361–376. [Google Scholar] [CrossRef] [PubMed]
  45. Zhang, Y.J.; Gan, R.Y.; Li, S.; Zhou, Y.; Li, A.N.; Xu, D.P.; Li, H.B. Antioxidant Phytochemicals for the Prevention and Treatment of Chronic Diseases. Molecules 2015, 20, 21138–21156. [Google Scholar] [CrossRef]
  46. Lodes, D.E.; Zhu, J.; Tsai, N.-P. E3 ubiquitin ligase Nedd4-2 exerts neuroprotective effects during endoplasmic reticulum stress. J. Neurochem. 2022, 160, 613–624. [Google Scholar] [CrossRef]
  47. Peng, Y.; Jin, H.; Xue, Y.H.; Chen, Q.; Yao, S.Y.; Du, M.Q.; Liu, S. Current and future therapeutic strategies for Alzheimer’s disease: An overview of drug development bottlenecks. Front. Aging Neurosci. 2023, 15, 1206572. [Google Scholar] [CrossRef]
  48. Lane, R.M.; Potkin, S.G.; Enz, A. Targeting acetylcholinesterase and butyrylcholinesterase in dementia. Int. J. Neuropsychopharmacol. 2006, 9, 101–124. [Google Scholar] [CrossRef]
  49. García-Ayllón, M.S.; Small, D.H.; Avila, J.; Sáez-Valero, J. Revisiting the Role of Acetylcholinesterase in Alzheimer’s Disease: Cross-Talk with P-tau and β-Amyloid. Front. Mol. Neurosci. 2011, 4, 22. [Google Scholar] [CrossRef]
  50. Chaurasiya, N.D.; Leon, F.; Muhammad, I.; Tekwani, B.L. Natural products inhibitors of monoamine oxidases—Potential new drug leads for neuroprotection, neurological disorders, and neuroblastoma. Molecules 2022, 27, 4297. [Google Scholar] [CrossRef]
  51. Martiz, R.M.; Patil, S.M.; Abdulaziz, M.; Babalghith, A.; Al-Areefi, M.; Al-Ghorbani, M.; Mallappa Kumar, J.; Prasad, A.; Mysore Nagalingaswamy, N.P.; Ramu, R. Defining the role of isoeugenol from Ocimum tenuiflorum against diabetes mellitus-linked Alzheimer’s disease through network pharmacology and computational methods. Molecules 2022, 27, 2398. [Google Scholar] [CrossRef]
  52. Gyebi, G.A.; Ogunyemi, O.M.; Adefolalu, A.A.; Rodríguez-Martinez, A.; López-Pastor, J.F.; Banegas-Luna, A.J.; Pérez-Sánchez, H.; Adegunloye, A.P.; Ogunro, O.B.; Afolabi, S.O. African-derived phytocompounds may interfere with SARS-CoV-2 RNA capping machinery via inhibition of 2′-O-ribose methyltransferase: An in silico perspective. J. Mol. Struct. 2022, 1262, 133019. [Google Scholar] [CrossRef] [PubMed]
  53. Gyebi, A.G.; Ogunyemi, O.M.; Ibrahim, I.M.; Afolabi, S.O.; Ojo, R.J.; Ejike, U.D.I.; Adebayo, J.O. Inhibitory potential of phytocompounds from Ocimum gratissimum against antiapoptotic BCL-2 proteins associated with cancer: An integrated computational study. Egypt. J. Basic Appl. Sci. 2022, 9, 588–608. [Google Scholar]
  54. Kollman, P.A.; Massova, I.; Reyes, C.; Kuhn, B.; Huo, S.; Chong, L.; Lee, M.; Lee, T.; Duan, Y.; Wang, W. The structures and free energies of complex molecules are calculated as follows: Combining molecular mechanics and continuum models. Acc. Chem. Res. 2000, 33, 889–897. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Effect of flavonoid-rich crude drug formulations (A,B) on acetylcholinesterase activity Legend: The mean ± SD (n = 3) is used to describe the data; t tests show that p < 0.0001. Standard medication: galantamine. ** is significant at p < 0.01 when Form A is compared to form B; **** is significant at p < 0.0001 when Form A and B is compared to galantamine.
Figure 1. Effect of flavonoid-rich crude drug formulations (A,B) on acetylcholinesterase activity Legend: The mean ± SD (n = 3) is used to describe the data; t tests show that p < 0.0001. Standard medication: galantamine. ** is significant at p < 0.01 when Form A is compared to form B; **** is significant at p < 0.0001 when Form A and B is compared to galantamine.
Life 14 01222 g001
Figure 2. Effect of flavonoid-rich crude drug formulations (A,B) on butyrylcholinesterase activity Legend: The mean ± SD (n = 3) is used to describe the data; t tests show that p < 0.001. Standard medication: galantamine. *** is significant at p < 0.001 when Form A and B is compared to galantamine.
Figure 2. Effect of flavonoid-rich crude drug formulations (A,B) on butyrylcholinesterase activity Legend: The mean ± SD (n = 3) is used to describe the data; t tests show that p < 0.001. Standard medication: galantamine. *** is significant at p < 0.001 when Form A and B is compared to galantamine.
Life 14 01222 g002
Figure 3. Effect of different concentrations of flavonoid-rich crude drug formulations A and B on the activity of rat brain monoamine oxidase (MAO) ex vivo. There was a significant decrease in monoamine oxidase activity in the range of extracts rich in flavonoids from crude drug formulations A and B when compared to the group induced exclusively with FeSO4. Legend: the mean ± SD (n = 3) is used to describe the data.
Figure 3. Effect of different concentrations of flavonoid-rich crude drug formulations A and B on the activity of rat brain monoamine oxidase (MAO) ex vivo. There was a significant decrease in monoamine oxidase activity in the range of extracts rich in flavonoids from crude drug formulations A and B when compared to the group induced exclusively with FeSO4. Legend: the mean ± SD (n = 3) is used to describe the data.
Life 14 01222 g003
Figure 4. Venn diagram showing the genes linked to AD and the bioactive chemicals from flavonoid-rich crude drug formulation A- and B-associated targets. Forty-seven genes related to both Alzheimer’s disease (AD) and C. sativa were identified.
Figure 4. Venn diagram showing the genes linked to AD and the bioactive chemicals from flavonoid-rich crude drug formulation A- and B-associated targets. Forty-seven genes related to both Alzheimer’s disease (AD) and C. sativa were identified.
Life 14 01222 g004
Figure 5. PPI network of target and AD-related genes created in Cytoscape utilizing (A,B) the STRING database and (B) the core key subnetwork of the top 22 nodes examined by CytoHubba.
Figure 5. PPI network of target and AD-related genes created in Cytoscape utilizing (A,B) the STRING database and (B) the core key subnetwork of the top 22 nodes examined by CytoHubba.
Life 14 01222 g005
Figure 6. Representation of functional annotation and enriched pathways. (A) GO biological processes. (B) Cellular components of GO. (C) Molecular function of GO. (D) KEGG pathway analysis.
Figure 6. Representation of functional annotation and enriched pathways. (A) GO biological processes. (B) Cellular components of GO. (C) Molecular function of GO. (D) KEGG pathway analysis.
Life 14 01222 g006
Figure 7. Superimposition of selected conformers from the docking analysis of the native ligand on the extracted conformation of (a) donepezil and (b) safinamide. The blue compounds represent the docked conformer, while the red conformation represents the conformation extracted from 4EY7 and 2V5Z.
Figure 7. Superimposition of selected conformers from the docking analysis of the native ligand on the extracted conformation of (a) donepezil and (b) safinamide. The blue compounds represent the docked conformer, while the red conformation represents the conformation extracted from 4EY7 and 2V5Z.
Life 14 01222 g007
Figure 8. Binding energies of bioactive compounds from flavonoid-rich extracts of crude drug formulations A and B that were identified by HPLC against target proteins.
Figure 8. Binding energies of bioactive compounds from flavonoid-rich extracts of crude drug formulations A and B that were identified by HPLC against target proteins.
Life 14 01222 g008
Figure 9. Interactions of top-ranked compounds and the reference inhibitor (donepezil) with amino acids in the 4EY7 active site. Sticks are used to represent the ligands. (ia,ib,ic) Three dimensional interactions, (ii) two dimensional interactions, and (S) cartoon representation showing ligands in the active site (a) Donepezil, (b) myricetin, and (c) quercetin.
Figure 9. Interactions of top-ranked compounds and the reference inhibitor (donepezil) with amino acids in the 4EY7 active site. Sticks are used to represent the ligands. (ia,ib,ic) Three dimensional interactions, (ii) two dimensional interactions, and (S) cartoon representation showing ligands in the active site (a) Donepezil, (b) myricetin, and (c) quercetin.
Life 14 01222 g009
Figure 10. Interactions of top-ranked compounds and reference inhibitor (decamethonium) with amino acids in the 6EP4 active site. Sticks are used to represent the ligands. (ia,ib,ic) Three dimensional interactions, (ii) two dimensional interactions, and (S) cartoon representation showing ligands in the active site (a) Decamethonium, (b) rutin, and (c) quercetin.
Figure 10. Interactions of top-ranked compounds and reference inhibitor (decamethonium) with amino acids in the 6EP4 active site. Sticks are used to represent the ligands. (ia,ib,ic) Three dimensional interactions, (ii) two dimensional interactions, and (S) cartoon representation showing ligands in the active site (a) Decamethonium, (b) rutin, and (c) quercetin.
Life 14 01222 g010
Figure 11. Interactions of top-ranked compounds and reference inhibitor (Safinamide) with amino acids in the 2V5Z active site. Sticks are used to represent the ligands. (ia,ib,ic) Three dimensional interactions, (ii) two dimensional interactions, and (S) cartoon representation showing ligands in the active site. (a) Safinamide, (b) kaempferol, and (c) quercetin.
Figure 11. Interactions of top-ranked compounds and reference inhibitor (Safinamide) with amino acids in the 2V5Z active site. Sticks are used to represent the ligands. (ia,ib,ic) Three dimensional interactions, (ii) two dimensional interactions, and (S) cartoon representation showing ligands in the active site. (a) Safinamide, (b) kaempferol, and (c) quercetin.
Life 14 01222 g011
Figure 12. Plots of thermodynamic parameters computed from the analysis of the MD trajectories of 4EY7 complex systems: (a) Backbone-Root Mean Square Deviation (RMSD), (b) per residue Root Mean Square Fluctuations (RMSF), (c) radius of gyration, (d) Surface Accessible Surface Area (SASA), and (e) number of hydrogen atoms.
Figure 12. Plots of thermodynamic parameters computed from the analysis of the MD trajectories of 4EY7 complex systems: (a) Backbone-Root Mean Square Deviation (RMSD), (b) per residue Root Mean Square Fluctuations (RMSF), (c) radius of gyration, (d) Surface Accessible Surface Area (SASA), and (e) number of hydrogen atoms.
Life 14 01222 g012
Figure 13. Plot of the Molecular Mechanics Generalized Born Surface Area (MM-GBSA) showing the amino acid residues of hAChE that contribute to the total binding free energy of (a) donepezil, (b) quercetin, and (c) myricetin.
Figure 13. Plot of the Molecular Mechanics Generalized Born Surface Area (MM-GBSA) showing the amino acid residues of hAChE that contribute to the total binding free energy of (a) donepezil, (b) quercetin, and (c) myricetin.
Life 14 01222 g013aLife 14 01222 g013b
Table 1. Ratio of crude extracts high in flavonoids utilized during formulation.
Table 1. Ratio of crude extracts high in flavonoids utilized during formulation.
Ratio of Crude Extracts High in Flavonoids
PlantsForm AForm B
B. vulgaris leaf2.5*
Persea americana seeds*2.5
Beta vulgaris root*2.5
Syzygium aromaticum2.5*
* Not part of formulation.
Table 2. The target enzyme binding site coordinates.
Table 2. The target enzyme binding site coordinates.
Dimensions6ep4 (Å)2v5z (Å)4ey7 (Å)
center_x−19.8351.92−13.35
center_y−41.39155.81−43.63
center_z47.1028.6727.41
Size x15.8113.9619.75
Size y17.4214.5313.23
Size z14.5914.5916.84
Table 3. Bioactive compounds identified in crude drug formulations.
Table 3. Bioactive compounds identified in crude drug formulations.
Compounds
Form AConc (mg/mL)Form BConc (mg/mL)
Gallic acid6.68Myricetin4.63
Caffeic acid8.92p-coumaric acid5.24
Syringic acid2.18Gallic acid6.72
Rutin4.14Caffeic acid8.26
Kaempferol2.47Quercetin1.67
Quercetin1.18Methyl gallate6.26
Rutin5.48
Table 4. Interaction of the phytochemicals selected from the docking analysis with the target proteins.
Table 4. Interaction of the phytochemicals selected from the docking analysis with the target proteins.
CompoundsProtein
Targets
Hydrogen BondsHydrophobic Interaction
Interacting ResiduesInteracting Residues
Donepezil4EY7Phe295Trp86 Phe338 His447 Tyr337 Tyr341 Trp286
MyricetinTyr133 Ser203 Glu202 Gly126 Trp86 Ser125 Gln71 Tyr72 Asp74 Tyr341 Tyr337Trp86 Tyr124
QuercetinGlu202 His447 Tyr341 Tyr337 Ser125 Asp74 Tyr72 Gln71 Trp86 Gly126 Tyr133Tyr124 Trp86
Decamethonium6EP4 His438 Trp82 Asp70 yr332 Thr120 Asn83
RutinTrp82 Tyr128 Thr122 Glu197 Ser198 Ser287 Thr120 Gly116 Gln119 Ala328 Asn68 His438 Asn83 Trp430 Gly78 Asp70Phe329 Trp430 Ala328
QuercetinTyr440 Trp430 Gly78 Trp82 His438 Tyr128 Glu197 Thr120 Gln119Trp82 Gly115
Safinamide2V5ZGln206Cys172 Ile316 Ile199 Tyr326 Tyr398 Leu171
KaempferolPro120 Ile199 Ty326 Gln206 Cys172 Ile198 Tyr4335 Leu164Ile316 Leu171 Tyr326 Ile198 Ile199 Cys172
QuercetinTyr435 Cys172 Pro102 Tyr326Ile199 Leu171 Cys172 Tyr326 Tyry398
Table 5. The means and standard deviations of several metrics that were examined from the MDS trajectories of the top-docked in complex with their corresponding targets.
Table 5. The means and standard deviations of several metrics that were examined from the MDS trajectories of the top-docked in complex with their corresponding targets.
ComplexesThermodynamic Parameters
RMSD (Å)RMSF (Å)SASA (Å2)RoG (Å)H-Bonds
4EY7_Donepezil1.68 ± 0.240.87 ± 0.7123,047.9 ± 465.9923.20 ± 0.10114.39 ± 9.33
4EY7_Quercetin1.71 ± 0.210.85± 0.4123,048.4 ± 49.5323.19 ± 0.09114.90 ± 8.41
4EY7_Myricetin1.72 ± 0.220.84 ± 0.5323,147.3 ± 487.8723.22 ± 0.09114.53 ± 9.41
Table 6. The means and standard deviations of several energy components that contribute to the free energy of binding between top-docked phytochemicals and target proteins.
Table 6. The means and standard deviations of several energy components that contribute to the free energy of binding between top-docked phytochemicals and target proteins.
SYSTEMΔVDWAALSΔEELΔEGBΔESURFΔGGASΔGSOLVΔTOTAL
4EY7_Donepezil−42.77 ± 3.09−9.72 ± 11.5438.57 ± 10.55–5.99 ± 0.417−52.59 ± 12.3432.7 ± 11.5−19.92 ± 3.62
4EY7_
quercetin
−31.26 ± 4.20−15.10 ± 7.5228.27 ± 5.27−4.10 ± 0.54−46.38 ± 9.3524.20 ± 5.10−22.20 ± 5.19
4EY7_
Myricetin
−38.23 ± 3.74−34.178± 10.2551.67 ± 6.76−5.10 ± 0.31−72.37 ± 9.5646.59 ± 6.75−25.78 ± 4.04
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ojo, O.A.; Ajayi-Odoko, O.A.; Gyebi, G.A.; Ayokunle, D.I.; Ogunlakin, A.D.; Ezenabor, E.H.; Olanrewaju, A.A.; Agbeye, O.D.; Ogunwale, E.T.; Rotimi, D.E.; et al. Network Pharmacology, Molecular Dynamics and In Vitro Assessments of Indigenous Herbal Formulations for Alzheimer’s Therapy. Life 2024, 14, 1222. https://doi.org/10.3390/life14101222

AMA Style

Ojo OA, Ajayi-Odoko OA, Gyebi GA, Ayokunle DI, Ogunlakin AD, Ezenabor EH, Olanrewaju AA, Agbeye OD, Ogunwale ET, Rotimi DE, et al. Network Pharmacology, Molecular Dynamics and In Vitro Assessments of Indigenous Herbal Formulations for Alzheimer’s Therapy. Life. 2024; 14(10):1222. https://doi.org/10.3390/life14101222

Chicago/Turabian Style

Ojo, Oluwafemi Adeleke, Omolola Adenike Ajayi-Odoko, Gideon Ampoma Gyebi, Damilare IyinKristi Ayokunle, Akingbolabo Daniel Ogunlakin, Emmanuel Henry Ezenabor, Adesoji Alani Olanrewaju, Oluwatobi Deborah Agbeye, Emmanuel Tope Ogunwale, Damilare Emmanuel Rotimi, and et al. 2024. "Network Pharmacology, Molecular Dynamics and In Vitro Assessments of Indigenous Herbal Formulations for Alzheimer’s Therapy" Life 14, no. 10: 1222. https://doi.org/10.3390/life14101222

APA Style

Ojo, O. A., Ajayi-Odoko, O. A., Gyebi, G. A., Ayokunle, D. I., Ogunlakin, A. D., Ezenabor, E. H., Olanrewaju, A. A., Agbeye, O. D., Ogunwale, E. T., Rotimi, D. E., Fouad, D., Batiha, G. E. -S., & Adeyemi, O. S. (2024). Network Pharmacology, Molecular Dynamics and In Vitro Assessments of Indigenous Herbal Formulations for Alzheimer’s Therapy. Life, 14(10), 1222. https://doi.org/10.3390/life14101222

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop