==== Front Molecules Molecules molecules Molecules 1420-3049 MDPI 33265992 10.3390/molecules25235636 molecules-25-05636 Article Metabolite Characterization and Correlations with Antioxidant and Wound Healing Properties of Oil Palm (Elaeis guineensis Jacq.) Leaflets via 1H-NMR-Based Metabolomics Approach https://orcid.org/0000-0002-5704-5053Che Zain Mohamad Shazeli 1 https://orcid.org/0000-0002-2997-7910Lee Soo Yee 1 Nasir Nadiah Mad 2 Fakurazi Sharida 3 https://orcid.org/0000-0002-7054-0649Shaari Khozirah 12* Sikorska Ewa Academic Editor 1 Natural Medicines and Products Research Laboratory (NaturMeds), Institute of Bioscience, Universiti Putra Malaysia, Serdang 43400, Malaysia; shazelizain@gmail.com (M.S.C.Z.); daphne.leesooyee@gmail.com (S.Y.L.) 2 Department of Chemistry, Faculty of Science, Universiti Putra Malaysia, Serdang 43400, Malaysia; nadiahmadnasir@upm.edu.my 3 Laboratory of Vaccines and Biomolecules (VacBio), Institute of Bioscience, Universiti Putra Malaysia, Serdang 43400, Malaysia; sharida@upm.edu.my * Correspondence: khozirah@yahoo.com.my; Tel.: +60-3-8942148 30 11 2020 12 2020 25 23 563614 11 2020 27 11 2020 © 2020 by the authors.2020Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).Oil palm (Elaeis guineensis Jacq.) leaflets (OPLs) are one of the major agricultural by-products generated from the massive cultivation of Malaysian palm oil. This biomass is also reported to be of potential value based on its health-improving effects. By employing proton nuclear magnetic resonance (1H-NMR) spectroscopy combined with multivariate data analysis (MVDA), the metabolite profile of OPLs was characterized and correlated with their antioxidant and wound healing properties. Principal component analysis (PCA) classified four varieties of extracts, prepared using solvents ranging from polar to medium polarity, into three distinct clusters. Cumulatively, six flavonoids, eight organic acids, four carbohydrates, and an amine were identified from the solvent extracts. The more polar extracts, such as, the ethyl acetate-methanol, absolute methanol, and methanol-water, were richer in phytochemicals. Based on partial least square (PLS) analysis, the constituents in these extracts, such as (+)-catechin, (−)-epicatechin, orientin, isoorientin, vitexin, and isovitexin, were strongly correlated with the measured antioxidant activities, comprising ferric reducing antioxidant power (FRAP), 2,2-diphenyl-1-picrylhydrazyl (DPPH), and nitric oxide (NO) free radical scavenging activities, as well as with cell proliferation and migration activities. This study has provided crucial evidence on the importance of these natural antioxidant compounds on the wound healing properties of OPL. Elaeis guineensis Jacq.oil palm leafmetabolite profilesolvent extractantioxidantwound healing1H-NMR metabolomics ==== Body 1. Introduction Oil palm (Elaeis guineensis Jacq.) leaflets (OPLs) are considered as an underutilized agricultural waste, though approximately 21.7 million tonnes of OPLs are produced annually from the Malaysian oil palm industry [1]. Studies have revealed the potential benefits of this biomass, and have eventually led to the use of OPLs as adsorbents [2], composites [3], soil fertility enhancers [4], fuels [5], papermaking materials [6], and ruminant feedstock [7]. Recently, studies have also reported OPL to contain phenolic compounds, and hence an alternative source of bioactive constituents for nutraceutical and pharmaceutical industries [8,9]. Solvent selection is one of the crucial contributors for effective extraction of phytoconstituents from plants [10]. It is necessary to consider the solvent’s polarity, reactivity, and ease of removal after extraction [11], as these factors contribute directly to the types of compounds that can be extracted, the extractable amount, and, hence, the biological activity of the extracts [12,13]. Therefore, before selecting the most suitable solvent for extracting target compounds of interest, it is necessary to compare the compositions and the bioactivity of the extractable constituents obtained from the different solvents [14,15]. Metabolomics is a robust technique for the study of small-molecule metabolite profiles of a plant matrix [16]. It is one of the fastest developing research areas, utilizing mostly high-field proton nuclear magnetic resonance spectroscopy (1H-NMR) or liquid chromatography mass spectrometry (LC-MS) as analytical platforms. With the use of multivariate statistical analysis (MVDA), samples can be differentiated based on their metabolite compositions. Through the interpretation of spectral data, and further comparison with databases and literature data, the important differentiating metabolites can be identified [17]. Thus, this technique is a highly promising tool for quality control and standardization of natural products [18]. The results obtained are highly reproducible and, more interestingly, MVDA in metabolomics studies also allows the identification of potential bioactive metabolites through correlations between the detected metabolites with the specific biological activity [19]. The utilization of 1H-NMR-based metabolomics for untargeted analysis to evaluate the correlation of known metabolites in plant extracts to their biological activities have been reported by several studies in the past such as on Lysilomo latisiliquum [20], Neptunia oleracea [21], and Aquilaria malaccensis Lamk [22]. Some aspects of metabolomics analysis of OPL have been previously reported by Vargas et al., [23] where they employed ultra-high-performance liquid chromatography-tandem mass spectrometry (UHPLC-MS/MS), as the analytical platform, to understand how transportation/storage conditions and extraction solvent affected its metabolite profiles [23]. However, to the best of our knowledge, there have been no studies conducted investigating the correlations between the multitude of metabolites in the OPL extract with any biological activity. In the present work, using 1H-NMR-based metabolomics approach, we performed a discriminative analysis of OPL extracted with solvents of different polarity (absolute methanol, aqueous methanol, ethyl acetate-methanol and ethyl acetate), identified the differentiating metabolites, and further examined the correlations between the OPL metabolites to their antioxidant and wound healing activities. The predicted bioactive metabolites were also relatively quantified. 2. Results 2.1. Visual Inspection of 1H-NMR Spectra of OPL Extracts and Metabolite Identification The representative 1H-NMR spectra of OPL samples extracted by different polarities of solvents is displayed in Figure 1. Previous studies on OPL have reported primary and secondary metabolites, which are organic acids, amino acids, simple carbohydrates, phenolic acids, luteolin, and apigenin derivatives [8,9,23,24,25,26]. Nevertheless, the data with regard to the secondary metabolites present in OPL are still considered limited. The peak assignment for signals in the 1D spectra were achieved based on comparisons with literature data [21,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45] and freely available online databases, including the Nuclear Magnetic Resonance Database (NMRDB) [46], Human Metabolome Database (HMDB) [47], and PubChem [48]. Basically, 3 regions of chemical importance were visually analyzed in the 1H-NMR spectra of OPL extracts prepared at different extraction solvent, viz., high-field (0.9–3.50 ppm), middle-field (3.50–5.45 ppm) and low-field regions (6.0–8.2 ppm). As expected, overlapping peaks make peak assignment of the compounds very challenging. Nevertheless, signals for 19 metabolites could be discerned from careful visual inspection of the spectra. Additionally, 2D experiments (J-resolved) provided additional information on signal splitting (Figure S1 in Supplementary Materials). The 1H-NMR spectra revealed the presence of metabolites from different classes including amino acids, organic acids, carbohydrates, and flavonoids. Collectively, 6 flavonoids, 8 organic acids, 1 amino acid, and 4 carbohydrates were recognized. The characteristic 1H-NMR chemical shifts corresponding to these recognized metabolites are presented in Table 1. On the whole, the efficiency of the solvents when extracting OPL metabolites improved with increasing polarity. The highly polar solvents were able to extract most of the carbohydrate metabolites, which could be assigned from the signals observed in the chemical shift region of 3.0–5.5 ppm [53]. The chemical shifts 5.42, 4.18, 5.18 and 4.62 ppm were recognized as sucrose, fructose, α-glucose, and β-glucose, respectively. The presence of these carbohydrates were in agreement with previous studies on oil palm leaf chemical composition [25,54]. The signals for these carbohydrates were also discernable in ethyl acetate, a medium polar solvent, albeit in lower quantities. The medium and highly polar solvents were also able to extract organic acids and the nutrient, choline. Except for fumaric acid, which was identified from the characteristic olefinic proton singlet at 6.54 ppm, the chemical shifts of the organic acids (citric acid, palmitic acid, succinic acid, aconitic acid, fatty acid, α-linolenic acid, and acetic acid) were identified based on their signals observed in the aliphatic region (0.90–3.50 ppm). These organic acids were detected in all the solvent extracts, except fumaric acid which was absent in ethyl acetate extracts. Extracts prepared from methanol, either on its own (absolute) or as admixtures with ethyl acetate or water, showed proton signals characteristic of flavonoids, specifically those in the downfield region between 6.00 and 8.00 ppm. Among the major flavonoids detected in these extracts were apigenin, luteolin, and catechin derivatives, in accordance with our previous findings using UHPLC-MS/MS analysis [25,26]. The flavonoids were detected as isomeric pairs—i.e., vitexin with isovitexin, orientin with isoorientin, and catechin with epicatechin (Figure 2). The chemical shifts 8.02, 6.89, and 6.77 ppm were assigned to vitexin while signals at 7.79, 6.88, and 6.52 ppm were assigned to isovitexin. The signals resonating at 7.54, 6.86, and 6.63 ppm were assigned to orientin while signals detected at 7.40, 6.89, 6.67, and 6.48 ppm were assigned to isoorientin. The proton signals of orientin and isoorientin were lower than those exhibited by the protons of vitexin and isovitexin at the same positions. Compared to vitexin and isovitexin, orientin and isoorientin experienced a shielding effect due to the additional hydroxyl group on ring B (Figure 2). This additional substituent acted as an electron donating group which increased the electron density of the nucleus, resulting in a slight shielding effect and hence lower chemical shifts. Intense peaks between 3.00 to 4.00 ppm were observed which were assigned to the sugar moiety attached to the aglycone of identified flavonoids. The peak assignments were based on comparison with literature data (Table 1) and confirmed with chemical shift data of commercial standards (Figure S2). Meanwhile, the catechin/epicatechin isomeric pair was detected based on the chemical shifts occurring at 7.04, 7.01, and 6.95 ppm. The peaks at 7.04 and 6.95 ppm were assigned to (−)-epicatechin, while the peak at 7.01 ppm was assigned to (+)-catechin. These signals assignment were in agreement with the chemical shifts reported previously [50]. The difference between (+)-catechin and (−)-epicatechin is the different configurations of C-2 and C-3. (+)-Catechin is the isomer with a trans-configuration while (−)-epicatechin has a cis-configuration (Figure 2). The trans-related H-2 and H-3 protons of epicatechin were slightly more deshielded in comparison to the cis-related protons of catechin. 2.2. Classification of OPL Extracts by Principal Component Analysis (PCA) Variation in metabolites in the different solvent extracts was examined by multivariate analysis (MVDA). The 1H-NMR datasets of the aqueous methanol, absolute methanol, ethyl acetate-methanol and ethyl acetate extracts were subjected to PCA to primarily distinguish between the extracts, and to identify the discriminating factors. The PCA model was obtained with good fitness (R2X) and high predictability (Q2) values of 0.995 and 0.960, respectively. There were no notable outliers. The first principal component (PC1) provided 79.7% of the variance, while the second principal component (PC2), gave 7.23%, cumulatively describing a total variance of 86.93%. The PCA score plot (Figure 3A) showed that different solvent extracts were grouped into three clusters with no notable outliers. The ethyl acetate extract was well separated from the methanolic extracts (aqueous methanol, absolute methanol, and ethyl acetate-methanol) by PC1. The absolute methanol and ethyl acetate-methanol extracts were also projected closer to each other, rather than with the aqueous methanol extract. This indicated that the phytochemical constituents in the absolute methanol and ethyl acetate-methanol extracts were very similar. The same metabolites may have been extracted by the two solvent systems but in varying quantities. The discriminating metabolites between the various clusters were determined from the loading plot (Figure 3B). The ethyl acetate extract contained higher amounts of α- and β-glucose, acetic acid, and fatty acids, notably, the essential fatty acid α-linolenic acid. Meanwhile, the methanolic solvent systems were able to extract a large number of compounds, comprising moderately polar to polar molecules. The methanolic solvent systems extracted both primary (viz. succinic acid, citric acid, aconitic acid, fumaric acid, choline, sucrose, and fructose) and secondary metabolites (viz. apigenin, luteolin, and catechin derivatives). However, although the ethyl acetate-methanol and absolute methanol extracts appeared to be in a separate cluster from the aqueous methanol extract in the score plot, there were no obvious discriminants between these three solvent extracts. The loading plot showed a congregation of most of the 1H-NMR chemical shifts close to PC2. Overall, as expected, the PCA results revealed that the methanolic solvent systems were capable of liberating both primary and secondary metabolites from OPL, while ethyl acetate mainly extracted some primary metabolites. Relative quantification of the identified metabolites (Figure 4) showed no significant difference in their extractability when methanolic solvent systems were used, while ethyl acetate, representing a less polar system, was mainly efficient for extraction of the fatty acids, particularly α-linolenic acid. The results indicated that the polar solvent systems were the better solvents for extraction or recovery of most OPL metabolites, especially the flavans and the four major OPL flavonoids. 2.3. Influence of OPL Extracts on Polyphenolic Contents, Antioxidant, and Wound Healing Properties As shown in Table 2, the phenolic contents, antioxidant and wound healing properties are influenced by the choice of solvent used for extraction. For total phenolic content (TPC), there was a positive relationship with solvent polarity, where the TPC increased with increasing solvent polarity. Aqueous methanol extract had the highest TPC, with a value of 393.61 mg GAE/g extract, followed by absolute methanol and ethyl acetate-methanol extracts, respectively. On the other hand, the absolute methanol extract showed the highest Total Flavonoid Content (TFC) with a value of 135.40 mg QCE/g extract, followed by aqueous methanol, and ethyl acetate-methanol extracts. Meanwhile, the ethyl acetate extract had the lowest TPC and TFC values, with 121.71 mg GAE/g and 5.94 mg QCE/g extract, respectively. Statistically, there were no significant differences (p < 0.05) in the TFC extracted by the three methanolic solvents. These results suggested that the methanolic solvent systems are more efficient for the extraction of polyphenols in OPL. The results of the antioxidant assays on the solvent extracts were similar to the trend of the polyphenolic contents, where the methanolic extracts exhibited significantly greater activities in comparison with ethyl acetate extracts. The methanolic extracts exhibited high ferric antioxidant power (FRAP) with values between 71.62 and 101.48 mg AAE/g extract. There was also no significant difference between the values. For the 2,2-diphenyl-1-picrylhydrazyl (DPPH) and nitric oxide (NO) free radical scavenging activities, the low IC50 values revealed the potency of the methanolic extracts. The aqueous and absolute methanol extracts showed the most significant scavenging of DPPH free radicals with IC50 values of 3.53 and 6.54 µg/mL, respectively. Meanwhile, the aqueous methanol extract recorded the lowest IC50 value of 18.77 µg/mL in the NO free radical scavenging, followed by absolute methanol with IC50 values of 67.64 µg/mL. These findings emphasized the efficiency of methanol for the extraction and recovery of antioxidant principles from OPL and are in agreement with Tahir et al. [8]. The various solvent extracts of OPL were also evaluated for their wound healing properties by assessing their efficacy in enhancing the proliferation and migration of 3T3 cells. Prior to this, to find the optimum concentration of extracts for cell proliferation and migration, various concentrations, ranging from 1.56 to 100 µg/mL, were preliminarily screened. The cells proliferated and migrated rapidly at a low concentration of 1.56 µg/mL. The results of cell proliferation and migration activities at this extract concentration are presented in Table 2. With the high cell proliferation activity (>95%), the results indicated that all the OPL extracts did not affect the cell viability of 3T3 cells, besides promoting cell growth. Meanwhile, a wound scratch test assay was conducted to assess the migration of cells. The results revealed that methanolic extracts showed high cell migration activity, enhancing the 3T3 cells’ migration at percentages ranging from 85.83 to 93.34%, with no significant difference between them. Overall, methanolic extracts were found to possess good cell proliferation and migration activities. Pearson’s correlation was performed to assess the relationship between the polyphenolic contents and the bioactivities. The results are presented in Table 3. There were positive correlations between the phenolic contents and the bioactivities, with correlations between TPC and cell proliferation, between TFC and DPPH free radical scavenging activity, and between TFC and cell migrations as the most significant ones. This strong positive relationship revealed the contribution of phenolic constituents, particularly the flavonoids, to their respective bioactivities. The promotion of cell proliferation by phenolic compounds and DPPH scavenging and cell migration enhancement by flavonoids have been reported previously [26,55,56]. In addition, the Pearson’s correlation also revealed a strong correlation between antioxidant and wound healing activities. These results could provide information on the role of free radical scavenging as one of the important events in wound healing [57]. Although Pearson’s correlation showed that the phenolic constituents contributed to the antioxidant and wound healing activities of the OPL extracts, the individual metabolites that are responsible for the bioactivities require further investigation. 2.4. Correlation between Identified Metabolites and Antioxidant and Wound Healing Activities of OPL Extracts To comprehend the relationship between the individual metabolites and the bioactivities, the correlation was further extended by computing the supervised MVDA. Partial least squares (PLSs) projection to latent structures was applied to determine the chemical constituents that are correlated to the bioactivities. PLS can relate the data of independent variables (X-matrix, 1H-NMR spectral data) to the data of dependent variables (antioxidant and wound healing activities), hence the potential biomarkers can be predicted. The PLS model obtained from the NMR spectral and bioactivity data of OPL extracts was validated by means of internal cross validation, permutation test, and regression analysis. Autofit of the PLS model in SIMCA revealed a good fitness (R2Y = 0.868) and predictive ability (Q2 = 0.60) of the model (Figure S3). The permutation plots presented in Figure S4 revealed that there was no overfitting of the PLS model. The Y-intercepts of R2 and Q2 were less than 0.3 and 0.05, respectively. These data indicated that the PLS model showed a satisfactory level of validity [41,42]. Similar to the PCA score plot, in the PLS biplot, the OPL extracts were separated into three clusters, as shown in Figure 5. The methanolic extracts were clustered on the positive side of PC1 while the ethyl acetate extract was on the negative side. All the Y variables (DPPH, NO, FRAP, cell proliferation, and cell migration activities) were projected on the same side as the methanolic extracts, demonstrating their better activities compared to the ethyl acetate extract. These results were consistent with the bioactivitiy results discussed previously. The phenolic constituents, particularly the flavans and flavonoids namely (+)-catechin, (−)-epicatechin, orientin, isoorientin, vitexin, and isovitexin, were located close to the Y variables, revealing their significant contribution towards the antioxidant and wound healing activities of OPL. These results supported the Pearson’s correlation results. Previous studies have reported that different plant matrices containing catechin, apigenin, and luteolin derivatives exhibited arrays of bioactivities comprising antioxidant, anti-inflammatory, and wound healing activities [8,25,26,55,58]. The primary metabolites comprising carbohydrates (sucrose and fructose) and organic acids (citric acid, succinic acid, fumaric acid, and aconitic acid) were also shown to be correlated with the antioxidant and wound healing activities. These findings are in accordance with other studies that reported primary metabolites, such as organic acids and carbohydrates, also contribute to antioxidant and wound healing activities of a plant extract [59,60]. 3. Discussion Overall, the findings revealed the methanolic extracts were rich in secondary metabolites, namely flavonoids, as displayed in proton NMR spectra. Both Pearson and PLS correlation studies demonstrated these groups of compounds were strongly correlated with studied antioxidant and wound healing properties. The role of these flavonoids as natural antioxidants often relates to their potential scavenging activity in combating excessive oxidative stress during wound healing. In normal metabolism, the level of oxidants and antioxidants is maintained; however, overproduction of oxidants can lead to oxidative damage that may spread all over the targeted cells such as DNA, lipids, and proteins. Hence, the presence of antioxidant molecules is able to shield these targeted cells at the risk of damage by free radicals or known as reactive oxygen species (ROS) such as super anion radical (O2•), singlet oxygen (O•), nitric oxide (NO•), hydroxyl radical (HO•), peroxynitrite (ONOO•), and hydrogen peroxide (H2O2) [61]. Studies reported that wounds cause the overproduction of ROS during wound repair processes as a defence mechanism against invading bacteria resulting in severe tissue damage at the targeted cells and even lead to neoplastic transformation which reducing the wound healing process [62,63]. Conversely, the antioxidant defence system works efficiently for neutralizing these free radicals by undergoing three main phases. Firstly, metal chelating proteins and enzymatic antioxidants including peroxidases will stop the production of free radicals followed by scavenging the free radicals by non-enzymatic or radical scavenging antioxidants to ensure no formation of oxidation chain. Lastly, enzymes such as lipases, proteases, and DNA repair enzymes will have roles in repairing and replacing the damages of the targeted cells and to reconstitute the damaged membranes [61]. The antioxidant activities expressed by flavonoids could be related to their structure-antioxidant activity relationships. The presence of one or more phenolic hydroxyl groups (-OH) in their structure plays a significant role in these activities [59,64,65,66]. The contribution of flavonoids, especially catechin, luteolin, and apigenin derivatives, in antioxidant action can be explained by direct scavenging of ROS via hydrogen atom donation, which is one of the common antioxidant action mechanisms of flavonoids [67]. To participate in the mechanism, the flavonoid structure should have an ortho-dihydroxy substitution in the B ring and C2-C3 double bond, or C4 carbonyl group in the C ring [67]. Flavanols (catechin isomers) and flavones (vitexin and orientin isomers) contain hydroxyl groups attached to carbons C4′ and/or C3′ of ring B (Figure 2). The hydroxyl groups attached to these positions could actively participate in the hydrogen transfer mechanism. For example, in the case of (+)-catechin (A), the hydrogen atom of the C4′-OH group was donated to the radical moiety (RO•) (Figure 6). The resulting unstable oxygen atom on C4′ then formed a hydrogen bond with hydroxyl group on C3′ forming a stabilized catechin radical intermediate (B). Then, the catechin radical stabilized a second radical moiety (RO•) by donating the hydrogen atom of the C3′-OH, forming the stable quinone structure C. 4. Materials and Methods 4.1. Reagents, Solvents, and Commercial Standards Analytical grade absolute methanol and ethyl acetate, deuterated methanol-d4 (CH3OH-d4), non-deuterated KH2PO4, sodium deuterium oxide (NaOD), trimethylsilyl propionic acid-d4 sodium salt (TSP) and deuterium oxide (D2O) were supplied by Merck (Darmstadt, Germany). Water was purified by a Milli-Q system (Millipore Lab, Bedford, MA, USA). The polyphenolic contents and antioxidant assays reagents (DPPH powder, Folin-Ciocalteu, sodium acetate, 2,4,6-tris(2-pyridyl)-S-triazine (TPTZ), ferric chloride, and sodium carbonate) and commercial standards (quercetin, gallic acid, and ascorbic acid) were purchased from Sigma (Aldrich, Darmstadt, Germany), while aluminum chloride was provided by HmbG Chemicals (Hamburg, Germany). Conventional organic solvents such as ethanol, methanol, ethyl acetate, and dimethyl sulfoxide (DMSO) were supplied by R&M Chemicals (Essex, UK). Reagents for wound healing assays such as trypsin and penicillin-streptomycin were bought from Nacalai Tesque (Kyoto, Japan), whereas phosphate buffer saline (PBS) and Dulbecco’s modified eagle medium (DMEM) were obtained from Solarbio (Beijing, China). Fetal bovine serum (FBS) and cell counting kit (CCK-8) were supplied by Tico Europe (Amstelveen, Netherlands) and Dojindo Lab (Kumamoto, Japan), respectively. The commercial standards of vitexin, orientin, isovitexin, and isoorientin with a purity of more than 98% were provided by ChemFaces (Wuhan, China). 4.2. Sampling and Preparation of Extracts The leaflets of Elaeis guineesis Jacq. were harvested from the University Agricultural Park, UPM, Malaysia. A piece of leaflet was sent for authentication to the Institute of Bioscience, UPM with the specimen number SK 3332/1. To ensure uniformity in samples, the sampling method routinely practiced by agronomists was followed [68]. The leaflets were obtained from the 17th frond from the apex of the tree. To minimize variation, the leaflets were harvested in the morning (8.00 am). A total of 24 uniform-sized leaflets were detached from the petiole of the oil palm frond and labelled accordingly, as shown in Figure 7. The leaflets were cleaned with water, their rachises were removed, and then were cut into small pieces. The leaflets were then oven-dried at 35 °C, until a constant weight was reached. The dried leaflets were ground and sieved to a fine powder with uniform size. The OPL powdered samples were extracted with four different solvent systems—i.e., ethyl acetate (100%), ethyl acetate-methanol (1:1), absolute methanol (100%), and methanol-water (4:1). For each of the solvent system, 0.1 g of OPL powder, placed in a 50 mL centrifuge tube, was mixed with 5 mL solvent, and vortexed for 0.5 min at the speed of at 3000 rpm. The extraction was performed using an ultrasonication method for 30 min at 25 °C in ultrasonic bath with frequency of 40 Hz. To obtain the crude dried extract, the mixture was centrifuged (4000× rpm, 15 min), and the supernatant was concentrated using vacuum evaporation and freeze-dried (Kansas City, MO, USA). A total of six replicates for each solvent system was prepared. 4.3. NMR Sample Preparation, Measurement, and Spectra Processing The sample preparation and NMR measurements were executed according to established protocol [40] with some modifications. For each sample, a 20 mg weight of the extract was placed into an Eppendorf tube, and mixed with an equal volume (300 µL) of deuterated methanol and prepared buffer (pH 6.0 and 0.01% TSP). Prior to analysis, the mixture was homogenized by an ultrasonication technique (15 min, 25 °C) and subsequently centrifuged for 15 min at a speed of 4000 rpm. For 1H-NMR spectral acquisition, 550 µL of filtered extract solution was pipetted into a clean NMR tube. The spectra were retrieved from a Varian INOVA NMR spectrometer (spectrometer (Varian Inc., Palo Alto, CA, USA)) with a frequency of 500 MHz. For each 1D spectrum, the conditions of the instrument were set at 24 °C with 64 scans, a 3.54 min acquisition time and 2 s relaxation delay. To suppress the water signal, a pre-saturation (PRESAT) pulse sequence was employed. The conditions produced a spectral width of −2 to 14 ppm. The 1H-NMR spectra of all extracts were processed using Chenomx software (Edmonton, AB, Canada). The spectral region (0.50–10.00 ppm) was bucketed into 243 integrated regions with δ 0.04 bin size. The signals of residual water (δ 4.70–4.88 ppm) and methanol (δ 3.24–3.33 ppm) were excluded during binning. To assist peak assignment, the 2D NMR experiment (J-resolved) was performed. For further confirmation, the assigned signals of orientin, vitexin, isoorientin, and isovitexin were verified with commercial standards (Figure S2). The data processing was carried out using MestRenova (Mestrelab Research, Santiago de Compostella, Spain). The identified primary and secondary metabolites were relatively quantified by plotting their respective mean peak area of the assigned signals observed in the 1H-NMR spectra. The discriminatory metabolites were selected based on the NMR signals which were not overlapping with other signals for the relative quantification, performed using SPSS software (SPSS Inc., Chicago, IL, USA). 4.4. Total Phenolic (TPC) and Flavonoid Content (TFC) Determination The TPC assay was executed using Folin-Ciocalteu (F-C) approach [69], with slight adjustments. Briefly, 0.1 mg/mL extract solution was prepared and an aliquot of 20 µL was mixed with 100 µL F-C reagent and 80 µL sodium carbonate solution (7.5%) in each well of a 96-well microtiter plate. The reaction was allowed to happen within 30 min incubation prior to the absorbance measurement at 750 nm using a microplate reader (Tecan Group Ltd., Männedorf, Switzerland). The TPC assay was carried out in six replicates for each extract and polyphenolic contents of the samples were evaluated from a calibration curve constructed for a series of gallic acid concentrations. The TPC values were expressed in milligrams of gallic acid equivalents per gram of extract (mg GAE/g extract). The TFC assay was executed using an aluminum chloride complex forming method [25]. Briefly, 0.1 mg/mL of extract solution was prepared and an aliquot of 125 µL was mixed with 375 µL 95% ethanol, 25 µL aluminum chloride solution (10%), 25 µL sodium acetate solution (1M), and 700 µL water in a 2 mL Eppendorf tube. Following this, 200 µL of the mixture was transferred to 96-well microtiter plate and incubated for 40 min at 25 °C. The absorbance measurement was performed using a microplate reader at 415 nm. The TFC assay was performed in six replicates for each extract and the flavonoid contents were calculated from a calibration curve constructed for a series of quercetin concentrations. The TFC values were stated in milligrams of quercetin equivalents per gram of extract (mg QCE/g extract). 4.5. In Vitro Antioxidant Assays The antioxidant activity of the solvent extracts were performed by ferric reducing antioxidant power (FRAP), 2,2-diphenyl-1-picrylhydrazyl (DPPH), and nitric oxide (NO) free radical scavenging assays using the following reagents: 2,4,6-tris(2-pyridyl)-1,3,5-triazine (TPTZ), DPPH reagent, N-(1-naphthyl)ethylenediamine dihydrochloride, and sulphanilamide. Measurements of absorbances were recorded on a microplate reader. Quercetin was employed as positive control in each assay. FRAP assay: The solvent extracts were measured for FRAP according to a previously described method [70], with slight adjustments. To prepare 25 mL FRAP reagent, 2.5 mL TPTZ solution (10 mM) and 2.5 mL ferric chloride solution (20 mM) were mixed, and the remaining volume was topped up with 30 mM acetate buffer with pH adjustment to 3.6. Prior to use, the mixture was incubated for 15 min at 37 °C. For the FRAP assay, 20 µL of extract was mixed with 180 µL of prepared reagent and the absorbance measurement was read at 593 nm after a 30 min incubation. The FRAP assay was performed in six replicates for each extract and the quantifications were calculated from a calibration curve constructed for a series of ascorbic acid concentrations. The results were stated in micrograms of ascorbic acid equivalents per milligram of extract (µg AAE/g extract). DPPH free radical scavenging assay: The assay on the OPL extracts was executed by following our recent published method [25]. Briefly, 0.1 mg/mL extract solution was prepared and serially diluted (100 to 0.7 µg/mL) inside 96-well microtiter plate. Subsequently, 100 µL of 0.059 mg/mL DPPH solution was added and incubated for 30 min. The absorbance measurement was read at 515 nm. The extract was carried out in six replicates and the radical inhibitions were presented as IC50 value in µg/mL. The scavenging activity (SA) was determined by using this formula: SA % = [(Ao − As)/Ao] × 100%(1) where Ao and As are absorbances of the solvent blank and extract solution, respectively. NO free radical scavenging assay: The assay was performed using a 96-well plate format as described previously [25]. To prepare Griess reagent, 0.1 g sulphanilamide and 0.01 g N-(1-naphthyl)-ethylenediamine dihydrochloride were weighed and mixed with 10 mL phosphoric acid (2.5%). For the assay, 1 mg/mL of extract solution was prepared and serially diluted (1000 to 0.98 µg/mL). Subsequently, 60 µL of extract solution was mixed with 60 µL sodium nitroprusside solution and incubated for 150 min. The absorbance measurement was read at 550 nm after the addition of 60 µL prepared Griess reagent. The extract was carried out in six replicates and the radical inhibitions were presented as an IC50 value in µg/mL. Calculation of the SA value was carried out according to the same formula as for the DPPH assay. 4.6. In Vitro Wound Healing Assays The wound healing properties of the solvent extracts were evaluated by cell proliferation and scratch (cell migration) assays, performed as described previously [26]. Allantoin was treated as the positive control for both assays. The wound healing assays were performed by using Swiss 3T3 albino mouse fibroblasts (Cell Line Service, Appelheim, Germany). A cell culture medium, comprised of high glucose DMEM, 10% FBS, and 1% penicillin-streptomycin, was used. The cells were maintained in a humidified 5% CO2 incubator at 37 °C. Prior to seeding, the cells at 70–80% confluence were used for all the experiments. A total of 1.563 µg/mL of solvent extracts and allantoin were prepared by mixing with a cell culture medium and filtered through 0.22 µm syringe filter for sterilization purpose. Cell proliferation assay: The 3T3 cells were seeded in a 96-well microtiter plate at a density of 1 × 105 in 100 µL culture medium per well. After 24 h incubation, the medium was replaced with 100 µL of the extract solution and incubated for a further 48 h. After this period, 10 µL of CCK-8 reagent was transferred into the mixture and continued incubating for 4 h. The absorbance of the sample, positive and negative (cells without treatment), was recorded on a microplate reader at 450 nm, against a solvent blank. Each extract solution was performed in six replicates. Cell migration assay: The 3T3 cells were seeded in 24 well plates at a density of 3 × 105 cells/mL per well. By employing a scratch assay technique, after 24 h incubation, a linear line was scratched on the confluent monolayered cells with a sterile P200 pipette tip. Care was taken when making the scratch to obtain a uniform size and length on each well. After adequately rinsing the well to eradicate the debris, the extract solution was pipetted into the well as treatment. The migration of the cell was monitored for wound closure and photographed on using phase contrast microscope, with the magnification set to 5× after 0 and 48 h of incubation. Using Image-J software (version 1.45, the area of the cells covered on the scratch area was measured and expressed as a percentage (%) of wound closure. Each extract solution was performed in six replicates. 4.7. Statistical Analysis The data retrieved from polyphenolic contents (TPC and TFC), antioxidant (FRAP, DPPH, and NO scavenging activities), and wound healing properties (cell proliferation and cell migration) assays were analyzed using Minitab software (Minitab Inc, State College, PA, USA) and GraphPad Prism software (San Diego, CA, USA) where the final values were presented as mean ± standard deviation of six replicates. To determine the significant differences of each datum, a one-way analysis of variance (ANOVA) was employed. For MVDA, PCA and PLS were analyzed using SIMCA-P+ software (Umetrics AB, Umea, Sweden, version 13.0) with Pareto scaling where the proton NMR chemical shifts (metabolites) were used to represent X-variables, while the antioxidant (1/IC50) and wound healing results were treated as Y-variables. 5. Conclusions The variation of metabolites in OPL extracts prepared using aqueous methanol, absolute methanol, ethyl acetate-methanol, and ethyl acetate as the solvent systems for extraction, as well as the correlation of the identified metabolites with the antioxidant and wound healing properties of extracts, were evaluated using a 1H-NMR-based metabolomics approach. A total of 19 metabolites, comprising six flavonoids, eight organic acids, four carbohydrates, and an amine, were identified in the various solvent extracts. For extraction purposes, the methanolic solvent systems (ethyl acetate-methanol, absolute methanol, and aqueous methanol) were more effective for the extraction of both primary and secondary metabolites. Based on the strong correlations with the antioxidant and wound healing activities, it could be concluded that the flavonoids (+)-catechin, (−)-epicatechin, orientin, isoorientin, vitexin, and isovitexin are the main bioactive metabolites. Several organic acids and carbohydrates constituents also contributed to the bioactivities. The findings of this study have emphasized the importance of these natural antioxidant compounds to the wound healing potential of OPL. From the perspective of quality control and assurance of therapeutic efficacy, this information will be of great value for future exploitation of this agriculture biomass into plant-based products, specifically for wound treatment purposes. Acknowledgments The first author gratefully acknowledges Universiti Putra Malaysia for supporting his study under the Graduate Research Fellowship (GRF) scheme. Sample Availability: Samples of the compounds are available from the authors. Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. Supplementary Materials The following are available online. Figure S1: 2D NMR 1H (J-resolved) spectrum of methanolic extract of oil palm (Elaeis guineensis Jacq.) leaflet, Figure S2: 1H-NMR spectra of commercial standards of orientin (16), isoorientin (17), vitexin (18) and isovitexin (19). Signals appeared in the region between 6.0–8.0 and 3.0–5.5 ppm were assigned to protons in flavonoid aglycone and 6- or 8-C-glycoside, respectively. Letter a labeled on peak indicates proton on C-6′. Figure S3: PLS score plot constructed from 1H-NMR spectral data of four different solvent extracts of OPL. The plot ellipse represents 95% hoteling T2 confidence. Figure S4: Permutation test results (100 permutations) for validation of PLS models built for (A) DPPH free radical scavenging, (B) NO free radical scavenging, (C) FRAP, (D) cell proliferation, and (E) cell migration activities. Click here for additional data file. Author Contributions K.S. and M.S.C.Z. conceptualized and designed the experiments; K.S. provided the research materials; M.S.C.Z. performed the experiments and drafted the manuscript; K.S., S.F., and M.S.C.Z. analyzed the data; K.S., N.M.N., and S.Y.L. validated the data and edited the manuscript. All authors have read and agreed to the published version of the manuscript. Funding This research received no external funding. Conflicts of Interest The authors declare no conflict of interest to reveal. Figure 1 (A) Representative proton nuclear magnetic resonance (1H-NMR) spectra (500 MHz, methanol-d4) of different solvent extracts of oil palm leaflet (OPL). (B) Part of 1H-NMR spectra of ethyl acetate extract within range of 4.2 to 5.4 ppm. (C) Part of 1H-NMR spectra of aqueous methanolic extract within range of 6.5 to 8.1 ppm. Identified metabolites: 1, sucrose; 2, fructose; 3, α-glucose; 4, β-glucose; 5, choline; 6, citric acid; 7, fumaric acid; 8, palmitic acid; 9, succinic acid; 10, aconitic acid; 11, fatty acid; 12, α-linolenic acid; 13, acetic acid; 14, (−)-epicatechin; 15, (+)-catechin; 16, orientin; 17, isoorientin; 18, vitexin; 19, isovitexin. Figure 2 Major flavonoids detected in OPL extracts. Isomer 1, vitexin (18)/isovitexin (19); Isomer 2, orientin (16)/isoorientin (17), and Isomer 3, (+)-catechin (15)/(−)-epicatechin (14). Figure 3 Principal component analysis (PCA) score (A) and loading (B) plots constructed from 1H-NMR spectral data of four different solvent extracts of OPL. The plot ellipse represents 95% hoteling T2 confidence. Identified metabolites: 1, sucrose; 2, fructose; 3, α-glucose; 4, β-glucose; 5, choline; 6, citric acid; 7, fumaric acid; 8, palmitic acid; 9, succinic acid; 10, aconitic acid; 11, fatty acid; 12, α-linolenic acid; 13, acetic acid; 14, (−)-epicatechin; 15, (+)-catechin; 16, orientin; 17, isoorientin; 18, vitexin; 19, isovitexin. Figure 4 Relative quantification of identified compounds in different solvent extracts of OPL based on the mean peak areas of selected 1H-NMR signals. (A) Carbohydrates and amines; (B) organic acids; (C) flavonoids. The discriminating chemical shifts (ppm) employed for the relative quantification were sucrose (δ 3.82), fructose (δ 4.18), α-glucose (δ 5.18), β-glucose (δ 4.62), choline (δ 3.20), citric acid (δ 2.74), succinic acid (δ 2.54), aconitic acid (δ 3.42), α-linolenic acid (δ 1.30), acetic acid (δ 2.07), palmitic acid (δ 0.9), fatty acid (δ 1.34), (+)-catechin (δ 7.02), (−)-epicatechin (δ 6.94), orientin (δ 7.54), isoorientin (δ 7.40), vitexin (δ 8.02), and isovitexin (δ 7.79). Bars labelled with the different letters (a,b,c,d) are significantly different between solvent extracts for the same metabolites (p < 0.05). Figure 5 Partial least square (PLS) loading biplot of bioactivities represented by ferric reducing antioxidant power (FRAP), 2,2-diphenyl-1-picrylhydrazyl (DPPH), and nitric oxide (NO) free radical scavenging activities, cell proliferation, and cell migration. Identified metabolites: 1, sucrose; 2, fructose; 3, α-glucose; 4, β-glucose; 5, choline; 6, citric acid; 7, fumaric acid; 8, palmitic acid; 9, succinic acid; 10, aconitic acid; 11, fatty acid; 12, α-linolenic acid; 13, acetic acid; 14, (−)-epicatechin; 15, (+)-catechin; 16, orientin; 17, isoorientin; 18, vitexin; 19, isovitexin. Figure 6 Schematic diagram of hydrogen transfer mechanism in (+)-catechin. (A): (+)-Catechin; (B): Hydrogen bond stabilized radical; (C): Stabilized quinone structure. Figure 7 Illustration of sampling method of the 24 mature leaflets from 17th oil palm frond of Elaeis guineensis Jacq. Six leaflets were sampled for each solvent system used for extraction. molecules-25-05636-t001_Table 1Table 1 Metabolites identified in different solvent extracts of OPL and key 1H-NMR chemical shifts (500 MHz, methanol-d4). No Metabolites 1H-NMR Chemical Shift (Multiplicity) Aqueous Methanol Absolute Methanol Ethyl Acetate-Methanol Ethyl Acetate Reference 1 Sucrose 5.42 (d); 3.82 (m) + + + + [21,27,30,38,39] 2 Fructose 4.18 (d) + + + + [21,29,39] 3 α-glucose 5.18 (d) + + + + [21,27,30,39] 4 β-glucose 4.62 (d) + + + + [21,29,30,39] 5 Choline 3.20 (s) + + + + [21,31,34,39] 6 Citric acid 2.74 (d) + + + + [29,30,39] 7 Fumaric acid 6.54 (s) + + + - [30,31,32] 8 Palmitic acid 1.66 (m); 0.90 (t) + + + + [27] 9 Succinic acid 2.54 (s) + + + + [30,38] 10 Aconitic acid 3.42 (s) + + + + [33] 11 Fatty acid 1.34 (m) + + + + [21,39,49] 12 α-linolenic acid 1.30 (brs); 2.35 (t) + + + + [27,29,30] 13 Acetic acid 2.07 (s) + + + + [35,36,37] 14 (−)-epicatechin 7.04 (s); 6.95 (s) + + + + [9,25,50,51] 15 (+)-catechin 7.01 (d) + + + - [9,25,50,51] 16 Orientin 7.54 (s); 6.86 (s); 6.63 (s); 3.00–4.00 (m) + + + - [8,25,26,32,39,52] 17 Isoorientin 7.40 (d); 6.89 (d); 6.67 (s); 6.48 (s); 3.00–4.00 (m); 1.96 (s) + + + - [8,25,26,32,39,52] 18 Vitexin 8.02 (d); 6.89 (d); 6.77 (s); 3.00–4.00 (m); 1.17 (t) + + + - [8,25,26,32,39,52] 19 Isovitexin 7.79 (d); 6.88 (d); 6.52 (s); 3.00–4.00 (m); 1.17 (t) + + + - [8,25,26,32,39,52] (+)—presence, (-)—absence of an identified metabolite in the OPL extract. molecules-25-05636-t002_Table 2Table 2 Phenolic contents, antioxidant, and wound healing properties of different OPL extracts. Solvent Systems Polyphenolic Contents Antioxidant Activities Wound Healing Properties TPC (mg GAE/g) TFC (mg QCE/g) FRAP (mg AAE/g) DPPH IC50 (µg/mL) NO IC50 (µg/mL) Proliferation (%) Migration (%) Ethyl acetate 121.71 ± 32.78 a 5.94 ± 3.58 a 16.26 ± 6.65 a 67.92 ± 14.16 a 213.34 ± 58.14 a 94.97 ± 2.12 a 64.34 ± 1.68 a Ethyl acetate-methanol 174.19 ± 32.40 b 121.48 ± 6.67 b 94.00 ± 23.94 b 10.74 ± 3.51 b 109.08 ± 29.61 b 98.41 ± 4.07 a 85.83 ± 3.25 b Absolute methanol 213.08 ± 41.61 b 135.40 ± 9.76 b 71.62 ± 21.05 b 6.54 ± 3.31 b,c 67.64 ± 21.97 c 100.13 ± 1.07 a 88.56 ± 6.94 b Aqueous methanol 393.61 ± 36.11 c 129.72 ± 8.7 b 101.48 ± 16.67 b 3.53 ± 1.30 c 18.77 ± 3.37 d 107.7 ± 13.11 a 93.34 ± 4.35 b Values labelled with the different letters (a,b,c,d) are significantly different between solvent extracts in the same tested parameters (p < 0.05). molecules-25-05636-t003_Table 3Table 3 Pearson’s coefficient of antioxidant and wound healing activities of different OPL extracts. Activity TPC TFC DPPH NO FRAP Proliferation Antioxidant TFC 0.605 DPPH 0.654 0.997 * NO 0.869 0.914 0.933 FRAP 0.705 0.921 0.945 0.887 Wound healing Proliferation 0.995 * 0.677 0.723 0.911 0.762 Migration 0.761 0.976 * 0.988 * 0.976 * 0.947 0.820 Values with asterisk (*) were significant at α = 0.05 based on Pearson’s correlation coefficient’s method. ==== Refs References 1. Malaysian Palm Oil Board (MPOB) Available online: http://www.mpob.gov.my/ (accessed on 25 April 2018) 2. Setiabudi H.D. Jusoh R. Suhaimi S.F.R.M. Masrur S.F. Adsorption of methylene blue onto oil palm (Elaeis guineensis ) leaves: Process optimization, isotherm, kinetics and thermodynamic studies J. Taiwan Inst. Chem. Eng. 2016 63 363 370 10.1016/j.jtice.2016.03.035 3. Rasat M.S.M. Wahab R. Yunus A.A.M. Moktar J. Ramle S.F.M. Kari Z.A. Yusoff M. Physical And Mechanical Properties of Bio-Composite Board from Compressed Oil Palm Fronds Adv. Nat. Appl. Sci. 2013 7 572 582 4. Haron K. Brookes P.C. Anderson J.M. Zakaria Z.Z. Microbial biomass and soil organic matter dynamics in oil palm (Elaeis guineensis jacq.) plantations, West Malaysia Soil. Biol. Biochem. 1998 30 547 552 10.1016/S0038-0717(97)00217-4 5. Khosravihaftkhany S. Morad N. Abdullah A.Z. Teng T.T. Ismail N. Biosorption of Pb (ii) and Fe (iii) from aqueous co-solutions using chemically pretreated oil palm fronds RSC Adv. 2015 5 106498 106508 10.1039/C5RA15325E 6. Mohd Kassim A.S. Aripin A.M. Ishak N. Zainulabidin M.H. Abang Zaidel D.N.F. Oil palm leaf fibre and its suitability for paper-based products ARPN J. Eng. Appl. Sci. 2016 11 7364 7369 7. Ghani A.A.A. Rusli N.D. Shahudin M.S. Goh Y.M. Zamri-Saad M. Hafandi A. Hassim H.A. Utilisation of oil palm fronds as ruminant feed and its effect on fatty acid metabolism Pertanika J. Trop. Agric. Sci. 2017 40 215 224 8. Tahir N.I. Shaari K. Abas F. Parveez G.K.A. Ishak Z. Ramli U.S. Characterization of apigenin and luteolin derivatives from oil palm (Elaeis guineensis Jacq.) Leaf using LC-ESI-MS/MS J. Agric. Food Chem. 2012 60 11201 11210 10.1021/jf303267e 23116142 9. Tahir N.I. Shaari K. Abas F. Ahmad Parveez G.K. Ahmad Tarmizi H. Ramli U.S. Identification of oil palm (Elaeis guineensis ) spear leaf metabolites using mass spectrometry and neutral loss analysis J. Oil Palm Res. 2013 25 72 83 10. Mushtaq M.Y. Choi Y.H. Verpoorte R. Wilson E.G. Extraction for metabolomics: Access to the metabolome Phytochem. Anal. 2014 25 291 306 10.1002/pca.2505 24523261 11. Bucar F. Wube A. Schmid M. Natural product isolation—How to get from biological material to pure compounds Nat. Prod. Rep. 2013 30 525 545 10.1039/c3np20106f 23396532 12. Hoste H. Torres-Acosta J.F.J. Sandoval-Castro C.A. Mueller-Harvey I. Sotiraki S. Louvandini H. Thamsborg S.M. Terrill T.H. Tannin containing legumes as a model for nutraceuticals against digestive parasites in livestock Vet. Parasitol. 2015 212 5 17 10.1016/j.vetpar.2015.06.026 26190131 13. Katritzky A.R. Fara D.C. Yang H. Tämm K. Tamm T. Karelson M. Quantitative Measures of Solvent Polarity Chem. Rev. 2004 104 175 198 10.1021/cr020750m 14719974 14. Azmir J. Zaidul I.S.M. Rahman M.M. Sharif K.M. Mohamed A. Sahena F. Jahurul M.H.A. Ghafoor K. Norulaini N.A.N. Omar A.K.M. Techniques for extraction of bioactive compounds from plant materials: A review J. Food Eng. 2013 117 426 436 10.1016/j.jfoodeng.2013.01.014 15. Jones W.P. Kinghorn A.D. Extraction of Plant Secondary Metabolites Natural Products Isolation 2nd ed. Sarker S.D. Latif Z. Human Press Inc. Totowa, NJ, USA 2006 Volume 20 323 351 16. Leiss K.A. Choi Y.H. Verpoorte R. Klinkhamer P.G.L. An overview of NMR-based metabolomics to identify secondary plant compounds involved in host plant resistance Phytochem. Rev. 2011 10 205 216 10.1007/s11101-010-9175-z 21765818 17. Verpoorte R. Choi Y.H. Kim H.K. Metabolomics: What’s new? Flavour Fragr. J. 2010 25 128 131 10.1002/ffj.1982 18. Heyman H.M. Meyer J.J.M. NMR-based metabolomics as a quality control tool for herbal products S. Afr. J. Bot. 2012 82 21 32 10.1016/j.sajb.2012.04.001 19. Martin A.C. Pawlus A.D. Jewett E.M. Wyse D.L. Angerhofer C.K. Hegeman A.D. Evaluating solvent extraction systems using metabolomics approaches RSC Adv. 2014 4 26325 26334 10.1039/C4RA02731K 20. Hernández-Bolio G.I. Kutzner E. Eisenreich W. de Jesús Torres-Acosta J.F. Peña-Rodríguez L.M. The use of 1 H-NMR Metabolomics to Optimise the Extraction and Preliminary Identification of Anthelmintic Products from the Leaves of Lysiloma latisiliquum Phytochem. Anal. 2017 29 413 420 10.1002/pca.2724 28895238 21. Lee S.Y. Abas F. Khatib A. Ismail I.S. Shaari K. Zawawi N. Metabolite profiling of Neptunia oleracea and correlation with antioxidant and α-glucosidase inhibitory activities using 1 H-NMR-based metabolomics Phytochem. Lett. 2016 16 23 33 10.1016/j.phytol.2016.02.014 22. Ismail S.N. Maulidiani M. Akhtar M.T. Abas F. Ismail I.S. Khatib A. Ali N.A.M. Shaari K. Discriminative Analysis of Different Grades of Gaharu (Aquilaria malaccensis Lamk.) via 1 H-NMR-Based Metabolomics Using PLS-DA and Random Forests Classification Models Molecules 2017 22 1612 10.3390/molecules22101612 23. Vargas L.H.G. Neto J.C.R. de Aquino Ribeiro J.A. Ricci-Silva M.E. Souza M.T. Rodrigues C.M. de Oliveira A.E. Abdelnur P.V. Metabolomics analysis of oil palm (Elaeis guineensis ) leaf: Evaluation of sample preparation steps using UHPLC-MS/MS Metabolomics 2016 12 153 10.1007/s11306-016-1100-z 24. Rozali N.L. Yarmo M.A. Idris A.S. Kushairi A. Ramli U.S. Metabolomics differentiation of oil palm (Elaeis guineensis Jacq.) spear leaf with contrasting susceptibility to Ganoderma boninense Plant Omics. 2017 10 45 52 10.21475/poj.10.02.17.pne364 25. Zain M.S.C. Lee S.Y. Teo C.Y. Shaari K. Adsorption and Desorption Properties of Total Flavonoids from Oil Palm (Elaeis guineensis Jacq.) Mature Leaf on Macroporous Adsorption Resins Molecules 2020 25 778 10.3390/molecules25040778 32054137 26. Zain M.S.C. Lee S.Y. Sarian M.N. Fakurazi S. Shaari K. In vitro wound healing potential of flavonoid c-glycosides from oil palm (Elaeis guineensis Jacq.) leaves on 3t3 fibroblast cells Antioxidants 2020 9 326 10.3390/antiox9040326 32316665 27. Azizan A. Ahamad Bustamam M.S. Maulidiani M. Shaari K. Ismail I.S. Nagao N. Abas F. Metabolite profiling of the microalgal diatom chaetoceros calcitrans and correlation with antioxidant and nitric oxide inhibitory Activities via1H NMR-Based Metabolomics Mar. Drugs 2018 16 154 10.3390/md16050154 29735927 28. Pariyani R. Ismail I.S. Ahmad Azam A. Abas F. Shaari K. Identification of the compositional changes in Orthosiphon stamineus leaves triggered by different drying techniques using 1H NMR metabolomics J. Sci. Food Agric. 2017 97 4169 4179 10.1002/jsfa.8288 28233369 29. Azizan A. Xin L.A. Abdul Hamid N.A. Maulidiani M. Mediani A. Abdul Ghafar S.Z. Zulaikha Zolkeflee N.K. Abas F. Potentially bioactive metabolites from pineapple waste extracts and their antioxidant and α-glucosidase inhibitory activities by 1 H-NMR Foods 2020 9 173 10.3390/foods9020173 32053982 30. Verpoorte R. Choi Y.H. Kim H.K. NMR-based metabolomics at work in phytochemistry Phytochem. Rev. 2007 6 3 14 10.1007/s11101-006-9031-3 31. Georgiev M.I. Ali K. Alipieva K. Verpoorte R. Choi Y.H. Metabolic differentiations and classification of Verbascum species by NMR-based metabolomics Phytochemistry 2011 72 2045 2051 10.1016/j.phytochem.2011.07.005 21807390 32. Khoo L.W. Kow A.S.F. Maulidiani M. Ang M.Y. Chew W.Y. Lee M.T. Tan C.P. Shaari K. Tham C.L. Abas F. 1 H-NMR metabolomics for evaluating the protective effect of Clinacanthus nutans (Burm. f) Lindau water extract against nitric oxide production in LPS-IFN-γ activated RAW 264.7 macrophages Phytochem. Anal. 2019 30 46 61 10.1002/pca.2789 30183131 33. Wishart D.S. Knox C. Guo A.C. Eisner R. Young N. Gautam B. Hau D.D. Psychogios N. Dong E. Bouatra S. HMDB: A knowledgebase for the human metabolome Nucleic Acids Res. 2009 37 D603 D610 10.1093/nar/gkn810 18953024 34. Mediani A. Abas F. Khatib A. Maulidiani H. Shaari K. Choi Y.H. Lajis N.H. 1 H-NMR-based metabolomics approach to understanding the drying effects on the phytochemicals in Cosmos caudatus Food Res. Int. 2012 49 763 770 10.1016/j.foodres.2012.09.022 35. Li Y. He Q. Geng Z. Du S. Deng Z. Hasi E. NMR-based metabolomic profiling of Peganum harmala L. reveals dynamic variations between different growth stages R. Soc. Open Sci. 2018 5 171722 10.1098/rsos.171722 30109044 36. Lu Y. Hu F. Miyakawa T. Tanokura M. Complex mixture analysis of organic compounds in yogurt by NMR spectroscopy Metabolites 2016 6 19 10.3390/metabo6020019 27322339 37. Fan S. Zhong Q. Fauhl-Hassek C. Pfister M.K.H. Horn B. Huang Z. Classification of Chinese wine varieties using 1H NMR spectroscopy combined with multivariate statistical analysis Food Control 2018 88 113 122 10.1016/j.foodcont.2017.11.002 38. Sulaiman F. Ahmad Azam A. Ahamad Bustamam M.S. Fakurazi S. Abas F. Lee Y.X. Ismail A.A. Mohd Faudzi S.M. Ismail I.S. Metabolite Profiles of Red and Yellow Watermelon (Citrullus lanatus ) Cultivars Using a 1 H-NMR Metabolomics Approach Molecules 2020 25 3235 10.3390/molecules25143235 39. Khoo L.W. Mediani A. Zolkeflee N.K.Z. Leong S.W. Ismail I.S. Khatib A. Shaari K. Abas F. Phytochemical diversity of Clinacanthus nutans extracts and their bioactivity correlations elucidated by NMR based metabolomics Phytochem. Lett. 2015 14 123 133 10.1016/j.phytol.2015.09.015 40. Kim H.K. Choi Y.H. Verpoorte R. NMR-based metabolomic analysis of plants Nat. Protoc. 2010 5 536 549 10.1038/nprot.2009.237 20203669 41. Razali M.T.A. Zainal Z.A. Maulidiani M. Shaari K. Zamri Z. Mohd Idrus M.Z. Khatib A. Abas F. Ling Y.S. Rui L.L. Classification of raw stingless bee honeys by bee species origins using the NMR- and LC-MS-based metabolomics approach Molecules 2018 23 2160 10.3390/molecules23092160 42. Pramai P. Abdul Hamid N.A. Mediani A. Maulidiani M. Abas F. Jiamyangyuen S. Metabolite profiling, antioxidant, and α-glucosidase inhibitory activities of germinated rice: Nuclear magnetic resonance based metabolomics study J. Food Drug Anal. 2016 26 47 57 10.1016/j.jfda.2016.11.023 43. Choi Y.H. Kim H.K. Linthorst H.J. Hollander J.G. Lefeber A.W. Erkelens C. Nuzillard J.M. Verpoorte R. NMR metabolomics to revisit the tobacco mosaic virus infection in Nicotiana tabacum leaves J. Nat. Prod. 2006 69 742 748 10.1021/np050535b 16724833 44. Ganbaatar C. Gruner M. Mishig D. Duger R. Schmidt A.W. Knölker H.-J. Flavonoid glycosides from the aerial parts of Polygonatum odoratum (Mill.) Druce growing in Mongolia Open Nat. Prod. J. 2015 8 1 7 10.2174/1874848101508010001 45. Yoon H. Eom S.L. Hyun J.Y. Jo G.H. Hwang D.S. Lee S.H. Yong Y.J. Lee Y.H. Lim Y.H. Park J.C. 1 H- and 13 C-NMR data on hydroxy/methoxy flavonoids and the effects of substituents on chemical shifts Bull Korean Chem. Soc. 2011 32 2101 2104 10.5012/bkcs.2011.32.6.2101 46. NMR Predict Available online: http://www.nmrdb.org (accessed on 20 October 2020) 47. Human Metabolome Database Available online: http://www.hmdb.ca/ (accessed on 20 October 2020) 48. PubChem Available online: https://pubchem.ncbi.nlm.nih.gov/ (accessed on 20 October 2020) 49. Wan-Nadilah W.A. Akhtar M.T. Shaari K. Khatib A. Hamid A.A. Hamid M. Variation in the metabolites and α-glucosidase inhibitory activity of Cosmos caudatus at different growth stages BMC Complement. Altern. Med. 2019 19 245 10.1186/s12906-019-2655-9 31488132 50. Berregi I. Santos J.I. Del Campo G. Miranda J.I. Quantitative determination of (-)-epicatechin in cider apple juices by 1H NMR Talanta 2003 61 139 145 10.1016/S0039-9140(03)00236-4 18969172 51. Abdul-Hamid N.A. Abas F. Ismail I.S. Shaari K. Lajis N.H. Influence of Different Drying Treatments and Extraction Solvents on the Metabolite Profile and Nitric Oxide Inhibitory Activity of Ajwa Dates J. Food Sci. 2015 80 H2603 H2611 10.1111/1750-3841.13084 26457883 52. Cuc N.T. Nhiem N.X. Yen P.H. Anh H.L.T. Van Minh C. Van Kiem P. Flavonoid glygosides from Antidesma ghaesembilla Vietnam J. Chem. 2015 53 94 97 53. Grasso D. Chin Y.P. Weber W.J. Jr. Structural and behavioral characterictiscs of a commercial humic acid and natural dissolved aquatic organic matter Chemosphere 1990 21 1181 1197 10.1016/0045-6535(90)90139-K 54. Ofori-Boateng C. Lee K.T. Sustainable utilization of oil palm wastes for bioactive phytochemicals for the benefit of the oil palm and nutraceutical industries Phytochem. Rev. 2013 12 173 190 10.1007/s11101-013-9270-z 55. Gothai S. Arulselvan P. Tan W. Fakurazi S. Wound healing properties of ethyl acetate fraction of Moringa oleifera in normal human dermal fibroblasts J. Intercult. Ethnopharmacol. 2016 5 1 6 10.5455/jice.20160201055629 27069722 56. Muhammad A.A. Aimi N. Pauzi S. Arulselvan P. Abas F. Fakurazi S. In Vitro Wound Healing Potential and Identification of Bioactive Compounds from Moringa oleifera Lam BioMed Res. Int. 2013 2013 974580 10.1155/2013/974580 24490175 57. Süntar I. Akkol E.K. Nahar L. Sarker S.D. Wound healing and antioxidant properties: Do they coexist in plants? Free Radic. Antioxid. 2012 2 1 7 10.5530/ax.2012.2.2.1 58. Abdul-Hamid N.A. Abas F. Ismail I.S. Tham C.L. Maulidiani M. Mediani A. Swarup S. Umashankar S. Zolkeflee N.K.Z. Metabolites and biological activities of Phoenix dactylifera L. pulp and seeds: A comparative MS and NMR based metabolomics approach Phytochem. Lett. 2019 31 20 32 10.1016/j.phytol.2019.03.004 59. Zhang B. Xia T. Duan W. Zhang Z. Li Y. Fang B. Xia M. Wang M. Effects of organic acids, amino acids and phenolic compounds on antioxidant characteristic of Zhenjiang aromatic vinegar Molecules 2019 24 3799 10.3390/molecules24203799 31652563 60. Zafar M.S. Quarta A. Marradi M. Ragusa A. Recent developments in the reduction of oxidative stress through antioxidant polymeric formulations Pharmaceutics 2019 11 505 10.3390/pharmaceutics11100505 61. Rezaeizadeh A. Some Mechanisms for Wound Healing Enhancing Effects of Bitter Gourd (Momordica charantia L.) Extract in Diabetic Rats Ph.D. Thesis Universiti Putra Malaysia Serdang, Malaysia 7 2013 62. Jorge M.P. Madjarof C. Ruiz A.L.T.G. Fernandes A.T. Rodrigues R.A.F. de Oliveira Sousa I.M. Foglio M.A. de Carvalho J.E. Evaluation of wound healing properties of Arrabidaea chica Verlot extract J. Ethnopharmacol. 2008 118 361 366 10.1016/j.jep.2008.04.024 18573628 63. Schäfer M. Werner S. Oxidative stress in normal and impaired wound repair Pharmacol. Res. 2008 58 165 171 10.1016/j.phrs.2008.06.004 18617006 64. Treml J. Šmejkal K. Flavonoids as Potent Scavengers of Hydroxyl Radicals Compr. Rev. Food Sci. Food Saf. 2016 15 720 738 10.1111/1541-4337.12204 65. Glevitzky I. Dumitrel G.A. Glevitzky M. Pasca B. Otrisal P. Bungau S. Cioca G. Pantis C. Popa M. Statistical Analysis of the Relationship Between Antioxidant Activity and the Structure of Flavonoid Compounds Rev. Chim. 2019 70 3103 3107 10.37358/RC.19.9.7497 66. Rice-Evans C.A. Miller N.J. Paganga G. Structure-antioxidant activity relationships of flavonoids and phenolic acids Free Radic Biol. Med. 1996 20 933 956 10.1016/0891-5849(95)02227-9 8743980 67. Amic D. Davidovic-Amic D. Beslo D. Rastija V. Lucic B. Trinajstic N. SAR and QSAR of the Antioxidant Activity of Flavonoids Curr. Med. Chem. 2007 14 827 845 10.2174/092986707780090954 17346166 68. Corley R.H.V. Tinker P.B.H. The Oil Palm 4th ed. Black Science Bath, UK 2003 0470750367 69. Lee S.Y. Mediani A. Maulidiani M. Khatib A. Ismail I.S. Zawawi N. Abas F. Comparison of partial least squares and random forests for evaluating relationship between phenolics and bioactivities of Neptunia oleracea J. Sci. Food Agric. 2018 98 240 252 10.1002/jsfa.8462 28580581 70. Benzie I.F.F. Strain J.J. The ferric reducing ability of plasma (FRAP) as a measure of “antioxidant power”: The FRAP assay Anal. Biochem. 1996 239 70 76 10.1006/abio.1996.0292 8660627