Next Article in Journal
Formulation of the Polysaccharide FucoPol into Novel Emulsified Creams with Improved Physicochemical Properties
Next Article in Special Issue
Effects of Sorafenib and Quercetin Alone or in Combination in Treating Hepatocellular Carcinoma: In Vitro and In Vivo Approaches
Previous Article in Journal
Bud-Poplar-Extract-Embedded Chitosan Films as Multifunctional Wound Healing Dressing
Previous Article in Special Issue
Recent Updates on the Functional Impact of Kahweol and Cafestol on Cancer
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Uncovering the Mechanisms of Active Components from Toad Venom against Hepatocellular Carcinoma Using Untargeted Metabolomics

1
National Traditional Chinese Medicine Clinical Research Base, Drug Research Center of Integrated Traditional Chinese and Western Medicine, The Affiliated Traditional Chinese Medicine Hospital of Southwest Medical University, Luzhou 646000, China
2
State Key Laboratories for Quality Research in Chinese Medicines, Faculty of Chinese Medicine, Macau University of Science and Technology, Macau 853, China
3
Institute of Integrated Chinese and Western Medicine, Southwest Medical University, Luzhou 646000, China
*
Authors to whom correspondence should be addressed.
Molecules 2022, 27(22), 7758; https://doi.org/10.3390/molecules27227758
Submission received: 5 October 2022 / Revised: 7 November 2022 / Accepted: 8 November 2022 / Published: 10 November 2022
(This article belongs to the Special Issue Natural Products as Antitumor Agents)

Abstract

:
Toad venom, a dried product of secretion from Bufo bufo gargarizans Cantor or Bufo melanostictus Schneider, has had the therapeutic effects of hepatocellular carcinoma confirmed. Bufalin and cinobufagin were considered as the two most representative antitumor active components in toad venom. However, the underlying mechanisms of this antitumor effect have not been fully implemented, especially the changes in endogenous small molecules after treatment. Therefore, this study was designed to explore the intrinsic mechanism on hepatocellular carcinoma after the cotreatment of bufalin and cinobufagin based on untargeted tumor metabolomics. Ultraperformance liquid chromatography with tandem mass spectrometry (UHPLC-MS/MS) was performed to identify the absorbed components of toad venom in rat plasma. In vitro experiments were determined to evaluate the therapeutic effects of bufalin and cinobufagin and screen the optimal ratio between them. An in vivo HepG2 tumor-bearing nude mice model was established, and a series of pharmacodynamic indicators were determined, including the body weight of mice, tumor volume, tumor weight, and histopathological examination of tumor. Further, the entire metabolic alterations in tumor after treating with bufalin and cinobufagin were also profiled by UHPLC-MS/MS. Twenty-seven active components from toad venom were absorbed in rat plasma. We found that the cotreatment of bufalin and cinobufagin exerted significant antitumor effects both in vitro and in vivo, which were reflected in inhibiting proliferation and inducing apoptosis of HepG2 cells and thereby causing cell necrosis. After cotherapy of bufalin and cinobufagin for twenty days, compared with the normal group, fifty-six endogenous metabolites were obviously changed on HepG2 tumor-bearing nude mice. Meanwhile, the abundance of α-linolenic acid and phenethylamine after the bufalin and cinobufagin intervention was significantly upregulated, which involved phenylalanine metabolism and α-linolenic acid metabolism. Furthermore, we noticed that amino acid metabolites were also altered in HepG2 tumor after drug intervention, such as norvaline and Leu-Ala. Taken together, the cotreatment of bufalin and cinobufagin has significant antitumor effects on HepG2 tumor-bearing nude mice. Our work demonstrated that the in-depth mechanism of antitumor activity was mainly through the regulation of phenylalanine metabolism and α-Linolenic acid metabolism.

Graphical Abstract

1. Introduction

Hepatocellular carcinoma (HCC), the sixth most common cancer and the third leading cause of cancer death in the world in 2020, remains a major challenge to global public health [1]. Surgical resection is still the primary treatment for HCC patients diagnosed at the early stage, but it is not ideal for HCC patients diagnosed at the advanced stage [2]. Currently, despite the use of sorafenib, which can prolong the survival rate, the concomitant problems also greatly affect the prognosis of patients, such as severe side effects, short-term effects, and drug resistance [3]. Traditional Chinese medicine (TCM) has attracted increasing interest as a potential antitumor drug repository due to its beneficial effects and relatively low toxicity [4].
Toad venom, also named as Chansu, is a dried product of secretion from Bufo bufo gargarizans Cantor or Bufo melanostictus Schneider, which was firstly recorded in Yao Xing Lun, dated back to the Tang Dynasty, for the treatment of brain dystrophy [5]. According to the Chinese Pharmacopoeia (2020 edition), toad venom possesses a series of pharmacological activities on detoxification and analgesia and has a cardiotonic effect, which is mainly applied for the treatment of ulcer, sore throat, sunstroke, fever, abdominal pain, vomiting, and diarrhea. With the in-depth study of toad venom, it could be considered to be a promising anticancer candidate due to its powerful activity against various cancer cells [2]. Previous phytochemical research revealed that bufadienolides are the main active component of toad venom, such as bufalin and cinobufagin [6]. Bufadienolides, acting as Na+/K+-ATPase inhibitors, can trigger intracellular Na+/Ca2+ exchange, increase intracellular calcium levels, and subsequently lead to an increase in cardiac contractility. However, in recent years, bufadienolides have received increasing attention for their antitumor activity. It has been demonstrated that bufalin (Buf) can inhibit HepG2 cell proliferation, invasion, migration, and adhesion, which have a definite effect on HCC [7]. Moreover, the HepG2 cell apoptosis induced by Buf was associated with an increase in Fas, Bax, and Bid expression and a decrease in Bcl-2 expression [8]. Cinobufagin (Cino), another major component of toad venom, has been used in clinical treatment of heart failure, sores, and pain [9]. Previously, we have demonstrated that Cino could effectively inhibit the growth of breast cancer in vitro and in vivo [10]. In addition, Cino exerted effective therapeutic activity against various cancers, such as liver cancer, lung cancer, melanoma, and colon cancer [9,11,12]. Notably, our research has also found that both Buf and Cino can effectively interact with HepG2 cells and be detected intracellularly, providing the basis of its antitumor activity [13]. Jinghui Zhang et al. demonstrated that the combined action of Buf and Cino acted synergistically in inducing apoptosis and inhibiting growth on HepG2 tumor-bearing nude mice [14]. Therefore, it is of great value to study these two active components of toad venom and their pharmacological activities. So far, the underlying mechanisms of toad venom for treating HCC have not been fully implemented, especially the changes in endogenous small molecules after treatment.
Metabolomics provides an effective approach to identify endogenous metabolite alterations and elucidate potential therapeutic mechanisms so as to reflect the overall functional state of the body. This is consistent with the integrity of TCM theory and presents the whole state of human body as under the influence of multiple factors [15]. Currently, metabolomics has a great potential for identifying and evaluating biomarkers and therapeutic targets of various cancers in clinical practice, such as HCC, lung cancer, colorectal cancer, prostate cancer, ovarian cancer, and other diseases [16]. In the present study, we identified the major components of toad venom in rat plasma. Actually, the compounds absorbed into plasma were considered biologically active and deserved further investigation [17]. On the basis of our previous study, we focused specifically on elaborating the antitumor activity and the possible mechanisms of Buf and Cino on HepG2 tumor-bearing nude mice. Concretely, nontargeted tumor metabolomics analysis was performed on ultraperformance liquid chromatography with tandem mass spectrometry (UHPLC-MS/MS) before and after cotreatment of Buf and Cino. Finally, this study will hopefully provide a novel insight to the potential mechanism of toad venom against HCC.

2. Results

2.1. In Vivo Compound Identification of Toad Venom in Rat Plasma

The major compounds of toad venom in rat plasma were detected using the UHPLC-HR-Q-Exactive Orbitrap MS system. The compound identification was performed based on chromatographic peaks, retention time, potential molecular formula, fragment ion information, and reference data. Previously, a total of 93 compounds in toad venom were characterized in the positive ion mode of the UHPLC-MS/MS system [13]. In this study, 27 compounds of toad venom extract in rat plasma were identified, including 19 prototype compounds and 8 metabolites. The detected prototype compounds included 16 bufogenins, 2 amino acids, and 1 alkaloid. Total ion chromatograms in the positive ion mode of blank plasma, toad venom plasma, and toad venom exact are shown in Figure 1. The details of the 27 identified compounds, including the retention time, chemical formula, measured mass, MS2 fragments, and chemical structure, are presented in Table 1 and Table S1 and Figure 2.

2.2. In Vitro Antitumor Effect of Buf and Cino

In this part, Buf and Cino were selected as model drugs to evaluate the anti-HCC activity of toad venom. Firstly, the in vitro antiproliferation activity of Buf and Cino alone or in combination were evaluated via CCK8 assay. As illustrated in Figure 3A,B, the cell viability of HepG2 cells were decreased with the increased drug concentration, showing apparent dose-dependent effects. The half-maximal inhibitory concentration (IC50) of the Buf and Cino were 26.753 μM and 86.025 μM, respectively. Subsequently, the combined therapeutic effect of Buf and Cino was also investigated (Figure 3C). According to quantitative analysis of the Chou–Talalay method, the CI values of Buf/Cino at molar ratios of 2:1, 1:1, and 1:2 were 0.871, 0.760, and 0.845, respectively. These data intuitively indicated that Buf/Cino with a smaller CI value has a stronger antitumor effect in vitro. Compared with the individual drugs, the combined in vitro antitumor efficacy was significantly increased. Therefore, to achieve a maximized effect, the combination of Buf and Cino with a molar ratio of 1:1 was chosen as the optimal regimen for further study. Finally, Hoechst 33342 staining showed a decrease in the number of cells in the treatment group and confirmed the production of apoptotic bodies, which proved that Buf and Cino had proapoptotic activity on HepG2 cells (Figure 3D).

2.3. In Vivo Antitumor Effect of Buf and Cino

To verify the antitumor effect of Buf and Cino in vivo, HepG2 cells’ tumor-bearing nude mice model was established. The whole experimental design and protocol are illustrated in Figure 4A. On day 20, the tumor tissues and main organs (heart, liver, spleen, lung, and kidney) of nude mice were isolated and weighed. Tumor volume of the mice treated with NS increased rapidly, showing no inhibitory effects on tumor growth. Compared with the NS group, the tumor growth of nude mice injected with Buf and Cino was statistically inhibited over time (Figure 4B,E). Notably, the tumor inhibition rate of the Mix group was 88.49% (Figure 4D). Meanwhile, there was no obvious decrease in body weight of nude mice, suggesting a suitable dosage during the whole treatment (Figure 4C). The H&E staining of major organs further verified the safety of drug treatment. All the organs had intact cell morphology, without nuclear shrinkage, cytoplasmic lysis, and inflammatory cell infiltration (Figure 5). As shown in Figure 6, the tumor cells in the NS group had the complete morphology, with clearly visible nuclei and no nuclear pyknosis. However, after the administration of Buf and Cino, the tumor cells became scattered and exhibited varying degrees of necrosis, such as nuclear fragmentation, nuclear dissolution, and even the disappearance of the nucleus. The proliferation and apoptosis of tumor cells were also evaluated by Ki-67 staining and TUNEL staining, respectively. Results demonstrated that the decreased positive Ki-67 expression and increased apoptosis rate of tumor cells were obviously observed compared with the NS group (Figure 6).

2.4. Nontargeted Metabolomics Analysis of Tumor Tissues

2.4.1. Stability Test of Tumor Metabolomic Method

Firstly, 20 μL of each sample was mixed into QC samples and detected five times. On PCA analysis, the distribution of QC samples with positive and negative modes were closely concentrated, suggesting a good experimental repeatability and stable detecting system (Figure 7E and Figure S1E). Further, Pearson correlation analysis was performed on QC samples, as shown in Figure 7D and Figure S1D. The experimental results showed that the correlation coefficients between QC samples were above 0.9, indicating that the experimental repeatability is good. Furthermore, Hotelling’s T2 test, Multivariate Control Chart (MCC), and relative standard deviation of ion peak abundance in QC samples were performed, which also reflected the true and reliable data with the detecting system (Figure 7A–C and Figure S1A–C).

2.4.2. Identification of Differential Metabolites

The structures of metabolites in tumor tissue samples were identified by matching the retention time, molecular mass (molecular mass error within <10 ppm), second-order fragmentation profile, collision energy, and other information of the metabolites in the local database. Ultimately, 621 metabolites were identified, including 363 metabolites in the positive ion mode and 258 metabolites in the negative ion mode. Among them, the proportions of lipid and lipid-like molecules, organic acids, and benzenoids were about 30%, 18%, and 13%, respectively. All metabolites detected in positive and negative ion mode were differentially analyzed based on univariate analysis. The metabolites with FC > 1.5 or FC < 0.67 and p-value < 0.05 between two groups were visualized in the form of a volcano diagram (Figure 8A). Further, the unsupervised 3D-PCA was applied to assess the overall distribution among all samples and the stability of the entire analysis process. As shown in Figure 8B, the two groups of samples were well-separated, suggesting a large number of different types of metabolites that presented between Mix and NS. Subsequently, the supervised PLS-DA were established to verify the difference between the two groups. As illustrated in Figure 8C, there was the satisfactory separation on the PLS-DA model between the NS and Mix groups. The values of R2Y and Q2 in the PLS-DA model between Mix and NS were 0.994 and 0.57, respectively (Figure 8D). Additionally, the data performed by permutation test revealed that the obtained model exhibited no overfitting behavior. Differential metabolites were screened with the requirement of VIP > 1 and p < 0.05 in the test. Totally, 89 potential differential metabolites were obtained, including 20 differential metabolites between Buf and NS, 13 differential metabolites between Cino and NS, and 56 differential metabolites between Mix and NS. The detailed information of these differential metabolites is presented in Table S2.

2.4.3. Hierarchical Cluster Analysis of Differential Metabolites

The expression of differential metabolites was calculated by distance matrix and analyzed by hierarchical cluster to comprehensively and visually obtain the differences in the expression patterns of metabolites in different samples. The clustering heatmap analysis is shown in Figure 9B. Compared with the NS group, the expression of 50 species of metabolites in the Mix group were upregulated, while 6 species of metabolites were downregulated (see Table S2). The tumor metabolomics results revealed that the abundance of metabolite species in nude mice tumor changed significantly after the administration of Buf and Cino. When compared with the NS group, the differential metabolites in the Mix group intervened on by Buf and Cino showed an upregulated trend, including phenylethylamine, linolenic acid, norvaline, L-serine methyl ester, L-leucyl-l-leucine methyl ester, and methyltyrosinate; also, differential metabolites exhibited a downregulated trend, such as Carbobenzyloxy-l-norvalyl-l-norleucine, 1,2-dioleoyl-sn-glycero-3-phosphatidylcholine, 1-stearoyl-2-linoleoyl-sn-glycerol, 7-hydroxymitragynine, 3-hydroxybenzaldehyde, and Pg 42:11. Figure 9A displays the correlation analysis of differential metabolites between the Mix and NS groups. Red color indicates a positive correlation, while blue color indicates a negative correlation. The higher value represents the closer correlation between two metabolites.

2.4.4. Pathway Analysis of Differential Metabolites

Differential metabolites pathway enrichment analysis is of great significance for revealing the pathogenesis of diseases. In this study, to explore the metabolic pathways associated with the efficacy of Buf and Cino in the treatment of HCC, the differential metabolites were imported into a MetaboAnalyst 5.0 database, and pathways with impact value > 0.1 were marked as potential target pathways. As a result, compared with the NS group, differential metabolite pathways of Buf and Cino used individually were both involved in nicotinate and nicotinamide metabolism (Figure 10A). Meanwhile, phenylalanine metabolism and α-linolenic acid metabolism were closely related to the therapeutic effects of Buf and Cino against HCC. As illustrated in Figure 10B,C, the differential metabolic pathways were colored and displayed by using the KEGG pathway mapper function according to the upregulated and downregulated information. On the one hand, phenylalanine could regulate substance metabolism as well as be a mediating factor of apoptosis. It is essential for the production of the nonessential tyrosine, the process of which was firstly discovered by Womack and Ross in 1934 [18]. This conversion is catalyzed by phenylalanine hydrolase, which is primarily active in the liver. Tyrosine and phenylalanine levels were decreased in the serum of HCC patients [19]. Moustafa R.K. et al. found that the pathway that allows the conversion of phenylalanine to tyrosine may lead to the accumulation of phenylalanine, which further leads to mitochondria-mediated apoptosis and thereby causes cancer cell death [20]. In this paper, both in vitro and in vivo assays have demonstrated the HepG2 cell apoptosis induced by the cotreatment of Buf and Cino. Therefore, we may speculate that the increased level of phenylalanine in HCC tumor after Buf and Cino intervention may be closely related to the decreased activity of phenylalanine hydrolase. On the other hand, in the fatty acid metabolic pathway, fatty acid synthesis is an important process that converts nutrients into metabolic intermediates for cell biosynthesis, energy storage, and signal transduction [21]. Cell proliferation is a general feature of all cancerous cells and requires fatty acids to synthesize signal transduction molecules and tumor cell membranes [22]. During the fatty acid synthesis process, fatty acid synthase plays a key role. The expression level of fatty acid synthase in cancer cells was significantly higher than that in normal cells [23]. A study had shown that α-linolenic acid, a polyunsaturated fatty acid, can inhibit the proliferation and metastasis of osteosarcoma cells by inhibiting the expression of fatty acid synthase [24]. Our results showed that α-linolenic acid was lowly expressed in HCC tumors, while Buf and Cino intervention increased its levels. Results of in vitro and in vivo experiments indicated that Buf and Cino have a notable effect on inhibiting HepG2 cell proliferation and inducing HepG2 cell apoptosis. Taken together, the anti-HCC effect of Buf and Cino may be attributed to the downregulation of fatty acid synthase activity and disorder of fatty acid synthesis.

3. Discussion

HCC is one of the severe malignancies that causes death in humans [25]. Screening high-efficiency and low-toxicity anti-HCC drugs from TCM has always been a research hotspot. TCM with multiple components, multiple efficacies, and multiple targets regulates the state of the body as a whole so as to play a therapeutic role. Since most oral TCMs have the first-pass effect, only components that are absorbed into the blood can function effectively [26]. Toad venom has a long history in the treatment of HCC, the main antitumor components of which are bufadienolides [5]. According to previous studies, Buf and Cino were considered as the main antitumor components in toad venom, the respective anticancer effects of which were closely linked with the downregulation of prosurvival proteins (Bcl-2) and upregulation of the proapoptotic proteins (Fas and Bax) [14,27]. Up to now, the anticancer mechanisms of Buf and Cino were evaluated individually, which were insufficient to guide their combined clinical use. Therefore, the mechanism of the combined effect of Buf and Cino should be further emphasized.
During this study, the absorbed components of toad venom in rat plasma were identified by using the UHPLC-HR-Q-Exactive Orbitrap MS system to better screen bioactive components. As mentioned previously, these two components played an important role in the treatment of HCC by toad venom. In our current study, the optimal ratio of Buf and Cino was selected as 1:1 via CCK-8 assay. Results of in vitro studies found that the combination therapy of Buf and Cino showed effective action on the inhibition of HepG2 cell proliferation and increase in HepG2 cell apoptosis, exhibiting an enhanced in vitro antitumor activity compared with individual treatment. Therefore, we explored the perturbation of endogenous metabolites in tumor tissue after the cotreatment of Buf and Cino, hoping to explain the in-depth mechanism of toad venom against HCC. After the intervention of Buf and Cino, the tumor growth of the Mix group was significantly inhibited. No obvious toxicity was observed in H&E staining results of major organs. An in vivo antitumor mechanism of Buf and Cino showed that Buf and Cino could inhibit cell proliferation and induce cell apoptosis on HepG2 tumor-bearing nude mice. Further, the tumor metabolomics results found that the abundance of many metabolite species in nude mice tumor changed significantly after the administration of Buf and Cino. The pathway enrichment of differential metabolites between the Mix group and the NS group were mainly involved in phenylalanine metabolism and α-linolenic acid metabolism.
It has been shown that the phenylalanine metabolism pathway was altered in gastric cancer and prostate cancer [28,29]. According to our studies, the biomarker that participated in phenylalanine metabolism in HCC was phenethylamine, which is an aromatic amino acid affecting the pathology of hepatitis B [14]. For example, phenylalanine metabolism efficiency decreased in the gut of patients with early chronic hepatitis B [30]. Moreover, phenylalanine could regulate substance metabolism as well as be a mediating factor of apoptosis. The induction of apoptosis by phenethylamine may be involved in three mechanisms: the FasR-mediated cell death receptor pathway, the Rho/ROCK pathway to activate mitochondria mediated apoptosis, and LAMP2 in graninase B signaling-mediated apoptosis [31]. Lamin B1 and PAK1 were found in the Fas/Fas ligand death receptor pathway, which have previously been shown to be involved in apoptosis [20]. Lamin B1, as the main component of the subnuclear layer, plays an important role in maintaining the integrity of the nuclear membrane. Destruction of nuclear membrane integrity is a sign of apoptosis. During the apoptosis process, Lamin B1 mRNA levels were decreased, which may be due to the induction of p53 or pRB tumor suppressor pathways [32]. References also indicated that downstream effector PAK1 of Fas/Fas ligand complex prevents apoptosis by restricting the expression of proapoptotic proteins or regulating post-translational modification of effectors [20]. Moustafa R. K et al. further confirmed that the altered phenylalanine metabolism activated several apoptosis pathways associated with key proteins such as HADHA and ACAT1 [20]. In view of the upregulation of phenethylamine in HCC with Buf and Cino cotreatment, an induction on cell apoptosis capacity can be speculated, which is consistent with the apoptotic phenotype observed in Buf and Cino-treated cells and tumor tissues.
α-linolenic acid metabolism is an extremely vital metabolic pathway that interferes with tumor proliferation and metabolism [33,34]. α-linolenic acid, a polyunsaturated fatty acid, has a variety of benefits, including improved immunity and anticancer and anti-inflammatory effects [35]. Due to its anti-inflammatory and antitumor properties, α-linolenic acid has been reported to be a potential adjuvant therapy candidate for breast, pancreatic, and colon cancer [36]. Currently, the main attention on the possible mechanism of α-linolenic acid inhibiting the growth of different breast cancer cell lines is increasing lipid peroxidation, altering the expression and activation of cell membrane receptors, changing transcription factors, and regulating tumor suppressor genes such as phosphatase and tensin homologous genes [37]. A study has shown that α-linolenic acid regulates the growth of cancer cells in breast cancer (MCF-7 and MDA-MB-231) and cervical cancer (SiHa and HeLa) by reducing NO release and inducing lipid peroxidation [38]. Julie K. et al. found that α-linolenic acid can reduce HER2-overexpressing breast cancer growth and demonstrated for the first time that docosahexaenoic acid (DHA) is responsible for the effects of α-linolenic acid-rich diets on HER2 signaling pathways [39]. Wang et al. found that α-linolenic acid could inhibit the migration of human triple-negative breast cancer cells by reducing the expression of Twist1 and inhibiting Twist1-mediated epithelial mesenchymal transformation [35]. Li et al. reported that dietary supplementation with α-linolenic acid induced n-3 LCPUFAs conversion and reduced prostate cancer growth in a mouse model [33]. Interestingly, the content of α-linolenic acid in HCC was lower than that of normal tissues [40]. In HepG2 cells, α-linolenic acid suppressed the proliferation, migration, and invasion. An evaluation of the underlying mechanisms suggested that α-linolenic acid enhances FXR expression and thus inhibits the Wnt/β-catenin signaling pathway, thereby hinting HCC progression [40]. Furthermore, Cis9 and trans11 conjugated linoleic acid induce HepG2 cell apoptosis by activating PPAR-γ signaling pathway [41]. Vecchini et al. reported that dietary intake of α-linolenic acid could decrease the expression of COX-2 and increase the apoptosis of liver cancer cells [42]. In this study, the abundance of α-linolenic acid detected in HepG2 tumor treated with Buf and Cino was significantly increased compared with the NS group, indicating their potential contribution to the progression of HCC.
Meanwhile, we noticed that amino acid metabolites were also altered in HepG2 tumor-bearing nude mice after the cotreatment of Buf and Cino, such as norvaline and Leu-Ala. Among them, norvaline, an isomer of valine, was involved in the cytotoxic activity of macrophages against breast tumor cells [43]. Furthermore, previous reports showed that norvaline could interact with gut microbiota to affect the occurrence of colorectal cancer [44]. A clinical study has confirmed that the levels of norvaline and Leu-Ala in plasma were significantly lower in patients with liver cancer than in healthy people. Compared with the NS group, the levels of these two amino acids were recovered after cotreatment of Buf and Cino, which is consistent with previous studies. Additionally, it is worth considering that some metabolites without ionization cannot be detected by UHPLC-MS/MS. So, in the future, nuclear magnetic resonance technology will be applied to make up for the deficiency of mass spectrometer detection [3]. Due to the limitation of the KEGG database and the MetaboAnalyst 5.0 database, several identified differential metabolites could not be matched with the corresponding KEGG ID, so they failed to participate in pathway enrichment, resulting in a small number of enriched pathways. Summarily, using Buf and Cino as model drugs, this study confirmed the antitumor activity of toad venom in HepG2 tumor-bearing nude mice and revealed its antitumor mechanism at a metabolic level. In our future work, transcriptomics and proteomics will be involved in deeper mechanism research of toad venom against HCC.

4. Materials and Methods

4.1. Materials and Reagents

Dried toad venom secretion originated from Bufo melanostictus Schneider was provided by Beijing newborn toad breeding center (Beijing, China) and authenticated as toad venom by Professor Qingrong Pu from the affiliated traditional Chinese medicine hospital of Southwest Medical University (Luzhou, China). Bufalin and cinobufagin were provided by Chengdu Ruifensi Biotechnology Co., LTD. Methanol, acetonitrile, and formic acid were purchased from Thermo Fisher Scientific (Shanghai, China) Co., Ltd. Cell-Counting-Kit-8 tests were obtained from Beyotime Biotechnology Co., Ltd. (Shanghai, China). Other reagents used in this study were of analytical reagent grade.

4.2. Cells and Animals

HepG2 cells were provided by Professor Jie Qing from the affiliated traditional Chinese medicine hospital of Southwest Medical University (Luzhou, China) and cultured in a DMEM high-glucose medium at 37 °C in an atmosphere containing 5% CO2. The female BALB/c nude mice (five weeks old) were purchased from Chongqing Tengxin Biotechnology Co., LTD (Chongqing, China). The nude mice were housed in a specific-pathogen-free (SPF) cleanliness level barrier system with a room temperature of 25 ± 2 °C and a relative humidity of 50 ± 10%. Animal welfare and experimental procedures were strictly performed in accordance with the Guide for the Care and Use of Laboratory Animals and the related ethics regulations of Southwest Medical University. Animal protocols were reviewed and approved by the Internal Animal Care and Use Committee of Southwest Medical University (Approval Number: 20211115-010).

4.3. In Vivo Compound Identification of Toad Venom in Rat Plasma

The sample preparation and UHPLC-MS/MS detection conditions are presented in the Supplementary Material. Based on the standard substances database established in our previous research, the prototypical compounds of toad venom in rat plasma were identified. Combined with reported references, metabolites of toad venom were further identified according to the first and second mass spectra of prototypes. The compounds absorbed in plasma were considered biologically active and deserved further study [45].

4.4. Cell Counting Kit-8 (CCK-8) Assay

The CCK-8 assay was used to study the antiproliferation of Buf and Cino. Briefly, HepG2 cells were grown in 96-well plates at a density of 1 × 105 per well prior to the treatment with various concentrations of Buf (1.563, 3.125, 6.25, 12.5, 25, 50, 100, 200, 400 nM), Cino (12.5, 25, 50, 100, 200, 400, 800, 1600 nM), and Buf/Cino (2:1, 1:1, 1:2) for 24 h. Subsequently, the culture medium was removed and replaced by 10 μL CCK-8 reagent for an additional 2 h incubation. Finally, the absorbance (Abs) of each well was measured at 450 nm using a microplate reader (Synergy2, Biotek, Highland Park, IL, USA). The cell viability was calculated. The combination index (CI) was calculated according to the Chou–Talalay method, used for quantitative analysis of combination therapy in vitro [46].
Cell viability = [(AM − AB)/(AC − AB)]
AM, AC, and AB are defined as the absorbances of the drug-treated group, control group, and blank group, respectively.
CI = ( D 50 ) A D A + ( D 50 ) B D B
where (D) A and (D) B are the concentration of A and B used in combined therapy to achieve 50% inhibitory effects; (D50) A and (D50) B are the concentration of A and B to achieve the same effects individually.

4.5. Hoechst 33342 Staining Assay

HepG2 cells were grown in 6-well plates at a density of 2.5 × 105 per well for 24 h. The culture medium was removed and replaced by Buf (12.5 nM), Cino (12.5 nM), and Buf/Cino (12.5/12.5 nM), respectively. After 24 h incubation, each well was added with 4% formaldehyde solution for the 30 min of fixation at −20 °C. After being washed twice with phosphate buffer, HepG2 cells were stained with Hoechst33342 solution for 10 min, then discarded and washed three times. Finally, a fluorescence microscope (Evos, Life Technologies, Carlsbad, CA, USA) was used to observe cell apoptosis.

4.6. In Vivo Antitumor Efficacy

To obtain tumor-bearing nude mice, the female BALB/c nude mice were injected with HepG2 cells (about 1 × 107 cells) into the right armpit. When the volume of tumors reached about 100 mm3, nude mice were divided randomly into four groups (seven mice per group) as follows, (1) Normal Saline (NS, ip every two days, same dosing volume as Buf group), (2) Buf (ip every two days, 1 mg/kg), (3) Cino (ip every two days, 1.14 mg/kg), and (4) Mix (Buf (ip, 1 mg/kg) plus Cino (ip, 1.14 mg/kg)). Tumor volume and body weight of nude mice were monitored every two days during 20-day treatment. At day 20, the nude mice were sacrificed by cervical dislocation. The tumor volume (V) and the tumor inhibition rate (TIR) were calculated, respectively. The excised organs were used to evaluate the safety of the therapy. Pathological examinations, including H&E staining, Ki-67 staining, and TUNEL staining, were performed on the tumor tissues. The detailed procedures are shown in the Supplementary Material.
V = 0.5 × length × wide2
TIR = [(WC − WE) / WC] × 100%
where WE and WC are defined as the tumor weight of the experiment group and control group.

4.7. Nontargeted Metabolomics Analysis of Tumor Tissues

4.7.1. Metabolites Extraction

Firstly, 80 mg tumor tissue was weighed and slowly thawed at 4 °C. After that, the sample was immersed in 200 μL water to homogenize and treated with 800 μL precooled methanol/acetonitrile solution (1:1, v/v) to sonicate for 30 min at 4 °C. Then, it was placed at −20 °C for 2 h to precipitate protein. Subsequently, after centrifuging for 20 min at 4 °C at the rate of 14,000 rpm, the supernatant was taken and freeze-dried. Finally, 100 μL acetonitrile/water solution (1:1, v/v) was added to redissolve and vortex for 1 min. After centrifugation again at 14,000 gpm for 15 min at 4 °C, the supernatant was taken for UHPLC-Q-Exactive Orbitrap MS analysis.

4.7.2. UHPLC-Q-Exactive Orbitrap MS Analysis

In this study, UHPLC-Q-Exactive Orbitrap MS analysis was conducted with the application of the ThermoFisher Scientific Vanquish Ultra-High Performance Liquid Chromatography system associated with the ThermoFisher Scientific Q-Exactive Orbitrap mass spectrometer. The chromatographic separation was performed on the Xbridge BEH HILIC column (2.1 × 100 mm, 2.5 μm). The temperatures of the sample disk and column oven were set at 4 °C and 25 °C, respectively. The mobile phase consisted of solvent A (water containing 25 mM ammonium acetate and 25 mM ammonia solution) and B (acetonitrile). At a 0.3 mL/min flow rate, the injection volume was kept at 2 μL. The gradient elution procedure was as follows: 0–1.5 min, 98% B, 1.5–12 min, 98–2% B, 12–14 min, 2% B, 14–14.1 min, 2–98% B, and 14.1–17 min, 98% B. In order to avoid the influence of instrument detection signal fluctuation, random sequence was used for continuous analysis of samples. Quality control (QC) samples were inserted into the sample queue to monitor and evaluate the stability of the system and the reliability of the experimental data.
Further, the samples were analyzed by Q-Exactive Orbitrap mass spectrometer. Both the primary and secondary spectrograms of the samples were collected by electrospray ionization (ESI) in positive and negative ion modes. ESI source and mass spectrum setting parameters were as follows: atomizing gas auxiliary heating 1 (Gas1); 60, auxiliary heating 2 (Gas2); 60, curtain gas (CUR); 30 psi, ion source temperature; 600 °C, spray voltage (ISVF); ±5500 V (positive and negative two modes). Primary mass charge ratio detection range: 80–1200 Da; resolution: 60,000; scanning accumulation time: 100 ms. The secondary stage adopted the segmental acquisition method, with a: scanning range: 70–1200 Da; secondary resolution, 30,000; scanning accumulation time: 50 ms; dynamic exclusion time: 4 s.

4.7.3. MS Data Processing and Analysis

The original data were converted into MzXML format by ProteoWizard software. XCMS software was used for peak alignment, retention time correction, and peak area extraction. PCA and PLS-DA were performed via SIMCA-14.1 software. The screening condition of differential metabolites were VIP > 1 and p < 0.05 in the t-test, which was chosen as the in-depth metabolic pathway analysis. Eventually, pathway enrichment was obtained from the MetaboAnalyst 5.0 (https://www.metaboanalyst.ca, 10 September 2022) platform. All the data were expressed as mean ± standard deviation (SD). The independent sample’s t-test of two groups was performed by SPSS 26.0 (SPSS Inc., Chicago, IL, USA).

5. Conclusions

The current work identified the active components of toad venom in rat plasma. Next, we conducted in vitro and in vivo experiments to explore the underlying mechanism of Buf and Cino against HCC. The antitumor effect after the cotreatment of Buf and Cino was mainly through regulating the phenylalanine metabolism and α-Linolenic acid metabolism. It is suggested that amino acid metabolites such as norvaline and Leu-Ala may play the important role regarding toad venom in treating HCC. Taken together, these results suggest the use of toad venom as a promising candidate for further clinical studies against HCC.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/molecules27227758/s1, Table S1. MS2 fragment of 27 compounds from toad venom extract detected in rat plasma. Table S2. Details of differential metabolites of HepG2 cell xenograft tumor among different groups compared with control group. Figure S1. Results of quality control. A Hotelling’s T2 test for the total sample in negative mode. B Multivariate control chart of QC samples in negative mode. C Relative standard deviation of QC samples in negative mode. D Pearson correlation analysis of QC samples in negative mode. E Principal component analysis of QC samples in negative mode. Figure S2. A and E Volcanograms showing up-and down-regulated metabolites between NS and Buf, NS and Cino group, respectively. Blue color indicates downregulated metabolites, while red color represents upregulated metabolites. B and F Score plot of 3D principal component analysis between NS and Buf, NS and Cino group, respectively. C and G Score plot of partial least-squares discriminant analysis between NS and Buf, NS and Cino group, respectively. D and H Corresponding validation plot between NS and Buf, NS and Cino group, respectively. Figure S3. A and C Correlation analysis of differential metabolites between NS and Buf, NS and Cino group, respectively. Red color indicates positive correlation and blue color indicates negative correlation. The higher the value, the stronger the correlation between the two metabolites. B and D Heatmap analysis of differential metabolites between NS and Buf, NS and Cino group, respectively. Red color indicates higher concentration of differential metabolites, blue color indicates concentration of differential metabolites. Figure S4. A-B Metabolic pathway of differential metabolites between NS and Buf, NS and Cino group enriched by MetaboAnalyst 5.0.

Author Contributions

Investigation, validation, data analysis, curation, and writing—original draft and visualization, P.L.; investigation, validation, and data analysis, Y.M.; data analysis, L.Y. and L.M.; visualization, C.S., Z.L. and M.M.; conceptualization and methodology, Q.S., S.Y. and W.R. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Special Funding for Traditional Chinese Medicine Scientific Research from Sichuan Administration of Traditional Chinese Medicine (Number: 2020JC0150), Projects of the Affiliated Traditional Chinese Medicine Hospital of Southwest Medical University (Number: [2019] 119 and [2020] 33), University-level project of Southwest Medical University (Number: 2021ZKQN118), Innovation Team and Talents Cultivation Program of National Administration of Traditional Chinese Medicine (Number: ZYYCXTD-C-202207), and Innovation Team of Sichuan Provincial Administration of Traditional Chinese Medicine (Number: 2022C007).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The Supplementary File shows all data produced in the course of this study.

Acknowledgments

Many thanks to the Jie Qing from the affiliated traditional Chinese medicine hospital of Southwest Medical University for the provision of HepG2 cells.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

ACAT1, Mitochondrial acetyl-CoA acetyltransferase; Buf, bufalin; Bcl-2, B-cell lymphoma-2; Bax, BCL2-Associated X; Cino, cinobufagin; CCK-8, cell counting Kit-8; CI, combination index; ESI, electrospray ionization; Fas, recombinant factor-related apoptosis; FasR, recombinant factor-related apoptosis receptor; FXR, Farnesoid X receptor; HCC: hepatocellular carcinoma; HADHA, hydroxyl-coenzyme A dehydrogenase/3-ketoacyl-coenzyme A thiolase; LAMP2, lysosomal-associated membrane protein 2; MCC, multivariate control chart; PCA, principal component analysis; PLS-DA, partial least squares discriminant analysis; QC, quality control; Rho/ROCK, Rho-associated kinase; TCM, traditional Chinese medicine.

References

  1. Wei, L.; Wang, Z.; Jing, N.; Lu, Y.; Yang, J.; Xiao, H.; Guo, H.; Sun, S.; Li, M.; Zhao, D.; et al. Frontier progress of the combination of modern medicine and traditional Chinese medicine in the treatment of hepatocellular carcinoma. Chin. Med. 2022, 17, 90. [Google Scholar] [CrossRef] [PubMed]
  2. Huang, J.; Chen, F.; Zhong, Z.; Tan, H.Y.; Wang, N.; Liu, Y.; Fang, X.; Yang, T.; Feng, Y. Interpreting the Pharmacological Mechanisms of Huachansu Capsules on Hepatocellular Carcinoma Through Combining Network Pharmacology and Experimental Evaluation. Front. Pharmacol. 2020, 11, 414. [Google Scholar] [CrossRef] [PubMed]
  3. Li, Z.; Bao, H. Anti-tumor effect of Inonotus hispidus petroleum ether extract in H22 tumor-bearing mice and analysis its mechanism by untargeted metabonomic. J. Ethnopharmacol. 2022, 285, 114898. [Google Scholar] [CrossRef] [PubMed]
  4. Fan, Y.; Ma, Z.; Zhao, L.; Wang, W.; Gao, M.; Jia, X.; Ouyang, H.; He, J. Anti-tumor activities and mechanisms of Traditional Chinese medicines formulas: A review. Biomed. Pharmacother. 2020, 132, 110820. [Google Scholar] [CrossRef] [PubMed]
  5. Wei, W.L.; Hou, J.J.; Wang, X.; Yu, Y.; Li, H.J.; Li, Z.W.; Feng, Z.J.; Qu, H.; Wu, W.Y.; Guo, D.A. Venenum bufonis: An overview of its traditional use, natural product chemistry, pharmacology, pharmacokinetics and toxicology. J. Ethnopharmacol. 2019, 237, 215–235. [Google Scholar] [CrossRef]
  6. Zhang, D.M.; Liu, J.S.; Tang, M.K.; Yiu, A.; Cao, H.H.; Jiang, L.; Chan, J.Y.; Tian, H.Y.; Fung, K.P.; Ye, W.C. Bufotalin from Venenum Bufonis inhibits growth of multidrug resistant HepG2 cells through G2/M cell cycle arrest and apoptosis. Eur. J. Pharmacol. 2012, 692, 19–28. [Google Scholar] [CrossRef]
  7. Qiu, D.Z.; Zhang, Z.J.; Wu, W.Z.; Yang, Y.K. Bufalin, a component in Chansu, inhibits proliferation and invasion of hepatocellular carcinoma cells. BMC Complement. Altern. Med. 2013, 12, 185. [Google Scholar] [CrossRef] [Green Version]
  8. Qi, F.; Inagaki, Y.; Gao, B.; Cui, X.; Xu, H.; Kokudo, N.; Li, A.; Tang, W. Bufalin and cinobufagin induce apoptosis of human hepatocellular carcinoma cells via Fas- and mitochondria mediated pathways. Cancer Sci. 2011, 102, 951–958. [Google Scholar] [CrossRef]
  9. Zhang, L.; Liang, B.; Xu, H.; Gong, Y.; Hu, W.; Jin, Z.; Wu, X.; Chen, X.; Li, M.; Shi, L.; et al. Cinobufagin induces FOXO1-regulated apoptosis, proliferation, migration, and invasion by inhibiting G9a in non-small-cell lung cancer A549 cells. J. Ethnopharmacol. 2022, 291, 115095. [Google Scholar] [CrossRef]
  10. Ren, W.; Chen, S.; Liao, Y.; Li, S.; Ge, J.; Tao, F.; Huo, Q.; Zhang, Y.; Zhao, Z. Near-infrared fluorescent carbon dots encapsulated liposomes as multifunctional nano-carrier and tracer of the anticancer agent cinobufagin in vivo and in vitro. Colloids Surf. B Biointerfaces 2019, 174, 384–392. [Google Scholar] [CrossRef]
  11. Yang, A.L.; Wu, Q.; Hu, Z.D.; Wang, S.P.; Tao, Y.F.; Wang, A.M.; Sun, Y.X.; Li, X.L.; Dai, L.; Zhang, J. A network pharmacology approach to investigate the anticancer mechanism of cinobufagin against hepatocellular carcinoma via downregulation of EGFR-CDK2 signaling. Toxicol. Appl. Pharmacol. 2021, 431, 115739. [Google Scholar] [CrossRef]
  12. Kim, G.H.; Fang, X.Q.; Lim, W.J.; Park, J.; Kang, T.B.; Kim, J.H.; Lim, J.H. Cinobufagin Suppresses Melanoma Cell Growth by Inhibiting LEF1. Int. J. Mol. Sci. 2020, 21, 6706. [Google Scholar] [CrossRef]
  13. Ren, W.; Han, L.; Luo, M.; Bian, B.; Guan, M.; Yang, H.; Han, C.; Li, N.; Li, T.; Li, S.; et al. Multi-component identification and target cell-based screening of potential bioactive compounds in toad venom by UPLC coupled with high-resolution LTQ-Orbitrap MS and high-sensitivity Qtrap MS. Anal. Bioanal. Chem. 2018, 410, 4419–4435. [Google Scholar] [CrossRef]
  14. Zhang, J.; Hong, Y.; Xie, P.; Chen, Y.; Jiang, L.; Yang, Z.; Cao, G.; Chen, Z.; Liu, X.; Chen, Y.; et al. Spatial Lipidomics Reveals Anticancer Mechanisms of Bufalin in Combination with Cinobufagin in Tumor-Bearing Mice. Front. Pharmacol. 2020, 11, 593815. [Google Scholar] [CrossRef]
  15. Hou, Z.; Song, F.; Xing, J.; Zheng, Z.; Liu, S.; Liu, Z. Comprehensive fecal metabolomics and gut microbiota for the evaluation of the mechanism of Panax Ginseng in the treatment of Qi-deficiency liver cancer. J. Ethnopharmacol. 2022, 292, 115222. [Google Scholar] [CrossRef]
  16. Zhao, Y.Y.; Cheng, X.L.; Vaziri, N.D.; Liu, S.; Lin, R.C. UPLC-based metabonomic applications for discovering biomarkers of diseases in clinical chemistry. Clin. Biochem. 2014, 47, 16–26. [Google Scholar] [CrossRef] [Green Version]
  17. Tian, D.; Gao, Q.; Lin, J.; Chang, Z.; Wang, Y.; Shi, Y.; Su, R.; Han, Z.; Ma, D. Uncovering the mechanism of the Shenzhi Jiannao formula against vascular dementia using a combined network pharmacology approach and molecular biology. Phytomedicine 2021, 90, 153637. [Google Scholar] [CrossRef]
  18. Ploder, M.; Neurauter, G.; Spittler, A.; Schroecksnadel, K.; Roth, E.; Fuchs, D. Serum phenylalanine in patients post trauma and with sepsis correlate to neopterin concentrations. Amino Acids 2007, 35, 303–307. [Google Scholar] [CrossRef]
  19. Wiggins, T.; Kumar, S.; Markar, S.R.; Antonowicz, S.; Hanna, G.B. Tyrosine, phenylalanine, and tryptophan in gastroesophageal malignancy: A systematic review. Cancer Epidemiol. Biomarkers Prev. 2015, 24, 32–38. [Google Scholar] [CrossRef] [Green Version]
  20. Ali, M.R.; Wu, Y.; Han, T.; Zang, X.; Xiao, H.; Tang, Y.; Wu, R.; Fernandez, F.M.; El-Sayed, M.A. Simultaneous Time-Dependent Surface-Enhanced Raman Spectroscopy, Metabolomics, and Proteomics Reveal Cancer Cell Death Mechanisms Associated with Gold Nanorod Photothermal Therapy. J. Am. Chem. Soc. 2016, 138, 15434–15442. [Google Scholar] [CrossRef]
  21. Rohrig, F.; Schulze, A. The multifaceted roles of fatty acid synthesis in cancer. Nat. Rev. Cancer 2016, 16, 732–749. [Google Scholar] [CrossRef]
  22. Wu, H.; Cheng, H.; Luo, S.; Peng, C.; Zhou, A.; Chen, Z.; Wu, H.; Li, Q. Use of cellular metabolomics and lipidomics to decipher the mechanism of Huachansu injection-based intervention against human hepatocellular carcinoma cells. J. Pharm. Biomed. Anal. 2022, 212, 114654. [Google Scholar] [CrossRef]
  23. Crous-Maso, J.; Palomeras, S.; Relat, J.; Camo, C.; Martinez-Garza, U.; Planas, M.; Feliu, L.; Puig, T. (-)-Epigallocatechin 3-Gallate Synthetic Analogues Inhibit Fatty Acid Synthase and Show Anticancer Activity in Triple Negative Breast Cancer. Molecules 2018, 23, 1160. [Google Scholar] [CrossRef] [Green Version]
  24. Fan, H.; Huang, W.; Guo, Y.; Ma, X.; Yang, J. alpha-Linolenic Acid Suppresses Proliferation and Invasion in Osteosarcoma Cells via Inhibiting Fatty Acid Synthase. Molecules 2022, 27, 2741. [Google Scholar] [CrossRef]
  25. Kong, W.S.; Shen, F.X.; Xie, R.F.; Zhou, G.; Feng, Y.M.; Zhou, X. Bufothionine induces autophagy in H22 hepatoma-bearing mice by inhibiting JAK2/STAT3 pathway, a possible anti-cancer mechanism of cinobufacini. J. Ethnopharmacol. 2021, 270, 113848. [Google Scholar] [CrossRef]
  26. Chen, J.; Liu, Z.; Deng, F.; Liang, J.; Fan, B.; Zhen, X.; Tao, R.; Sun, L.; Zhang, S.; Cong, Z.; et al. Mechanisms of Lian-Gui-Ning-Xin-Tang in the treatment of arrhythmia: Integrated pharmacology and in vivo pharmacological assessment. Phytomedicine 2022, 99, 153989. [Google Scholar] [CrossRef] [PubMed]
  27. Li, F.J.; Hu, J.H.; Ren, X.; Zhou, C.M.; Liu, Q.; Zhang, Y.Q. Toad venom: A comprehensive review of chemical constituents, anticancer activities, and mechanisms. Arch. Pharm. 2021, 354, e2100060. [Google Scholar] [CrossRef] [PubMed]
  28. Zhao, Y.; Lv, H.; Qiu, S.; Gao, L.; Ai, H. Plasma metabolic profiling and novel metabolite biomarkers for diagnosing prostate cancer. RSC Adv. 2017, 7, 30060–30069. [Google Scholar] [CrossRef] [Green Version]
  29. Lario, S.; Ramirez-Lazaro, M.J.; Sanjuan-Herraez, D.; Brunet-Vega, A.; Pericay, C.; Gombau, L.; Junquera, F.; Quintas, G.; Calvet, X. Plasma sample based analysis of gastric cancer progression using targeted metabolomics. Sci. Rep. 2017, 7, 17774. [Google Scholar] [CrossRef] [Green Version]
  30. Yao, X.; Li, Z.; Gong, X.; Fu, X.; Xiao, X.; He, M.; Huang, B.; Xu, Z. Total saponins extracted from Abrus cantoniensis Hance suppress hepatitis B virus replication in vitro and in rAAV8-1.3HBV transfected mice. J. Ethnopharmacol. 2020, 249, 112366. [Google Scholar] [CrossRef]
  31. Qi, G.; Wang, B.; Zhang, Y.; Li, H.; Li, C.; Xu, W.; Jin, Y. Living-Cell Imaging of Mitochondrial Membrane Potential Oscillation and Phenylalanine Metabolism Modulation during Periodic Electrostimulus. Anal. Chem. 2019, 91, 9571–9579. [Google Scholar] [CrossRef]
  32. Freund, A.; Laberge, R.M.; Demaria, M.; Campisi, J. Lamin B1 loss is a senescence-associated biomarker. Mol. Biol. Cell 2012, 23, 2066–2075. [Google Scholar] [CrossRef]
  33. Fang, Z.; Huang, H.; Wang, L.; Lin, Z. Identification of the alpha linolenic acid metabolism-related signature associated with prognosis and the immune microenvironment in nasopharyngeal carcinoma. Front. Endocrinol. 2022, 13, 968984. [Google Scholar] [CrossRef]
  34. Gu, Z.; Suburu, J.; Chen, H.; Chen, Y.Q. Mechanisms of omega-3 polyunsaturated fatty acids in prostate cancer prevention. Biomed. Res. Int. 2013, 2013, 824563. [Google Scholar] [CrossRef] [Green Version]
  35. Wang, S.C.; Sun, H.L.; Hsu, Y.H.; Liu, S.H.; Lii, C.K.; Tsai, C.H.; Liu, K.L.; Huang, C.S.; Li, C.C. alpha-Linolenic acid inhibits the migration of human triple-negative breast cancer cells by attenuating Twist1 expression and suppressing Twist1-mediated epithelial-mesenchymal transition. Biochem. Pharmacol. 2020, 180, 114152. [Google Scholar] [CrossRef]
  36. Li, Q.; Liu, X.; Wang, X.; Qiu, S.; Byambasuren, K.; Dang, L.; Wang, Z. Antiproliferative Ability and Fluorescence Tracking of alpha-Linolenic Acid-Loaded Microemulsion as Label-Free Delivery Carriers in MDA-MB-231 Cells. J. Agric. Food Chem. 2019, 67, 11518–11526. [Google Scholar] [CrossRef]
  37. Wiggins, A.K.; Mason, J.K.; Thompson, L.U. Growth and gene expression differ over time in alpha-linolenic acid treated breast cancer cells. Exp. Cell Res. 2015, 333, 147–154. [Google Scholar] [CrossRef]
  38. Kim, K.B.; Nam, Y.A.; Kim, H.S.; Hayes, A.W.; Lee, B.M. alpha-Linolenic acid: Nutraceutical, pharmacological and toxicological evaluation. Food Chem. Toxicol. 2014, 70, 163–178. [Google Scholar] [CrossRef]
  39. Mason, J.K.; Klaire, S.; Kharotia, S.; Wiggins, A.K.; Thompson, L.U. alpha-linolenic acid and docosahexaenoic acid, alone and combined with trastuzumab, reduce HER2-overexpressing breast cancer cell growth but differentially regulate HER2 signaling pathways. Lipids Health Dis. 2015, 14, 91. [Google Scholar] [CrossRef] [Green Version]
  40. Feng, S.; Xie, X.; Chen, C.; Zuo, S.; Zhao, X.; Li, H. Alpha-linolenic acid inhibits hepatocellular carcinoma cell growth through Farnesoid X receptor/beta-catenin signaling pathway. Nutr. Metab. 2022, 19, 57. [Google Scholar] [CrossRef]
  41. Lu, G.; Zhang, G.; Zheng, X.; Zeng, Y.; Xu, Z.; Zeng, W.; Wang, K. c9, t11- conjugated linoleic acid induces HCC cell apoptosis and correlation with PPAR-γ signaling pathway. Am. J. Transl. Res. 2015, 7, 2752–2763. [Google Scholar] [PubMed]
  42. Vecchini, A.; Ceccarelli, V.; Susta, F.; Caligiana, P.; Orvietani, P.; Binaglia, L.; Nocentini, G.; Riccardi, C.; Calviello, G.; Palozza, P.; et al. Dietary alpha-linolenic acid reduces COX-2 expression and induces apoptosis of hepatoma cells. J. Lipid Res. 2004, 45, 308–316. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Zhu, Q.; Zai, H.; Zhang, K.; Zhang, X.; Luo, N.; Li, X.; Hu, Y.; Wu, Y. L-norvaline affects the proliferation of breast cancer cells based on the microbiome and metabolome analysis. J. Appl. Microbiol. 2022, 133, 1014–1026. [Google Scholar] [CrossRef]
  44. Coker, O.O.; Liu, C.; Wu, W.K.K.; Wong, S.H.; Jia, W.; Sung, J.J.Y.; Yu, J. Altered gut metabolites and microbiota interactions are implicated in colorectal carcinogenesis and can be non-invasive diagnostic biomarkers. Microbiome 2022, 10, 35. [Google Scholar] [CrossRef]
  45. Qu, J.; Ke, F.; Liu, Z.; Yang, X.; Li, X.; Xu, H.; Li, Q.; Bi, K. Uncovering the mechanisms of dandelion against triple-negative breast cancer using a combined network pharmacology, molecular pharmacology and metabolomics approach. Phytomedicine 2022, 99, 153986. [Google Scholar] [CrossRef] [PubMed]
  46. Wei, M.M.; Wu, P.; Gao, J.N.; Wang, Y.J.; Zhou, J.L.; Xie, T.; Dong, X.M.; Fang, C.; Song, G.Q.; Huang, Z.H. A combination index and glycoproteomics-based approach revealed synergistic anticancer effects of curcuminoids of turmeric against prostate cancer PC3 cells. J. Ethnopharmacol. 2021, 267, 113467. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Identification of active components of toad venom in rat blood. (A) Blank blood. (B) Toad venom in rat blood. (C) Toad venom extract.
Figure 1. Identification of active components of toad venom in rat blood. (A) Blank blood. (B) Toad venom in rat blood. (C) Toad venom extract.
Molecules 27 07758 g001
Figure 2. The detailed chemical structure of identified components of blood in rat plasma.
Figure 2. The detailed chemical structure of identified components of blood in rat plasma.
Molecules 27 07758 g002
Figure 3. The HepG2 cells viability were cultured in different concentrations of drugs for 24 h. (A) Buf, (B) Cino, (C) used in combination. (D) Hoechst 33342 staining was performed by the treatment of Buf, Cino, and the combination, respectively. “*” means p < 0.05, and “**” means p < 0.01 between two groups.
Figure 3. The HepG2 cells viability were cultured in different concentrations of drugs for 24 h. (A) Buf, (B) Cino, (C) used in combination. (D) Hoechst 33342 staining was performed by the treatment of Buf, Cino, and the combination, respectively. “*” means p < 0.05, and “**” means p < 0.01 between two groups.
Molecules 27 07758 g003
Figure 4. The cotreatment of bufalin and cinobufagin inhibited the HepG2 tumor growth. (A) The designed experimental scheme on in vivo nude mice model. (BD) Changes in tumor volume, body weight, and tumor weight during 20 days’ drug intervention. (E) The excised tumor images with the treatment of NS, Buf, Cino, and in combination for 20 days. Values were presented as the mean ± S.D., n = 6. “**” means p < 0.01 between two groups.
Figure 4. The cotreatment of bufalin and cinobufagin inhibited the HepG2 tumor growth. (A) The designed experimental scheme on in vivo nude mice model. (BD) Changes in tumor volume, body weight, and tumor weight during 20 days’ drug intervention. (E) The excised tumor images with the treatment of NS, Buf, Cino, and in combination for 20 days. Values were presented as the mean ± S.D., n = 6. “**” means p < 0.01 between two groups.
Molecules 27 07758 g004
Figure 5. H&E staining analysis of major organs (heart, liver, spleen, lung, kidney) of HepG2 tumor-bearing nude mice model after 20 days’ drug intervention (original magnification: 200×).
Figure 5. H&E staining analysis of major organs (heart, liver, spleen, lung, kidney) of HepG2 tumor-bearing nude mice model after 20 days’ drug intervention (original magnification: 200×).
Molecules 27 07758 g005
Figure 6. Pathological analysis of tumor with different treatment groups. H&E staining, pale pink areas suggest necrosis area. Ki-67 staining, Brown color represents the nucleus of vigorously proliferating cells. TUNEL staining, green fluorescence represents the nucleus of apoptotic cells. (Original magnification: 200× and 400×).
Figure 6. Pathological analysis of tumor with different treatment groups. H&E staining, pale pink areas suggest necrosis area. Ki-67 staining, Brown color represents the nucleus of vigorously proliferating cells. TUNEL staining, green fluorescence represents the nucleus of apoptotic cells. (Original magnification: 200× and 400×).
Molecules 27 07758 g006
Figure 7. Results of quality control. (A) Hotelling’s T2 test for the total sample in positive mode. (B) Multivariate control chart of QC samples in positive mode. (C) Relative standard deviation of QC samples in positive mode. (D) Pearson correlation analysis of QC samples in positive mode. (E) Principal component analysis of QC samples in positive mode.
Figure 7. Results of quality control. (A) Hotelling’s T2 test for the total sample in positive mode. (B) Multivariate control chart of QC samples in positive mode. (C) Relative standard deviation of QC samples in positive mode. (D) Pearson correlation analysis of QC samples in positive mode. (E) Principal component analysis of QC samples in positive mode.
Molecules 27 07758 g007
Figure 8. (A) Volcanograms showing up- and downregulated metabolites between NS and Mix group. Blue color indicates downregulated metabolites, while red color represents upregulated metabolites. (B) Score plot of 3D principal component analysis between NS and Mix group. (C) Score plot of partial least-squares discriminant analysis between NS and Mix group. (D) Corresponding validation plot from NS and Mix group.
Figure 8. (A) Volcanograms showing up- and downregulated metabolites between NS and Mix group. Blue color indicates downregulated metabolites, while red color represents upregulated metabolites. (B) Score plot of 3D principal component analysis between NS and Mix group. (C) Score plot of partial least-squares discriminant analysis between NS and Mix group. (D) Corresponding validation plot from NS and Mix group.
Molecules 27 07758 g008
Figure 9. (A) Correlation analysis of differential metabolites between NS and Mix group. Red color indicates positive correlation and blue color indicates negative correlation. The higher the value, the stronger the correlation between the two metabolites. (B) Heatmap analysis of differential metabolites between NS and Mix group. Red color indicates higher concentration of differential metabolites, blue color indicates concentration of differential metabolites.
Figure 9. (A) Correlation analysis of differential metabolites between NS and Mix group. Red color indicates positive correlation and blue color indicates negative correlation. The higher the value, the stronger the correlation between the two metabolites. (B) Heatmap analysis of differential metabolites between NS and Mix group. Red color indicates higher concentration of differential metabolites, blue color indicates concentration of differential metabolites.
Molecules 27 07758 g009
Figure 10. (A) Metabolic pathway of differential metabolites between NS and Mix group enriched by MetaboAnalyst 5.0. (B,C) KEGG enrichment pathway of α-linolenic acid metabolism and phenylalanine metabolism, respectively.
Figure 10. (A) Metabolic pathway of differential metabolites between NS and Mix group enriched by MetaboAnalyst 5.0. (B,C) KEGG enrichment pathway of α-linolenic acid metabolism and phenylalanine metabolism, respectively.
Molecules 27 07758 g010
Table 1. UHPLC-HR-Q-Exactive Orbitrap MS data of 27 compounds from toad venom extract detected in rat plasma.
Table 1. UHPLC-HR-Q-Exactive Orbitrap MS data of 27 compounds from toad venom extract detected in rat plasma.
NoCompoundTypeRT
(Min)
Chemical FormulaError (Ppm)[M+H]+
Measured Mass
[M+H]+
Calculated Mass
1Bufotenidine-5-O-
glucuronide
M2.54C19H27O7N2−0.399395.18112395.1812
2BufotenidineP2.62C13H19ON2−0.020219.14917219.1492
3Pitneloyl arginineP5.28C13H25O5N4−1.060317.18185317.1820
4Suberoyl arginineP7.43C14H27O5N4−0.533331.19702331.1976
51-hydroxylarenobufaginM9.15C24H33O7−1.546433.22128433.2148
65-hydroxylarenobufaginM9.65C24H33O7−1.477433.22131433.2148
716-hydroxylarenobufaginM10.01C24H33O7−1.200433.22153433.2148
8HydratalbufarenoginM10.06C24H35O7−0.490435.23709435.2377
911α-hydroxytelocinobufaginP10.52C24H35O6−1.491419.24234419.2428
10Ψ-BufarenoginP10.75C24H33O6−1.522417.22638417.2272
11GamabufotalinP11.34C24H35O5−1.837403.24756403.2479
12BufarenoginP11.94C24H33O60.539417.22726417.2272
1319-hydroxylbufalinP12.43C24H35O5−0.101403.24799403.2479
14ArenobufaginP12.88C24H33O6−1.522417.22733417.2272
15DesacetylbufotalinP13.16C24H35O5−0.324403.24774403.2479
16HellebrigenolP13.38C24H35O60.035419.24289419.2428
17HellebrigeninP13.76C24H33O6−0.132417.22711417.2272
18BufotalininP15.97C24H31O6−1.241415.21106415.2115
19TelocinobufaginP22.47C24H35O5−0.101403.24774403.2479
20CinobufaginolP24.43C26H35O7−0.784459.23801459.2377
21DesacetylcinobufaginP26.84C24H33O5−0.901401.23135401.2323
22BufalinP26.95C24H35O4−1.178387.25244387.2530
23ResibufogeninP28.69C24H33O4−0.950385.23697385.2373
24CinobufaginP28.89C26H35O6−0.265443.24225443.2428
25UnknownM29.76C17H27O2−1.013263.20032263.2084
26Cinobufagin-3-O-
glucuronide
M31.19C31H41O12−0.076605.45453605.4545
27UnknownM31.57C20H33O2−1.824305.24695305.2553
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liang, P.; Ma, Y.; Yang, L.; Mao, L.; Sun, Q.; Sun, C.; Liu, Z.; Mazhar, M.; Yang, S.; Ren, W. Uncovering the Mechanisms of Active Components from Toad Venom against Hepatocellular Carcinoma Using Untargeted Metabolomics. Molecules 2022, 27, 7758. https://doi.org/10.3390/molecules27227758

AMA Style

Liang P, Ma Y, Yang L, Mao L, Sun Q, Sun C, Liu Z, Mazhar M, Yang S, Ren W. Uncovering the Mechanisms of Active Components from Toad Venom against Hepatocellular Carcinoma Using Untargeted Metabolomics. Molecules. 2022; 27(22):7758. https://doi.org/10.3390/molecules27227758

Chicago/Turabian Style

Liang, Pan, Yining Ma, Luyin Yang, Linshen Mao, Qin Sun, Changzhen Sun, Zengjin Liu, Maryam Mazhar, Sijin Yang, and Wei Ren. 2022. "Uncovering the Mechanisms of Active Components from Toad Venom against Hepatocellular Carcinoma Using Untargeted Metabolomics" Molecules 27, no. 22: 7758. https://doi.org/10.3390/molecules27227758

APA Style

Liang, P., Ma, Y., Yang, L., Mao, L., Sun, Q., Sun, C., Liu, Z., Mazhar, M., Yang, S., & Ren, W. (2022). Uncovering the Mechanisms of Active Components from Toad Venom against Hepatocellular Carcinoma Using Untargeted Metabolomics. Molecules, 27(22), 7758. https://doi.org/10.3390/molecules27227758

Article Metrics

Back to TopTop