Phytochemical Analysis, Network Pharmacology and in Silico Investigations on Anacamptis pyramidalis Tuber Extracts

Anacamptis pyramidalis (L.) Rich. forms part of the Orchidaceae family that is highly valued for its horticultural as well as therapeutic benefits. The present study set out to investigate the inhibitory activity of A. pyramidalis tubers against key biological targets for the management of type 2 diabetes, Alzheimer disease, and skin hyperpigmentation. In addition, the antioxidant potential of the extracts was also assessed using multiple methods. The detailed phytochemical profiles of the extracts were determined using high-performance liquid chromatography. Based on qualitative phytochemical fingerprint, a network pharmacology analysis was conducted as well. Parishin was identified from the water extract only, whereas gastrodin and caffeic acid derivatives were present in the methanol extract. The methanol extract exhibited high inhibitory activity against tyrosinase (69.69 mg kojic acid equivalent/g extract), α-amylase (15.76 mg acarbose equivalent/g extract), and α-glucosidase (20.07 mg acarbose equivalent/g extract). Similarly, the methanol extract showed highest antioxidant potential (22.12, 44.23, 45.56, and 29.38 mg Trolox equivalent/g extract, for 2,2-diphenyl-1-picrylhydrazyl (DPPH), 2,2′-azino-bis(3-ethylbenzothiazoline-6-sulfonic acid) (ABTS), CUPric Reducing Antioxidant Capacity (CUPRAC), and Ferric Reducing Antioxidant Power (FRAP) assays, respectively). Finally, the results of network pharmacology analysis, besides corroborating traditional uses of plant extracts in the management of cold and flu, confirmed a direct involvement of identified phytochemicals in the observed enzyme inhibitory effects, especially against tyrosinase, α-amylase, and α-glucosidase. Furthermore, based on the results of both colorimetric assays and network pharmacology analysis related to the activity of A. pyramidalis extracts and identified phytocompounds on enzymes involved in type 2 diabetes, a docking study was conducted in order to investigate the putative interactions of oxo-dihydroxy octadecenoic acid trihydroxy octadecenoic acid against aldose reductase, peroxisome proliferator-activated receptor (PPAR)-α, dipeptidyl peptidase (DPP)-IV, and α-glucosidase. Docking analysis suggested the inhibitory activity of these compounds against the aforementioned enzymes, with a better inhibitory profile shown by oxo-dihydroxy octadecenoic acid. Overall, the present findings supported the rationale for the use of A. pyramidalis as source of bioactive metabolites and highlight, today more than ever, for the strong necessity of linkage strategy between wild resource valorization and conservation policy.

aim to elucidate the putative target proteins underlying the observed biological effects, whereas a docking analysis was conducted to confirm the interactions between selected phytochemicals and enzymes underlying the observed biological effects. It is expected that data generated from this investigation will provide valuable insights on the pharmacological potential of A. pyramidalis tubers, thus validating the traditional uses of this plant and improving the local product chain, also in terms of sustainability.

Phytochemical Profile
The standard assessment of bioactive compounds, in terms of phenolic and flavonoid contents is presented in Table 1. From data presented, A. pyramidalis tuber possessed low concentration of flavonoids. The phenolic contents of methanol and water extracts were not significantly (p > 0.05) different, thereby suggesting that solvents used were equally potent in extracting phenols from A. pyramidalis tuber. The major shortcoming of standard spectrophotometric methods in the determination of bioactive contents is the interference of other nonphenolic compounds, such as vitamins, which might yield inaccurate data. The use of high-performance liquid chromatography coupled to mass spectrometry detection in the field of phytochemistry provides explicit and detailed profiles of plant extracts, which might be crucial in understanding interaction with specific proteins. Values expressed are means ± S.D. of three parallel measurements. GAE: Gallic acid equivalent; RE: Rutin equivalent. Different letters ( a and b ) indicate significant differences in the extracts (p < 0.05).
To our best knowledge, this is the first report on the phytochemical profile of A. pyramidalis. The most abundant compounds were parishin and gastrodin derivatives, as well as a caffeic acid derivative. Similar fragmentation profiles have been reported in Gastrodia elata, also from the Orchidaceae family. Identification of compounds was performed by mass spectrometry (MS n )in negative mode. The identification was performed by comparison with analytical standards, when available, and by comparison of mass fragmentation profiles (MS/MS) with data from scientific literature. Compounds were numbered according to their retention times, keeping the same numeration in methanol and water extracts. Table 2 shows the characterization of compounds, whereas chromatograms are shown in Figure 1.
The fragmentation of compound 1 corresponded to a disaccharide formed by two hexose units: Loss of 162 daltons (Da) and fragment ions at m/z 179 and 161. It was present in both aqueous and methanolic extracts.
Compound 2 was identified as citric acid due to its [M − H] − at m/z 191 and fragment ions at m/z 173 and 111 (comparison with an analytical standard). Several citric acid glycosides were also characterized in A. pyramidalis extracts. Compounds 3 and 6 exhibited deprotonated molecular ions at m/z 459 and fragment ions at m/z 173 and 111 (typical of citric acid). This fragmentation has been previously reported for parishin G, 2-[4-O-(β-d-glucopyranosyl)benzyl] citrate, in G. elata [12]. Compounds 10 and 13 exhibited the same fragmentation patterns; they were characterized as parishin B and parishin C, respectively, taking into account their elution order [13]. Compound 16, with [M − H] − at m/z 995, was identified as parishin by comparison with bibliographic data [13].   [12]. Compounds 10 and 13 exhibited the same fragmentation patterns; they were characterized as parishin B and parishin C, respectively, taking into account their elution order [13]. Compound 16, with [M − H] − at m/z 995, was identified as parishin by comparison with bibliographic data [13].
Compound 8 corresponded to the formate adduct of roseoside (vomifoliol glucoside or drovomifoliol-O-β-D-glucopyranoside) [14]. This compound has not been previously reported in species of the Orchidaceae family.
Compounds 15, 19, and 20 were tentatively characterized as dihydroxybenzoic acid derivatives due to the 153→109 fragmentation. Compound 22 was characterized as a caffeic acid derivative due to the fragment ions at m/z 179 and 135 (comparison with an analytical standard).
Two acacetin derivatives were found in both methanol and aqueous extracts: Compounds 23 and 24. Both of them presented the aglycone acacetin at m/z 283, which was identified by its fragment ion at m/z 268.
Compounds 25 and 26, only present in the aqueous extract, were characterized as the oxylipins oxo-dihydroxy-octadecenoic acid and trihidroxy-octadecenoic acid [15], bioactive compounds that are produced in the oxidative metabolism of polyunsaturated fatty acids. Compound 8 corresponded to the formate adduct of roseoside (vomifoliol glucoside or drovomifoliol-O-β-d-glucopyranoside) [14]. This compound has not been previously reported in species of the Orchidaceae family.
Compounds 15, 19, and 20 were tentatively characterized as dihydroxybenzoic acid derivatives due to the 153→109 fragmentation. Compound 22 was characterized as a caffeic acid derivative due to the fragment ions at m/z 179 and 135 (comparison with an analytical standard).
Two acacetin derivatives were found in both methanol and aqueous extracts: Compounds 23 and 24. Both of them presented the aglycone acacetin at m/z 283, which was identified by its fragment ion at m/z 268.
Compounds 25 and 26, only present in the aqueous extract, were characterized as the oxylipins oxo-dihydroxy-octadecenoic acid and trihidroxy-octadecenoic acid [15], bioactive compounds that are produced in the oxidative metabolism of polyunsaturated fatty acids.
A relative semi-quantification was performed to check for the most abundant compounds. It was done by measuring peak areas of each compound in MS mode using the extracted ion chromatograms, with precursor ion [M − H] − . The relative percentage of each compound was calculated and is shown in Table 3, in which the heat map highlights the most abundant compounds (the darker the color, the higher the concentration). It can be observed that gastrodin derivatives (38.5%) and a caffeic acid derivative (compound 22, 16%) were the most abundant compounds in methanol extract, whereas the aqueous extract was also rich in gastrodin derivatives (35%), followed by caffeic acid derivative (16%), parishin G isomers (12%), and parishin B (7%). In terms of the total content of bioactive compounds, the methanol extract presented approximately two-fold concentration compared to the aqueous extract.

Enzyme Inhibition
Enzyme inhibitors have received due interest in the management of several diseases due to their role in pathophysiological mechanisms [16]. In fact, enzymes drive biological processes and as such have become the main strategy in drug design. Type 2 diabetes is a chronic disease, characterized by hyperglycemia. To control the progression of these diseases, dietary recommendations such as enzyme inhibitors could play a pivotal role. Acarbose, approved by the Food and Drug Administration for the treatment of type 2 diabetes in adults, competitively and reversibly inhibits pancreatic α-amylase and membrane-bound intestinal α-glucosidase. However, side effects associated to the use of acarbose, along with escalating prevalence of type 2 diabetes worldwide, are urging the scientific community to find safer alternatives. From Table 4, it is shown that the methanol extract of A. pyramidalis tubers was a more potent inhibitor of α-amylase (15.76 mg acarbose equivalent (ACAE)/g extract) and α-glucosidase (20.07 mg ACAE/g extract) compared with the water extract (5.23 mg ACAE/g extract, for both enzymes). The comorbidity of type 2 diabetes and Alzheimer disease has been supported by epidemiological evidences [17]. While referring to Alzheimer disease, researchers coined the term "type 3 diabetes" to underline the shared molecular and cellular features associated with insulin resistance, cognitive decline, and memory deficit [18]. Cholinesterase enzymes, responsible for the hydrolysis of the neurotransmitter acetylcholine, have been targeted in the management of Alzheimer disease. It has been observed that in the early stage of the disease a considerable increase in the activity of acetylcholinesterase was noted while butyrylcholinesterase activity, which has been often overlooked, shoots up to 90% in the late stages of Alzheimer disease, exacerbating the condition of the patient [19]. In the present study, the inhibitory capabilities of A. pyramidalis tubers' extracts were assessed on both acetylcholinesterase and butyrylcholinesterase. As shown in Table 4, the water extract exhibited the lowest inhibitory action on acetylcholinesterase and butyrylcholinesterase. A higher inhibition was observed for the methanol extract; this could be attributed to the presence of caffeic acid or acacetin as well as the other unknown compounds [20,21]. Also, we observed a high correlation (R > 0.9) total bioactive compounds' and enzyme inhibitory properties ( Figure 2). The comorbidity of type 2 diabetes and Alzheimer disease has been supported by epidemiological evidences [17]. While referring to Alzheimer disease, researchers coined the term "type 3 diabetes" to underline the shared molecular and cellular features associated with insulin resistance, cognitive decline, and memory deficit [18]. Cholinesterase enzymes, responsible for the hydrolysis of the neurotransmitter acetylcholine, have been targeted in the management of Alzheimer disease. It has been observed that in the early stage of the disease a considerable increase in the activity of acetylcholinesterase was noted while butyrylcholinesterase activity, which has been often overlooked, shoots up to 90% in the late stages of Alzheimer disease, exacerbating the condition of the patient [19]. In the present study, the inhibitory capabilities of A. pyramidalis tubers' extracts were assessed on both acetylcholinesterase and butyrylcholinesterase. As shown in Table 4, the water extract exhibited the lowest inhibitory action on acetylcholinesterase and butyrylcholinesterase. A higher inhibition was observed for the methanol extract; this could be attributed to the presence of caffeic acid or acacetin as well as the other unknown compounds [20,21]. Also, we observed a high correlation (R > 0.9) total bioactive compounds' and enzyme inhibitory properties (Figure 2).  The increased demand for depigmenting products derived from natural products has been fueled by increased public concern towards synthetic agents. Tyrosinase inhibitors play a key role in the formulation of depigmenting products. In fact, the inhibition of tyrosinase, a copper-containing, rate-limiting enzyme responsible for the biosynthesis of melanin, has been directly related to the skin-lightening effect of dermatological and cosmetic products used for the management of epidermal hyperpigmentation conditions. The methanol extract (69.69 mg·kainic acid equivalents (KAE)/g extract) of A. pyramidalis tubers showed higher inhibition against tyrosinase compared with the water extract (11.09 mg·KAE/g extract). In silico molecular docking studies conducted by Pei and colleagues [22] reported that gastrodin, identified in a significant amount in the methanol extract of A. pyramidalis tubers, interacted primarily with histidine residues of tyrosinase active site. It was suggested that acacetin, identified Agastache rugosa Kuntze and A. pyramidalis tubers (Table 3), could be related to the observed tyrosinase inhibition [23].

Antioxidant Activity
Assessing antioxidant properties of plant extracts is crucial in the evaluation of plants' bioactivity and plants' ability to prevent and/or mitigate health problems. Investigations have demonstrated that intake of antioxidants could prevent or delay the onset/progress of several human ailments [24][25][26]. In order to obtain a comprehensive understanding of the antioxidant potential of A. pyramidalis tubers' extracts, multiple antioxidant assays were employed. Results of the total antioxidant capacity, estimated using the phosphomolybdenum assay, are presented in Table 4. The methanol extract showed highest activity (0.73 mmol·trolox equivalents (TE)/g) for the phosphomolybdenum assay. A similar trend was observed for radical scavenging properties. Data shown in Table 5 revealed higher Trolox equivalent values for methanol extract, as regards the 2,2-diphenyl-1-picrylhydrazyl(DPPH) and 2,2 -azino-bis(3-ethylbenzothiazoline-6-sulfonic acid) (ABTS) radical scavenging assays. Likewise, a higher reducing activity was observed for the methanol extract in the CUPRAC and FRAP assays (45.56 and 29.38 mg·TE/g extract). The metal chelating properties of the extracts was also assessed. It was observed that the water extract exhibited high metal chelating potential. According to correlation analysis, a strong correlation between total phenolic, flavonoid content, and antioxidant properties (except for metal chelating) was observed ( Figure 2). From this point, observed metal chelating activity could attribute to nonphenolic metal chelators in the extracts.

Prediction of Pharmacologic Targets and Pharmacokinetic Profile
Based on qualitative phytochemical fingerprint assessment, a network pharmacology approach was conducted on selected extracts' phytochemicals, in order to elucidate the putative target proteins underlying both the observed biological effects and the traditional uses. According to the pharmacokinetic predictions of gastrointestinal adsorption (results depicted in Supplementary Materials: Pharmacokinetic folder), carried out through SwissADME bioinformatic platform, acacetin, caffeic acid, dihydroxy-benzoic acid, gastrodin, oxo-dihydroxy-octadecenoic acid, parishin A, parishin B, parishin C, roseoside, and trihydroxy-octadecenoic acid were assayed for the identification of putative target proteins (SwissTargetPrediction bioinformatic platform). As expected, this last bioinformatic platform yielded a wide plethora of targets interacting with the selected secondary metabolites. Details about the bioinformatic analysis were reported in Supplementary Materials (target proteins folder), whereas component-target analysis, conducted with Cytoscape software, is illustrated in Figure 3. Specifically, component-target analysis indicated the capability of both roseoside and gastrodin to interact with tyrosinase. The putative interaction between gastrodin and tyrosinase was consistent with previous docking studies [23]. However, the sole roseoside could interact with glucosidase, whereas the interaction of the sole gastrodin with maltase glucoamylase was predicted, as well. Overall, these predicted interactions were consistent, at least in part, with the enzyme inhibitory effects induced by A. pyramidalis extracts, which could represent promising sources of bioactive compounds for the treatment of skin hyperpigmentation and type 2 diabetes. Regarding the putative antidiabetic activity of the tested extracts, it is also of noteworthy interest to highlight the interaction of oxo-dihydroxy-octadecenoic and trihydroxy-octadecenoic acid with dipeptidyl peptidase (DPP)-IV and peroxisome proliferator-activated receptor (PPAR)-α, which represent key targets of anti-diabetic therapy [27]. These compounds were also predicted to interact with multiple prostanoid receptors and enzymes, which could be at the basis of the efficacy of A. pyramidalis infusion against flu and cold [5]. Regarding the traditional uses of this plant against neurological and psychiatric diseases [4], some concerns arise from the network pharmacology approach on selected phytocompounds. On one side, caffeic acid, gastrodin, and acacetin were predicted to interact with target proteins involved in neurotransmitter pathways (i.e., dopamine, norepinephrine, serotonin, adenosine). On the other hand, the lack of any blood-brain barrier permeant property shown by selected compounds (depicted in Supplementary Materials: Pharmacokinetic folder) rule out a direct effect of A. pyramidalis-derived phytochemicals in central circuitries involved in neurological and neuropsychiatric effects. Nevertheless, considering the putative interactions with dipeptidyl peptidase IV, an enzyme deeply involved in the gut-brain axis mediated by neuropeptides [28], the present findings cannot exclude indirect effects induced by A. pyramidalis extracts on brain functions.

Docking Results
Based on the results of enzyme inhibition assays and network pharmacology analysis, a docking study was conducted to explore putative interactions between selected A. pyramidalis phytochemicals, namely oxo-dihydroxy octadecenoic acid and trihydroxy octadecenoic acid, and key enzymes involved in type 2 diabetes (i.e., aldose reductase, PPAR-α, DPP-IV, and αglucosidase). The calculated binding free energy, inhibition constant (Ki), and the nonbonding interactions represent the outcomes of any docking study. These parameters are essential to rank, compare, and design inhibitors for any drug design study. The results of docking study are listed in Table 6. Interestingly, oxo-dihydroxy octadecenoic acid showed higher binding affinity than trihydroxy octadecenoic acid toward the selected enzymes. In addition, both compounds showed

Docking Results
Based on the results of enzyme inhibition assays and network pharmacology analysis, a docking study was conducted to explore putative interactions between selected A. pyramidalis phytochemicals, namely oxo-dihydroxy octadecenoic acid and trihydroxy octadecenoic acid, and key enzymes involved in type 2 diabetes (i.e., aldose reductase, PPAR-α, DPP-IV, and α-glucosidase). The calculated binding free energy, inhibition constant (K i ), and the nonbonding interactions represent the outcomes of any docking study. These parameters are essential to rank, compare, and design inhibitors for any drug design study. The results of docking study are listed in Table 6. Interestingly, oxo-dihydroxy octadecenoic acid showed higher binding affinity than trihydroxy octadecenoic acid toward the selected enzymes. In addition, both compounds showed the same trend in which the binding affinity went from aldose reductase, PPAR-α, DPP-IV, and α-glucosidase enzymes. Obviously, hydrogen bonds are the abundant interaction of these compounds, however, the number of these interactions are different for each compound and each enzyme. The hydrogen bond interactions and other hydrophobic interactions of the studied compounds are shown in Figure 4. Table 6. The calculated binding free energy, ∆G, in kcal/mol, inhibition constant, K i , the key residues and the number of hydrogen atoms of the docked compounds.

Plant Material and Preparation of Extracts
The plant materials were collected in Turkey (Kastamonu) in 2019 from wild, flowering plants. Taxonomical identification was performed by co-author, the botanist Dr. Ismail Senkardes. A voucher specimen is conserved at Marmara University, Faculty of Pharmacy Herbarium. Collection was performed to obtain representative samples without damaging the consistence of wild populations. The tubers were manually separated, cleaned, roughly sectioned, and taken to dryness in a ventilated oven (temperature 45 • C) until reaching a constant weight. The dry plant material was powdered using a laboratory mill and directly used for extraction or stored in a vacuo-plastic bag and stored in the dark until used.
Extraction by maceration was used to obtain methanol extract. Briefly, pulverized tuber (5 g) was stirred with 100 mL of methanol for 24 h at room temperature. The solution was then filtered and the solvent was evaporated by using rotary evaporator. Infusion was selected for water extract. Briefly, the plant material (5 g) was kept in the boiled water (100 mL) for 15 min. After cooling, the extract was filtered and lyophilized. The obtained dry extracts were stored at +4 • C in the dark.

Profile of Bioactive Compounds
To obtain total levels of phenolic and flavonoid content in the extracts, colorimetric assays were used as described in our previous paper [29]. Gallic acid ((GAE) for total phenolic (TPC)) and rutin ((RE) for total flavonoid (TFC)).

Chromatographic Analysis
For HPLC-MS analysis, 5 mg of dried extract (DE) were re-dissolved in 1 mL of methanol, filtered through 0.45 µm polytetrafluoroethylene (PTFE) membrane filters, and 10 µL of the solution was injected. The phytochemical profile was obtained with an Agilent Series 1100 liquid chromatograph with a G1315B diode array detector and an ion trap mass spectrometer (Esquire 6000, Bruker Daltonics, (Billerica, Massachusetts, USA) with an electrospray interface. A Luna Omega Polar C 18 analytical column of 150 × 3.0 mm and 5 µm particle size (Phenomenex) was used. The method was adapted from a previous article of our research groups [30]. Separation was performed with a mobile phase of water-formic acid (100:0.1, v/v) and CH 3 CN. The following program was used: (1) Initial mobile phase, 10% CH 3 CN; (2) linear increase from 10% to 25% CH 3 CN (0-25 min); (3) 25% CH 3 CN (25-30 min); (4) linear increase from 25% to 50% CH 3 CN (30-40 min); (5) linear increase from 50% to 100% CH 3 CN (40-42 min); and (6) 100% CH 3 CN (42-47 min). Then, CH 3 CN percentage was returned to the initial mobile phase, with a 7 min stabilization time. The flow rate was 0.4 mL min −1 . The scan range was m/z 100-1200 with a speed of 13,000 Da/s. The electrospray ionization (ESI) conditions were: Drying gas (N 2 ) flow rate and temperature, 10 mL/min and 365 • C; nebulizer gas (N 2 ) pressure, 50 psi; capillary voltage, 4500 V; capillary exit voltage, −117.3 V. We used the auto MS n mode for the acquisition of MS n data, with isolation width of 4.0 m/z, and fragmentation amplitude of 0.6 V (MS n up to MS 4 ).

Determination of Antioxidant and Enzyme Inhibitory Effects
The antioxidant potential of the extracts was evaluated by phosphomolybdenum, antiradical (DPPH and ABTS), reducing power (FRAP and CUPRAC) and ferrous chelating assays, as described by Grochowski et al. [31]. Trolox equivalents were used for expression of antioxidant activities. Ethylenediaminetetraacetic acid was employed as a reference compound for the metal chelating assay.

Prediction of Putative Targets and Pharmacokinetics
Putative targets were identified through the bioinformatic method recently described by Gu and colleagues [32]. Briefly, chemical structures were prepared and converted in canonical "Simplified Molecular Input Line Entry System" (SMILES) using ChemSketch software (12.0 version). The SMILES were then processed by the SwissTargetPrediction (http://www.swisstargetprediction.ch/) and SwissADME (http: //www.swissadme.ch/index.php) platforms, for predicting putative targets and pharmacokinetic profile, respectively. The name of identified targets were normalized according to UniProt database (https: //www.uniprot.org/). Finally, Cytoscape software (3.7.2 version) was used to create a component-target illustration network.

Docking Calculations
In order to investigate the binding affinity and the nonbonding interactions of the compounds, oxo-dihydroxy octadecenoic acid and trihydroxy octadecenoic acid against aldose reductase, DPP-IV, PPAR-α, and α-glucosidase enzymes, docking calculations were performed using Autodock4 software (version 4.2) (Molinspiration Database) [33]. The crystal structures of the target enzymes were downloaded from Protein Data Bank (PDB). PDBID:4GCA was used to get the crystal structure of aldose reductase enzyme in which the enzyme was crystalized with IDD 1219 inhibitor, PDBID:4N8D and PDBID:1K7L were used in the case of DPP-IV and PPAR-α enzymes, respectively, and PDBID:5NN5 for the crystal structure of α-glucosidase enzyme. All the co-crystalized molecules, such as the inhibitor and water molecules, were removed and all the proteins were neutralized by adding polar hydrogen atoms and Kollman united atom charges. The initial 3D structures of the two compounds were optimized using AM1 semi-empirical method [34] and the structures were saved in mol2 format. Autogrid 4 was used to design a grid box with 60 × 60 × 60 dimensions with 0.375 Å distance between points. Lamarckian genetic algorithm was used to scan 250 conformations for each inhibitor. The values of the calculated binding affinity were clustered and ranked in the output file. Discovery studio 5.0 visualizer was used to view the results and study the enzyme-inhibitor nonbonding interactions.

Statistical Analysis
The antioxidant and enzyme inhibitory activities results are given as mean ± standard deviation (S.D.). The results were statistically evaluated using the student t-test (α = 0.05). Pearson correlation coefficients were also calculated for total bioactive compounds and biological activities. Statistical analysis was carried out using SPSS v. 14.0 program (SPSS Inc, Chicago, IL, USA).

Conclusions
The current investigation highlighted the pharmacological potential of orchids, with a particular focus on traditionally used A. pyramidalis tubers. This study also supported the use of cutting-edge technologies to unravel the phytochemical profile of herbal extracts. The identification of phytochemicals could then sustain observed biological activities, with particular regards to enzyme inhibition. Parishin G isomer-2, gastrodin, and caffeic acid derivatives were the main compounds in the tested extracts and the extracts exhibited potent enzyme inhibitory, radical scavenging, and reducing properties. In addition, the outcomes of network pharmacology and docking calculations showed the affinity of oxo-dihydroxy octadecenoic acid and trihydroxy octadecenoic acid against aldose reductase, PPAR-α, DPP-IV, and α-glucosidase enzymes, which suggests the pharmacological potential of A. pyramidalis phytochemicals against type 2 diabetes. Funding: Technical and human support provided by Centro de Instrumentación Científico Técnica (CICT) of Universidad de Jaén (UJA, MINECO, Junta de Andalucía, FEDER) are gratefully acknowledged. The study was supported by the Italian Ministry of University funds (FAR 2019) granted to Giustino Orlando, Annalisa Chiavaroli, and Claudio Ferrante.

Conflicts of Interest:
The authors declare no conflict of interest.