Next Article in Journal
Effect of Potato Pulp Pectic Polysaccharide on the Stability of Acidified Milk Drinks
Next Article in Special Issue
Development of Advanced Chemometric-Assisted Spectrophotometric Methods for the Determination of Cromolyn Sodium and Its Alkaline Degradation Products
Previous Article in Journal
New Freeze-Dried Andean Blueberry Juice Powders for Potential Application as Functional Food Ingredients: Effect of Maltodextrin on Bioactive and Morphological Features
Previous Article in Special Issue
Comprehensive Review on Application of FTIR Spectroscopy Coupled with Chemometrics for Authentication Analysis of Fats and Oils in the Food Products
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Molecules
Volume 25
Issue 23
10.3390/molecules25235636
molecules-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Metabolite Characterization and Correlations with Antioxidant and Wound Healing Properties of Oil Palm (Elaeis guineensis Jacq.) Leaflets via 1H-NMR-Based Metabolomics Approach

by
Mohamad Shazeli Che Zain
1,
Soo Yee Lee
1,
Nadiah Mad Nasir
2,
Sharida Fakurazi
3 and
Khozirah Shaari
1,2,*
1
Natural Medicines and Products Research Laboratory (NaturMeds), Institute of Bioscience, Universiti Putra Malaysia, Serdang 43400, Malaysia
2
Department of Chemistry, Faculty of Science, Universiti Putra Malaysia, Serdang 43400, Malaysia
3
Laboratory of Vaccines and Biomolecules (VacBio), Institute of Bioscience, Universiti Putra Malaysia, Serdang 43400, Malaysia
*
Author to whom correspondence should be addressed.
Molecules 2020, 25(23), 5636; https://doi.org/10.3390/molecules25235636
Submission received: 14 November 2020 / Revised: 25 November 2020 / Accepted: 27 November 2020 / Published: 30 November 2020
(This article belongs to the Special Issue Chemometrics Tools Used in Analytical Chemistry)
Browse Figures
Versions Notes

Abstract

:
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.

Graphical Abstract

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%
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.

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.

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.

Acknowledgments

The first author gratefully acknowledges Universiti Putra Malaysia for supporting his study under the Graduate Research Fellowship (GRF) scheme.

Conflicts of Interest

The authors declare no conflict of interest to reveal.

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. [Google Scholar] [CrossRef] [Green Version]
  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. [Google Scholar]
  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. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef]
  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. [Google Scholar]
  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. [Google Scholar]
  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. [Google Scholar] [CrossRef]
  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. [Google Scholar]
  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. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef] [Green Version]
  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. [Google Scholar] [CrossRef] [PubMed]
  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. [Google Scholar] [CrossRef] [PubMed]
  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. [Google Scholar] [CrossRef]
  15. Jones, W.P.; Kinghorn, A.D. Extraction of Plant Secondary Metabolites. In Natural Products Isolation, 2nd ed.; Sarker, S.D., Latif, Z., Eds.; Human Press Inc.: Totowa, NJ, USA, 2006; Volume 20, pp. 323–351. [Google Scholar]
  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. [Google Scholar] [CrossRef] [Green Version]
  17. Verpoorte, R.; Choi, Y.H.; Kim, H.K. Metabolomics: What’s new? Flavour Fragr. J. 2010, 25, 128–131. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef] [Green Version]
  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. [Google Scholar] [CrossRef]
  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 1H-NMR Metabolomics to Optimise the Extraction and Preliminary Identification of Anthelmintic Products from the Leaves of Lysiloma latisiliquum. Phytochem. Anal. 2017, 29, 413–420. [Google Scholar] [CrossRef]
  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 1H-NMR-based metabolomics. Phytochem. Lett. 2016, 16, 23–33. [Google Scholar] [CrossRef]
  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 1H-NMR-Based Metabolomics Using PLS-DA and Random Forests Classification Models. Molecules 2017, 22, 1612. [Google Scholar] [CrossRef] [Green Version]
  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. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  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. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  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. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  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. [Google Scholar] [CrossRef] [PubMed]
  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 1H-NMR. Foods 2020, 9, 173. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Verpoorte, R.; Choi, Y.H.; Kim, H.K. NMR-based metabolomics at work in phytochemistry. Phytochem. Rev. 2007, 6, 3–14. [Google Scholar] [CrossRef] [Green Version]
  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. [Google Scholar] [CrossRef]
  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. 1H-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. [Google Scholar] [CrossRef] [Green Version]
  33. Wishart, D.S.; Knox, C.; Guo, A.C.; Eisner, R.; Young, N.; Gautam, B.; Hau, D.D.; Psychogios, N.; Dong, E.; Bouatra, S.; et al. HMDB: A knowledgebase for the human metabolome. Nucleic Acids Res. 2009, 37, D603–D610. [Google Scholar] [CrossRef] [PubMed]
  34. Mediani, A.; Abas, F.; Khatib, A.; Maulidiani, H.; Shaari, K.; Choi, Y.H.; Lajis, N.H. 1H-NMR-based metabolomics approach to understanding the drying effects on the phytochemicals in Cosmos caudatus. Food Res. Int. 2012, 49, 763–770. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Lu, Y.; Hu, F.; Miyakawa, T.; Tanokura, M. Complex mixture analysis of organic compounds in yogurt by NMR spectroscopy. Metabolites 2016, 6, 19. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  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. [Google Scholar] [CrossRef]
  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 1H-NMR Metabolomics Approach. Molecules 2020, 25, 3235. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef]
  40. Kim, H.K.; Choi, Y.H.; Verpoorte, R. NMR-based metabolomic analysis of plants. Nat. Protoc. 2010, 5, 536–549. [Google Scholar] [CrossRef]
  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.; et al. Classification of raw stingless bee honeys by bee species origins using the NMR- and LC-MS-based metabolomics approach. Molecules 2018, 23, 2160. [Google Scholar] [CrossRef] [Green Version]
  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. [Google Scholar] [CrossRef] [Green Version]
  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. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef] [Green Version]
  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. 1H- and 13C-NMR data on hydroxy/methoxy flavonoids and the effects of substituents on chemical shifts. Bull Korean Chem. Soc. 2011, 32, 2101–2104. [Google Scholar] [CrossRef] [Green Version]
  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. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  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. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef]
  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. [Google Scholar]
  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. [Google Scholar] [CrossRef] [Green Version]
  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. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  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. [Google Scholar] [CrossRef] [Green Version]
  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. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  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. [Google Scholar] [CrossRef] [Green Version]
  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, July 2013. [Google Scholar]
  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. [Google Scholar] [CrossRef]
  63. Schäfer, M.; Werner, S. Oxidative stress in normal and impaired wound repair. Pharmacol. Res. 2008, 58, 165–171. [Google Scholar] [CrossRef]
  64. Treml, J.; Šmejkal, K. Flavonoids as Potent Scavengers of Hydroxyl Radicals. Compr. Rev. Food Sci. Food Saf. 2016, 15, 720–738. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef] [PubMed]
  68. Corley, R.H.V.; Tinker, P.B.H. The Oil Palm, 4th ed.; Black Science: Bath, UK, 2003; ISBN 0470750367. [Google Scholar]
  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. [Google Scholar] [CrossRef] [PubMed]
  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. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Molecules 25 05636 g001 550
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 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.
Molecules 25 05636 g001
Molecules 25 05636 g002 550
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 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).
Molecules 25 05636 g002
Molecules 25 05636 g003 550
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 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.
Molecules 25 05636 g003
Molecules 25 05636 g004 550
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 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).
Molecules 25 05636 g004
Molecules 25 05636 g005 550
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 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.
Molecules 25 05636 g005
Molecules 25 05636 g006 550
Figure 6. Schematic diagram of hydrogen transfer mechanism in (+)-catechin. (A): (+)-Catechin; (B): Hydrogen bond stabilized radical; (C): Stabilized quinone structure.
Figure 6. Schematic diagram of hydrogen transfer mechanism in (+)-catechin. (A): (+)-Catechin; (B): Hydrogen bond stabilized radical; (C): Stabilized quinone structure.
Molecules 25 05636 g006
Molecules 25 05636 g007 550
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.
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 g007
Table
Table 1. Metabolites identified in different solvent extracts of OPL and key 1H-NMR chemical shifts (500 MHz, methanol-d4).
Table 1. Metabolites identified in different solvent extracts of OPL and key 1H-NMR chemical shifts (500 MHz, methanol-d4).
NoMetabolites1H-NMR Chemical Shift (Multiplicity)Aqueous MethanolAbsolute MethanolEthyl Acetate-MethanolEthyl AcetateReference
1Sucrose5.42 (d); 3.82 (m)++++[21,27,30,38,39]
2Fructose4.18 (d)++++[21,29,39]
3α-glucose5.18 (d)++++[21,27,30,39]
4β-glucose4.62 (d)++++[21,29,30,39]
5Choline3.20 (s)++++[21,31,34,39]
6Citric acid2.74 (d)++++[29,30,39]
7Fumaric acid6.54 (s)+++-[30,31,32]
8Palmitic acid1.66 (m); 0.90 (t)++++[27]
9Succinic acid2.54 (s)++++[30,38]
10Aconitic acid3.42 (s)++++[33]
11Fatty acid1.34 (m)++++[21,39,49]
12α-linolenic acid1.30 (brs); 2.35 (t)++++[27,29,30]
13Acetic acid2.07 (s)++++[35,36,37]
14(−)-epicatechin7.04 (s); 6.95 (s)++++[9,25,50,51]
15(+)-catechin7.01 (d)+++-[9,25,50,51]
16Orientin7.54 (s); 6.86 (s); 6.63 (s); 3.00–4.00 (m)+++-[8,25,26,32,39,52]
17Isoorientin7.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]
18Vitexin8.02 (d); 6.89 (d); 6.77 (s); 3.00–4.00 (m); 1.17 (t)+++-[8,25,26,32,39,52]
19Isovitexin7.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.
Table
Table 2. Phenolic contents, antioxidant, and wound healing properties of different OPL extracts.
Table 2. Phenolic contents, antioxidant, and wound healing properties of different OPL extracts.
Solvent SystemsPolyphenolic ContentsAntioxidant ActivitiesWound 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 acetate121.71 ± 32.78 a5.94 ± 3.58 a16.26 ± 6.65 a67.92 ± 14.16 a213.34 ± 58.14 a94.97 ± 2.12 a64.34 ± 1.68 a
Ethyl acetate-methanol174.19 ± 32.40 b121.48 ± 6.67 b94.00 ± 23.94 b10.74 ± 3.51 b109.08 ± 29.61 b98.41 ± 4.07 a85.83 ± 3.25 b
Absolute methanol213.08 ± 41.61 b135.40 ± 9.76 b71.62 ± 21.05 b6.54 ± 3.31 b,c67.64 ± 21.97 c100.13 ± 1.07 a88.56 ± 6.94 b
Aqueous methanol393.61 ± 36.11 c129.72 ± 8.7 b101.48 ± 16.67 b3.53 ± 1.30 c18.77 ± 3.37 d107.7 ± 13.11 a93.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).
Table
Table 3. Pearson’s coefficient of antioxidant and wound healing activities of different OPL extracts.
Table 3. Pearson’s coefficient of antioxidant and wound healing activities of different OPL extracts.
Activity TPCTFCDPPHNOFRAPProliferation
AntioxidantTFC0.605
DPPH0.6540.997 *
NO0.8690.9140.933
FRAP0.7050.9210.9450.887
Wound healingProliferation0.995 *0.6770.7230.9110.762
Migration0.7610.976 *0.988 *0.976 *0.9470.820
Values with asterisk (*) were significant at α = 0.05 based on Pearson’s correlation coefficient’s method.
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.

Share and Cite

MDPI and ACS Style

Che Zain, M.S.; Lee, S.Y.; Nasir, N.M.; Fakurazi, S.; Shaari, K. Metabolite Characterization and Correlations with Antioxidant and Wound Healing Properties of Oil Palm (Elaeis guineensis Jacq.) Leaflets via 1H-NMR-Based Metabolomics Approach. Molecules 2020, 25, 5636. https://doi.org/10.3390/molecules25235636

AMA Style

Che Zain MS, Lee SY, Nasir NM, Fakurazi S, Shaari K. Metabolite Characterization and Correlations with Antioxidant and Wound Healing Properties of Oil Palm (Elaeis guineensis Jacq.) Leaflets via 1H-NMR-Based Metabolomics Approach. Molecules. 2020; 25(23):5636. https://doi.org/10.3390/molecules25235636

Chicago/Turabian Style

Che Zain, Mohamad Shazeli, Soo Yee Lee, Nadiah Mad Nasir, Sharida Fakurazi, and Khozirah Shaari. 2020. "Metabolite Characterization and Correlations with Antioxidant and Wound Healing Properties of Oil Palm (Elaeis guineensis Jacq.) Leaflets via 1H-NMR-Based Metabolomics Approach" Molecules 25, no. 23: 5636. https://doi.org/10.3390/molecules25235636

APA Style

Che Zain, M. S., Lee, S. Y., Nasir, N. M., Fakurazi, S., & Shaari, K. (2020). Metabolite Characterization and Correlations with Antioxidant and Wound Healing Properties of Oil Palm (Elaeis guineensis Jacq.) Leaflets via 1H-NMR-Based Metabolomics Approach. Molecules, 25(23), 5636. https://doi.org/10.3390/molecules25235636

Article Metrics

Back to TopTop