- 1Department of Rehabilitation Medicine, Teaching Base of Huludao Central Hospital Affiliated to Jinzhou Medical University, Huludao, Liaoning, China
- 2Department of Family Medicine, Shengjing Hospital of China Medical University, Shenyang, Liaoning, China
- 3Department of Science and Education, Teaching Base of Huludao Central Hospital Affiliated to Jinzhou Medical University, Huludao, Liaoning, China
- 4Department of Pharmacy, Teaching Base of Huludao Central Hospital Affiliated to Jinzhou Medical University, Huludao, Liaoning, China
- 5Clinical Laboratory,Teaching Base of Huludao Central Hospital Affiliated to Jinzhou Medical University, Huludao, Liaoning, China
Background: Lung adenocarcinoma (LUAD) causes millions of deaths annually worldwide. Parthenolide (PTL), extracted from traditional Chinese herbal medicines, has various biological activities. In this study, we investigated the effects of PTL on amino acid metabolism and oxidative stress in LUAD cells.
Methods: This study identified differential proteins and potential mechanisms of action of PTL in LUAD cells through label-free quantitative proteomics and protein-protein interaction networks. Combined with targeted amino acid metabolomics, we confirmed the results of GO and KEGG analyses. On this basis, the potential targets of PTL in LUAD were identified through network pharmacology, molecular simulation docking, and molecular dynamics simulations. Finally, the effects of PTL on amino acid metabolism and oxidative stress in LUAD were verified using in vivo and in vitro experiments.
Results: PTL treatment of LUAD cells resulted in significant changes in expression of 157 proteins. GO and KEGG enrichment analyses showed that these proteins were involved in amino acid metabolism and oxidative stress response. Targeted amino acid metabolomics further confirmed that PTL affected amino acid metabolism in LUAD. Network pharmacology, molecular docking, and molecular dynamics simulations identified GCTG as a potential target of PTL in LUAD. Meanwhile, in vitro and in vivo experimental results indicated that PTL targeting GCTG affected the proliferation, amino acid metabolism, and oxidative stress levels of LUAD cells.
Conclusion: PTL affects proliferation, amino acid metabolism, and oxidative stress levels of LUAD cells via targeting GCTG. Therefore, our study provides new insights into the prevention and treatment of LUAD with PTL, which may lay the foundation for future research directions
1 Background
Lung cancer is one of the most common type of cancers worldwide and a leading cause of cancer-related deaths (1). The World Health Organization (WHO) estimates that there are approximately 2 million new cases and 1.76 million deaths each year, and these numbers are constantly increasing (2). Medical experts have deepened their understanding of the biological behavior and pathogenesis of lung cancer over the past few decades and there have been breakthroughs in its treatment (3, 4). However, the 5-year survival rate of lung cancer is still low, and the recurrence rate of the disease and poor quality of life of patients have not been effectively addressed (5). In China, these issues have become even more prominent, becoming a major public health issue that imposes a heavy burden on society and families.
According to the histopathological classification of lung cancer, adenocarcinoma accounts for approximately half of its incidence rate (6, 7). Therefore, diagnosis and treatment of lung adenocarcinoma (LUAD) is an important part of lung cancer research. Although the proliferation rate of LUAD is lower than that of small-cell lung cancer, the occurrence of malignant biological behaviors is relatively late (8, 9). However, owing to the large population of China, early cancer screening is not yet universal. Most patients reach an advanced stage when they visit the hospital, and their 5-year survival rate is not optimistic (10). Currently, surgery, radiotherapy, and chemotherapy are the cornerstones of LUAD treatment; moreover, emerging interventional therapies and immunotherapies in recent years have provided patients with more choices (11, 12). However, the uncertainty of LUAD prognosis and various toxic side effects not only bring more difficulties to the treatment itself, but also seriously affect the quality of life of the patients (13). Therefore, the identification of safer and more effective treatments is urgently warranted in LUAD research.
Traditional Chinese medicine (TCM) has been used in China for thousands of years. Modern medicine has been proven to have significant therapeutic effects on various diseases, including tumors (14). In the treatment of tumors, TCM is mainly used to alleviate adverse reactions after surgery, radiotherapy, and chemotherapy as well as to improve the quality of life of the patients (15, 16). However, each TCM decoction contains many types of herbs, and each herb contains multiple compounds. Therefore, it is difficult to explain the effects of TCM on tumors using specific mechanisms, which greatly hinders the research progress in this field (17). Therefore, there is a critical need to extract effective single components from various herbs and explore their effects on tumors.
Parthenolide (PTL) is a sesquiterpene lactone derived primarily from Tanacetum parthenium (18). The nucleophilic properties of the methylene gamma-lactone ring and epoxy group enable them to rapidly interact with biological targets, induce cellular oxidative stress, and exert biological effects (19, 20). PTL is widely used to treat high fever, headache, toothache, rheumatoid arthritis, and other inflammatory diseases (21). In recent years, the antitumor effect of PTL has been observed in many tumors, such as hepatocellular carcinoma and breast cancer, and it does not affect normal cell function.
Several studies have reported the effects of PTL on LUAD. They mostly focus on the action of PTL on the NF-kB pathway and STAT3 pathway, leading to cell apoptosis (22, 23). However, no studies have explored the impact of PTL on LUAD from the perspectives of metabolic reprogramming and oxidative stress. As is well known, the malignant growth of tumors requires continuous acquisition of energy, leading to significant metabolic heterogeneity between tumor tissue and normal tissue. Since the proposal of the Warburg effect, tumor metabolic reprogramming has been the focus of many scientists. Recent research has demonstrated that in addition to sugar metabolism, which provides energy for tumor cell proliferation, amino acid metabolism also plays an irreplaceable role in the occurrence and development of tumors. The occurrence and progression of tumors are often accompanied by an imbalance between the oxidative and antioxidant systems, leading to intracellular DNA damage, protein oxidation, and lipid peroxidation (24). The metabolism of amino acids such as glutamine and serine may cause changes in oxidative stress levels within tumor cells (25, 26). Therefore, we investigated whether PTL affects amino acid metabolism and oxidative stress levels in LUAD cells.
2 Methods
2.1 Reagents
PTL was purchased from GLPBIO (Montclair, USA) and the human LUAD cell line, H1975, was purchased from the Shanghai Institutes for Biological Sciences, China. The RPMI 1640 medium was purchased from Thermo Fisher Scientific (Waltham, MA, USA). Fetal bovine serum was purchased from Oumashi Biotechnology Co., Ltd. (Shanghai, China). All amino acids and stable isotope-labelled standards were obtained from Sigma-Aldrich (St. Louis, MO, USA). Analytical grade ammonium acetate was obtained from Sigma-Aldrich (St. Louis, MO, USA). Methanol (Optima LC-MS), acetonitrile (Optima LC-MS), and formic acid (Optima LC-MS) were purchased from Thermo Fisher Scientific (FairLawn, NJ, USA). Ultrapure water was purchased from Millipore (Billerica, MA). TCEP (tris (2-carboxyethyl) phosphine) and IAA (iodoacetamide) were purchased from Sigma-Aldrich (St. Louis, MO, USA). Sequence-grade trypsin (sequence grade) was purchased from Promega (Madison, WI, USA). NP-40 (Nonidet P 40) and acetone were purchased from Sangon Biotech (Shanghai, China). The protease inhibitor cocktail was obtained from Kangchen Biotech (Shanghai, China). RNAkey Total RNA Extraction Kit was purchased from Beijing Saiwen Innovation Biotechnology Co., Ltd. The qPCR quantitative analysis and reverse transcription kits were purchased from Novozan Biotechnology(Nanjing,China). The poly (vinylidene fluoride) membrane and ECL reagent were purchased from Sigma (USA). The relevant antibodies were provided by Abclonal (Beijing, China). An amino acid extraction kit was purchased from Solarbio (Beijing, China). The NADPH/NADP detection kit was obtained from Wuhan Saiyu Biotechnology Co., Ltd. The GSH/GSSG detection kit was provided by Nanjing Jiancheng Bioengineering Research Institute.
2.2 Cell culture and treatment
H1975 cells were cultured in complete RPMI 1640 medium containing 10% fetal bovine serum, 100 U/mL streptomycin, and 100 U/mL penicillin. H1975 cells were incubated in an environment of 5% CO2 and 37°C. According to the appropriate concentration (inhibitory concentration 50 [IC50]) standardized in the initial experiment, PTL was dissolved in 0.05% DMSO, diluted to 10 μM with PBS, added to the treatment group, and an equal amount of complete culture medium or culture medium lacking amino acids was added to the control group.
2.3 Label-free quantitative proteomics
The process was divided into nine steps as follows: 1) sample lysis; 2) BCA(Bicinchoninic Acid Assay) quantification; 3) acetone precipitation; 4) protein reconstitution, reduction, alkylation, and enzymatic hydrolysis; 5) removal of SDC(sodium deoxycholate); 6) peptide desalination; 7) Nano HPLC(High Performance Liquid Chromatography) separation; 8) LC-MS/MS(Liquid Chromatography-Tandem Mass Spectrometry); and 9) MaxQuant analysis and LFQ(label-free quantification).
2.4 Protein-protein interaction network construction and GO and KEGG analysis
We constructed a PPI network of differentially expressed proteins using the String data platform (https://cn.string-db.org/), and the final graph was implemented using the R package”igraph” (1.2.6). The biological processes, molecular functions, and related signaling pathways of the differentially expressed proteins were analyzed using MetaScape(http://metascape.org/gp/index.html#/main/step1) software. The parameters related to the enrichment analysis were set to a minimum overlap of 3, p-value cutoff of 0.01, and minimum enrichment of 1.5. Using R language - ggplot2 [3.4.4] package was used for GSEA (Gene Set Enrichment Analysis) Analysis.
2.5 Standard solution preparation and extraction of metabolites
Stock solutions of the individual amino acids were mixed and prepared in an amino acid-free matrix to obtain a series of amino acid calibrators. Certain concentrations of L-Alanine-d4 and Phenylanine-d2 were mixed as internal standards (IS). The stock solution of all of these and the working solution were stored in -20°C. The samples were added to 450 μL water and 2 mL of acetonitrile/methanol (1:1), which contained mixed internal standards, and vortexed for mixing. Next, the samples were placed in liquid nitrogen and thawed at room temperature. Subsequently, the samples underwent ultrasonic treatment for 4 min, was placed in liquid nitrogen again for 5 min, and thawed at room temperature. The samples were centrifuged at 12,000 rpm for 20 min at 4°C. Finally, the supernatant was injected into a LC-MS/MS system for analysis.
2.6 LC-MS method
An ultra-high performance liquid chromatography coupled to tandem mass spectrometry (UHPLC-MS/MS) system (ExionLC™ AD UHPLC-QTRAP 6500+, AB SCIEX Corp., Boston, MA, USA) was used to quantitate amino acids in Novogene Co., Ltd. (Beijing, China). Separation was performed on an ACQUITY UPLC BEH Amide column (2.1×100mm, 1.7μm), which was maintained at 50°C. The mobile phase, consisting of 0.1% formic acid in 5 mM ammonium acetate (solvent A) and 0.1% formic acid in acetonitrile (solvent B), was used at a flow rate of 0.30 mL/min. The solvent gradient was set as follows: initial 80% B, 0.5 min; 80–70% B, 2 min; 70–45% B, 4 min; 45–80% B, 6.01min; 80% B, 9 min.
The mass spectrometer was operated in positive multiple reaction mode (MRM). Parameters were as follows: IonSpray Voltage (5500 V), Curtain Gas (35 psi), Ion Source Temp (550°C), and Ion Source Gas of 1 and 2 (50 and 60 psi).
2.7 Construction of pharmacological networks
To explore the relevant targets of the disease, the keyword “lung adenocarcinoma” was searched using OMIM (http://omim.org/), GeneCards (http://www.genecards.org/), and DisGeNET (https://disgenet.com/) databases. Relevant targets of parthenolide were extracted from the SEA search server (https://sea.bkslab.org/), PharmMapper server (https://lilab-ecust.cn/pharmmapper/index.html), and SwissTargetPrediction(http://swisstargetprediction.ch/). The intersection of disease targets, drug targets, and proteomic results was obtained. The intersection of these targets was considered as a potential target for the treatment of LUAD with parthenolide. Subsequently, the standard names of the targets were retrieved from UniProt KB (http://www.uniprot.org/). Relevant targets were obtained from differentially expressed amino acid metabolites using MetScape in Cytoscape 3.8.0. Finally, a compound–target–amino acid network diagram was constructed that included the relationships between parthenolide, related target genes, and amino acids.
2.8 Molecular docking
The PDB files of the targets were downloaded from the PDB database (RCSB PDB, http://www.rcsb.org/pdb/). They were preprocessed using Pymol software, and protein hydrogenation and charge calculations were performed using Autodock software and converted into pdbqt files. The SDF file for the compound was downloaded from PubChem (https://pubchem.ncbi.nlm.nih.gov/). Chem3Dpro software was used to calculate its minimum energy, convert it into a mol file, and use autodock to convert it into pdbqt; then Vina molecule docking was used, with the ligand position in the PDB file as the docking interface; the CONF file of each folder was the docking location; Vs_desult was the docking result; The binding energy between components and targets, as well as the number of hydrogen bonds formed, were important criteria for evaluating molecular docking results. The lower the binding energy and the more hydrogen bonds there are, the more stable the binding becomes, and greater the possibility of target molecule interactions. Pymol (open-source) was used to draw 3D and 2D graphs, and Discovery Studio 4.5 Client was used to display docking sites.
2.9 Molecular dynamic simulation
A 100 ns MD simulation was performed on the composite using Gromacs 2023. The protein adopted CHARMM 36 force field parameters, and the ligand topology was constructed using GAFF2 force field parameters. The protein-ligand complex was placed in a cubic box under periodic boundary conditions. The box was filled with water molecules using the TIP3P water model. The particle grid Ewald (PME) and Verlet algorithms were used to separate the electrostatic interactions. Subsequently, 100000 steps of the isothermal isostatic ensemble equilibrium and isothermal isobaric ensemble equilibrium were performed using a coupling constant of 0.1 ps and a duration of 100 ps. The van der Waals and Coulomb interactions were calculated using a cut-off value of 1.0 nm. Finally, the system was subjected to molecular dynamics simulations using Gromacs 2023 at a constant temperature (300 K) and pressure (1 bar) for a total duration of 100 ns.
2.10 RT-qPCR(Reverse Transcription Quantitative Polymerase Chain Reaction)
Cells were divided into three groups: control, PTL, and PTL + OE-GCTG (also known as GGCT). The concentration of PTL was 10 μM. qPCR experiments were conducted 48 h after cell transfection. Plasmids were prepared using GenePharma (Shanghai, China). qPCR of extracted mRNA was performed according to the manufacturer’s instructions. The primer sequence was as follows: GCTG (upstream): 5’- TCGGCGGCGTTCTTCTGTG-3’; GCTG (downstream):
5’- GCTATCCCTCCATGCCAAGTTTG-3’; GAPDH (upstream):
5’- GGAGTCCACTGGCGTCTTCAC-3’; GAPDH (downstream):
5’- ATTGCTGATGATCTTGAGGCTGTTG-3’.
In addition, RT-qPCR experiments involving human tissues were performed after obtaining informed consent from the patients and were approved by the Ethics Committee of Shengjing Hospital (Ethics Code: 2023PS1232K). For basic experimental analysis. 32 pairs of cancer and adjacent normal tissue samples were collected from patients who underwent the same surgery at Shengjing Hospital between January 2019 and December 2021.
2.11 Western blot
Proteins were extracted from the three groups of cells mentioned earlier, quantified, and denatured. Sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) was used to separate the same amount of protein, which was then transferred to a polyvinylidene fluoride membrane with a pore diameter of 0.45 μM. Subsequently, the membrane was sealed with a rapid blocking solution and incubated overnight with primary antibody at 4°C. After washing, the membrane was incubated with a secondary antibody at room temperature for 1 h. Finally, ECL reagent was used to visualize the protein bands on a gel imaging system (GE Healthcare). In addition, the grayscale values of each band were analyzed using ImageJ software.
2.12 Cell clone formation experiment
In each group, 800 cells/well were inoculated in a 6-well culture plate, which consisted of complete culture medium of 2 ml per well, with three wells in each group. The inoculated cells were shaken well and incubated on a clean bench for 20 minutes, then cultured in the incubator. The medium was changed every 3 days, the cell status observed, the supernatant discarded, and the cells washed once with PBS until the number of cells in most individual clones in the well exceeded 50. Next, 1 mL of 4% paraformaldehyde was added to each well and the cells were fixed at 4°C for 60 min and washed once with PBS. Crystal violet staining solution (1 mL) was added to each well and the cells were stained for 15 min. The cells were then washed several times with PBS, air-dried, photographed, and counted.
2.13 Amino acid content determination experiment
Cells were incubated using RPMI 1640 medium without amino acids. Cells were collected into centrifuge tubes and the supernatant was discarded after centrifugation. According to the instructions, 1 mL of the corresponding reagent was added to every 5 million cells, and the cells were sonicated (power 20%, sonication for 3 seconds, interval of 10 seconds, repeated 30 times). Subsequently, the liquid was centrifuged at 10000 rpm and 4°C for 10 minutes, and the supernatant was collected. The enzyme-linked immunosorbent assay (ELISA) reader was preheated for more than 30 min, and the wavelength was adjusted to 570 nm. Amino acid content in the supernatant was measured and calculated according to the manufacturer’s instructions.
2.14 NADPH/NADP detection
Eight ml NADPH probe buffer (component B-II) and 2 ml NADPH probe (component B-I) was mixed together. Two μl of storage solution was diluted with PBS to obtain 2 μM working solution. The cells were seeded in a 6-well plate, counted to 1x106, and treated. The cells were resuspended in 50 μl of component D and repeatedly freeze thawed thrice in liquid nitrogen. A 96-well plate was prepared and 50 μl of PBS solution was added to blank and standard wells 1–6, respectively; 50 μl of NADPH working solution was added to standard well 7. Fifty microliters of the processed sample suspension was added to each experimental well; next, 50 μl of the mix was added to all wells and incubated at room temperature in the dark for 1 h. A ELISA reader was used to detect the OD value at 460 nm.
2.15 GSH/GSSG detection
Cells were plated in a 6-well plate and treated before counting. The number of cells in the sample were adjusted to 1x106, and the cells were resuspended using 110 μl of reagent D. Configuration of 50 μM GSH standard and 50 μM GSSG standard was performed. Next, 10 μl PBS was added to the blank well of the 96-well plate, 10 μl GSH standard to the standard well, and 10 μl of the test cell suspension to the experimental well. Then, 100 μl of reagent A and 10 μl of reagent B was added to each well in sequence, and incubated at room temperature for 2 minutes. Then 50 μl of reagent C was added to each well. After 30 s, the OD value at 450 nm was measured using the ELISA reader and the reading was recorded (A1); after 10 minutes, OD value was measured again under the same conditions and the reading was recorded (A2). Result analysis: GSH (μM)=[Measurement well (A2-A1)/standard (A2-A1)] x standard concentration (50 μM). GSSG was determined as described above.
2.16 Reactive oxygen species detection
5×105 cells were seeded in a 6-well plate and subjected to the corresponding treatment. Next, 10 μM ROS probe reagent was added to each well (leaving one blank well). After incubation at 37°C for 30 minutes, the treated cells were digested with trypsin and resuspended. After adjusting the parameters, the average fluorescence intensity of the blank well cells was set to 10 and the average fluorescence intensity of the other cells was sequentially detected by adding ROS probes.
2.17 Pull-down assay
Biotinylated drug (10 μM) was incubated with streptavidin beads (10 μL) in binding buffer (200 μL total volume) for 1 h at 4°C with rotation.Beads were washed 3× with ice-cold binding buffer to remove unbound drug.Target protein lysate (100 μg) was added to the drug-bound beads and incubated for 2 h at 4°C with rotation.Beads were pelleted using a magnetic rack and washed 3× with binding buffer.
2.18 Animal experiments
All experimental procedures involving animals were approved by the Animal Ethics Committee of the Shengjing Hospital (Ethics Code: 2023PS1263K). Twelve female nude mice (BALB/C, 4–5 weeks old, weight 13–15 g) were purchased from Beijing Huafukang Biotechnology Co., Ltd. All nude mice were kept under SPF conditions (22 ± 0.5°C, 50 ± 10% relative humidity, free access to food and water). Nude mice were housed in a SPF-grade animal breeding center for 1 week. Subsequently, 12 nude mice were randomly divided into control, PTL, and PTL+OE-GCTG groups (n=4). Lentivirus preparation was performed by Genepharma. Normal H1975 or OE-GCTG H1975 cells were digested, centrifuged, resuspended in serum-free medium, and counted to a final concentration of 1x107/ml. Cells were randomly inoculated into the posterior lower part of the right forelimb of nude mice using a sterile syringe, with an injection volume of 120 μl. Normal H1975 cells were inoculated in the control and PTL groups, whereas OE-GCTG H1975 cells were inoculated in the PTL + OE-GCTG group. PTL was dissolved in 0.05% DMSO and diluted in PBS. The PTL and PTL+OE-GCTG groups were intraperitoneally injected with PTL [(20 mg/kg, dose selection based on the results of our previous study (unpublished)]; results showed that 20 mg/kg PTL significantly inhibited the volume of LUAD tumors in nude mice. Therefore, this dose was selected for subsequent analyses. The control group was administered an equal volume of DMSO+PBS without any other treatment. The drugs were injected three times a week for four consecutive weeks, and the body weight and tumor volume of the nude mice were measured twice a week. After 4 weeks, the mice were euthanized using the cervical dislocation method after anesthesia. Photographs of the tumor were taken, and relevant experiments were conducted.
2.19 Statistical analysis
The extracted proteomic and amino acid metabolomic data were analyzed using univariate statistical analysis, multivariate statistical analysis, hierarchical clustering, and correlation analysis. Univariate statistical analysis was conducted using the Student’s t-test and multiple variation analysis. When the p-value was less than 0.05, the distribution of the measured data showed a decrease of less than 0.67 times or an increase of more than 1.5 times, which was considered a significant difference between cells treated with estradiol. The multivariate statistical analysis included unsupervised principal component analysis (PCA), supervised partial least squares discriminant analysis (PLSDA), and orthogonal partial least squares discriminant analysis (OPLSDA). The discrimination of the tested substances was determined by the variable importance (VIP>1) and p-value (p<0.05) of the projection (VIP) parameter. Gel Pro Analyzer and Image J software were used for grayscale analysis, while SPSS 22.0 and GraphPad Prism7.0 software were used for statistical analysis. Quantitative data were presented in the form of mean ± standard deviation. After testing whether the data in each group followed a normal distribution, a t-test or one-way ANOVA was used to detect intergroup differences. P<0.05 indicated statistical significance.
3 Results
3.1 Quantitative proteomics analysis shows significant differences in protein expression
As shown in Figure 1, the flowchart illustrates the study process. A total of 157 differentially expressed proteins were screened by mass spectrometry-based parallel reaction monitoring (PRM), including 60 upregulated and 97 downregulated proteins (Figure 2A, Supplementary Tables S1, S2). As shown in Figure 2B, the significantly altered proteins formed a cluster, indicating that they have similar expression patterns and may be closely related to metabolic processes. As shown in Figure 2C, the PPI network revealed significant differential protein interactions in this group, which may indicate that it plays a special role in certain molecular processes or cellular functions.

Figure 1. Flowchart of the experimental process. The effects of PTL on amino acid metabolism and oxidative stress in lung adenocarcinoma were revealed through quantitative proteomics, targeted amino acid metabolomics, network pharmacology, and experimental verification. PTL, Parthenolide.

Figure 2. Quantitative proteomics reveals 157 significantly differentially expressed proteins. (A) The volcano map contains 60 up-regulated proteins and 97 down regulated proteins; (B) Cluster analysis shows that differentially expressed proteins with similar expression patterns form a cluster; (C) Protein-protein interaction (PPI) network diagram reveals differential protein interactions. The concentration of PTL was 10 μM. The experiment included three biological replicates.
3.2 GO and KEGG analysis results indicate that differential proteins converge in amino acid metabolism processes and anti-oxidative stress pathways
As shown in Figure 3A and Table 1, the differentially expressed proteins mainly focus on “oxidoreductase activity,” “protein homodimerization activity,” and “nucleotide transferase activity” in terms of molecular function. As shown in Figure 3B and Table 2, the differential proteins mainly focus on “response to extracellular stimuli,” “glutamine family amino acid metabolism process,” and “amide metabolism process” in terms of cellular biology processes. Subsequently, these proteins were found to play biological roles in pathways such as “glutathione metabolism,” “alanine, aspartate, and glutamate metabolism,” and “biosynthesis of amino acids” (Figure 3C, Table 3). Similarly, GSEA indicated that the functions of the differentially expressed proteins were mainly focused on amino acid metabolism and antioxidative stress (Figure 3D, Table 4). As shown in Figure 3E, the three most distinct biological effects were demonstrated by GSEA enrichment analysis. Therefore, the results of the GO and KEGG analyses indicate that the proteins affected by PTL may alter the amino acid metabolism and oxidative stress status in LUAD.

Figure 3. Enrichment terms for differential proteins. (A) Molecular functions in GO analysis; (B) Biology processes in GO analysis; (C) KEGG pathway analysis; (D) GSEA enrichment analysis; (E) The top three biological processes ranked. KEGG,Kyoto Encyclopedia of Genes and Genomes; GO,Gene Ontology; GSEA,Gene Set Enrichment Analysis.
3.3 Targeted amino acid metabolomics shows that the amino acid content of LUAD cells was influenced by PTL
As shown in Figure 4A, 23 amino acids were detected in this study; we found significant changes in the content of 12 amino acids after PTL administration (Supplementary Tables S3). Interestingly, all 12 amino acids significantly decreased after PTL administration. As shown in Figure 4B, these significantly altered amino acids formed a cluster, indicating that they had similar metabolic patterns and may be closely related to energy metabolism. Correlation analysis was performed to determine the degree of association between significant amino acids, with the aim of gaining a deeper understanding of the relationships between amino acid types in biological processes. Correlation analysis showed that each amino acid was significantly positively correlated with the other 11 amino acids (P<0.05,r>0.7) (Figure 4C). These results indicate that PTL reduces amino acid production in LUAD, thereby affecting energy metabolism.

Figure 4. Quantitative results of 23 amino acids after PTL action on lung adenocarcinoma cells. (A) Significant differences were observed in 12 amino acids after PTL administration; (B) Cluster analysis shows that significantly altered amino acids have similar metabolic patterns; (C) Correlation analysis shows a positive correlation between 12 amino acids. Ns stands for p ≥ 0.05, ∗p < 0.05, ∗∗p < 0.01. The concentration of PTL was 10 μM. PTL, Parthenolide; Ns, non-significant. The experiment included six biological replicates.
3.4 The results of multivariate statistical analysis validate the significant differences in metabolites after administration
From the PCA scoring chart, it was concluded that there was good repeatability within the groups, the sample data were very similar, and there was good discrimination between the groups (Figure 5A). Similarly, the PLS-DA and OPLS-DA models (Figures 5B, C) showed significant differences between the two groups. External permutation test was performed to avoid overfitting, with a R2 value of 0.429, a Q2 value of -0.279 and with a R2 value of 0.227, a Q2 value of -0.786, respectively (Figures 5D, E). The Variable Importance in Projection (VIP) scores showed that 11 significantly differentially expressed amino acid metabolites had scores greater than 1, indicating that these metabolites enhanced the explanatory power of the model (Figure 5F, Table 5). The volcano map consistently showed similar results (Figure 5G).

Figure 5. Multivariate statistical analysis confirm significant differences in metabolites after administration (A) PCA scoring chart; (B) PLS-DA scoring chart; (C) OPLS-DA scoring chart; (D, E) External permutation test; (F) The Variable Importance in Projection scores; (G) The volcano map indicates that these metabolites enhance the explanatory power of the model. OPLS-DA,orthogonal partial least squares discriminant analysis; PLS-DA,partial least squares discriminant analysis; PCA, principal component analysis.

Table 5. Eleven amino acid components significantly altered by treatment of LUAD cells with parthenolide.
3.5 Network pharmacology analysis to screen the targets of PTL on LUAD
Based on the SwissTarget Prediction (STP), SEA search server, and PharmMapper databases, 148 target genes were associated with PTL (Figure 6A). Based on DisGeNET, GeneCards Search, and OMIM databases, 70 target proteins were associated with LUAD (Figure 6B). When these selected genes were crossed with differentially expressed genes using quantitative proteomics, 29 potential targets of PTL against LUAD were identified (Figure 6C).

Figure 6. Screening of PTL targets through network pharmacology analysis (A) Screening PTL targets based on public databases; (B) Screening targets for lung adenocarcinoma based on public databases; (C) Venn diagram shows 29 potential targets of PTL acting on lung adenocarcinoma. PTL, Parthenolide.
Subsequently, in Cytoscape 3.9.1, 11 different amino acid metabolites regulated by PTL were introduced into MetScape and 41 amino acid metabolite targets were obtained (Figure 7A). Finally, a compound-target-metabolite network consisting of 82 nodes (PTL, 70 related target genes, and 11 amino acid metabolites) and 153 edges was constructed (Figure 7B). Based on parameters such as degree, betweenness centrality, proximity, and stress in the network, GCTG (also known as GGCT), DAO, and BAAT were identified as key targets (Table 6).

Figure 7. Screening of 11 amino acid metabolite targets regulated by PTL using Cytoscape 3.9.1 software. (A) Forty-one amino acid metabolite targets were obtained using MetScape software; (B) A compound-target-metabolites network diagram was constructed.
3.6 Molecular docking analysis of PTL binding to potential targets
We conducted molecular docking studies using AutoDock Vina to further elucidate the interactions between PTL and potential targets. Docking analysis of GCTG showed that PTL formed hydrogen bonds with the SER A: 24 and TRP A: 64 residues at the active site. The van der Waals interactions included TYR A: 22, ASN A:25, ALA A: 69, and TYR A:139 residues. The binding energy between GCTG and PTL was -10.5 (Figure 8A). However, in the molecular docking analysis of DAO and BAAT with PTL, there was no interaction with hydrogen bonds, and the binding energies of the two were -8.5 and -7, respectively (Figures 8B, C). In summary, GCTG had the lowest binding energy for molecular docking with PTL and hydrogen-bonding interactions. Therefore, GCTG was selected as a key target for subsequent experiments.

Figure 8. Molecular docking analysis of PTL binding to potential targets (A–C) Molecular docking analysis demonstrated the binding of GCTG, DAO, and BAAT to PTL, respectively. GCTG, Gamma-Glutamylcyclotransferase; DAO, D-Amino Acid Oxidase; BAAT, Bile Acid-CoA : Amino Acid N-Acyltransferase; PTL, Parthenolide.
3.7 Molecular dynamics analysis shows that there is stable binding between PTL and GCTG
First, the balance of the simulation system was evaluated using root mean square deviation (RMSD). Generally speaking, an RMSD value less than 2 Å and reaching stability within 10–20 ns indicates high conformational stability of the system. As shown in Figure 9A, the composite system reached equilibrium at around 20ns and eventually fluctuated at approximately 1.8 Å. Therefore, PTLs exhibit high stability when bound to a target protein. Further analysis revealed that the radius of gyration (Rg) and solvent-accessible surface area (SASA) of the composite system fluctuated stably during motion, indicating that the composite system remained stable and compact throughout the simulation (Figures 9B, C). The number of hydrogen bonds between PTL and the target protein during the dynamic process is shown in Figure 9D, and the number of hydrogen bonds between the complex system ranged from 0 to 8, indicating that the complex system had good hydrogen bond interactions. As shown in Figure 9E, the root mean square fluctuation (RMSF) of the composite was relatively low (mostly below 2 Å), resulting in low flexibility and high stability. In summary, the binding of the complex was stable, and the complex exhibited good hydrogen bonding properties. Therefore, PTL bind well to its target proteins.

Figure 9. Molecular dynamics analysis confirms the stable binding between PTL and GCTG (A) Root mean square deviation (RMSD); (B) Radius of gyration (Rg); (C) Solvent accessible surface area; (D) Hydrogen bond numbers; (E) Root mean square fluctuation (RMSF). PTL, Parthenolide.
3.8 PTL targeting GCTG affects the proliferation, amino acid metabolism, and oxidative stress levels in LUAD
As shown in Figures 10A–C, bioinformatics analysis revealed that GCTG was significantly upregulated in tumor tissues and could serve as a potential biomarker for predicting prognosis. High expression was significantly associated with poor prognosis. RT-qPCR results of 32 clinical samples also showed that GCTG was significantly overexpressed in cancer tissues (Figure 10D). As shown in Table 7, GCTG expression significantly correlated with tumor size and staging. Therefore, both dry and wet experiments indicated that GCTG was significantly upregulated and had a pro-cancer effect in LUAD. As shown in (Figure 10E), compared to that in the control group, PTL significantly inhibited the expression of GCTG protein, whereas the rescue experiment partially reversed the expression of this protein. As shown in Figures 10F, G, PTL significantly inhibited the proliferation and amino acid metabolism of LUAD cells, whereas rescue experiments partially reversed these biological processes. Correspondingly, as shown in Figures 10H–J, PTL targeting GCTG enhanced ROS expression while inhibiting the GSH/GSSG and NADPH/NADP ratios. To further confirm whether PTL has a direct mechanism of action on GCTG, a pull-down assay was performed, and it was found that PTL directly regulates GCTG (Figure 10K)

Figure 10. In vitro experiment on the effect of PTL targeting GCTG in lung adenocarcinoma. (A) High expression of GCTG in tumor tissues in TCGA database; (B) High expression of GCTG is significantly associated with poor prognosis; (C) ROC curve analysis for GCTG expression in LUAD and paracancerous tissues; (D) GCTG is highly expressed in tumor tissues of 32 pairs of samples; (E) PTL inhibits the expression of GCTG protein; (F) PTL targets GCTG to inhibit the total amino acid content of lung adenocarcinoma cells; (G) PTL targets GCTG to inhibit the proliferation of lung adenocarcinoma cells; (H) PTL targets GCTG to promote the increase of ROS in lung adenocarcinoma cells; (I) PTL targets GCTG to inhibit GSH/GSSG ratio; (J) PTL targets GCTG to inhibit NADPH/NADP ratio. (K) Pull down assay confirms that PTL directly acts on GCTG. Ns stands for p ≥ 0.05, ∗p < 0.05, ∗∗p < 0.01, ∗∗∗p < 0.001. The concentration of PTL was 10 μM. PTL, Parthenolide. All the above in vitro experiments were conducted with three biological replicates.

Table 7. Relationship between GCTG and clinical pathological characteristics of lung adenocarcinoma patients.
As shown in Figure 11A, nude mice were euthanized and their tumors were removed for photography. PTL significantly decreased tumor weight and volume (Figures 11B, C). Immunohistochemical analysis revealed that PTL inhibited GCTG expression (Figure 11D). The amino acid content determination experiment also showed that PTL inhibited amino acid metabolic reprogramming in LUAD tissues (Figure 11E). The aforementioned rescue experiments partially reversed this phenomenon.

Figure 11. In vivo experiment of PTL targeting GCTG to regulate lung adenocarcinoma. (A–C) PTL targets GCTG to inhibit the proliferation of lung adenocarcinoma; (D) Immunohistochemical experiments show that PTL inhibits the expression of GCTG; (E) PTL targeting GCTG inhibits the content of total amino acids in lung adenocarcinoma tissue. Ns stands for p ≥ 0.05, ∗∗p < 0.01, ∗∗∗p < 0.001. The intraperitoneal injection dose of PTL was 20 mg/kg. PTL, Parthenolide. In vivo experiments included four biological replicates.
4 Discussion
PTL is a natural compound, extracted from the small white chrysanthemum (Chrysanthemus morifolium) and has received widespread attention in recent years owing to its various biological activities (27, 28). It has been confirmed that pyrethrum lactone has many effects, such as antioxidant, anti-inflammatory, and immune regulation effects, and has significant therapeutic effects on pulmonary hypertension, diabetes, and certain types of infectious diseases (29–31). PTL also exhibits promising anti-tumor properties, particularly through its interactions with critical signaling pathways such as NF-kB and STAT3, which regulate cell survival and apoptosis (32–36). However, studies on the effects of PTL on LUAD have predominantly focused on these pathways, neglecting the potential role of metabolic reprogramming and oxidative stress in tumor progression. This study aimed to bridge this gap by investigating the metabolic alterations and oxidative stress imbalance induced by PTL in LUAD cells.
To the best our knowledge, this is the first study to investigate the differential protein expression in LUAD cells treated with PTL using PRM technology. Subsequently, through GO and KEGG analyses, differentially expressed proteins were successfully identified in biological processes related to amino acid metabolism. Based on the above proteomic results, we found that 12 of the 23 amino acids were significantly differentially expressed in PTL-treated cells. The combination of the multiple omics experimental methods mentioned above can not only shorten the research time but also accurately position our research direction towards amino acid metabolism reprogramming and oxidative stress.
To further investigate the molecular mechanisms underlying the effects of PTL on amino acid metabolism and oxidative stress in LUAD, we used network pharmacology and molecular docking to explore its downstream targets. Results identified GCTG as a potential site of action for the alteration of amino acid metabolism and oxidative stress in PTL-induced LUAD. γ-glutamyltransferase (GCTG, also known as GGCT) is an important enzyme in glutathione metabolism, catalyzing the reaction of gamma glutamyl peptide to produce 5-hydroxyproline and free amino acids (37, 38). GCTG is highly conserved in many species including bacteria, plants, and other higher organisms (39). Reportedly, GCTG is highly expressed in a variety of cancers, including breast cancer, lung cancer, bladder cancer, colon cancer, osteosarcoma and glioma, suggesting that it may be the oncogene in these tumors (40–43). Moreover, GCTG also promotes the malignant behavior of tumors (44–46). Because GCTG is a crucial catalytic enzyme in amino acid metabolism and its high expression in LUAD tissue was confirmed in clinical specimens, we speculate that PTL inhibits GCTG and thus affects the amino acid reprogramming of cancer cells. In addition, GSEA analysis suggested that PTL can cause changes in oxidative stress-related pathways in tumor cells; therefore, experiments detecting ROS expression levels were also included in the study.
To confirm this hypothesis, both in vivo and in vitro experiments were conducted. Our results indicate that PTL delays tumor growth and amino acid metabolism and disrupts the oxidative stress balance by inhibiting GCTG expression. Amino acid metabolism, an important component of energy metabolic reprogramming, has attracted much attention in the field of cancer. Physapubescin I, a natural compound extracted from edible herbaceous fruits, increases intracellular glutamine levels and correspondingly reduces glutamate and its downstream metabolites, thereby inhibiting cell proliferation and inducing apoptosis (47). Ophiopogon is a tetracyclic diterpenoid derived from plants belonging to the Lamiaceae and Camellia genera. Reportedly, oridonin reverses the resistance of pancreatic cancer cells, PANC-1/Gem, to gemcitabine by inhibiting glutathione S-transferase (GST pi) (48). In addition, oxidative stress is also one of the hotspots in the field of cancer. Baicalin induces ROS production in bladder cancer cells, thereby enhancing the expression of apoptosis-related proteins (49). Resveratrol has a two-way regulatory effect on tumors: low dose activates SIRT1 to maintain oxidative stability, whereas high dose induces mitochondrial ROS bursts. This results in the reversal of chemotherapeutic resistance in stem cells (50). Overall, there has been relatively little research on the effects of traditional Chinese medicine monomers on amino acid metabolism in tumors. However, there is currently no study on the impact of PTL on energy metabolism and oxidative stress in LUAD. Therefore, this study explored the effects of PTL on amino acid metabolism, oxidative stress, and downstream molecular mechanisms in LUAD using through multi omics/network, pharmacology/molecular simulation docking, and experimental validation methods. These findings provide valuable insights into the therapeutic potential of PTL, not only as a standalone treatment but also in combination with existing therapies, thereby enhancing the overall treatment efficacy and improving patient outcomes. The significance of this study lies in its potential to contribute to a more comprehensive understanding of LUAD treatment, paving the way for innovative approaches that address both the biological and metabolic aspects of tumor growth. Interestingly, PTL was previously considered an antioxidant in non-tumor studies; however, in this study, it was found to be an oxidant. This may be due to changes in the nature of PTL caused by different types of diseases.
This study has some limitations. First, the three basic forms of energy metabolism often occur simultaneously in tumors. Therefore, in addition to amino acid metabolism, future studies should investigate the effects of PTL on lipid and glucose metabolisms. Second, to make our results more specific, further research is required to determine which specific amino acids are affected by PTL targeting GCTG. Finally, the relationship between amino acid metabolism and oxidative stress requires further investigation.
5 Conclusion
This is the first study to report the impact of PTL on amino acid metabolic reprogramming and oxidative stress in LUAD. We also explored the effects of PTL on the downstream target, GCTG, using multiple omics, network pharmacology, and wet experiments. This study provides a potential new mechanism by which PTL delays the proliferation of LUAD; however, further research is warranted.
Data availability statement
The data presented in the study are deposited in the iProx (https://www.iprox.cn/page/home.html), accession numbers IPX0013067000 and PXD067441.
Ethics statement
The studies involving humans were approved by Ethics Committee of Shengjing Hospital. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study. The animal study was approved by the Animal Ethics Committee of the Shengjing Hospital. The study was conducted in accordance with the local legislation and institutional requirements.
Author contributions
JiyL: Writing – original draft. YL: Data curation, Writing – review & editing. JiaL: Writing – review & editing, Software. SS: Writing – review & editing.
Funding
The author(s) declare financial support was received for the research and/or publication of this article. This study was supported by the 345 Talent Project of Shengjing Hospital.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2025.1642866/full#supplementary-material
Abbreviations
LUAD, Lung adenocarcinoma; PTL, Parthenolide; TCM, Traditional Chinese medicine; BCA, Bicinchoninic Acid Assay; SDC, sodium deoxycholate; HPLC, High Performance Liquid Chromatography; PCA, principal component analysis; OPLSDA, orthogonal partial least squares discriminant analysis; PLSDA, partial least squares discriminant analysis.
References
1. Bade BC and Dela Cruz CS. Lung cancer 2020: epidemiology, etiology, and prevention. Clin Chest Med. (2020) 41:1–24. doi: 10.1016/j.ccm.2019.10.001
2. Nasim F, Sabath BF, and Eapen GA. Lung cancer. Med Clin North Am. (2019) 103:463–73. doi: 10.1016/j.mcna.2018.12.006
3. Li Y, Yan B, and He S. Advances and challenges in the treatment of lung cancer. BioMed Pharmacother. (2023) 169:115891. doi: 10.1016/j.biopha.2023.115891
4. Zhang Y, Deng C, Zheng Q, Li Y, and Chen H. The grading system of lung adenocarcinoma: ever-evolving concepts. J Thorac Oncol. (2022) 17:e30. doi: 10.1016/j.jtho.2021.11.015
5. Hardardottir H, Jonsson S, Gunnarsson O, Hilmarsdottir B, Asmundsson J, Gudmundsdottir I, et al. Advances in lung cancer diagnosis and treatment - a review. Laeknabladid. (2022) 108:17–29. doi: 10.17992/lbl.2022.01.671
6. Hutchinson BD, Shroff GS, Truong MT, and Ko JP. Spectrum of lung adenocarcinoma. Semin Ultrasound CT MR. (2019) 40:255–64. doi: 10.1053/j.sult.2018.11.009
7. Denisenko TV, Budkevich IN, and Zhivotovsky B. Cell death-based treatment of lung adenocarcinoma. Cell Death Dis. (2018) 9:117. doi: 10.1038/s41419-017-0063-y
8. Zeng J, Ding X, Ding J, and Wang X. Histological transformation into SCLC: An important resistance mechanism of NSCLC upon immunotherapy. Front Immunol. (2023) 14:1275957. doi: 10.3389/fimmu.2023.1275957
9. Pascoe HM, Knipe HC, Pascoe D, and Heinze SB. The many faces of lung adenocarcinoma: A pictorial essay. J Med Imaging Radiat Oncol. (2018) 62:654–61. doi: 10.1111/1754-9485.12779
10. Zhang Y, Fu F, Zhang Q, Li L, Liu H, Deng C, et al. Evolutionary proteogenomic landscape from pre-invasive to invasive lung adenocarcinoma. Cell Rep Med. (2024) 5:101358. doi: 10.1016/j.xcrm.2023.101358
11. Ruiz-Cordero R and Devine WP. Targeted therapy and checkpoint immunotherapy in lung cancer. Surg Pathol Clin. (2020) 13:17–33. doi: 10.1016/j.path.2019.11.002
12. Wang Q, Li M, Yang M, Yang Y, Song F, Zhang W, et al. Analysis of immune-related signatures of lung adenocarcinoma identified two distinct subtypes: implications for immune checkpoint blockade therapy. Aging (Albany NY). (2020) 12:3312–39. doi: 10.18632/aging.102814
13. Kuhn E, Morbini P, Cancellieri A, Damiani S, Cavazza A, Comin CE, et al. Adenocarcinoma classification: patterns and prognosis. Pathologica. (2018) 110:5–11.
14. Zhang X, Qiu H, Li C, Cai P, and Qi F. The positive role of traditional Chinese medicine as an adjunctive therapy for cancer. Biosci Trends. (2021) 15:283–98. doi: 10.5582/bst.2021.01318
15. Li Z, Feiyue Z, and Gaofeng L. Traditional Chinese medicine and lung cancer–From theory to practice. BioMed Pharmacother. (2021) 137:111381. doi: 10.1016/j.biopha.2021.111381
16. Liu Y, Yang S, Wang K, Lu J, Bao X, Wang R, et al. Cellular senescence and cancer: Focusing on traditional Chinese medicine and natural products. Cell Prolif. (2020) 53:e12894. doi: 10.1111/cpr.12894
17. He J, Yin P, and Xu K. Effect and molecular mechanisms of traditional chinese medicine on tumor targeting tumor-associated macrophages. Drug Des Devel Ther. (2020) 14:907–19. doi: 10.2147/DDDT.S223646
18. Freund RRA, Gobrecht P, Fischer D, and Arndt HD. Advances in chemistry and bioactivity of parthenolide. Nat Prod Rep. (2020) 37:541–65. doi: 10.1039/C9NP00049F
19. Araujo TG, Vecchi L, Lima PMAP, Ferreira EA, Campos IM, Brandao DC, et al. Parthenolide and its analogues: A new potential strategy for the treatment of triple-negative breast tumors. Curr Med Chem. (2020) 27:6628–42. doi: 10.2174/0929867326666190816230121
20. Liu X and Wang X. Recent advances on the structural modification of parthenolide and its derivatives as anticancer agents. Chin J Nat Med. (2022) 20:814–29. doi: 10.1016/S1875-5364(22)60238-3
21. He Y, Yu Q, Ma X, Lv D, Wang H, Qiu W, et al. A metabolomics approach reveals metabolic disturbance of human cholangiocarcinoma cells after parthenolide treatment. J Ethnopharmacol. (2024) 328:118075. doi: 10.1016/j.jep.2024.118075
22. Hayashi S, Hatashita M, Hayashi A, Matsumoto H, Shioura H, Kitai R, et al. Thermosensitization by parthenolide in human lung adenocarcinoma A549 cells and p53- and hsp72-independent apoptosis induction via the nuclear factor-kappaB signal pathway. Int J Mol Med. (2008) 21:585–92. doi: 10.3892/ijmm.21.5.585
23. Song JM, Qian X, Upadhyayya P, Hong KH, and Kassie F. Dimethylaminoparthenolide, a water soluble parthenolide, suppresses lung tumorigenesis through down-regulating the STAT3 signaling pathway. Curr Cancer Drug Targets. (2014) 14:59–69. doi: 10.2174/15680096113136660108
24. Jelic MD, Mandic AD, Maricic SM, and Srdjenovic BU. Oxidative stress and its role in cancer. J Cancer Res Ther. (2021) 17:22–8. doi: 10.4103/jcrt.JCRT_862_16
25. Abulikemu A, Zhao X, Xu H, Li Y, Ma R, Yao Q, et al. Silica nanoparticles aggravated the metabolic associated fatty liver disease through disturbed amino acid and lipid metabolisms-mediated oxidative stress. Redox Biol. (2023) 59:102569. doi: 10.1016/j.redox.2022.102569
26. Ruocco C, Segala A, Valerio A, and Nisoli E. Essential amino acid formulations to prevent mitochondrial dysfunction and oxidative stress. Curr Opin Clin Nutr Metab Care. (2021) 24:88–95. doi: 10.1097/MCO.0000000000000704
27. Wang M and Li Q. Parthenolide could become a promising and stable drug with anti-inflammatory effects. Nat Prod Res. (2015) 29:1092–101. doi: 10.1080/14786419.2014.981541
28. Li N, Wang Y, Wang X, Sun N, and Gong YH. Pathway network of pyroptosis and its potential inhibitors in acute kidney injury. Pharmacol Res. (2022) 175:106033. doi: 10.1016/j.phrs.2021.106033
29. Yao Y, Song L, Zuo Z, Chen Z, Wang Y, Cai H, et al. Parthenolide attenuates hypoxia-induced pulmonary hypertension through inhibiting STAT3 signaling. Phytomedicine. (2024) 134:155976. doi: 10.1016/j.phymed.2024.155976
30. Liu L, Feng L, Gao J, Hu J, Li A, Zhu Y, et al. Parthenolide targets NLRP3 to treat inflammasome-related diseases. Int Immunopharmacol. (2023) 119:110229. doi: 10.1016/j.intimp.2023.110229
31. Hao QF, Wang BB, Zhang W, Qiu W, Liu QL, Li XM, et al. NF-κB inhibitor parthenolide promotes renal tubules albumin uptake in type 2 diabetic nephropathy. Chin Med Sci J. (2020) 35:31–42. doi: 10.24920/003573
32. Sztiller-Sikorska M and Czyz M. Parthenolide as cooperating agent for anti-cancer treatment of various Malignancies. Pharm (Basel). (2020) 13. doi: 10.3390/ph13080194
33. Liu X, Wang C, Li S, Qu L, Yin F, Lu D, et al. Parthenolide derivatives as PKM2 activators showing potential in colorectal cancer. J Med Chem. (2021) 64:17304–25. doi: 10.1021/acs.jmedchem.1c01380
34. Kim SL, Liu YC, Seo SY, Kim SH, Kim IH, Lee SO, et al. Parthenolide induces apoptosis in colitis-associated colon cancer, inhibiting NF-κB signaling. Oncol Lett. (2015) 9:2135–42. doi: 10.3892/ol.2015.3017
35. Cai Z, Gao L, Hu K, and Wang QM. Parthenolide enhances the metronomic chemotherapy effect of cyclophosphamide in lung cancer by inhibiting the NF-kB signaling pathway. World J Clin Oncol. (2024) 15:895–907. doi: 10.5306/wjco.v15.i7.895
36. Zhu S, Sun P, Bennett S, Charlesworth O, Tan R, Peng X, et al. The therapeutic effect and mechanism of parthenolide in skeletal disease, cancers, and cytokine storm. Front Pharmacol. (2023) 14:1111218. doi: 10.3389/fphar.2023.1111218
37. Kageyama S, Ii H, Taniguchi K, Kubota S, Yoshida T, Isono T, et al. Mechanisms of tumor growth inhibition by depletion of γ-glutamylcyclotransferase (GGCT): A novel molecular target for anticancer therapy. Int J Mol Sci. (2018) 19. doi: 10.3390/ijms19072054
38. Kageyama S, Hanada E, Ii H, Tomita K, Yoshiki T, Kawauchi A, et al. Gamma-glutamylcyclotransferase: A novel target molecule for cancer diagnosis and treatment. BioMed Res Int. (2015) 2015:345219. doi: 10.1155/2015/345219
39. He Z, Sun X, Wang S, Bai D, Zhao X, Han Y, et al. Ggct (γ-glutamyl cyclotransferase) plays an important role in erythrocyte antioxidant defense and red blood cell survival. Br J Haematol. (2021) 195:267–75. doi: 10.1111/bjh.17775
40. Lin Z, Xiong L, Zhou J, Wang J, Li Z, Hu H, et al. γ-glutamylcyclotransferase knockdown inhibits growth of lung cancer cells through G0/G1 phase arrest. Cancer Biother Radiopharm. (2015) 30:211–6. doi: 10.1089/cbr.2014.1807
41. Ii H, Nohara Y, Yoshiya T, Masuda S, Tsuda S, Oishi S, et al. Identification of U83836E as a γ-glutamylcyclotransferase inhibitor that suppresses MCF7 breast cancer xenograft growth. Biochem Biophys Res Commun. (2021) 549:128–34. doi: 10.1016/j.bbrc.2021.02.103
42. Hanada E, Kageyama S, Murai R, Kubota S, Ii H, Nakata S, et al. Pro-GA, a novel inhibitor of γ-glutamylcyclotransferase, suppresses human bladder cancer cell growth. Anticancer Res. (2019) 39:1893–8. doi: 10.21873/anticanres.13297
43. Zhong X, Ni J, Jia Z, Yan H, Zhang Y, Liu Y, et al. CBX3 is associated with metastasis and glutathione/glycosphingolipid metabolism in colon adenocarcinoma. J Gastrointest Oncol. (2022) 13:246–55. doi: 10.21037/jgo-22-97
44. Li HN, Zhang HM, Li XR, Wang J, Xu T, Li SY, et al. MiR-205-5p/GGCT attenuates growth and metastasis of papillary thyroid cancer by regulating CD44. Endocrinology. (2022) 163. doi: 10.1210/endocr/bqac022
45. Taniguchi K, Ii H, Kageyama S, Takagi H, Chano T, Kawauchi A, et al. Depletion of gamma-glutamylcyclotransferase inhibits cancer cell growth by activating the AMPK-FOXO3a-p21 axis. Biochem Biophys Res Commun. (2019) 517:238–43. doi: 10.1016/j.bbrc.2019.07.049
46. Shen SH, Chen SF, Guo JF, and Wang ZX. The GGCT and REST positive feedback loop promotes tumor growth in Glioma. Transl Oncol. (2024) 49:102083. doi: 10.1016/j.tranon.2024.102083
47. Yang KY, Wu CR, Zheng MZ, Tang RT, Li XZ, Chen LX, et al. Physapubescin I from husk tomato suppresses SW1990 cancer cell growth by targeting kidney-type glutaminase. Bioorg Chem. (2019) 92:103186. doi: 10.1016/j.bioorg.2019.103186
48. Wang B, Shen C, Li Y, Zhang T, Huang H, Ren J, et al. Oridonin overcomes the gemcitabine resistant PANC-1/Gem cells by regulating GST pi and LRP/1 ERK/JNK signalling. Onco Targets Ther. (2019) 12:5751–65. doi: 10.2147/OTT.S208924
49. Choi EO, Park C, Hwang HJ, Hong SH, Kim GY, Cho EJ, et al. Baicalein induces apoptosis via ROS-dependent activation of caspases in human bladder cancer 5637 cells. Int J Oncol. (2016) 49:1009–18. doi: 10.3892/ijo.2016.3606
Keywords: lung adenocarcinoma, parthenolide, amino acid metabolism, oxidative stress, tumor
Citation: Liu J, Liu Y, Li J and Shen S (2025) Effects of parthenolide on amino acid metabolism and oxidative stress in lung adenocarcinoma based on quantitative proteomic analysis, targeted amino acid metabolomics, network pharmacology, and experimental validation. Front. Oncol. 15:1642866. doi: 10.3389/fonc.2025.1642866
Received: 07 June 2025; Accepted: 31 July 2025;
Published: 01 September 2025.
Edited by:
Yinu Wang, Northwestern University, United StatesReviewed by:
Davinder Singh, The Ohio State University, United StatesWang Weifang, The Third Affiliated Hospital to Changchun University of Chinese Medicine, China
Copyright © 2025 Liu, Liu, Li and Shen. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Shuang Shen, U3MwMDQ1MDY3OEAxNjMuY29t
†These authors have contributed equally to this work