Skip to main content

Metabolomics combined with network pharmacology reveals the potential development value of Campanumoea javanica Bl. and its metabolite differences with Codonopsis Radix

Abstract

Campanumoea javanica Bl. (CJ) traditionally used in Southwestern China, is now widely consumed as a health food across the nation. Due to its similar efficacy to Codonopsis Radix (CR) and their shared botanical family, CJ is often used as a substitute for CR. According to the Chinese Pharmacopoeia, Codonopsis pilosula var. modesta (Nannf.) L.T. Shen (CPM), Codonopsis pilosula (Franch.) Nannf. (CP), and Codonopsis tangshen Oliv. (CT) are the primary sources of CR. However, details on the differences in composition, effectiveness, and compositional between CJ and CR are still limited. Besides, there is little evidence to support the application of CJ as a drug. In this study, we employed widely targeted metabolomics, network pharmacology analysis, and molecular docking to explore the disparities in metabolite profiles between CJ and CR and to predict the pharmacological mechanisms of the dominant differential metabolites of CJ and their potential medicinal applications. The widely targeted metabolomics results indicated that 1,076, 1,102, 1,102, and 1,093 compounds, most phenolic acids, lipids, amino acids, and flavonoids, were characterized in CJ, CPM, CP, and CT, respectively. There were an average of 1061 shared compounds in CJ and CRs, with 95.07% similarity in metabolic profiles. Most of the metabolites in CJ were previously unreported. Twelve of the seventeen dominant metabolites found in CJ were directly associated with treating cancer and lactation, similar to the traditional medicinal efficacy. The molecular docking results showed that the dominant metabolites of CJ had good docking activity with the core targets PIK3R1, PIK3CA, ESR1, HSP90AA1, EGFR, and AKT1. This study provides a scientific basis for understanding the similarities and differences between CJ and CR at the metabolome level, offering a theoretical foundation for developing innovative medications from CJ. Additionally, it significantly enhances the metabolite databases for both CJ and CR.

Peer Review reports

Introduction

Campanumoea javanica Bl. (CJ), is known as “Tudangshen” in Chinese and has been utilized in the national medicine system of China for hundreds of years. CJ is a perennial herbaceous dicotyledonous plant widely distributed in the tropical and subtropical regions of eastern Asia. It belongs to the genus Campanumoea, which is part of the family Campanulaceae. CJ is not only a commonly used herb medicine in the Han nationality but also nine minority habitual drugs including Yi, Miao, and Yao. CJ was first recorded in the Qing Dynasty’s An Illustrated Book of Plants [1] and later recorded in the Chinese Pharmacopoeia [2], Guizhou Quality Standards for Traditional Chinese Medicines and Ethnic Herbs [3], Guizhou Standards for the Preparation of Traditional Chinese Medicinal Pieces [4], and Standards for the Preparation of Traditional Chinese Medicinal Pieces of the Guangxi Zhuang Autonomous Region [5]. CJ has been used in Chinese medicine to treat internal injuries due to pulmonary asthenia, spleen deficiency diarrhea, lactation, and infantile enuresis [6]. To date, CJ has been produced into a variety of Chinese patent drugs for national use, such as Benglu Pill, Baidai Pill, Qinghuoyangyuan Capsules, Jinwugutong Capsules, Hulisan Capsules. The Hulisan Capsules is a famous bloodstasis type of Chinese patent drug in Yunnan, selling more than 158 million yuan in 2014 [7,8,9,10]. In addition, CJ is also an important edible plant and is often used as a medicinal dietary ingredient [11]. These biological activities were attributed to its high content of polysaccharides [12]. CJ is naturally distributed widely in Southern China, including Guizhou, Sichuan, Yunnan, and Guangxi provinces [13].

Codonopsis Radix (CR) is a famous traditional Chinese medicine with a long history of use in pharmaceuticals and edibles. It is one of the most popular herbs of traditional Chinese medicine today and is well regarded by the public as a nourishing drug. The CJ and CR demonstrate a certain level of similarity in efficacy and belong to the same family. Consequently, the CJ is often used as a substitute for CR. In fact, the genus of CJ is Campanumoea within the family Campanulaceae, while CR belongs to the genus Codonopsis in the same family. Moreover, CJ and CR are each unique in terms of their efficacy; for example, there is a significant difference between them in the treatment of pediatric enuresis and lactation disorders. The confused application of the two might result in dangerous effects, including medical mishaps. Unfortunately, although a few studies have identified them, effective methods to distinguish them are still lacking. It is difficult to accurately identify its varieties in the process of market circulation. CJ and CR are often mixed into one thing, which causes unnecessary difficulties in the treatment of clinical diseases in Chinese medicine. Up until now, there have been only a few phytochemical studies on CJ. Zhang et al. isolated 6 compounds from CJ, including two lignans (syringin and dangshen glycosides I) and two flavonoids (5-hydroxy-4’, 6, 7-trimethoxyflavone, and 5-hydroxy-4’, 7-dimethoxyflavone) [14]. In addition, zanthocapensol, dangshen glycosides II, and 2 anti-angiogenic active ingredients (4E, 8E, 12E-triene-10-yne-1, 6, 7-tetradecanetriol, and 9-(2-tetrahydropyran)-8E-en-4, 6-diyne-3-nonanol) were isolated from CJ [15]. According to the available literature, more than 50 compounds including glycosides, alkynes, flavonoids, and phenylpropanoid compounds have been isolated from CJ [16,17,18].

Currently, little is known about the phytochemical distinctions between CJ and CR, as well as the potential medical development value of CJ. In this study, ultra-performance liquid chromatography-electrospray ionization-tandem mass spectrometry (UPLC-ESI-MS/MS) integrating chemometrics was applied for metabolite profiling and discrimination of CJ and CR. Targeted metabolomics analysis was used to reveal metabolic differences between CJ and CR in metabolism and to screen for dominant differential metabolites of CJ. The network pharmacology was employed to predict the pharmacological effects of dominant differential metabolites and explore the potential medicinal value of CJ. Furthermore, we used molecular docking techniques to dock their core targets with the active ingredients to confirm the reliability of the network pharmacological predictions. This is a comprehensive investigation of CJ phytochemistry and compares the chemical findings of CJ and CR, providing a better understanding of the differences between CJ and CR, and giving references for the quality control of herbs and the medicinal development of CJ.

Materials and methods

Sample collection

From August 2021 to October 2021, a total of twelve samples were collected in the production regions. CJ samples were collected from Guiding County, Guizhou province, while Codonopsis pilosula (Franch.) Nannf. (CP), Codonopsis pilosula var. modesta (Nannf.) L.T. Shen (CPM), and Codonopsis tangshen Oliv. (CT) were obtained from Weiyuan County (Gansu province), Wen County (Gansu province), and Wuxi County (Chongqing municipality), respectively (Fig. 1). Professor Faming Wu, College of Pharmacy, Zunyi Medical University, identified CJ, CT, CP, and CPM as the roots of Campanumoea javanica Bl., Codonopsis tangshen Oliv., Codonopsis pilosula (Franch.) Nannf., and Codonopsis pilosula var. modesta (Nannf.) L.T. Shen, respectively. Each sample weighed between 1,000 and 1,500 g, and three samples was collected from each region. The samples were promptly brought to the laboratory, washed, and separated into aboveground and subsurface sections that were packed in self-sealing bags. For metabolite extraction, the subterranean section was held at -80 °C. All samples are stored in the laboratory of Zunyi Medical University.

Fig. 1
figure 1

The regional distribution of samples

Sample preparation

Sample preparation was conducted with the method of Zhang KaiXian [19]. In a Scientz-100 F freeze-dryer (Scientz, China), fresh CJ and CR samples were vacuum freeze-dried. Subsequently, they were ground into a powder using a mixer (MM 400, Retsch, Germany) at 30 Hz for 1.5-2 min. Approximately 50 mg of the sample was dissolved in 1.2 mL of 70% methanol extract (Merck, Germany), which was then vortexed (IKA VORTEX GENIUS 3, VG 3 S25) every 30 min for 30 s for a total of 6 times. The sample was placed at 4℃ overnight. Afterward, it underwent a 10-minute centrifugation at 12,000 rpm (Allegra 64R, Beckman, USA). Following the centrifugation, the supernatant was filtered using a 0.22-µm microporous membrane and then transferred to a vial for analysis using ultra-performance liquid chromatography-electrospray ionization-tandem mass spectrometry (UPLC-ESI-MS/MS) analysis.

UPLC-ESI-MS/MS analysis

Broadly targeted metabolomics employed AB sciex TripleTOF660 high resolution mass spectrometry for qualitative detection of mixed samples, followed by AB sciex6500 QTRAP quantification, combining the benefits of untargeted and targeted metabolomics to achieve precise qualification via high resolution mass spectrometry. Using UPLC-ESI-MS/MS for broadly targeted metabolomics, the sample extracts were analysed [20]. Chromatographic separation was performed on a Nexera X2 UPLC machine (Shimadzu, Kyoto, Japan) utilising an Agilent SB-C18 column (1.8 μm, 2.1 mm, 100 mm; Agilent, Foster City, CA, USA). The mobile phase consisted of acetonitrile combined with 0.1% formic acid (solvent B) and purified water containing 0.1% formic acid (solvent A). 0–9 min (95%→5% A), 9–10 min (5% A), 10–11.1 min (5%→95% A), and 11.1–14 min (95% A) comprised the gradient elution program. The injection volume was 2 µL, the flow rate was 0.35 mL/min, and the temperature setting for the column oven was 40 °C.

The alternative technique for injecting the effluent was a triple quadrupole-linear ion trap mass spectrometer (6500 Q-TRAP UPLC-MS/MS; Applied Biosystems, Framingham, MA, USA) equipped with an ESI Turbo Ion-Spray interface. The ESI ion source was operated with the following parameters: turbo spray ion source; 550 °C temperature; 5,500 V spray voltage (positive ion mode)/-4,500 V spray voltage (negative ion mode); gas I and gas II pressures of 50 and 60 psi respectively; curtain gas pressure of 25.0 psi; and collision-activated dissociation set to high. Triple quadrupole and linear ion trap modes were utilised for the calibration of the instrument and mass, using solutions of 10 and 100 µmol/L polypropylene glycol, respectively. Triple quadrupole pictures were obtained using the multiple reaction monitoring (MRM) assays with the collision gas (nitrogen) set to medium. For each MRM transition, the collision energy and delustering potential were tuned to maximise efficiency. Based on the metabolites that eluted at that time, a predetermined set of MRM transitions was looked for throughout each session. To gather and examine mass spectra, AB Sciex, Foster City, CA, USA, provided Analyst v1.6.3.

The metabolites in the project samples were subjected to qualitative and quantitative mass spectrometry analysis utilizing the Metware Biotechnology Co., Ltd. (MWDB), KEGG compound database, and MRM. Metabolite identification was conducted utilizing retention time (RT), MS2 fragmentation, metabolite precise mass, and MS2 fragment isotopic distribution. The secondary spectrum and RT of metabolites in the project samples are intelligently matched one by one with the secondary spectrum and RT of the company’s database through the company’s self-developed intelligent secondary spectrum matching method with the MS and MS2 tolerances set to 20 ppm and 20 ppm, respectively, and the RT tolerance of 0.2 min.

To ensure the precision of the metabolite annotations, interference signals were eliminated at the outset of each study. These signals included the repeated signals of K+, Na+, and NH4+ ions, fragment ions, and the isotope signal. We used Multi Quant v3.0.2 (AB Sciex) for MRM-based measurement of metabolites. In MRM mode, interference is initially eliminated by screening the precursor ions (parent ions) of the target substances through the four-stage rod and excluding the ions corresponding to other molecular weight substances. The collision chamber induces the precursor ions to ionize and break, resulting in the formation of numerous fragment ions. These fragment ions are subsequently filtered through the triple four-stage rod in order to select a characteristic fragment ion necessary for excluding the interference of non-target substances. Once the metabolite spectra of various samples were acquired, the peak area of each peak in the mass spectra was integrated. Subsequently, the mass spectra of identical metabolites in different samples were adjusted to account for the integration of the peaks [21,22,23].

Network pharmacology

Using the PubChem database (https://pubchem.ncbi.nlm.nih.gov/), the canonical SMILES of the prominent differential metabolites were searched. The relevant target proteins were then predicted using the Swiss Target Prediction database (http://www.swisstargetprediction.ch/). Once the pertinent target proteins were sorted out and the repeated targets were removed, the targets of the dominating differential metabolites were found [24].

A network was constructed and visualized with default parameters and similarity thresholds using STRING (http://string-db.Org, high confidence > 0.9) and Cytoscape v3.9.1 [25, 26]. The Kyoto Encyclopaedia of Genes and Genomes (KEGG, https://www.genome.jp/kegg/), Gene Ontology (GO, http://geneontology.org/) functional enrichment, and disease enrichment analyses were carried out using the David database (https://david.ncifcrf.gov/) (P < 0.01).

Molecular docking

The PPI results were imported into Cytoscape 3.9.1 software in TSV text format, and the topological attributes of the results were analyzed by using the “Network Analyzer” function to obtain four topological parameters (Degree centrality, Betweenness centrality, Closeness centrality, Eigenvector centrality) [27].

The active ingredients screened from CJ were validated through molecular docking with core targets using Autodock software. The 3D structures of the ligand active ingredients were obtained in mol 2 format from Zinc (http://zinc.docking.org/substances/home/) and PubChem databases, while the three-dimensional structures of target receptor proteins were retrieved from the protein database (RCSB) (http://www.rcsb.org/), then PDB format files were downloaded. Subsequently, water molecules and excess inactive ligands were removed from the proteins using PyMOL 2.4.1 software. The processed results were imported into AutoDock Tools 1.5.6 software for hydrogenation, charge addition, and output in pdbqt format [28].

The receptor and ligand pdbqt formats were imported into AutoDock to define the molecular docking range. By setting the target protein as the center of the grid, the center coordinates (center 3.48/2.398/11.641) and box size (size 66/52/62) parameters were adjusted to ensure complete coverage of the protein by the docking box.

In AutoDock, molecular docking involves identifying protein macromolecules, inserting small drug molecules, and setting up operational methods and docking parameters. The pdbqt format is utilized for calculating the minimum binding energy and visualizing the results with Origin software. OpenBabelGUI is used to convert the combined pdbqt format to PDB format, and the docking results are visualized using PyMOL 2.4.1 software [29].

Statistical analysis

The raw data was arranged using Excel 365 for Microsoft Office. Principal component analysis, Pearson correlation analysis, hierarchical cluster analysis, K-means cluster analysis, orthogonal partial least squares-discriminant analysis (OPLS-DA). Metabolites with variable importance in projection (VIP, OPLS-DA) and |log2fold change|>1 were considered differential metabolites (DMs) [30].

Results and analysis

Widely targeted metabolomics analysis

The differences in metabolite accumulation between Campanumoea javanica bl. And Codonopsis Radix

To learn more about the variations in metabolite accumulation between CJ and CR, we analyzed the metabolites using UPLC-ESI-MS/MS. In CJ and CR, about 1,116 metabolites were discovered and identified, categorized into ten distinct chemical types (Fig. 2A). The number of metabolites was 1,102 (CPM) = 1,102 (CP) > 1,093 (CT) > 1,076 (CJ). Phenolic acids were the most prevalent class (214–216). The remaining classes were lipids (145–150), flavonoids (104–119), amino acids and their derivatives (117–119), organic acids (110–113), alkaloids (111–112), derivatives (75), terpenoids (10–13), lignans and coumarins (55–60), others (128–131).

A Venn diagram was applied to analyze the differences in metabolite classes between CJ and CR aiming better to understand the variations in metabolites (Fig. 2B). The results showed that 1038 metabolites were shared between CJ and three CRs, including phenolic acids, lipids, amino acids and derivatives, alkaloids, organic acids, flavonoids, nucleotides and derivatives, lignans, coumarins, and terpenoids. The numbers of metabolites in each category were 204, 146, 116, 124, 108, 105, 98, 75, 52, and 10, respectively. Compared with CJ, 17 or 26 more compounds were identified in CPM, CT, and CP, mainly flavonoids and phenolic acids.

CJ shared 1057, 1063, and 1063 metabolites with CP, CT, and CPM, respectively (mean 1061), and the average similarity of metabolic profiles was 95.07%, and no unique compounds were found in CJ. The shared metabolites among three CRs were 1082–1090 (mean 1085), and the average similarity of metabolic profiles was 97.19%. The results indicated that CJ and the three CRs had remarkably similar metabolic profiles, thereby providing material basis supporting the potential use of CJ as a substitute for CR (Fig. 2C).

Fig. 2
figure 2

Analysis of differences in metabolite accumulation between Codonopsis Radix and Campanumoea javanica Bl. (A) Type distribution of metabolites; (B) Venn diagram of metabolites; (C) The quantity of metabolites shared by different samples

We performed principal component analysis, Pearson correlation analysis, and hierarchical cluster analysis to further assess fundamental characteristics and differences of metabolites between CJ and CR. According to the hierarchical cluster analysis results, CP, CPM, and CT metabolites were grouped in one category, whereas CJ metabolites were classified separately. It suggests that CP, CPM, and CT metabolites exhibit distinct accumulation characteristics compared to CJ (Fig. 3A). Pearson correlation analysis showed that the correlation coefficient of samples in the group was greater than that of species. So, the reliability of the test data was high, and the results of subsequent analyses were credible (Fig. 3B). Dim1 and Dim2 accounted for 28.6% and 17.9% of the overall variance, respectively, according to principal component analysis. Dim1 significantly separated CJ and CT (CP, CPM), and Dim2 significantly separated CT from the other three groups (Fig. 3C). Thus, the results of the analysis indicated that the metabolic profiles of CJ are distinct from those of CP, CT, and CPM, with some differences in content profiles.

Fig. 3
figure 3

Multivariate statistical analysis of different metabolites. (A) Hierarchical cluster analysis of Codonopsis Radix and Campanumoea javanica Bl. metabolites; (B) Pearson correlation analysis of Codonopsis Radix and Campanumoea javanica Bl. metabolites; (C) Principal component analysis of Codonopsis Radix and Campanumoea javanica Bl. metabolites

K-means cluster analysis was applied to a better understanding of the relative content trend of metabolites in distinct species. Approximately 1,116 metabolites were divided into four subclasses (Fig. 4). In sub-class 1, 218 compounds (19.53%) were relatively high in CP, mainly Lipids (17.43%), Phenolic acids (15.60%), and Flavonoids (12.84%). In sub-class 2, 251 compounds (22.49%) were relatively high in CJ, the majority of which were phenolic acids (26.29%), lipids (19.12%), and flavonoids (13.15%). In sub-class 3, 277 compounds (24.82%) were relatively high in CT, mainly phenolic acids (20.94%), amino acids and derivatives (14.08%), and alkaloids (13.36%). In sub-class 4, 370 compounds (33.15%) were relatively high in CPM, the majority of which were phenolic acids (17.13%), organic acids (16.57%), and saccharides (15.17%). Metabolic pathway annotation was performed by the KEGG database for 251 compounds enriched in CJ, and a total of 67 metabolic pathways were enriched, which were mainly enriched in Metabolic pathways (ko01100), Biosynthesis of secondary metabolites (ko01110), ABC Transporters (ko02010), and Biosynthesis of amino acids (ko01230).

Fig. 4
figure 4

K-means clustering analysis of the patterns of 1,116 differential metabolites in Codonopsis Radix and Campanumoea javanica Bl

Screening for differential metabolites

To screen for differential metabolites, we modeled the comparison of CJ with CPM, CT, and CP by OPLS-DA, and the Q2 values were all greater than 0.90, which is an excellent model (Fig. 5A and B, and 5C). The screening of differential metabolites were conducted using criteria of | log2(fold change) |>1 and VIP ≥ 1. CJ was compared with CP, CPM, and CT in a two-by-two comparison to obtain the following three comparison groups of differential metabolites: CP vs. CJ, CPM vs. CJ, and CT vs. CJ. The specific expressions of DMs are illustrated in Fig. 6.

Fig. 5
figure 5

OPLS-DA score chart. A: CP vs. CJ; B: CT vs. CJ; C: CPM vs. CJ

Compared with CP, CT, and CPM, 114,112 and 84 DMs were up-regulated in CJ (Fig. 6A), but 314, 262 and 356 DMs were down-regulated in CJ (Fig. 6B). We discovered that there were 35, 11, and 41 DMs up-regulated only in CP vs. CJ, CPM vs. CJ, and CT vs. CJ, respectively, and 32, 72, and 56 DMs were down-regulated, respectively. 45 DMs were all up-regulated in CP vs. CJ, CPM vs. CJ, and CT vs. CJ, including phenolic acids (20 DMs), flavonoids (10 DMs), alkaloids (4 DMs), others (3 DMs), amino acids and derivatives (2 DMs), lipids (2 DMs), terpenoids (2 DMs), lignans and coumarins (1 DM), and organic acids (1 DM).

Fig. 6
figure 6

The specific expressions of DMs (A) Up-regulated differential metabolites in comparison groups; (B) Down-regulated differential metabolites in comparison groups

Screening for dominant differential metabolites

Based on analysis of metabolite difference, the dominant differential metabolites in CJ were identified by screening for up-regulated DMs and those with a significance level (P < 0.05) in CJ. Total of 17 dominant differential metabolites were found in CJ (Fig. 7). Among them, 12 components’ Canonical SMILES were retrieved in the PubChem database, including 1-caffeoylquinic acid*, 1,3-O-dicaffeoylglycerol, tachioside, 4-O-caffeoylquinic acid*, radicamine A, stearamide, 10- formyltetrahydrofolic acid, oleamide (9-octadecenamide), 2,3-dimethylmaleic anhydride, isololiolide, sec-O-glucosylhamaudol, and cyclohexane-1,3-dione.

Fig. 7
figure 7

Hierarchical clustering heat map of Campanumoea javanica Bl. dominant differential metabolites

Network pharmacological analysis of CJ dominant differential metabolites

The corresponding target proteins for the twelve dominant differential metabolites in CJ were estimated, leading to the identification of 229 targets. The protein-protein interaction (PPI) network was constructed with 229 targets using Cytoscape 3.9.1 (Fig. 8A). The targets with ≥ 17 degrees were PIK3R1, PIK3CA, PIK3CD, ESR1, HSP90AA1, AKT1, EGFR, suggesting a more crucial function for the targets in the PPI network map.

On 229 targets, GO functional enrichment analysis, KEGG pathway enrichment analysis, and disease enrichment analysis were performed utilizing the David database (6.8, https://david.ncifcrf.gov/). The statistical significance threshold of P < 0.01 was employed for the GO functional analysis, KEGG pathway enrichment analysis, and disease enrichment analysis. Approximately 873 items in total, including 589 biological processes (BP), 97 cell compositions (CC), 187 molecular functions (MF), 147 KEGG pathway analysis entries, and six disease enrichment analysis entries, were derived from the GO functional analysis. Visualization of the top 20 entries for each category was performed via the bioinformatics website. (Figure 8C and D, http://www.bioinformatics.com.cn/). The composition of BP primarily consists of xenobiotic stimuli, a one-carbon metabolic process, and positive regulation of MAP kinase activity. CC primarily included plasma membrane, cytosol, and membrane raft. Enzyme binding, endopeptidase activity, and carbonate dehydratase activity were MF’s primary inclusions in the GO functional analysis. Nitrogen metabolism, cancer (prostate cancer, proteoglycans in cancer, bladder cancer), and signaling pathway (HIF-1 signaling pathway, Fc epsilon RI signaling pathway, thyroid hormone signaling pathway) were among the major topics covered by KEGG pathway analysis. Using Cytoscape software (3.7.0), a “drug-component-target-network” map of CJ was generated based on the evaluation of 12 components, 229 targets, and 20 pathways (Fig. 8B).

Fig. 8
figure 8

(A) protein-protein interaction network diagram of Campanumoea javanica Bl.; (B) Campanumoea javanica Bl. dominant chemical composition-target-pathways network; (C) KEGG analysis of signaling pathways enriched by candidate genes; (D) GO analysis of candidate targets, including molecular function, cellular component, and biological processes

Molecular docking of core targets and dominant differential metabolites

Based on the criteria of top 5 Degree values and topological parameters greater than the median (HSP90AA1, AKT1, EGFR had the Same Degree value and 4 topological parameters of PIK3CD were below the median), 6 target proteins were finally selected (PIK3R1, PIK3CA, ESR1, AKT1, HSP90AA1, EGFR) and molecularly docked with 12 dominant metabolites. 72 receptor-ligand docking results were finally obtained, as shown in Fig. 9.

Fig. 9
figure 9

Molecular docking heat map of the dominant differential metabolites of Campanumoea javanica Bl. with key targets

The binding energies of the 72 receptor-ligand pairs were all ≤ -5.0 kcal/mol, with 7 pairs having energies ≤ -9.0 kcal/mol, accounting for 9.7%. This indicates a strong binding affinity between the predicted active ingredients and key targets in this study. Finally, the top 7 docking results based on molecular docking scores were visualized using PyMol 2.4.1 software, as shown in Fig. 10.

Fig. 10
figure 10

Visualization diagram of molecular docking. (A: PIK3CA ~ sec-O-glucosylhamaudol; B: AKT1 ~ sec-O-glucosylhamaudol; C: HSP90AA1 ~ sec-O-glucosylhamaudol; D: HSP90AA1 ~ 1,3-O-dicaffeoylglycerol; E: AKT1 ~ 10-formyltetrahydrofolic acid; F: AKT1 ~ 1,3-O-dicaffeoylglycerol; G: AKT1 ~ 4-O-caffeoylquinic acid*)

Discussion

Differences between Campanumoea javanica bl. and Codonopsis Radix

The metabolite analysis results showed that 1038 (93.01%) shared metabolites were identified in CJ, CPM, CP, and CT, and 251 compounds (22.49%) in CJ had content dominance, indicating that CJ and CR metabolite species were highly similar, but with significant differences in content. In addition, CPM, CP, and CT had 1073 (96.15%) shared metabolites, and the cluster analysis clustered them into one group, indicating their metabolite profiles were almost identical. This experiment greatly enhanced the metabolic database of CR compared to Tang’s study [31]. The above evidence confirmed the correctness of the Chinese Pharmacopoeia in treating CPM, CP, and CT all as sources of CR. It was consistent with the “similar species, similar effects” hypothesis advanced by Xie [32]. Subsequently, on the basis of this experiment, in-depth research can be conducted to verify the scientific basis of the folk use of CJ as a CR substitute. There are many successful cases have been obtained by applying the same family or genus of plants to discover new drug sources. For example, reserpine, an active ingredient in the treatment of hypertension, was originally only obtained from the Rauvolfia serpentine. In recent years, researchers have isolated and identified it in other plant from the same genus [33].

Great potential for food exploitation of Campanumoea javanica bl

Our research found that CJ was rich in terpenoids, flavonoids, and lipids, with the relative content of 251 compounds including 66 phenolic acids, 48 lipids, and 33 flavonoids was higher than that of CR. Commonly employed as biological preservatives in food, phenolic acids function primarily to prevent oxidative damage by acting as antioxidants. Additionally, they predominantly contribute to the enhancement of overall health through their antimicrobial and anti-inflammatory properties [34]. Terpenoids are the most prevalent class of chemicals in natural products and extensively dispersed in nature, as a significant class of secondary metabolites in plants. They are widely used as raw materials in pharmaceuticals, food, and cosmetics, as well as taking part in physiological processes like plant growth and development and environmental response. In the future, they may also be used as natural additive, such as essence, sweeteners, nutrition enhancers [35]. Flavonoids are necessary nutrients that are crucial to health care but cannot be produced by the human body and must instead be consumed from outside sources. Daily consumption of fruits and vegetables high in flavonoids can significantly lower the risk of developing cancer. Foods that are high in these CJ components will be developed into healthy congee, savory dishes, cold dishes, etc., not only to meet the needs of the public’s healthcare and broaden the range of uses for CJ but also to increase the income of the farmers in the producing areas [36, 37]. One of the nutrients in a living thing is called lipid, which also acts as a carbon and energy reserve for the body as well as a mediator of cellular signaling pathways and a stress regulator, among other functions. Lipids are also an essential part of the biofilm. The flavor and nutritional value of foods that are lipid-rich are directly correlated with their lipid content and composition [38]. To conduct application and processing research on these CJ plant resources, CJ can thus draw on the existing application forms and processing methods of CR. This helps to increase the success rate of new product development for CJ resources, lower the risk of product development, and thus promote the in-depth development and application of CJ plant resources.

Campanumoea javanica bl. Has the potential for new drug exploitation

Polyethylenimine can be modified by the dominant metabolite 2,3-dimethylmaleic anhydride from CJ to improve the antitumor effects of conventional therapy [39], whereas isooliolide changes the expression of proteins that are crucial for the apoptotic cascade and have antitumor properties [40]. The 12 major dominant differential metabolite components of CJ were also predicted as targets using the network pharmacology method, and the results revealed that they primarily involved 229 targets, including PIK3CA, PIK3R1, HSP90AA1, EGFR, MAPK1, etc., and the targets were enriched in 147 pathways, including Nitrogen metabolism, Prostate cancer, HIF-1 signaling pathway, and proteoglycans in cancer. Sun et al. discovered that PIK3CA is responsible for encoding the catalytic subunit p110a protein and that its release and activation can keep PI3K in a continuously active state. Additionally, changes in PIK3CA can induce other downstream PI3K kinases, which can activate the entire pathway and produce a pro-inflammatory effect [41]. Wu et al. discovered that PIK3R1 was a negative regulator of PI3K/Akt/mTOR signaling and that lipopolysaccharide can downregulate its expression and contribute to the progression of inflammation [42]. Wang et al. discovered that HSP90AA1 is a related gene that promotes the progression of squamous carcinoma in the lung [43]. Jin et al. discovered that EGFR has rapid endocytosis activity and is expressed in numerous solid tumors, and Li et al. discovered that prolactin targets are MAPK1 and EGFR, etc. [44, 45]. As a result, CJ is thought to perform prominently in treating inflammation, cancer, breastfeeding. Based on network pharmacology and molecular docking, computer simulations were used to predict and validate the potential mechanism of action of the major metabolites of CJ. Subsequent ex vivo and in vivo efficacy experiments are needed to validate and provide theoretical support for the clinical application of CJ.

Conclusions

In this study, we employed metabolomics, network pharmacology, and molecular docking verification to delve into the metabolites of CJ and CR and predict the pharmacodynamics of the dominant differential metabolites of CJ. The results of the metabolomic analysis showed an average of 1061 (95.07%) shared metabolites in CJ and CRs but with remarkably different metabolite profiles. CJ had a higher relative concentration of 251 components, including 66 phenolic acids, 48 lipids, and 33 flavonoids compared to CRs. These components indicated that CJ possessed good edible characteristics and food development value. Furthermore, the dominant differential metabolites of CJ have shown medical utility in treating cancer, inflammation, and lactation, indicating its value for new drug development. However, further exploration is required to understand the specific mechanism of the different components and pharmacological effects due to the lack of consideration of the relative content differences and interactions between the components. This research has provided initial insights into the metabolite differentiation between CJ and CR and the dominant metabolites from the perspective of plant metabolomics. It also provides a conceptual framework for CJ’s advancement of novel pharmaceuticals and significantly improves the metabolite databases on CJ and CR.

Data availability

The data supporting the fndings of this study are available from the corresponding author, W.F. (wufaming@zmu.edu.cn), upon request.

Abbreviations

CJ:

Campanumoea javanica Bl

CPM:

Codonopsis pilosula var. modesta (Nannf.) L.T. Shen

CP:

Codonopsis pilosula (Franch.) Nannf

CT:

Codonopsis tangshen Oliv

CR:

Codonopsis Radix

RT:

Retention time

GO:

Gene Ontology

BP:

Biological processes

CC:

Cell compositions

MF:

Molecular functions

ESI:

Electrospray ionization

PPI:

Protein-protein interaction

MRM:

Multiple reaction monitoring

DMs:

Differential metabolites

VIP:

Variable importance in projection

MWDB:

Metware database

KEGG:

Kyoto Encyclopedia of Genes and Genomes

OPLS-DA:

Orthogonal partial least squares-discriminant analysis

UPLC-ESI-MS/MS:

Ultra-performance liquid chromatography-electrospray ionization-tandem mass spectrometry

References

  1. Wu QR. An Illustrated Book of plants. Shanghai: Commercial; 1957.

    Google Scholar 

  2. Pharmacopoeia Committee of the Ministry. of Health of the People’s Republic of China. Pharmacopoeia of the people’s republic of China. People’s health publishing house; 1977.

  3. Guizhou Provincial Drug Administration. Quality standards of Chinese medicinal materials and ethnic medicinal materials in Guizhou Province. Guiyang: Guizhou Science and Technology; 2003.

    Google Scholar 

  4. Guizhou Food and Drug Administration. Processing specification of Chinese herbal medicine decoction pieces in Guizhou Province. Guiyang: Guizhou Science and Technology Publishing House; 2005.

    Google Scholar 

  5. Guangxi Zhuang Autonomous Region Food and Drug Administration. Processing specification of traditional Chinese medicine decoction pieces in Guangxi Zhuang Autonomous. Guangxi: Region Guangxi Science and Technology; 2007.

    Google Scholar 

  6. Tu CL, Liao SW. Advances in Chemical constituents and pharmacological action from roots of Campanumoea Javanica. Dru Eva. 2022;19(16):1022–4. https://doi.org/10.19939/j.cnki.1672-2809.2022.16.18.

    Article  Google Scholar 

  7. Tian J, Xu YL. Thin layer chromatography was used to simultaneously detect Campanumoea javanica bl. In Qinghuo Yangyuan Capsules and Jinwu Gutong Capsules. Guizhou Med J. 2011;35(9):836–7. https://doi.org/10.3969/j.issn.1000-744X.

    Article  Google Scholar 

  8. Jia MR, Li XW. Summary of Chinese ethnic Medicine annals. Beijing: China pharmaceutical science and technology; 2005.

    Google Scholar 

  9. Cheng TY, Liu YY, Hu HB. Study on specific Chromatogram and multi-components determination of alkaloids in Hulisan tablets by HPLC. Chin J Mod App Pha. 2021;38:1594–600. https://doi.org/10.13748/j.cnki.issn1007-7693.2021.13.010.

    Article  CAS  Google Scholar 

  10. Yunnan Provincial People’s Government Portal. http://www.hh.gov.cn/zmhd/jytabljg/2018njytabljg/2018nrdjybljg/201809/t20180926_302038.html. Accessed 20 Jan 2024.

  11. Zhou MM, Liu L, Zhou ZH. Reviews of Campanumoea javanica bl. In different ethnic medicine in clinical application. Chin J Ethn Ethn. 2018;27:37–41. https://doi.org/10.3969/j.issn.1007-8517.2018.24.zgmzmjyyzz201824012.

    Article  Google Scholar 

  12. Peng M, Zhang ZD, Yang J. Inhibitory effects of 14 polysaccharides on Sarcocarcinoma S180 of mice. J Mou Agr Bio. 2011;30(1):56–9. https://doi.org/10.15958/j.cnki.sdnyswxb.2011.01.006.

    Article  Google Scholar 

  13. Editorial Committee of Flora of China. Chinese Academy of Sciences. Beijing: Flora of China science; 1993.

    Google Scholar 

  14. Zhang ZJ, Yang XS, Zhu WS. Campanumoea javanica bl. Study on chemical constituents. Chin Traditional Herb Drugs. 2005:28–30.

  15. Yang DS, Li ZL, Wang X, Yang YP, Peng WB, Liu KC, et al. Chemical constituents from roots of Campanumoea Javanica and their antiangiogeneic activities. Chin Tradit Herb Drugs. 2015;46:470–5. https://doi.org/10.7501/j.issn.0253-2670.2015.04.002.

    Article  CAS  Google Scholar 

  16. Wu LF, Jing YS, Zhang ZD, Yang J. A comparative study of different extraction methods for polysaccharides from Campanumoea Javanica roots. Foo Sci. 2012;33:45–8.

    Google Scholar 

  17. Shen Y. Bio-Pharmacy and Quality Standard Research of Radix. Yunnan university of traditional Chinese medicine. 2013:79.

  18. Pan L, Zhou LG, Deng Y, Wang X, Zhou GP, Yuan MM. Chemical constituents of Ethyl acetate Extract Part Campanumoea Javanica subsp. Javanica. Chin Med Mat. 2018;41:1119–21. https://doi.org/10.13863/j.issn1001-4454.2018.05.022.

    Article  Google Scholar 

  19. Zhang KX, Zhang DL, Yang QF, Long LT, Xie J, Wang Y, et al. Integrated widely targeted metabolomics and network pharmacology revealed quality disparities between Guizhou and conventional producing areas of Codonopsis Radix. Front Nutr. 2023;10:1271817. https://doi.org/10.3389/fnut.2023.1271817.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Zhao YS, Sun CY, Wang SZ, Zhang ML, Li YL, Xue QH, et al. Widely targeted metabolomic, transcriptomic, and metagenomic profiling reveal microbe-plant-metabolic reprogramming patterns mediated by Streptomyces pactum Act12 enhance the fruit quality of Capsicum annuum L. Food Res Int. 2023;166:112587. https://doi.org/10.1016/j.foodres.2023.112587.

    Article  CAS  PubMed  Google Scholar 

  21. Huang X, Wang W, Gong T, Wickell D, Kuo LY, Zhang X, et al. The flying spider-monkey tree fern genome provides insights into fern evolution and arborescence. Nat Plants. 2022;8(5):500–12. https://doi.org/10.1016/j.msec.2019.04.081.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Ning K, Hou C, Wei X, Zhou Y, Zhang S, Chen Y, et al. Metabolomics Analysis Revealed the Characteristic Metabolites of Hemp Seeds Varieties and metabolites responsible for antioxidant Properties. Front Plant sci. 2022;13:904163. https://doi.org/10.3389/fpls.2022.904163.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Wang Y, Zhang H, Ri HC, An Z, Wang X, Zhou JN, et al. Deletion and tandem duplications of biosynthetic genes drive the diversity of triterpenoids in Aralia elata. Nat Commun. 2022;13:2224. https://doi.org/10.1038/s41467-022-29908-y.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Song XH, He Y, Liu M, Yang Y, Yuan Y, Yan J, et al. Mechanism underlying Polygonum capitatum effect on Helicobacter pylori-associated gastritis based on network pharmacology. Bioorg Chem. 2021;114:105044. https://doi.org/10.1016/j.bioorg.2021.105044.

    Article  CAS  PubMed  Google Scholar 

  25. Qu J, Ke F, Liu Z, Yang X, Li X, Xu H, et al. Uncovering the mechanisms of dandelion against triple-negative breast cancer using a combined network pharmacology, molecular pharmacology, and metabolomics approach. Phytomedicine. 2022;99:153986. https://doi.org/10.1016/j.phymed.2022.153986.

    Article  CAS  PubMed  Google Scholar 

  26. Sherman BT, Hao M, Qiu J, Jiao X, Baseler MW, Lane HC, et al. DAVID: a web server for functional enrichment analysis and functional annotation of gene lists (2021 update). NAR. 2022;50(W1):W216–21. https://doi.org/10.1093/nar/gkac194.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Fan ZX, Li XN, Jin M, Wang ZL, Gao ZX. Mechanism of Mongolian medicine Sendeng-9 decoction in the treatment of brucellosis based on network pharmacology based on network pharmacology. J Baotou Med Coll. 2024;40(05):72–8. https://doi.org/10.16833/j.cnki.jbmc.2024.05.014.

    Article  Google Scholar 

  28. Shu YX, Guo JT, Xu WF, Zhang YP, Zhang KF, Shi XY, et al. Action mechanism of periploca forrestii in against alzheimer’s disease based on network pharmacology and molecular docking. Chin J Gerontol. 2023;43(22):5484–92. https://doi.org/10.3969/j.issn.1005-9202.2023.22.031.

    Article  CAS  Google Scholar 

  29. Wang Y, Zhang YW, Wang YJ, Shu XY, Lu CR, Shao SL, et al. Using Network Pharmacology and Molecular Docking to explore the mechanism of Shan Ci Gu (Cremastra appendiculata) against Non-small Cell Lung Cancer. Front Chem. 2021;9:682862. https://doi.org/10.3389/fchem.2021.682862.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Jiao YL. Transcriptomic and metabolomic analyses reveal the flavor of bitterness in the tip shoots of Bambusa oldhamii Munro. Sci Rep. 2023;13(1):14853. https://doi.org/10.1038/s41598-023-40918-8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Tang XY, Zeng JX, Wang XX, Xu WY, Zhao PC, Fan CL, et al. Chemical and metabolic profiling of Codonopsis Radix extract with an integrated strategy using ultra-high-performance liquid chromatography coupled with mass spectrometry. J Sep Sci. 2023;46(2):e2200723. https://doi.org/10.1002/jssc.202200723.

    Article  CAS  PubMed  Google Scholar 

  32. Xie ZW. Study on the theory of traditional Chinese medicine varieties. China J Chin Materia Med. 1991:195–9.

  33. Chen XY, Guo XY, Li C, Li Y, Yao SC, Zhang ZJ, et al. Research progress on pharmacological properties and monoterpenoid indole alkaloids biosynthetic pathway in plants of genus Rauvolfia. Chin Tradit Herb Drugs. 2019;50:2004–12. https://doi.org/10.7501/j.issn.0253-2670.2019.08.035.

    Article  Google Scholar 

  34. Luo Y, Li YC, Wang M, Zhou L, Meng FB, Jiang LS. Effects of grafting methods and raw materials on the physicochemical properties and biological activities of phenolic acids grafted oat β-glucan. Food Res Int. 2023;173(Pt 1):113250. https://doi.org/10.1016/j.foodres.2023.113250.

    Article  CAS  PubMed  Google Scholar 

  35. Zhang MZ, Wang LN, Zhao W, Gao D, Zhang YX, Yin KD. Effects of salt/alkali-tolerant phosphate-solubilizing bacteria on lipid metabolome of soybean seedlings. Chin J Oil Crop Sci. 2023:1–11. https://doi.org/10.19802/j.issn.1007-9084.2022274.

  36. Wang ZZ, Zhao Y, Pan WW, Zhao YR, Zhang HX, Wang RF. Study on Antitumor Activity of flavonoids from Zizyphus Jujuba Mill. J Chin Ins Foo Sci Tec. 2023;23:90–8. https://doi.org/10.16429/j.1009-7848.2023.03.010.

    Article  Google Scholar 

  37. Zheng P, Wen M, Liu QY, Wang X, Zuo YJ. Optimization of extraction process and evaluation of antioxidant and antitumor activities of total flavonoids from Scutellaria Barbata. Sci Technol food Ind. 2023;1–16. https://doi.org/10.13386/j.issn1002-0306.2023020199.

  38. Zhang JH, Liu WQ, Luo HM. Research progress on the activity of terpenoids in medicinal plants. Wor Sci Tec Sci Tec Tra Chin Med. 2018;20:419–30. https://doi.org/10.11842/wst.2018.03.015.

    Article  Google Scholar 

  39. Huang X, Cao J, Zhan Y, Liu T, Yan H. Polyethylenimine modified with 2,3-dimethylmaleic anhydride potentiates the antitumor efficacy of conventional chemotherapy. Mater Sci Eng C Mater Biol Appl. 2019;102:558–68. https://doi.org/10.1016/j.msec.2019.04.081.

    Article  CAS  PubMed  Google Scholar 

  40. Vizetto-Duarte C, Custodio L, Gangadhar KN, Lago JH, Dias C, Matos AM, et al. Isololiolide, a carotenoid metabolite isolated from the brown alga Cystoseira tamariscifolia, is cytotoxic and able to induce apoptosis in hepatocarcinoma cells through caspase-3 activation, decreased Bcl-2 levels, increased p53 expression and PARP cleavage. Phytomedicine. 2016;23:550–7. https://doi.org/10.1016/j.phymed.2016.02.008.

    Article  CAS  PubMed  Google Scholar 

  41. Sun CA, Wang FY, Li SX. The influence of Shao Yao Tang to signal pathway expression of TLR2、PIK3CA and AKT in colon tissues of Uc rat models. Jiangsu J Tradit Chin Med. 2019;51(11):83–6. https://doi.org/10.3969/j.issn.1672-397X.2019.11.030.

    Article  Google Scholar 

  42. Wu T, Zhang GQ, Li C, Zhang LT. Effect of miR-155 on the PI3K/AKt/mTOR signaling pathway in rheumatoid arthritis rats by targeting PIK3R1. Chin J Osteoporos. 2022;28:1593–8. https://doi.org/10.3969/j.issn.1006-7108.

    Article  Google Scholar 

  43. Wang L, Zhao H, Zhang L, Luo H, Chen Q, Zuo X. HSP90AA1, ADRB2, TBL1XR1, and HSPB1 are chronic obstructive pulmonary disease-related genes that facilitate squamous cell lung cancer progression. Oncol Lett. 2020;19:2115–22. https://doi.org/10.3892/ol.2020.11318.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Jin MQ, Shuang CZ, Shen YL. Development and anti-tumor study of bispecific antibody drug conjugate targeting TPBG and EGFR. Chin Bio. 2023;43(5):11–23. https://doi.org/10.13523/j.cb.2212010.

    Article  Google Scholar 

  45. Li MT, Shi L, Hou C, Yu M, Li JK. Molecular mechanism of galactopoietic pill in hypogalactia treatment based on network pharmacology and molecular docking. J Shaanxi nor U. 2023;51(3):75–85. https://doi.org/10.15983/j.cnki.jsnu.2023206.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This research was funded by the National Natural Science Foundation of China (82060687, 81760679), the Science and Technology Support Plan Project of Guizhou province (Qiankehe Foundation [2021] General 248), and the State Key Laboratory Open Subjects (Qiankehe Platform Talent [2017] 5101).

Author information

Authors and Affiliations

Authors

Contributions

P.J. experimented, analyzed the data and wrote the article. L.S. experimented, analyzed the data and wrote the article. W.X. participated in part of the experiments. L.S. revised the manuscript. X.J. refined the ideas. W.Y. participated in part of the experiments. Y.Q. revised the manuscript. W.F. designed the study and revised the manuscript. Z.D. designed the study and revised the manuscript.

Corresponding authors

Correspondence to Faming Wu or Delin Zhang.

Ethics declarations

Ethics approval and consent to participate

We declare that all the experimental plants were collected with permission from local authorities of agricultural department and all experimental studies on plants were complied with relevant institutional, national, and international guidelines and legislation. All samples are conserved at the School of Pharmacy, Zunyi Medical University.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Peng, J., Liu, S., Wu, X. et al. Metabolomics combined with network pharmacology reveals the potential development value of Campanumoea javanica Bl. and its metabolite differences with Codonopsis Radix. BMC Plant Biol 24, 683 (2024). https://doi.org/10.1186/s12870-024-05401-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12870-024-05401-0

Keywords