Vibrational Spectroscopy-Based Chemometrics Analysis of Clinacanthus nutans Extracts after Postharvest Processing and Extract Effects on Cardiac C-Kit Cells

Chemical constituents in plants can be greatly affected by postharvest processing, and it is important to identify the factors that lead to significant changes in chemistry and bioactivity. In this study, attenuated total reflectance-Fourier transform infrared (ATR-FTIR) spectroscopy was used to analyze extracts of Clinacanthus nutan (C. nutans) leaves generated using different parameters (solvent polarities, solid-liquid ratios, ultrasonic durations, and cycles of extraction). In addition, the effects of these extracts on the viability of cardiac c-kit cells (CCs) were tested. The IR spectra were processed using SIMCA-P software. PCA results of all tested parameter sets were within acceptable values. Solvent polarity was identified as the most influential factor to observe the differences in chemical profile and activities of C. nutans extracts. Ideal extraction conditions were identified, for two sample groups (G1 and G2), as they showed optimal total phenolic content (TPC) yield of 44.66 ± 0.83 mg GAE/g dw and 45.99 ± 0.29 mg GAE/g dw and CC viability of 171.81 ± 4.06% and 147.53 ± 6.80%, respectively. Validation tools such as CV-ANOVA (p < 0.05) and permutation (R2 and Q2 plots were well intercepted to each other) have further affirmed the significance and reliability of the partial least square (PLS) model of solvent polarity employed in extraction. Hence, these approaches help optimize postharvest processes that encourage positive TPC and CCs results in C. nutans extracts.

Myocardial infarction remains a major cause of heart failure and contributes to high morbidity and mortality worldwide. Despite the availability of advanced therapeutics, it is still challenging to reverse heart remodeling caused by cardiomyocyte loss [25]. Studies have reported that the adult heart possesses endogenous cardiac c-kit cells that can regenerate dead myocardium [26,27], and they can be activated from the quiescent state by growth factors [28] or by factors and matrices secreted by mesenchymal stem cells [29][30][31]. However, growth factors are expensive and limited by challenges such as low plasma stability, short biological half-lives, and low specificity to target organs [32]. An alternative remedy sourced from the plant has been studied in Geum japonicum extract and was reported to be helpful for myocardial regeneration [33]. us, plantderived compounds could be another therapeutic regimen for cardiac repair. C. nutans, a widely grown local herb, may also prove to be helpful for treating heart disease. e herb can potentially prevent insulin resistance-induced cardiovascular diseases by reducing metabolic effects and transcriptional changes induced by a high fat and high cholesterol diet [34].
Plant metabolomic analysis involves the application of mathematical and statistical calculations of metabolite variations in plants of different species under different processing conditions [35]. Sample variation analyses are commonly performed with analytical tools such as Fourier transform infrared (FTIR) spectroscopy, liquid chromatography-mass spectroscopy (LCMS), and nuclear magnetic resonance (NMR). Multivariate data analysis (MVDA), including principal component analysis (PCA) and partial least squares to latent structures (PLS), helps reduce the dimensionality of big datasets to facilitate visualization of sample distribution and to correlate the chemical profile with the tested activity.
Studies have shown that extraction parameters may be the main factors in determining the polyphenol content and the bioactivity of constituents of plant extracts [36]. erefore, the goal of this study was to optimize postharvest processing of C. nutans leaves (i.e., extraction parameters) by identifying the best extraction solvent, solid-liquid ratio, duration of extraction, and the number of extraction cycles. We employed attenuated total reflectance-FTIR (ATR-F-TIR) spectroscopy to profile the metabolites present in the samples and chemometrics techniques to investigate the quality and safety of C. nutans extracts after postharvest processing. We also tested the effects of these extracts on the viability of mouse endogenous cardiac c-kit cells.  SK 1980/11). e dried leaves were pulverized into a fine powder using a grinder (Ultra Centrifugal Mill ZM200, Retsch, Haan, Germany). e powder was then sealed in a plastic bag and kept at room temperature until further use.
e mixtures went through a sonication process according to their respective conditions. Extracts were centrifuged at 6000 rpm for 15 min. Supernatants collected were vacuum filtered through Whatman No.1 filter paper (Little Chalfont, Bucks, UK) and lyophilized under centrifugal vacuum using a laboratory evaporator (Genevac EZ-2 Series, Ipswich, England). e lyophilized extracts were kept in a desiccator before analysis.

ATR-FTIR
Analysis. FTIR spectra were obtained using a ermo Nicolet iS10 FTIR spectrometer (Madison, WI, USA), equipped with a smart iTR accessory. A small amount of each extract (n � 3 technical replicates) was applied neatly on the ATR crystal and scanned in absorption mode. IR measurements were made at the mid-IR range between 4000 and 650 cm −1 , with a spectral resolution of 4 cm −1 .

Spectral
Processing. All spectra were converted into percentage transmittance form for spectrum preprocessing. Each IR spectrum was then subjected to baseline correction and smoothing using Omnic software (Version 7.3, Madison, WI, USA). Data exported from Omnic software were compiled into an MS Excel (Version 2016, Redmond, WA, USA) file and subjected to SIMCA-P software (Version 13.0, Umetrics, Umeå, Sweden) for MVDA analysis [40].

Total Phenolic Content.
e total phenolic content (TPC) of each extract was assessed by the FC method, as described previously [41,42] with some modifications. Gallic acid was used as the standard for this assay. For each 10 µL test sample or "standard," 600 µL of FC reagent (10-fold diluted) and 400 µL of sodium carbonate solution were added, mixed well, and kept in the dark for 60 min. Next, 150 µL of each reaction mixture was transferred to a 96-well plate. e absorbance value was acquired at 765 nm using a FLUOstar ® Omega microplate reader (BMG Labtech, Ortenberg, Germany). A gallic acid calibration curve was constructed and used to calculate the TPC of each sample. Results were expressed as mg gallic acid equivalent (GAE)/g dry weight (dw) of extract.

Cytotoxicity Assay.
Cardiac c-kit cells (CCs) were isolated from C57/BL mice, as reported previously. Passage 5 CCs were seeded into a 96-well plate with 10% FBS at a density of 1000 cells per well and incubated at 37°C and 5% CO 2 . Next, 5 mg of each C. nutans extract was prepared with 2.5 mL of serum-free DMEM/F12 to yield a stock concentration of 2 mg/mL. e extracts went through a sonication process to ensure solubility and were filtered across a 0.22 µm syringe filter. Extracts then were serially diluted to 50, 100, 500, and 1000 µg/mL. After one hour of incubation, C. nutans crude leave extracts were added to CCs and incubated for 96 h. All experiments were conducted in triplicate. Serum-free DMEM/F12 was served as the control of the study. Cytotoxicity of C. nutans was examined using the presto blue assay. Fluorescence readings were taken using the FLUOstar ® Omega microplate reader with fluorescence excitation and emission wavelengths of 540-570 nm and 580-610 nm, respectively. Cell viability in percentage was calculated using the following formula:

Statistical
Analysis. Data obtained from TPC and CCs cytotoxicity assays were processed in MS Excel and analyzed using SPSS (Version 22.0, IBM, Armonk, NY, USA). Oneway analysis of variance (ANOVA) was used to examine for significant differences. Unit variance (UV) scaling was performed as the pretreatment step in MVDA via the SIMCA-P software. PCA was used as the initial pattern recognition technique for visualization of the clustering of samples in score scatter plots. Correlation of samples and bioactivities was observed in the biplot of PLS. Hotelling's T2 was among the diagnostic tools used to identify strong outliers. R 2 (goodness of fit) and Q 2 (goodness of prediction) estimates were used to validate the PLS model. e predictive ability of the PLS model was also tested using crossvalidated analysis of variance (CV-ANOVA) and the 100 random permutation function. Variable importance in projection (VIP) values ≥1 served as the cut-off point for the identification of the top influential signals for the PLS model [43].

Effects of Extraction Parameters on TPC.
e parameters indicated a range of extraction yields for crude C. nutans extracts, as shown in Table S1. e TPC of crude extracts (Table 1) was impacted profoundly by solvent polarity, S-L ratio, extraction duration, and the number of extraction cycles. e 25% ethanol crude extract had the highest TPC (45.99 mg GAE/g dw), and the value was statistically significantly higher (p < 0.05) than that of all other solvent polarities.
Other studies reported a similar result [44,45] in which C. nutans extract produced by reflux extraction yielded the highest TPC at the same solvent polarity (25% ethanol). is higher TPC yield was because the higher water content in the 25% ethanol system caused the plant material to swell, increasing the contact between solvent and plant matrix [46]. Although polyphenols are mostly polar, a highly polar solvent system, such as water, commonly produces extracts with a high level of impurities such as sugars, soluble proteins, and organic acids and can hinder the extraction of phenolic compounds [47][48][49]. In the current study, higher ethanol concentrations (50-100% ethanol) decreased TPC value, and absolute ethanol (100%) reduced the phenolic yield in the extract. e results were expressed in mean ± standard error of the mean (SEM), where n � 3. Values marked with the same letter subscript represent comparison within the same extraction parameter, while values with the same superscripts portrayed significant differences (p < 0.05) within samples of the same parameter set. Statistical significance was analyzed using SPSS one-way ANOVA.
Based on these results, the 25% ethanol solvent system was suitable for extracting the optimal number of polyphenols. All samples differed significantly in TPC value (p < 0.05), except for the 50% and 75% ethanol extracts. erefore, when testing for different S-L ratios in the subsequent experiments, 25% ethanol was used.
When 25% ethanol extracts, with the five different S-L ratios, were tested, the sample extracted at 1 : 15 S-L ratio had the highest TPC value (47.58 mg GAE/g dw). e TPC value increased as the S-L ratio increased to 1 : 15. However, from that point, the value started to decrease.
is change may be due to the force driven by the concentration gradient between solid and solvent, as described by principles of mass transfer. A higher S-L ratio increases the concentration gradient and diffusion rate to Evidence-Based Complementary and Alternative Medicine allow greater extraction of solids (phenolic compounds) by solvent until equilibrium is reached [22]. e TPC value of the 1 : 15 S-L ratio sample was significantly greater (p < 0.05) than that of the other S-L ratio samples tested, except for the 1 : 10 S-L ratio sample. e selection of the 1 : 15 S-L ratio for subsequent testing was due to its higher TPC value.
When 25% ethanol extracts, with an S-L ratio of 1 :15, were tested through three extraction cycles, the TPC of the highest value (58.52 mg GAE/g dw) came from the 30 min duration extract. e TPC value increased significantly with increasing duration of extraction (p < 0.05) up to 30 min but declined thereafter. In addition to the overheating effect from extended extraction time, long-term exposure to light and oxygen in the extraction medium could cause phenolic oxidation [50]. Additionally, Fick's second law of diffusion could be at work in longer durations (i.e., solute concentrations in the plant matrix and solvent system could reach final equilibrium [51]). ese results suggest that prolonged extraction duration does not increase the yield of phenolic compounds.
Extraction usually is repeated several times, and extracts are then combined to enhance the extraction efficiency [52,53]. When testing the 25% ethanol extract with an S-L ratio of 1 : 15 and 30-minute extraction duration with different numbers of extraction cycles, multiple cycles of extraction were found to improve the extraction yield of polyphenols. ree extraction cycles produced the highest TPC value (45.49 mg GAE/g dw), which was significantly greater than that of the cycle numbers (p < 0.01). In another study, three extraction cycles showed greater yields of TPC and TFC than two and four extraction cycles [54]. Hence, three extraction cycles were optimal for extracting polyphenol compounds from C. nutans leaves.

Effects of Extraction Parameters on CC Activity.
CCs have the potential in salvaging the damaged heart's manipulation of the cells in vitro and boost its function in repairing damaged myocardium and can serve as an integrative strategy in increasing the therapeutic effects and outcome. Solvent polarity, S-L ratio, extraction duration, and the number of extraction cycles had profound impacts on the viability of CCs (Table 1). To examine the cytotoxicity of the C. nutans leave extracts, CCs were treated with extracts for 96 h (Figures 1(a)-1(d)). C. nutans was not cytotoxic to CCs. Instead, it promoted cell growth, even at its highest concentration (1000 µg/mL). Cells treated with the C. nutans aqueous extract (0% ethanol) at 100 µg/mL had significantly greater viability than cells treated with the other concentrations (p < 0.05). However, CC viability decreased with increasing ethanol concentrations, demonstrating a dose-dependent relationship (Figure 1(a)). e suggestion is that highly polar polyphenols in the aqueous extract were more favorable to the survival and growth of CCs. As the S-L ratio increased, the viability of CCs declined (Figure 1(b)). CCs treated with the 1 : 5 S-L ratio treatment exhibited the highest viability-significantly greater than that of the control (p < 0.05). is indicates that polyphenols, extracted at lower S-L ratios and lower ethanol concentrations, were more favorable to CC growth.  (Figures 1(c) and 1(d)). e highest CC viability was observed with C. nutans extract single extraction in 10 min. ese observations collectively suggest that extracts with the highest TPC do not necessarily contribute to the highest cell viability. Nonetheless, a shorter exposure time to the ethanol-containing extraction solvent results in less solvent contact with the leaves could also contribute to better CC viability. It appears that not only increased extraction cycles produce extracts with more compounds or impurities that could affect CC viability [55], the exposure time to the extraction solvent could also affect CC bioactivity.

FTIR Spectrum
Analysis. C. nutans leave extracts, scanned through ATR cells, had their spectra compared within their respective parameters to assess the functional group differences and similarities among the extracts. Based on the comparison among parameter 1 (solvent polarity) tested extracts, several peaks were similar for all samples. ese include peaks assigned to O-H stretching (3500-3200 cm -1 ), asymmetric stretching of C-H (2935-2915 cm -1 ), C�C aromatic stretch (∼1600 cm -1 ), O-H bending of phenol or tertiary alcohol groups (1410-1310 cm -1 ), primary alcohol groups (∼1050 cm -1 ), C-H bending of end methylene (∼900 cm -1 ), and C-H bending of benzene rings (∼820 and ∼770 cm -1 ). A profound difference was observed between the aqueous and ethanol C. nutans extracts (Table 2), with 10 and 11 known peaks, respectively. Symmetrical C-H stretching was absent in aqueous extract. Ester (υ C�O) peaked at 1726 cm −1 and differentiated the absolute ethanol extract from all other parameter 1 tested samples. Furthermore, the presence of C-O stretching from ester or glycoside near 1246 cm −1 made the absolute ethanol extract more distinctive than the others. IR analysis indicated the presence of aromatic compounds in all extracts attributed to the aromatic benzene ring stretching (C�C) at ∼1600 cm −1 and C-H bending (∼820 cm −1 and ∼770 cm -1 ). Nitro groups (asymmetrical N-O stretching) of aromatics at ∼1514 cm −1 were observed in both aqueous and aqueous-ethanolic extracts. Meanwhile, symmetrical N-O stretching of the aromatic nitro groups (1360-1290 cm -1 ), visible in the IR peak of the aqueousethanolic extracts (25-75%), became a prerequisite for identifying aqueous-ethanolic C. nutans extracts. Functional groups such as hydroxyl, carbonyl, and aromatic rings indicated in C. nutans extracts, characterized the presence of phenolic and flavonoids compounds.
In parameter 2 tests (solid-liquid ratio [S-L ratio] variation), there were six major peaks identified for all samples (Table S2) C-H stretching of alkene (3091 cm -1 ) was observed in the 1 : 25 S-L ratio extracts. However, secondary alcohol groups (∼1070 cm -1 ) and end methylene (∼900 cm -1 ) were not found in the 1 : 25 S-L ratio extract.
Extracts produced from different extraction times (parameter 3) did not vary significantly. Twelve peaks were assigned into their respective functional groups (Table S3). ere was no difference in chemical profiles found among the samples extracted from shorter periods (10, 20, and 30 min). However, extracts with a longer extraction time (50 min) lacked secondary alcohol groups (υ C-O) at ∼1070 cm −1 . is lack of secondary alcohol groups was probably due to the effect of overheating by sonication [56], signifying that prolonged sonication could potentially cause the decomposition of bioactive constituents. e number of extraction cycles (parameter 4) had slight effects on the FTIR fingerprint of the samples (Table S4). Extracts extracted at one and two-cycle(s) lacked benzene ring (δ C-H) at ∼770 cm −1 , which might be due to short sonication exposure. Bending C-H of the benzene ring (∼770 cm -1 ) was found in extracts where the number of extraction cycles increased to 3, 4, and 5. However, long-term exposure to ultrasonic cavitation can cause certain compounds to degrade. For instance, crude extract acquired from five extraction cycles eliminated O-H stretching (3500-3200 cm -1 ) of hydroxyl and C-O stretching (∼1070 cm -1 ) from secondary alcohol groups. 3cycled extraction was shown as the optimal condition to attain the most functional compound groups in C. nutans extract. A peak of C�O ester stretching at 1202 cm −1 was only found in the 3-cycled sonication extract, suggesting that the number of extraction cycles can alter the extract's quality.

Classifications of Samples of Each Parameter Set by PCA.
PCA was conducted using the processed ATR-FTIR data to illustrate the pattern of distribution of samples within each parameter set (Figures 2(a)-2(d)). e first two principal components in the first parameter set yielded a total variation of 87.6%, followed by the second parameter with 79.4% of the total variation. PCA for the third and fourth parameters achieved maximum variation of 82.3% and 89.8%, respectively. Based on PCA, all tested groups within each parameter set exhibited an acceptable value of goodness of fit (R 2 X [cum] � 0.876, 0.794, 0.822, and 0.898, respectively) and predictability (Q 2 [cum] � 0.826, 0.681, 0.662, and 0.709, respectively). e difference between R 2 X (cum) and Q 2 (cum) was within 0.2, confirming that the PCA models for each parameter set were reliable.
For the first parameter set (solvent polarity) in the score scatter plot, PC1 separated G1 (aqueous), G2 (25% ethanol), and G3 samples (50% ethanol) from G4 (75% ethanol) and G5 (100% ethanol) samples (Figure 2(a)). An assessment was made of Hotelling's T2 values to identify the presence of any strong outliers at two levels, specifically at 95% and 99% confidence  Evidence-Based Complementary and Alternative Medicine intervals. Among the tested samples, none were identified as outliers as they all fell within the interval of 99%. It was deducible that based on the polarity of polyphenols, PC1 made a differentiation because G1 to G3 samples may possess a similar majority of highly polarised polyphenols. e close clustering within G1 to G3 group samples was not distinctively different from the unsupervised PCA model. However, the G4 and G5 samples discriminated from the rest of the groups. For the second parameter set (S-L ratio), only G6 samples (1 : 5 S-L ratio) were distinctively different from the other group samples based on PC1. G2, G7, G8, and G9 samples (1 : 15, 1 : 10, 1 : 20, and 1 : 25, respectively) overlapped themselves in the PCA plot, which shows that these samples shared similar compounds to a significant degree. Hotelling's T2 showed no strong outliers present in these groups.
For samples from the third parameter set (extraction times), G2 (30 min extraction) was the only well-clustered group. PC1 divided the samples into two major parts with G10 (10 min extraction), G12 (40 min extraction), and G13 (50 min extraction) in the positive region and G11 (20 min extraction) and G13 in the negative region of the plot. One sample from G10 and G13 overlapped with G12 samples. ese could be described as moderate outliers, which may have arisen from technical error. Meanwhile, no strong outliers were observed.
Groups from the number of extraction cycles (parameter 4) were closely clustered with each other. G17 (5 extraction cycles) samples tended to overlap with G2 (3 extraction cycles) samples, and this could be a moderate outlier. A single strong outlier was detected in G16 (4 extraction cycles) samples, as it fell outside the 95% Hotelling's T2 ellipse, but the datum was accepted in the analysis because the confidence interval was within 99%.

Correlations of Sample Profile with Bioactivities via PLS.
All tested parameter sets were subjected to PLS analysis. However, only the first parameter set (solvent polarity) passed the internal R 2 Y and Q 2 cumulative cross-validation stage in SIMCA-P software. Validation of the PLS model for solvent polarity occurred, with acceptable goodness of fit (R 2 Y � 0.801) and good predictive power (Q 2 cumulative � 0.741). e biplot of the PLS model revealed the distribution of samples and their correlation to the TPC and CCs activities. Data showed that G1 and G2 clustered closely for both of the activities (Figure 3(a) (d) Figure 2: (a-d) PCA-X score scatter plots of parameters 1-4 (solvent polarity, solid-liquid ratio, extraction duration, and the cycle of extraction) at PC1 and PC2.
Evidence-Based Complementary and Alternative Medicine 7 total of 1786 signals with VIP ≥1. e first tenth of VIP values ≥1 was selected as the top influential signals contributing to the PLS models at PC1 and PC2. In addition to the evaluation of VIP, CV-ANOVA was employed to confirm the significance of the prediction model in the PLS model [57]. CV-ANOVA showed that the correlation between the chemical profile and TPC and CCs activity was significant (p < 0.05). e indication was that the solvent polarity used in extraction could significantly affect the chemical profile and bioactivities. e PLS biplot was affirmed with 100 random permutations to assess the reliability of the PLS model [57]. Based on the result, the Q 2 regression line that is colored green was negatively intercepted (Figures 3(b) and 3(c)). Furthermore, R 2 and Q 2 plots were well intercepted to each other in activities, particularly in TPC, and were in line with a previous study by Chew et al. [43]. erefore, the effect of solvent polarity on bioactivities presented in the PLS model is reliable for the study.

Conclusion
e application of the ATR-FTIR chemometrics approach revealed sample variation in each tested parameter set via their respective PCA models. PLS model demonstrated the relationship among samples from the first parameter set and their correlation to TPC and cytotoxicity results. From the chemometric-based optimization techniques, solvent polarity was the most influential factor affecting the observed chemical differences and bioactivities of C. nutans extracts. G1 (aqueous, S-L ratio of 1 : 15, and three cycles of extraction) and G2 (25% ethanol, S-L ratio of 1 : 15, and three cycles of extraction) were identified as the optimum conditions for extracting C. nutans extracts with the ideal TPC values and no toxicity to CCs. MVDA results were subjected to several validation tools such as R 2 and Q 2 , VIP, CV-ANOVA, and permutation plots to confirm the reliability of the PLS model. ese approaches proved effective  Evidence-Based Complementary and Alternative Medicine as a guide to optimize postharvest processing parameters and provide extracts with the greatest TPC bioactivities and CC viability.

Data Availability
All the data generated or analyzed in this study are included in this manuscript. e ATR-FTIR data for parameters 2, 3, and 4 used to support the findings of this study are included within the supplementary information file.

Conflicts of Interest
e authors declare that they have no conflicts of interest. Table S1: extraction yield (% w/w) of all the parameters. Table S2: list of functional groups reported in different solidliquid (S-L) ratios of C. nutans extracts along with their tentative identification (parameter 2). Table S3: major similarity peaks of functional groups reported in the duration of C. nutans extraction times along with their tentative identification (parameter 3).