Deuterated methanol-d₄ (CH₃OH-d₄), non-deuterated KH₂PO₄, sodium deuterated oxide (NaOD), trimethylsilyl propionic acid-d₄ sodium salt (TSP), deuterium oxide (D₂O) and methanol were purchased from Merck (Darmstadt, Germany). Liquid nitrogen was supplied by MOX-Linde Sdn. Bhd. (Petaling Jaya, Malaysia). Standard compounds such as andrographolide, neoandrographolide and deoxyandrographolide were previously isolated by Professor Dr. Johnson Stanslas from Pharmacotherapeutics Unit, Department of Medicine, Faculty of Medicine and Health Sciences, Universiti Putra Malaysia. RAW264.7 cells were purchased from American Type Culture Collection (ATCC, Rockville, USA) under the catalogue number of ATCC® TIB-71TM. Dulbecco’s Modified Eagle’s Medium (DMEM) with and without phenol red, fetal bovine serum (FBS), penicillin streptomycin (Pen Strep) and TrypLETM Express were purchased from Gibco by Life Technologies Inc. (Eggenstein, Germany). Phosphate buffer saline (PBS) was purchased from Invitrogen by Life Technologies. Lipopolysaccharides (LPS) from Escherichia coli 0111:B4 and interferon-gamma (IFN-γ) were purchased from Sigma Aldrich Co. (St. Louis, MO). DMSO was purchased from Merck (Darmstadt, Germany). PBS and 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) were purchased from Invitrogen by Life Technologies Inc. (Eggenstein, Germany) and DMSO was purchased from Merck (Darmstadt, Germany).

The A. paniculata seeds from three germplasm accessions of Malaysia (Accession 11,265, 11,341 and 11,248 which will be designated as H, P and T, respectively) were collected in July 2011 from Ladang Dua, Universiti Putra Malaysia (UPM). The germination and cultivation procedure in this study was performed according to Yusof et al. (2014) [12] . Leaf samples were harvested in January, February and March 2012 assigned as 120, 150 and 180 DAT at two different time sessions (morning: 7:30 - 10:30 am; evening: 2:30 - 5.30 pm). Each sample was quenched with liquid nitrogen and ground into powder form before stored at −80^°C prior to drying by freeze-dry system (Labconco, Texas, USA) to remove the water content.

The FTIR sample preparation and measurement in this study were performed according to Yusof et al. (2014) [12] . The freeze-dried leaf powder of A. paniculata (150 mg) was extracted by sonication in 5 mL commercial-grade methanol (30 min, 40˚C). The extraction step was repeated twice and the total combined supernatant was filtered through filter paper (Whatman, 125 mm) and evaporated using miVac Modular Concentrator series (Genevac Ltd, Ipswich, UK). The crude extracts were then kept at −80˚C until IR analysis. The crude samples were freeze-dried again before submitting to IR analysis to prevent absorbance interference caused by the water content. IR spectra were recorded using a Spectrum 100 FTIR spectrometer (PerkinElmer, Waltham, MA, USA) with a single reflectance horizontal ATR cell equipped with a zinc selenide crystal. The analyses were made in absorption mode on a small amount of crude sample deposited on the ATR crystal. The spectra of five biological replicates and pure compounds in three technical replicates were obtained in the frequency range 4000 - 280 cm⁻¹ with a resolution of 4 cm⁻¹ and a total accumulation of 16 scans.

Freeze-dried plant leaf (20 mg) was transferred to a 2.0 mL Eppendorf tube. A volume of 1.5 mL of a mixture of KH₂PO₄ buffer (pH 6.0) in D₂O containing 0.05% trimethylsilylpropionic acid sodium salt (TSP, w/w) and methanol-d4 (1:1) was added to the plant sample. The extraction procedure in this study was performed according to Kim et al. (2010), with a minor modification [13] . The Eppendorf tubes were vortexed for 1 minute at room temperature and ultrasonicated for 20 minutes. This step was followed by centrifuging for 10 minutes at 10,000 rpm to remove the insoluble components from which a 0.75 mL clear supernatant was transferred into the 5 mm NMR tube. ¹H NMR spectra were recorded at 26˚C on a 500 MHz Varian INOVA NMR spectrometer operating at a proton NMR frequency of 499.887 MHz. The acquisition time of each ¹H-NMR measurement was 3.53 minutes, which consisted of 64 scans. A presaturation sequence was used to suppress the residual H₂O signal. The resulting spectra were manually phased and baseline was corrected using the Chenomx software (v.5.1, Alberta, Canada) with a consistent setting for all spectra. Spectral intensities were binned by equal width (δ 0.04) corresponding to the region of δ 0.50 - 10.00. The region of δ 4.70 - 4.90 representing water, and δ 3.23 - 3.36 representing residual methanol were excluded. 2D J-resolved was recorded for additional support in the metabolite identification.

Andrographolide and neoandrographolide were selected from the ¹H NMR spectra and their concentrations are listed in Table 1. Deoxyandrographolide, α-glucose, β-glucose, sucrose, choline and analine were excluded from the quantification due to their overlapping in the 1D NMR spectra. Relative quantification of the identified metabolites

from the A. paniculata was based on the mean peak area of the ¹H NMR signal of interest after binning using Chenomx software. Data was analyzed using the Minitab Statistical Software 14.12.0, © 2004 Minitab Inc. PA 16801-9928, USA. Among the groups, significance test was performed using the one-way ANOVA at 5% significance level.

The anti-inflammation activity of A. paniculata crude extracts was evaluated through Nitric oxide (NO) suppression activity based on Lee’s method with some modifications [14] . Each A. paniculata crude extract of 2 mg was dissolved in 5 μL of DMSO/995 μL DMEM without phenol red. It was then sonicated for 10 minutes without heating, and centrifuged for 10 minutes at 1000 rpm. The solubility of 50% for the methanol crude extract was determined earlier as to estimate the actual concentration of the prepared crude extracts. From the resulted supernatant, 260 μL was pipetted out into 1040 μL DMEM without phenol red as to minimize the use of DMSO in the sample preparation. This is the stock extract which was further diluted into desired concentrations via serial dilutions. Two stock concentrations of 1000 and 10,000 μg/mL were subjected to the two-fold dilution in the screening stage. For the actual assay, the concentrations of A. paniculata extracts used were in the range of 100 μg/mL - 1 mg/L in a serial dilution.

DMSO of 0.1% in DMEM without phenol red was utilized as a control vehicle in the serial dilutions in order to maintain the final concentration of 0.1% DMSO in each well. These various concentrations were prepared separately in the aseptic microcentrifuge tubes before being dispensed into the wells.

The RAW264.7 murine macrophages were cultured in 75 cm² plastic culture flask in DMEM with phenol red containing HEPES, L-glutamine supplemented with 10% FBS and 1% penicillin/streptomycin in 5% CO₂ atmosphere at 37˚C. After the cells have reached the confluency of 80% - 90%, they were detached from the flask wall by rinsing off the flask with 5 ml of PBS solution followed by 2 ml of TrypLE^TM Express before incubated for 10 minutes. DMEM with phenol red of 4 mL was used to rinse out the detached cells which then centrifuged at 3000 rpm, 4˚C for 10 minutes by using Rotina 35R (Andreas Hettich GmbH & Co., Tuttlingen, Germany). The supernatant was then removed and the cells were re-suspended with fresh DMEM containing HEPES, and L-glutamine with 5% FBS but without phenol red. The cells were counted by using a standard trypan blue cell counting technique.

The cell concentration was adjusted to 1 × 10⁵ cells/mL in the same medium (DMEM without phenol red) before seeded in a 96-well flat bottom tissue culture plates at initial seeding densities of 5 × 10⁴ cells/50 μL per well. The ready plates with seeded cells were incubated at 37˚C (5% CO₂ atmosphere) overnight to allow anchoring of cells to the wells.

The treatment of the cells was done by replenishing the medium in each well with 50 μL DMEM without phenol red containing triggering agent of 1 ng/mL recombinant murine IFN-γ and 1 μg/mL LPS. For each dilution of the plant extracts, 50 μL was added into the appropriate wells in four replicates. The assay plate was incubated for 17 hours at 37˚C in 5% CO₂ atmosphere. There were five kinds of controls prepared in duplicates for each plate. The controls were media only without cells (blank control), cells in media only, cells in media containing triggering agent plus 0.1% DMSO (negative control), cells in media containing triggering agent only, and cells containing triggering agents plus curcumin (positive control).

Nitrite concentration was determined by Griess assay through the reaction of 50 μL of cell culture supernatant and 50 μL of Griess reagent (1% sulphanilamide and 0.1% N-(1-naphthyl)ethylenediamine dihydrochloride in 2.5% phosphoric acid) at room temperature. The optical density was measured at 550 nm after 5 minutes of incubation at room temperature by Spectramax Plus (Molecular Devices) UV/Vis microplate reader. The accumulated NO release was determined spectrophotometrically by the total nitrite in the medium. The nitrite concentration in the samples was determined by regression analysis using serial dilutions of sodium nitrite as a standard. The percentage inhibition was calculated based on the ability to inhibit nitrite below the level produced by cells cultured in the presence of triggering agents and DMSO, which was considered as 0% inhibition.

In order to avoid any false positive result of NO inhibition, the treated cells were also submitted for MTT assay as to measure any cytotoxic effect on the cells by using the Heras’s method [15] . Supernatant in all wells was removed followed by the addition of 100 μL DMEM with phenol red. Subsequently, 20 μL of 3-(4,5-dimethylthiazol-2-yl)- 2,5-diphenyltetrazolium bromide (MTT, 5 mg/ml) was added to each well and the plate was incubated in 5% CO₂ atmosphere at 37˚C for 4 h. The mixture of culture media and MTT in all wells were discarded and the purplish formazan crystals formed were dissolved in DMSO and further incubated for 5 minutes. The absorbance of colour intensity was then measured at 570 nm at the room temperature wherein the absorbance of formazan in the control of untreated cells was taken as 100% viability.

Data were analyzed by using GraphPad Prism software version 5.01 (GraphPad Software Inc., San Diego, USA). The graph of half maximum inhibitory concentration (IC₅₀) for each sample was obtained from logarithm graph plot. In which each graph showed the goodness of fit, R² > 0.95 and within the 95% of the confidence interval with p-values < 0.05 were considered significant. All results were expressed as the mean ± S.E.M of n, where n represented the number of wells replication. Differences between the mean values for different groups were analyzed by Analysis of variance (ANOVA) using SPSS 16.0 (Statistical Package for the Social Sciences).

Multivariate data analysis by partial least-squares analysis (PLS) was performed with the SIMCA 13.0 software using scaling based on Variance.

In the present work, metabolomics analysis was applied to discriminate the metabolite variation among the three accessions of P, T and H. Metabolites from this plant including the major andrographolide, deoxyandrographolide, neoandrographolide and carbohydrates (α-glucose, β-glucose and sucrose) as well as amino acids (alanine and choline) were successfully identified. Figure 1 shows the ¹H-NMR spectra of A. paniculata leaf extracts displaying both the qualitative and quantitative variations of some metabolites. The compounds were identified based on the 2D J-resolved (Figure 2) and comparison with the reported NMR signals of those isolated metabolites from this plant. Other compounds were assigned by comparison with the NMR spectra of the reference compounds measured under the same conditions as extracts. All of the spectra showed signals in the aliphatic (δ 0.5 - 3.0), aromatic (δ 5.5 - 8.0) and carbohydrate region (δ 3.0 - 5.5).

Table 1 lists the characterized metabolites which could be recognized by signals shown in those three main regions. The peaks in at δ 0.64 and δ 1.00 were typical signals of the neoandrographolide aliphatic protons, whereas andrographolide signals were detectable at δ 0.72, 1.16, 1.32 and 1.92. Some amino acids namely alanine and choline were also observed in the aliphatic region at 1.48 and 3.20, respectively. Alpha-glucose, β-glucose and sucrose were ascribed based on the signals displayed in the carbohydrate region of δ 3.0 - 5.5. Other andrographolide signals were observed at δ 4.09 and 4.54. As for deoxyandrographolide, the signal was detectable at δ 4.16. Andrographolide and

neoandrographolide signals were also observed in the aromatic region of δ 6.91 and 7.40 respectively.

Relative quantification of the identified major compounds which are andrographolide and neoandrographolide in A. paniculata extracts was based on the mean peak area of the ¹H NMR signals from the three different accessions (P, T and H), which were harvested at three different ages (120, 150 & 180 DAT) and two different time sessions of morning and evening. A. paniculata accessions from different ages were able to be differentiated by the intensities of andrographolide (binned at 0.74 ppm) and deoxyandrographolide (binned at 0.62 ppm). The significance (p < 0.05) difference between the samples at different age for each accession was conducted by One-Way ANOVA using MINITAB 14. Morning and evening session data for 120, 150 & 180 DAT of three different accessions (P, T and H) were combined in this relative quantification analysis.

The concentration of andrographolide at 180 DAT of P accession was significantly higher (p < 0.05) as compared to 120 and 150 DAT as shown in Figure 3(a). Similarly, T accession of the harvesting age at 180 DAT showed the highest concentration of andrographolide (Figure 3(b)), but with no significant difference when compared to 120 and 150 DAT. The highest concentration of andrographolide for H accession was observed at 150 DAT instead of 180 DAT with insignificant difference between the three

DATs (Figure 3(c)) (Table 1). The concentration of neoandrographolide as in Figures 3(a)-(c) was observed to be more abundant in 180 DAT harvesting age for H and T accessions, whereas 150 DAT harvesting age for P accession. However, the p value exhibited no significance between all DATs for each accession (Table 1).

Investigation on the in vitro anti-inflammatory effect of the A. paniculata accessions methanol crude extracts was carried out via nitric oxide assay. A. paniculata has been used traditionally for inflammation since ages [16] . Thus, based on the claims, this is an attempt made in identifying the possible metabolite changes due to the harvesting time and age variations in relation to the anti-inflammation activity. The correlation of the quantitative IC₅₀ values from NO assay with the data of FT-IR and FT-NMR has paved the way in identifying the possible biomarkers.

Since the A. paniculata leaf extracts have shown promising results at the screening stage, the IC₅₀ values were then further determined as summarized in Table 2 and Figure 4. The potential of the A. paniculata crude extracts as anti-inflammatory agents

Standard: IC₅₀ curcumin = 14.69 μM; P = Pasir Puteh, T = Tanjung Ipoh, H = Harapan.

were assessed from the inhibition of the NO production without the induction of cell death. Thus, another test was carried out in order to determine the cell viability of the induced macrophages.

It could be observed that all the three accessions yielded the highest NO activity at 180 DAT harvesting age of morning session for T and H, and evening session for P with IC₅₀ values of 13.99, 10.61 and 18.51 µg/mL, respectively. In relation to the NO inhibition, P and H accessions gave a significant difference between 120 and 150 DAT, but no significant difference between the morning and evening harvesting times. However, T accession exhibited a significant variation between morning and evening harvesting time for each of the harvesting age.

Generally, the small IC₅₀ value is the lowest concentration that could inhibit the production of NO by at least 50%. However, the statistical computation only recognizes big values in representing a significant activity. Thus, prior to any MVDA computation step, the IC₅₀ values were converted into their antiradical scavenging activity (1/IC₅₀) readings in order to able the statistically translation of the activity.

PLS was applied using a validation model with a degree of over fit between the variables and the responses [17] . PLS was performed to obtain a more precise view on the correlation between the NO assay and the phytochemical present with the effect of harvesting age variations of 120, 150 and 180 DAT. This approach was implemented because PLS has greater success in linking NO assay with the phytochemicals and in providing models for prediction of the highest bioactivity among the different harvesting ages [17] [18] .

Three clusters were obtained from the PLS bi-plot that combined the score and loading plot of the FTIR data (Figure 5). A differentiation analysis of the metabolites between the 180 DAT harvesting age samples from the other two different ages of the

three A. paniculata accessions (P, T and H) correlated to NO assay was accomplished (Figure 5). Among the three accessions, P showed that 180 DAT harvesting age was strongly correlated with the NO inhibition activity (Figure 5(a)).

The other two harvesting ages were projected to the negative side of the plot by PC1 which was opposite of the NO. The peak signals were characterized as those belonging to the major compounds of andrographolide and neoandrographolide. The same data distribution was observed over the ellipse for the T accession of different harvesting ages model (Figure 5(b)). The closest harvesting age to the NO inhibitory property was 180 DAT. This group was projected nearer to the NO in the positive quadrant of both PC1 and PC2.

There was a different observation for H accession, whereby all the three harvesting ages were well-separated in groups (Figure 5(c)). The 180 DAT was well discriminated from the other two harvesting ages in the positive side of PC2 which gave 180 DAT as the nearest to having NO inhibition property.

A. paniculata samples were analyzed with NMR spectrometer to support the data obtained from FTIR analysis. The PLS loading bi-plot data gave 180 DAT of P and T accessions as well separated from 120 and 150 DAT by PC1 (Figure 6(a) and Figure 6(b)), whereby the 180 DAT of P and T were strongly correlated to NO activity. The characteristic of 180 DAT was observed to be having higher level of the major compounds of andrographolide, neoandrographolide and deoxyandrographolide which could be correlated to the NO inhibition property. This is in agreement with the obtained anti-inflammatory activity data in which exhibited that the 180 DAT of P and T accessions gave the lowest IC₅₀ value for NO. Other metabolites that contributed to the clustering were carbohydrates and amino acids which are shown in the loading bi-plot of Figure 6(a) and Figure 6(b).

H accession gave a different observation wherein both 120 and 180 DAT were close to NO assay compared to 150 DAT as shown in Figure 6(c). Even though both of the DATs of 120 and 180 situated close to the NO activity, the clustering of these two could be established as 120 DAT was closest to NO by PC1 whereas 180 DAT by PC2 with 180 DAT having the highest NO inhibition. Metabolites that contributed to this separation were andrographolide, neoandrographolide and carbohydrates with higher content of neoandrographolide and amino acids in 120 DAT.

The correlation between anti-inflammatory and metabolite variations among the three different harvesting ages of P, T and H accession using NMR spectroscopy in this study is agreeing with the data obtained from the FTIR analysis. Both of the spectroscopic data are in congruent when correlated to the NO inhibition activity.

PLS model validation was performed by checking the Q² and R² accumulative means after cross validation and permutation test with 100 permutations (Figure 7). According to Eriksson et al. (2006), the permutation test will generate Q² value higher than 0.5, R² value must always be larger than the value of Q² for a good model [19] . Thus, with

the higher significance and lower variation of the NMR-adopted-permutation has helped each model to be validated as a good model. In the current study, the Q² and R² values were higher than 0.5 (P accession: R²Y = 0.75833, Q²Y = 0.71733; T accession: R²Y = 0.84973, Q²Y = 0.77553; H accession: R²Y = 0.75833, Q²Y = 0.57511) which suggested that all of the models meet the criteria of the validation and prediction performances.