Scrophulariae Radix (SR) is one of the oldest and most frequently used Chinese herbs for oriental medicine in China. Before clinical use, the SR should be processed using different methods after harvest, such as steaming, “sweating”, and traditional fire-drying. In order to investigate the difference in chemical constituents using different processing methods, the two-dimensional (2D) 1H-13C heteronuclear single quantum correlation (1H-13C HSQC)-based metabolomics approach was applied to extensively characterize the difference in the chemical components in the extracts of SR processed using different processing methods. In total, 20 compounds were identified as potential chemical markers that changed significantly with different steaming durations. Seven compounds can be used as potential chemical markers to differentiate processing by sweating, hot-air drying, and steaming for 4 h. These findings could elucidate the change of chemical constituents of the processed SR and provide a guide for the processing. In addition, our protocol may represent a general approach to characterizing chemical compounds of traditional Chinese medicine (TCM) and therefore might be considered as a promising approach to exploring the scientific basis of traditional processing of TCM.
Two-dimensional (2D) 1H–13C heteronuclear single quantum coherence (HSQC) has been increasingly applied to metabolomics studies because it can greatly improve the resolving capability compared with one-dimensional (1D) 1H NMR. However, preprocessing methods such as peak matching and alignment tools for 2D NMR-based metabolomics have lagged behind similar methods for 1D 1H NMR-based metabolomics. Correct matching and alignment of 2D NMR spectral features across multiple samples are particularly important for subsequent multivariate data analysis. Considering different intensity dynamic ranges of a variety of metabolites and the chemical shift variation across the spectra of multiple samples, here, we developed an efficient peak matching and alignment algorithm for 2D 1H–13C HSQC-based metabolomics, called global intensity-guided peak matching and alignment (GIPMA). In GIPMA, peaks identified in all spectra are pooled together and sorted by intensity. Chemical shift of a stronger peak is regarded to be more accurate and reliable than that of a weaker peak. The strongest undesignated peak is chosen as the reference of a new cluster if it is not located within the chemical shift tolerance of any existing peak cluster (PC), or otherwise it is matched to an existing PC and the aligned chemical shift of the PC is updated as the intensity-weighted average of the chemical shifts of all peaks in the cluster. Setting an optimum chemical shift tolerance (Δδo) is critical for the peak matching and alignment across multiple samples. GIPMA dynamically searches for and intelligently selects the Δδo for peak matching to maximize the number of valid peak clusters (vPC), that is, spectral features, among multiple samples. By GIPMA, fully automatic peakwise matching and alignment do not require any spectrum as initial reference, while the chemical shift of each PC is updated as the intensity-weighted average of the chemical shifts of all peaks in the same PC, which is warranted to be statistically more accurate. Accurate chemical shifts for each representative spectral feature will facilitate subsequent peak assignment and are essential for correct metabolite identification and result interpretation. The proposed method was demonstrated successfully on the spectra of six model mixtures consisting of seven typical metabolites, yielding correct matching of all known spectral features. The performance of GIPMA was also demonstrated on 2D 1H–13C HSQC spectra of 87 real extracts of 29 samples of five Dendrobium species. Hierarchical cluster analysis (HCA) and principal component analysis (PCA) of the 87 matched and aligned spectra by GIPMA generates correct classification of the 29 samples into five groups. In summary, the proposed algorithm of GIPMA provided a practical peak matching and alignment method to facilitate 2D NMR-based metabolomics studies.
Abstract We presented a strategy utilizing 2D NMR‐based metabolomic analysis of crude extracts, categorized by different pharmacological activities, to rapidly identify the primary bioactive components of TCM. It was applied to identify the potential bioactive components from Scutellaria crude extracts that exhibit anti‐non‐small cell lung cancer (anti‐NSCLC) activity. Four Scutellaria species were chosen as the study subjects because of their close phylogenetic relationship, but their crude extracts exhibit significantly different anti‐NSCLC activity. Cell proliferation assay was used to assess the anti‐NSCLC activity of four species of Scutellaria . 1 H− 13 C HSQC spectra were acquired for the chemical profiling of these crude extracts. Based on the pharmacological classification (PCA, OPLS‐DA and univariate hypothesis test) were performed to identify the bioactive constituents in Scutellaria associated with the anti‐NSCLC activity. As a result, three compounds, baicalein, wogonin and scutellarin were identified as bioactive compounds. The anti‐NSCLC activity of the three potential active compounds were further confirmed via cell proliferation assay. The mechanism of the anti‐NSCLC activity by these active constituents was further explored via flow cytometry and western blot analyses. This study demonstrated 2D NMR‐based metabolomic analysis of pharmacologically classified crude extracts to be an efficient approach to the identification of active components of herbal medicine.