- Research article
- Open Access
- Open Peer Review
Antioxidants and α-glucosidase inhibitors from Neptunia oleracea fractions using 1H NMR-based metabolomics approach and UHPLC-MS/MS analysis
BMC Complementary and Alternative Medicinevolume 19, Article number: 7 (2019)
Neptunia oleracea is a plant cultivated as vegetable in Southeast Asia. Previous works have revealed the potential of this plant as a source of natural antioxidants and α-glucosidase inhibitors. Continuing our interest on this plant, the present work is focused in identification of the bioactive compounds from different polarity fractions of N. oleracea, namely hexane (HF), chloroform (CF), ethyl acetate (EF) and methanol (MF).
The N. oleracea fractions were obtained using solid phase extraction (SPE). A metabolomics approach that coupled the use of proton nuclear magnetic resonance (1H NMR) with multivariate data analysis (MVDA) was applied to distinguish the metabolite variations among the N. oleracea fractions, as well as to assess the correlation between metabolite variation and the studied bioactivities (DPPH free radical scavenging and α-glucosidase inhibitory activities). The bioactive fractions were then subjected to ultra-high performance liquid chromatography tandem mass spectrometry (UHPLC–MS/MS) analysis to profile and identify the potential bioactive constituents.
The principal component analysis (PCA) discriminated EF and MF from the other fractions with the higher distributions of phenolics. Partial least squares (PLS) analysis revealed a strong correlation between the phenolics and the studied bioactivities in the EF and the MF. The UHPLC-MS/MS profiling of EF and MF had tentatively identified the phenolics present. Together with some non-phenolic metabolites, a total of 37 metabolites were tentatively assigned.
The findings of this work supported that N. oleracea is a rich source of phenolics that can be potential antioxidants and α-glucosidase inhibitors for the management of diabetes. To our knowledge, this study is the first report on the metabolite-bioactivity correlation and UHPLC–MS/MS analysis of N. oleracea fractions.
Diabetes mellitus is a metabolically complicated disease in which the patient experiences a high blood glucose level. Over the past few decades, the occurrence of diabetes mellitus has increased drastically. Globally, the number of diabetic cases in adult population has increased from 108 million in 1980 to 422 million in 2014 . Among diabetic patients, the majority suffer from type 2 diabetes, which is associated with insulin insensitivity in the body tissue. The impaired glucose metabolism leads to excessive production of reactive oxygen species (ROS), enhances oxidative stress in the body and ultimately leads to complications such as cardiovascular disease, retinopathy, nephropathy and impaired wound healing. One of the therapeutic approaches for diabetes management is to slow down the absorption of postprandial glucose by inhibiting carbohydrate hydrolyzing enzymes, such as α-glucosidase . However, the currently available α-glucosidase inhibitors including acarbose, voglibose and miglitol produce gastrointestinal side effects . Hence, efforts are being made to search for antioxidants and α-glucosidase inhibitors from natural resources for the management of diabetes and its complications.
Neptunia oleracea is a plant cultivated as vegetable in Southeast Asia. It has a pleasant flavor and can be eaten raw or cooked. Proximate analysis of N. oleracea has shown that this plant is a good source of crude protein, crude fiber and important minerals, such as potassium and calcium . Furthermore, this plant has been reported to contain a high level of pro-vitamin A carotenoids, which are important in preventing vitamin A deficiency . Studies have revealed that this plant exhibited several biological activities, including antiulcer , hepatoprotective , antiviral , 5α-reductase inhibitory , analgesic and antiinflammatory  activities. Besides, to identify antioxidants and α-glucosidase inhibitors from natural resources, N. oleracea has been studied by our research group for its free radical scavenging and α-glucosidase inhibitory properties. The results have shown that this plant is a prominent source of free radical scavengers and α-glucosidase inhibitors; with phenolics as the potential candidates for these active compounds [11, 12].
Most of the previous studies of N. oleracea focused on the crude extracts and the metabolites contributing to these pharmacological properties were remained unassigned [6,7,8,9]. Often, the metabolite profile of crude extracts is complicated due to the large amount of metabolites that could be present, making the identification process challenging. One of the approach to facilitate the analysis and make the identification of the important metabolites more effective is to fractionate the crude extract by different polarity. Hence, in this study, the crude extract of N. oleracea was fractionated according to polarity using solid phase extraction (SPE). The metabolite variation among these fractions as well as the correlation between the metabolites and the studied bioactivities (DPPH free radical scavenging and α-glucosidase inhibitory activities) were assessed using a 1H NMR-based metabolomics approach. The potential bioactive constituents present in the bioactive fractions were profiled and identified using ultra-high performance liquid chromatography tandem mass spectrometry (UHPLC–MS/MS). These procedures help to highlight the bioactive fractions of N. oleracea for the DPPH free radical scavenging and α-glucosidase inhibitory activities, and to reveal the bioactive compounds, as continuation of our previous efforts on the investigation of N. oleracea as a source of antioxidant and α-glucosidase inhibitors.
Chemicals and reagents
Phosphate buffer, glycine, α-glucosidase enzyme, p-nitrophenyl-α-D-glucopyranose (PNPG) and 2,2-diphenyl-1-picrylhydrazyl (DPPH) were purchased from Sigma-Aldrich (Hamburg, Germany). Analytical grade solvents including absolute ethanol, hexane, chloroform, ethyl acetate and methanol were supplied by Merck (Darmstadt, Germany). The NMR chemicals including deuterated dimethyl sulfoxide-d6 (DMSO-d6) and trimethylsilyl propionic acid-d4 sodium salt (TSP) were also supplied by Merck (Darmstadt, Germany). The LC/MS grade acetonitrile and formic acid were purchased from Merck (Darmstadt, Germany). Water was prepared using a Milli-Q purification system.
Neptunia oleracea was planted in Universiti Putra Malaysia Agricultural Park by spreading the stems onto a pond. The plant was authenticated by an in-house Botanist (Dr. Shamsul Khamis) of the Institute of Bioscience, Universiti Putra Malaysia and the voucher specimen (SK2516/14) was deposited in the Herbarium Biodiversity Unit of the institute.
Sample drying and extraction
The optimized conditions of drying and extraction (freeze drying and sonication with absolute ethanol) were selected based on previous works on N. oleracea [11, 12]. The three-months-old plants were harvested. Immediately after the harvest, the leaves were separated from the stems of the plants. The leaves were cleaned with water and then kept in freezer at − 80 °C. The frozen leaves were then lyophilized until constant weight. The dried leaves were then powdered and used for extraction. To produce the crude extract, 4 g of powdered samples were soaked in 100 mL of absolute ethanol and subjected to sonication for 1 h in an ultrasonic bath sonicator (Branson, 1418510E-MTH model, Danbury, USA). The mixtures were then centrifuged at 13, 000 rpm for 30 min to separate the supernatant and the precipitates. The supernatant was collected and concentrated. Six individual extractions were performed to obtain 6 replications of samples.
Fractionation of crude extract
The crude extracts were subjected to solid phase extraction (SPE) in order to obtain fractions with different polarities. A weight of 1 g crude extract was premixed with 2 g of silica powder and loaded onto the SPE column cartridge (Isolute®, Mid Glamorgan, UK). The crude extract was then eluted consecutively with hexane, chloroform, ethyl acetate and methanol to yield hexane (HF), chloroform (CF), ethyl acetate (EF) and methanol (MF) fractions. Each solvent was allowed to pass through the cartridge until the eluent obtained was colorless. Each fraction collected was concentrated using rotary evaporator and finally lyophilized and kept in chiller at 4 °C for further analysis. Six individual SPE were performed in order to obtain the 6 replicates.
Total phenolic content (TPC)
The TPC was determined using a previously reported procedure . Samples were prepared by dissolving the fractions in DMSO. A total of 20 μL of samples and 100 μL of Folin-Ciocalteu reagent were mixed in the 96-well plates and incubated for 5 min. Later, 80 μL of 7.5% sodium carbonate solution was added. The plate was then incubated in the dark for 30 min and checked for the absorbance at 765 nm using a microplate reader (SPECTRAmax PLUS, Sunnyvale, CA, USA). The analysis was performed in 3 determinations for every sample. A standard curve of gallic acid was obtained to calculate the TPC and the results were expressed in μg GAE/mg extract.
DPPH free radical scavenging assay
The DPPH free radical scavenging activity of the fractions was evaluated according to the method previously described , with some modifications. The samples were prepared in various concentrations by serial dilution of stock solutions (1.0 mg/mL) prepared in DMSO. The test samples (50 μL) was mixed with DPPH (100 μL) in the wells. A set of wells was designated as reagent blank where the test sample was replaced by equal volume of sample solvent. The absorbance was then checked at 515 nm using microplate reader after incubation of the plate in dark for 30 min. The percentage of inhibition was calculated as % inhibition = [(AB-AS)/AB] × 100, where AB and As are the absorbance of reagent blank and tested samples, respectively. The analysis was performed in 3 determinations for every sample. The results are expressed as IC50 value in μg/mL. Quercetin was used as positive control in this assay.
α-Glucosidase inhibition assay
The inhibitory activity on α-glucosidase was assayed according to the previously described procedure , with some modifications. Samples from each fraction were prepared using DMSO at 1 mg/mL concentration as stock and further diluted to 8 serial dilutions using DMSO containing buffer. A total of 10 μL test sample, 130 μL of 30 mM phosphate buffer and 10 μL of α-glucosidase enzyme (prepared in 50 mM phosphate) were pre-incubated in the 96-well plate for 5 min at room temperature. A set of wells was designated as negative control where the test sample was replaced by equal volume of sample solvent. Next, 50 μL of PNPG substrate were added and the plate was further incubated for 15 min. The reaction was discontinued by addition of 50 μL of 2 M glycine (pH 10) and the absorbance was then read at 405 nm using a microplate reader. The percentage of inhibition was calculated as % inhibition = [(An-As)/An] × 100%, where An and As are the absorbance values of the negative control and test samples, respectively. The analysis was performed in 3 determinations for every sample. Results were expressed as IC50 values in μg/mL. Quercetin was used as positive control in this assay.
The NMR measurements were carried out based the method reported by Mediani et al. , with some modifications. Five milligrams of samples were mixed with 500 μL DMSO-d6 (containing 0.1% TSP). The mixture was then sonicated for 15 min and centrifuged at 13,000 rpm for 10 min. The supernatant was transferred to NMR tubes. The 1H NMR spectra were acquired using a 500 MHz Varian INOVA NMR spectrometer running at a frequency of 499.887 MHz and a temperature of 26 °C. The time to acquire each spectrum with presat setting was 3.54 min, comprising of 64 scans with a width of 20 ppm. Phase and baseline of all the 1H-NMR spectra were corrected using Chenomx software v. 8.1 (Alberta, Canada).
Ultra-high performance liquid chromatography tandem mass spectrometry analysis (UHPLC-MS/MS) of bioactive fractions
The UHPLC-MS/MS analysis of the bioactive fractions was carried out according to Zilani et al. , with modifications. The UHPLC was performed using a reversed-phase UHPLC Thermo system fitted with a Hypersil Gold, Thermo C18 column (2.1 mm × 100 mm, 1.9 μm). The mobile phase was 0.1% formic acid in water (solvent A) and 0.1% formic acid in LC–MS grade acetonitrile (solvent B). The analysis time was 33 min, and the flow rate was 250 μL/min. The injection volume was set to 5 μL. The programmed gradient proceeded using the following sequence for solvent B: 5% at 0 min, 15% at 6 min, 20% at 15 min, 20% at 18 min, 25% at 28 min, 100% at 30 min, 100% at 33 min. The samples were prepared by dissolving 4 mg of fraction in 2 mL of HPLC grade methanol, which was then centrifuged and filtered through a 0.22 µm nylon membrane into a 2-mL screw-capped sample vial. Mass spectra identification was attained using the Thermo Finnigan model (San Jose, CA, USA) Thermo Scientific™ Q Exactive™ Hybrid Quadrupole-Orbitrap mass spectrometer with an ESI source coupled to a Surveyor UHPLC binary pump and auto-sampler. This LC–MS/MS system combines quadrupole precursor ion selection with high resolution accurate mass (HRAM) Orbitrap detection and it is useful for targeted and untargeted screening and the identification of unknown compounds. The collision-induced dissociation (CID) energy was at 35%. This system was monitored by Xcalibur 2.2 and Mass Frontier 7.0 software. The negative and positive ion mass spectra were obtained in full ion scan mode (200–1700 amu) at a scan rate of 0.5 Hz. The most abundant ions in each scan were selected and subjected to MS/MS analysis.
Minitab software version 16 (Minitab Inc., State College, PA, USA) and GraphPad InStat version 2.02 (San Diego, USA) were used for the analysis of TPC, DPPH and α-glucosidase inhibition assay result data. Results were expressed as mean ± SD of 6 replicates. For multivariate data analysis (MVDA), the 1H-NMR spectra were bucketed and converted to ASCII files using Chenomx software. A total of 246 integrated regions were obtained after binning the region δ 0.5–10.0 with a width of δ 0.04. The residual signal of DMSO was excluded at δ 2.42–2.62. The data file was then imported into SIMCA-P software version 13.0 (Umeå, Sweden) for MVDA. Principal component analysis (PCA) and partial least-square analysis (PLS) were performed with Pareto scaling method. Pearson’s correlation was performed using MetaboAnalyst 2.5, a free metabolomics data analytical tool available online (http://www.metaboanalyst.ca).
Visual inspection of 1H NMR spectra of Neptunia oleracea fractions
The N. oleracea fractions (HF, CF, EF and MF) were subjected to 1H NMR analysis and the 1H NMR spectra were examined to identify the metabolites present in each fraction. The representative 1H NMR spectra of the fractions are shown in Fig. 1. Several classes of metabolites were identified from the fractions, including amino acids, fatty acids, phytosterols, sugars, flavonoids, triterpenes and phenolic acids. The chemical shifts of the identified metabolites are presented in Table 1. Primary metabolites were identified using the Chenomx database, whereas secondary metabolites were identified by comparison with reported data [18,19,20,21,22,23,24].
In the aliphatic region (δ 0.5–3.0) of the HF and the CF, signals of fatty acids and some amino acids, including leucine, alanine and valine were detected. The methyl signals (δ 0.8–0.9) assigned to a triterpene, namely oleanolic acid was also found in these fractions. Besides, signals of phytosterol were also detected in the CF. In contrast, the sugar signals (δ 3.0–5.0) were more prominent in the EF and the MF. Sugars such as α-glucose, β-glucose, and fructose were detected in both fractions. In the aromatic region (δ 5.5–8.5), signals belonging to flavonoids were also detected in these polar fractions. These signals were assigned for catechin and derivatives of quercetin, kaempferol, myricetin, and apigenin. Furthermore, the signals of three phenolic acids (caffeic, gallic and 3,4-O-dimethylgallic acids) were also observed in these fractions. These flavonoids and phenolic acids are the metabolites that were identified in previous works to be potential and important antioxidants and α-glucosidase inhibitors in N. oleracea [11, 12]. In addition, the signals for the anomeric protons of the arabinosyl, rhamnosyl and glucosyl moieties of the flavonoid derivatives were more obvious in the EF and MF as compared to the crude extracts and were found at δ 5.27 (d, J = 6.8 Hz), δ 5.15 (d, J = 1.5 Hz) and δ 5.33 (d, J = 10.2 Hz), respectively.
Visual inspection of the 1H NMR spectra revealed both qualitative and quantitative variations among the different fractions. The metabolites, including phytosterol, triterpene, fatty acids and amino acids resided mainly in the HF and CF, while the sugars and phenolics were more likely to reside in the EF and MF. However, for unbiased interpretation of the results, the 1H NMR spectral data of these fractions were further subjected to MVDA. The unsupervised MVDA, PCA, was used to differentiate these fractions based on their metabolite contents and to highlight the metabolites that contributed to the differentiation.
Metabolic classification of Neptunia oleracea fractions
The PCA score and loading plots of the 1H NMR data from the N. oleracea fractions are shown in Fig. 2. The score plot shows a clear discrimination of the more polar fractions (MF and EF) from the less polar fractions (CF and HF) by PC1, which contributed to most of the variance (52.4%). The more polar fractions are clustered on the positive side of PC1, while the less polar fractions are located on the other side. In addition, the MF and EF samples are projected close to one another in the score plot. This finding indicates that the EF and MF are closely related in their phytochemical contents and might contain the same metabolites in different intensities. From the loading plot, it is obvious that the less polar fractions (HF and CF) were separated from the other fractions by their higher contents of phytosterol, oleanolic acid, amino acids and fatty acids. Meanwhile, the EF and MF were discriminated from the HF and CF by their higher content of sugars and identified phenolics. These findings support the observation obtained via the visual inspection of the 1H NMR spectra. The presence of the identified phenolics in the EF and MF might contribute to the potent DPPH free radical scavenging and α-glucosidase inhibitory activities of the EF and MF. This can be observed in the results of the TPC, DPPH free radical scavenging and α-glucosidase inhibitory activities of the fractions, which showed significant TPC and inhibitions of the DPPH free radical and the α-glucosidase enzyme by the EF and MF (Fig. 3). In order to confirm the contribution of the phenolics to the bioactivities of EF and MF, the correlation was assessed between the identified phenolics and the studied bioactivities.
Correlation between phenolics and bioactivities of Neptunia oleracea fractions
The PLS model was constructed for the correlation evaluation of phenolics with the bioactivities of EF and MF, where the NMR chemical signals were input as X variables and the 1/IC50 values of the DPPH radical scavenging and α-glucosidase inhibitory activities were the Y variables. The resulted PLS model was validated by goodness of fit, prediction of Y and a permutation test. Autofit of the PLS model revealed its good fitness (R2Y = 0.925) and predictive ability (Q2 value = 0.888). The permutation test presented in Fig. 4 indicated that the PLS model does not show overfitting, with the Y-intercepts of R2 and Q2 less than 0.3 and 0.05, respectively, and the R2 lines far from being horizontal . These findings showed that the PLS model is valid and can appropriately correlate the identified phenolics with the studied bioactivities.
The PLS biplot is a combination of score and loading plots. Distance of the X and Y variables to the sample clusters indicates the extent of their contribution to the characteristics of the respective cluster. As shown in Fig. 5, the more polar fractions were separated from the less polar ones by the PC1, as previously observed in the PCA. Furthermore, the two Y variables (DPPH radical scavenging and α-glucosidase inhibitory activities) were projected on the same side with the EF and MF, indicating their better activities compared to those of the less polar fractions. Comparing the more polar fractions, the MF was more active than the EF, as the MF was closer to the Y variables. These results were in agreement with the bioactivities findings presented earlier (Fig. 3). In addition, all the identified phenolics were on the side where the active fractions located. This shows their positive contributions to the DPPH free radical scavenging and α-glucosidase inhibitory activities of the EF and MF. Among these phenolics, flavonoids including derivatives of quercetin and kaempferol play important roles in the free radical scavenging and α-glucosidase inhibitory activities of the EF; while derivatives of apigenin and phenolic acids were observed to be the important phytochemical markers for the MF.
Pearson’s correlation was also performed to support the correlation of the metabolites with the DPPH and α-glucosidase inhibitory activities and the results are presented in Fig. 6. Pearson’s correlation showed positive correlations between the studied bioactivities and identified phenolics. This result demonstrates the relationship of these phenolics with the DPPH and α-glucosidase inhibitory activities of N. oleracea. Moreover, these metabolites had also high correlation with each other as well as with the sugars.
UHPLC–MS/MS profiling of Neptunia oleracea bioactive fractions
The EF and MF exhibited the most significant DPPH free radical scavenging and α-glucosidase inhibitory activities as the phenolic compounds were mainly distributed in these fractions. Hence, these fractions were subjected to UHPLC–MS/MS analysis to identify the derivatives of the phenolics present. The respective total ion chromatograms (TIC) of the EF and MF are presented in Fig. 7, while Table 2 summarizes the retention time (RT), MS/MS data and distribution of the identified metabolites. The TIC profiles showed that most of the prominent peaks were attributed to the presence of flavonoids. A total of 37 metabolites were tentatively identified based on the MS/MS data in comparison with the literature. These identified metabolites were classified into 7 groups, namely quercetin and its derivatives, apigenin derivatives, myricetin derivatives, kaempferol and its derivatives, other flavonoids, phenolic acid derivatives and other metabolites. Most of the identified metabolites were present in both the EF and MF. However, it was observed that the derivatives of phenolic acids as well as some flavonoid di- and triglycosides were only present in the MF, while the flavonoid aglycones were detected only in the EF.
In the present study, ethanolic crude extract of N. oleracea was subjected to SPE to yield 4 fractions with different polarities (HF, CF, EF and MF). Metabolites from different classes were identified in the fractions. Both scrutiny of the 1H NMR spectra and unsupervised MVDA (PCA) showed that the HF and CF contained higher contents of phytosterol, oleanolic acid, amino acids and fatty acids; while EF and MF were characterized by their higher content of sugars and identified phenolics. The principle of “like-dissolves-like” can explain the distribution of the identified metabolites in the different fractions. The SPE of N. oleracea crude extract involved separating the metabolites relative to their solubility in the eluting solvent. The long hydrocarbon chain of fatty acids and the cyclic hydrocarbon structures of triterpene and phytosterol made these metabolites have low polarities . Thus, these metabolites were more likely dissolved and eluted by the solvents with low polarities, such as hexane and chloroform. In contrast, the sugar attachment and the presence of a polar functional group, such as phenolic hydroxyl and carbonyl groups, make the phenolics polar and hence more soluble in solvents with relatively higher polarities, such as ethyl acetate and methanol. The better recovery of phenolics by these solvents compared to hexane and chloroform has also been reported previously [27, 28].
The EF and MF showed higher TPC and exhibited more potent DPPH free radical scavenging and α-glucosidase inhibitory activities compared to HF and CF (Fig. 3). The correlation findings revealed the contributions of phenolics towards the studied bioactivities. Via PLS analysis, derivatives of quercetin and kaempferol were observed to play important roles in the free radical scavenging and α-glucosidase inhibitory activities of the EF; while derivatives of apigenin and phenolic acids were the important bioactivity contributors of the MF. These results were in agreement with those in previous studies that revealed the significance of apigenin, quercetin and kaempferol derivatives as well as caffeic acid as the antioxidants and α-glucosidase inhibitors in N. oleracea [11, 12]. The potent antioxidant and α-glucosidase inhibitory activities of these metabolites can be explained by the presence and number of phenolic hydroxyl groups [29, 30]. Furthermore, ortho-dihydroxyl (catechol) substitution had been reported to impart high antioxidant and α-glucosidase inhibitory activities to molecules . Hence, the presence of catechol moiety may contribute to the potent activities of these metabolites, particularly the quercetin derivatives and caffeic acid. In addition, Pearson’s correlation suggested synergistic effects of the phenolics towards the studied bioactivities due to their strong correlations with each other. Besides, the result revealed that the high intensity of phenolics was associated with the sugar present. This positive relationship between carbon-rich secondary metabolites such as phenolics and the sugar has also been reported previously [32, 33].
In order to identify the derivatives of phenolics present in the bioactive fractions, the EF and MF were subjected to UHPLC-MS/MS analysis. The UHPLC coupled with tandem MS represents a robust tool to characterize the metabolites in complex mixtures, owing to the combination of separation capability of UHPLC and identification ability of tandem MS. Together with some non-phenolic metabolites, a total of 37 metabolites were tentatively identified in the EF and MF. Most of these metabolites were present in both the EF and MF, with the exception of some metabolites which present in either one of the fraction, as shown in Fig. 7 and Table 2. This difference might be due to the different polarities of the metabolites, which were affected by the presence and number of the hydroxyl groups as well as sugar substitutions . Identification of the metabolites by the UHPLC-MS/MS are discussed in detail.
Identification of quercetin and its derivatives
In the EF and MF, a total of 10 metabolites were identified as quercetin derivatives based on the presence of an aglycone fragment ion at m/z 301 and the characteristic fragment ions at m/z 271 and 151 in their MS/MS spectra . Metabolites 10 and 22 showed the same deprotonated molecular ion at m/z 609 but with a different pattern of fragmentation. Metabolite 10 was identified as quercetin 3-O-rhamnoside-7-O-glucoside based on the fragment ions at m/z 463, 447 and 301, which arose from the subsequent loss of rhamnose (m/z 146) and glucose (m/z 162) moieties. This fragmentation pattern was agreed with a previously reported pattern . Meanwhile the fragment ions at m/z 463 and 447 were absent in metabolite 22, showing less prominent fragmentation of this metabolite. This finding suggests that the rhamnose and glucose are probably attached together as a rutinose. This observation was similar to those reported by Sanchez-Rabaneda et al. . Hence, metabolite 22 was assigned as quercetin-3-O-rutinoside (rutin).
The transition of the deprotonated molecular ion of metabolite 14 (m/z 755) to the quercetin ion (m/z 301) showed the loss of 2 rhamnosyl and 1 glucosyl moieties. This metabolite was identified as quercetin 3-O-(2,6-di-O-rhamnosylglucoside) by comparing it with reported MS/MS data . Metabolite 15 gave a deprotonated molecular ion at m/z 625 and fragment ions at m/z 463 and 301, showing the loss of two glucose moieties. Hence, this metabolite was identified as quercetin-3,7-di-O-glucoside . Metabolites 26, 27, 32 resulted in a quercetin fragment ion from the loss of glucuronyl (m/z 176), glucosyl (m/z 162) and rhamnosyl (m/z 146) moieties, respectively. Consequently, these metabolites were identified as quercetin-3-O-glucoronide, quercetin-3-O-glucoside and quercetin-3-O-rhamnoside, respectively [38, 40]. Metabolites 28 and 30 had both the same deprotonated molecular ion at m/z 433 and the same fragmentation pattern. However, based on the elution order and the previously reported data , metabolite 28 was assigned as quercetin-3-O-xyloside; while the later eluting metabolite 30 was assigned as quercetin-3-O-arabinoside. A methylated type of quercetin derivative was also detected. With the presence of a deprotonated molecular ion at m/z 623 and a fragment ion at m/z 315 due to the loss of rutinoside (m/z 308), metabolite 29 was identified as isorhamnetin-3-O-rutinoside . Lastly, metabolite 36 was identified as quercetin based on the deprotonated molecular ion at 301 in its MS/MS spectrum .
Identification of apigenin derivatives
Thalang et al.  had reported the presence of apigenin derivatives in N. oleracea extracts. However, the identity of these derivatives had not yet been revealed. In this present study, four derivatives of apigenin were tentatively identified in the EF and MF of N. oleracea. These metabolites showed the loss of m/z 120 and 90 in their MS/MS spectra, which are the typical fragmentation of C-glycoside that corresponding to the cross-ring cleavages in the sugar unit . With the deprotonated molecular ion at m/z 593 and C-glycoside fragmentation pattern, metabolite 12 was identified as apigenin-6,8-di-C-glucoside (vicenin 2). This fragmentation behavior was consistent with the previous characterization . Both metabolite 23 and 24 showed the same deprotonated molecular ion at m/z 431 and common characteristic fragment ions at m/z 341 (loss of m/z 90) and 311 (loss of m/z 120). However, after comparison of the elution order and with the previously reported data , metabolite 23 was assigned as apigenin-6-C-glucoside (vitexin) while metabolite 24 was assigned as apigenin-8-C-glucoside (isovitexin). The presence of the fragment ion at m/z 353 also differentiated isovitexin from vitexin, as reported by Sanchez-Rabaneda et al. . Metabolite 21 showed a deprotonated molecular ion at m/z 577. The transition of m/z 577 to 413 showed the loss of a rhamnose unit, while the presence of the ion at m/z 457 (due to loss of m/z 120) revealed the cross-ring cleavage of the hexose moiety and the linkage between this hexose and the rhamnose is at the 1–2 position . With this spectral evidence, metabolite 21 was elucidated as vitexin-2″-O-rhamnoside.
Identification of myricetin derivatives
Metabolites 16, 18, 19 and 20 were identified as myricetin derivatives based on the presence of a fragment ion at m/z 316, corresponding to the myricetin aglycone fragment in the MS/MS spectra . Metabolite 16 was identified as myricetin-3-O-rutinoside based on the transition of m/z 625 to 316, which revealed the loss of a rutinose moiety . Metabolites 18, 19 and 20 were assigned as myricetin-3-O-glucoside, myricetin-3-O-rhamnoside and myricetin-3-O-arabinoside based on the deprotonated molecular ions at m/z 479, 463 and 449, respectively. The transition of these ions to m/z 316 revealed the losses of the respective sugar moieties [38, 45].
Identification of kaempferol and its derivatives
Four metabolites (25, 31, 34 and 35) were detected and assigned as kaempferol derivatives in the EF and MF. All these metabolites showed characteristic fragment ions at m/z 285 due to the loss of sugar moieties. The transition of m/z 593 to m/z 285 in metabolite 25 indicated the loss of the O-rhamanosyglucosyl moiety and the presence of the fragment ion at m/z 327 indicated the presence of a 1 → 2 isomer . Therefore, metabolite 25 was identified as kaempferol 7-O-(2″-rhamnosyl)-glucoside. Meanwhile, the MS/MS data of metabolites 31, 34 and 35 showed the losses of m/z 162, 132 and 146, which were characteristic of the presence of glucosyl, arabinosyl and rhamnosyl moieties, respectively. Hence, metabolites 31, 34 and 35 were identified as kaemperol-3-O-glucoside, kaempferol-3-O-arabinoside and kaempferol-3-O-rhamnoside, respectively. These assignments were consistent with previously reported data . Kaempferol aglycone (metabolite 37) was detected in the EF with the presence of fragment ions at m/z 217 and 151 .
Identification of other flavonoids
Besides the derivatives of quercetin, kaempferol, myricetin and apigenin, derivatives of other flavonoids were also detected in the EF and MF of N. oleracea. Metabolites 8 and 17 were detected and elucidated as C-glycosides of luteolin based on the losses of m/z 90 and 120. Comparison of the spectral data with those reported had identified metabolite 8 and 17 as luteolin-6,8-di-C-glucoside (leucenin 2) and luteolin-8-C-glucoside (orientin), respectively . Catechin and its derivative were also detected. Metabolite 9 showed a deprotonated molecular ion at m/z 289 and yielded fragment ions at m/z 245, 179 and 137, which matched with those reported for catechin . Hence, metabolite 9 was assigned as catechin. Meanwhile, metabolite 13 yielded fragment ions m/z 331 and 169, which was corresponding to the epigallocatechin and gallic acid moieties, respectively. Subsequently, metabolite 13 was identified as (−)-epigallocatechin-3-gallate . Metabolite 33 showed a deprotonated molecular ion at m/z 461. The MS/MS data showed the loss of a hexose moiety and the subsequent fragmentation pattern of gliricidin, based on comparison with previously described characterization . Hence, metabolite 33 was identified as gliricidin-O-hexoside.
Identification of phenolic acid derivatives
Derivatives of phenolic acids were also detected and identified as well. Analysis of the deprotonated molecular ions at m/z 331 and 483 (metabolites 3 and 7, respectively) showed the presence of galloylglucose. In the fragmentation of ion m/z 331, the loss of a glucose moiety yielded the fragment ion m/z 169, which originated from a gallic acid. Hence, metabolite 3 was identified as monogalloylglucose . Metabolite 7 was assigned as digalloylglucose based on comparison of the fragmentation pattern with those reported . Metabolite 2 and 5 yielded fragment ions at m/z 147 (cinnamic acid moiety) and m/z 197 (syringic acid moiety), and they were identified as unknown derivatives of cinnamic and syringic acids, respectively. Metabolite 6 was identified as protocatechuic acid-O-pentoside, by the ions m/z 153, 135 and 109, which originate from protocatechuic acid and its subsequent fragmentation . Metabolite 11, with a deprotonated molecular ion at m/z 367 and fragmentation data that corresponding to methylgallate, was identified as a dimer of methylgallate .
Identification of other metabolites
A derivative of phenylethanoid was also detected in the EF. Metabolite 4 was identified as hydroxytyrosol hexoside as the transition of the deprotonated molecular ion at m/z 315 to the fragment ion of hydroxytyrosol at m/z 153 showed the loss of a hexose moiety . In addition to all the aforementioned phenolic constituents, an oligosaccharide was identified. Metabolite 1, detected in the MF was identified as a hexose polymer, as the fragment ions at m/z 341 and 179 revealed the losses of hexose moieties. Besides, the fragment ions at m/z 143, 131, 119, 113, 89, and 71 are the typical fragments that arise from hexose , thus further confirming the presence of hexose.
Neptunia oleracea had been fractionated into 4 fractions with different polarities (HF, CF, EF and MF). The correlation study identified the EF and MF as the bioactive fractions of N. oleracea with regards to the DPPH scavenging and α-glucosidase inhibitory activities. Both activities showed strong correlation with the identified phenolics. Via the UHPLC-MS/MS profiling of the EF and MF, a total of 37 metabolites, with mostly phenolics were tentatively identified. Most of these identified metabolites were present in both the EF and MF. To the best of our knowledge, the differentiation, metabolite-bioactivity correlation and UHPLC-MS/MS profiling of bioactive fractions were performed for the first time on N. oleracea. The findings in this work are comprehensively supporting our previous work regarding the contribution of phenolics to the DPPH free radical scavenging and α-glucosidase inhibitory activities of this plant, with more detailed identification of these bioactive metabolites. In summary, N. oleracea is a rich source of phenolics that can be potential antioxidants and α-glucosidase inhibitors used in the management of diabetes.
Ethyl acetate fraction
Reactive oxygen species
Solid phase extraction
Total ion chromatograms
Total phenolic content
Ultra-high performance liquid chromatography tandem mass spectrometry analysis
World Health Organization. Diabetes (Fact sheets) http://www.who.int/en/news-room/fact-sheets/detail/diabetes. 2017; Accessed 11 Dec 2017.
Mohamed EA, Siddiqui MJ, Ang LF, Sadikun A, Chan SH, Tan SC, Asmawi MZ, Yam MF. Potent α-glucosidase and α-amylase inhibitory activities of standardized 50% ethanolic extracts and sinensetin from Orthosiphon stamineus Benth as anti-diabetic mechanism. BMC Complement Altern Med. 2012;12:176.
Javadi N, Abas F, Mediani A, Hamid AA, Khatib A, Simoh S, Shaari K. Effect of storage time on metabolite profile and alpha-glucosidase inhibitory activity of Cosmos caudatus leaves - GCMS based metabolomics approach. J Food Drug Anal. 2015;23:433–41.
Saupi N, Zakaria MH, Bujang JS, Arshad A. The proximate compositions and mineral contents of Neptuniaoleracea Loureiro, an aquatic plant from Malaysia. Emirates J Food Agric. 2015;27:266–74.
Tee E, Lim CL. Carotenoid composition and content of Malaysian vegetables and fruits by the AOAC and HPLC methods. Food Chem. 1991;41:309–39.
Bhoomannavar VS, Patil VP, Hugar S, Nanjappaiah HM, Kalyane N. Anti-ulcer activity of Neptunia oleracea lour. Pharmacologyonline. 2011;3:1015–20.
Bhoomannavar VS, Shivakumar SI, Hallikeri CS, Hatapakki BC. Hepatoprotective activity of leaves of Neptunia oleracea lour in carbon tetrachloride induced rats. Res J Pharm Biol Chem Sci. 2011;2:309–14.
Nakamura Y, Murakami A, Koshimizu K, Ohigashi H. Identification of pheophorbide a and its related compounds as possible anti-tumor promoters in the leaves of Neptunia oleracea. Biosci Biotechnol Biochem. 1996;60:1028–30.
Kumar N, Chaiyasut C. Health promotion potential of vegetables cultivated in northern Thailand: a preliminary screening of tannin and flavonoid contents, 5α-reductase inhibition, astringent activity, and antioxidant activities. J Evid Based Complementary Altern Med. 2017;22:573–9.
Paul S, Choudhury SN, De B. Structural elucidation of a bioactive compound from the leaves of Neptunia prostrate. Asian J Chem. 2012;24:1469–72.
Lee SY, Abas F, Khatib A, Ismail IS, 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.
Lee SY, Mediani A, Maulidiani M, Khatib A, Ismail IS, 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–52.
Zhang Q, Zhang J, Shen J, Silva A, Dennis DA, Barrow CJ. A simple 96-well microplate method for estimation of total polyphenol content in seaweeds. J Appl Phycol. 2006;18:445–50.
Wan C, Yuan T, Cirello AL, Seeram NP. Antioxidant and α-glucosidase inhibitory phenolics isolated from highbush blueberry flowers. Food Chem. 2012;135:1929–37.
Deutschländer MS, Van de Venter M, Roux S, Louw J, Lall N. Hypoglycaemic activity of four plant extracts traditionally used in South Africa for diabetes. J Ethnopharmacol. 2009;30:619–24.
Mediani A, Abas F, Khatib A, Maulidiani H, Shaari K, Choi YH, Lajis NH. 1H-NMR-based metabolomics approach to understanding the drying effects on the phytochemicals in Cosmos caudatus. Food Res Int. 2012;49:763–70.
Zilani MN, Sultana T, Rahman SA, Anisuzzman M, Islam MA, Shilpi JA, Hossain MG. Chemical composition and pharmacological activities of Pisum sativum. BMC Complement Altern Med 2017;1:171.
Benavides A, Montoro P, Bassarello C, Piacente S, Pizza C. Catechin derivatives in Jatropha macrantha stems: characterisation and LC/ESI/MS/MS quali-quantitative analysis. J Pharm Biomed Anal. 2006;40:639–47.
Fossen T, Larsen Å, Kiremire BT, Andersen ØM. Flavonoids from blue flowers of Nymphaea caerulea. Phytochemistry. 1999;51:1133–7.
Gogna N, Hamid N, Dorai K. Metabolomic profiling of the phytomedicinal constituents of Carica papaya L. leaves and seeds by 1H NMR spectroscopy and multivariate statistical analysis. J Pharm Biomed Anal. 2015;115:74–85.
Jin JL, Lee YY, Heo JE, Lee S, Kim JM, Yun-Choi HS. Anti-platelet pentacyclic triterpenoids from leaves of Campsis grandiflora. Arch Pharm Res. 2004;27:376–80.
Park BJ, Matsuta T, Kanazawa T, Park CH, Chang KJ, Onjo M. Phenolic compounds from the leaves of Psidium guajava II. Quercetin and its glycosides. Chem Nat Compd. 2012;48:477–9.
Wang DM, Pu WJ, Wang YH, Zhang YJ, Wang SS. A new isorhamnetin glycoside and other phenolic compounds from Callianthemum taipaicum. Molecules. 2012;17:4595–603.
Yaya S, Amian K, Benjamin B, Fanté B, Sorho S, Amadou TS, Jean-marie C. Flavonoids and gallic acid from leaves of Santaloides afzelii (connaraceae). Rasayan J Chem. 2012;5:332–7.
Maulidiani AF, Khatib A, Shitan M, Shaari K, Lajis NH. Comparison of partial least squares and artificial neural network for the prediction of antioxidant activity in extract of Pegaga (Centella) varieties from 1H nuclear magnetic resonance spectroscopy. Food Res Int. 2013;54:852–60.
Stoker HS. General, organic, and biological chemistry. 5th ed Nelson Education. 2012.
Ado MA, Abas F, Ismail IS, Ghazali HM, Shaari K. Chemical profile and antiacetylcholinesterase, antityrosinase, antioxidant and α-glucosidase inhibitory activity of Cynometra cauliflora L. leaves. J Sci Food Agric. 2015;95:635–42.
Barchan A, Bakkali M, Arakrak A, Pagán R, Laglaoui A. The effects of solvents polarity on the phenolic contents and antioxidant activity of three Mentha species extracts. Int J Curr Microbiol Appl Sci. 2014;3:399–412.
Kongpichitchoke T, Hsu JL, Huang TC. Number of hydroxyl groups on the B-ring of flavonoids affects their antioxidant activity and interaction with phorbol ester binding site of PKCδ C1B domain: in vitro and in silico studies. J Agric Food Chem. 2015;63:4580–6.
Ma K, Han J, Bao L, Wei T, Liu H. Two sarcoviolins with antioxidative and α-glucosidase inhibitory activity from the edible mushroom Sarcodon leucopus collected in Tibet. J Nat Prod. 2014;77:942–7.
Wang Y, Xiang L, Wang C, Tang C, He X. Antidiabetic and antioxidant effects and phytochemicals of mulberry fruit (Morus alba L.) polyphenol enhanced extract. PLoS One. 2013;8:e71144.
Kim JK, Park SY, Lim SH, Yeo Y, Cho HS, Ha SH. Comparative metabolic profiling of pigmented rice (Oryza sativa L.) cultivars reveals primary metabolites are correlated with secondary metabolites. J Cereal Sci. 2013;57:14–20.
Matt P, Krapp A, Haake V, Mock HP, Stitt M. Decreased rubisco activity leads to dramatic changes of nitrate metabolism, amino acid metabolism and the levels of phenylpropanoids and nicotine in tobacco antisense RBCS transformants. Plant J. 2002;30:663–77.
Welch CR, Wu Q, Simon JE. Recent advances in anthocyanin analysis and characterization. Curr Anal Chem. 2008;4:75–101.
Kumar S, Chandra P, Bajpai V, Singh A, Srivastava M, Mishra DK, Kumar B. Rapid qualitative and quantitative analysis of bioactive compounds from Phyllanthus amarus using LC/MS/MS techniques. Ind Crop Prod. 2015;69:143–52.
Materska M. Bioactive phenolics of fresh and freeze-dried sweet and semi-spicy pepper fruits (Capsicum annuum L.). J Funct Foods. 2014;7:269–77.
Sanchez-Rabaneda F, Jauregui O, Casals I, Andres-Lacueva C, Izquierdo-Pulido M, Lamuela-Raventos RM. Liquid chromatographic/electrospray ionization tandem mass spectrometric study of the phenolic composition of cocoa (Theobroma cacao). J Mass Spectrom. 2003;38:35–42.
Elsadig Karar MG, Quiet L, Rezk A, Jaiswal R, Rehders M, Ullrich MS, Brix K, Kuhnert N. Phenolic profile and in vitro assessment of cytotoxicity and antibacterial activity of Ziziphus spina-christi leaf extracts. Med Chem (Los Angeles). 2016;6:143–56.
Francescato LN, Debenedetti SL, Schwanz TG, Bassani VL, Henriques AT. Identification of phenolic compounds in Equisetum giganteum by LC-ESI-MS/MS and a new approach to total flavonoid quantification. Talanta. 2013;105:192–203.
Ibrahim RM, El-Halawany AM, Saleh DO, El NEMB, El-Shabrawy AE-RO, El-Hawary SS. HPLC-DAD-MS/MS profiling of phenolics from Securigera securidaca flowers and its anti-hyperglycemic and anti-hyperlipidemic activities. Rev Bras Farmacogn. 2015;25:134–41.
Li Z, Meng F, Zhang Y, Sun L, Yu L, Zhang Z, Peng S, Guo J. Simultaneous quantification of hyperin, reynoutrin and guaijaverin in mice plasma by LC-MS/MS: application to a pharmacokinetic study. Biomed Chromatogr. 2016;30:1124–30.
Brito A, Ramirez JE, Areche C, Sepúlveda B, Simirgiotis MJ. HPLC-UV-MS profiles of phenolic compounds and antioxidant activity of fruits from three citrus species consumed in northern Chile. Molecules. 2014;19:17400–21.
Thalang NV, Trakoontivakorn G, Nakahara K. Determination of antioxidant activity of some commonly consumed leafy vegetables in Thailand. JIRCAS J Sci Pap. 2001;9:39–46.
Li YQ, Zhou FC, Gao F, Bian JS, Shan F. Comparative evaluation of quercetin, isoquercetin and rutin as inhibitors of α-glucosidase. J Agric Food Chem. 2009;57:11463–8.
Wu Y, Jiang X, Zhang S, Dai X, Liu Y, Tan H, Gao L, Xia T. Quantification of flavonol glycosides in Camellia sinensis by MRM mode of UPLC-QQQ-MS/MS. J Chromatogr B. 2016;1017:10–7.
Es-Safi NE, Gómez-Cordovés C. Characterization of flavonoid glycosides from fenugreek (Trigonella foenum-graecum) crude seeds by HPLC-DAD-ESI/MS analysis. Int J Mol Sci. 2014;15:20668–85.
Zhao Y, Li X, Zeng X, Huang S, Hou S, Lai X. Characterization of phenolic constituents in Lithocarpus polystachyus. Anal Methods. 2014;6:1359–63.
Simirgiotis JM, Schmeda-Hirschmann G, Bórquez J, Kennelly JE. The Passiflora tripartita (Banana passion) fruit: a source of bioactive flavonoid C-glycosides isolated by HSCCC and characterized by HPLC–DAD–ESI/MS/MS. Molecules. 2013;18:1672–92.
Sun J, Liang F, Bin Y, Li P, Duan C. Screening non-colored phenolics in red wines using liquid chromatography/ultraviolet and mass spectrometry/mass spectrometry libraries. Molecules. 2007;12:679–93.
Savic I, Nikolic V, Savic I, Nikolic L, Jovic M, Jovic M. The qualitative analysis of the green tea extract using ESI-MS method. Savrem Tehnol. 2014;3:30–7.
Ye M, Yang WZ, Di LK, Qiao X, Li BJ, Cheng J, Feng J, Guo DA, Zhao YY. Characterization of flavonoids in Millettia nitida var. hirsutissima by HPLC/DAD/ESI-MSn. J Pharm Anal. 2012;2:35–42.
Wyrepkowski CC, Lu D, Gomes M, Sinhorin AP, Vilegas W, De Grandis RA, Resende FA, Varanda EA, Campaner L. Characterization and quantification of the compounds of the ethanolic extract from Caesalpinia ferrea stem bark and evaluation of their mutagenic activity. Molecules. 2014;19:16039–57.
Hofmann T, Nebehaj E, Albert L. Antioxidant properties and detailed polyphenol profiling of European hornbeam (Carpinus betulus L.) leaves by multiple antioxidant capacity assays and high-performance liquid chromatography/multistage electrospray mass spectrometry. Ind Crop Prod. 2016;87:340–9.
Eyles A, Jones W, Riedl K, Cipollini D, Schwartz S, Chan K, Herms DA, Bonello P. Comparative phloem chemistry of Manchurian (Fraxinus mandshurica) and two north American ash species (Fraxinus americana and Fraxinus pennsylvanica). J Chem Ecol. 2007;33:1430–48.
Llorent-Martinez EJ, Spinola V, Gouveia S, Castilho PC. HPLC-ESI-MSn characterization of phenolic compounds, terpenoid saponins, and other minor compounds in Bituminaria bituminosa. Ind Crop Prod. 2015;69:80–90.
This study was financially supported by Universiti Putra Malaysia under Research University Grant Scheme (RUGS) (9362700). The first author extends her appreciation to the Ministry of Higher Education Malaysia for the scholarship provided.
Availability of data and materials
The datasets used and/or analyzed during the current study available from the corresponding author on reasonable request.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.