Hui-Juan Cao1†
Wei Zhou2†Xiao-Le Xian1†Shu-Jun Sun3Pei-Jie Ding1Chun-Yu Tian1Fu-Ling Tian1Chun-Hua Jiang1Ting-Ting Fu1Shu Zhao1
Jian-Ye Dai2,4,5*- 1Traditional Chinese Medicine College, North China University of Science and Technology, Tangshan, China
- 2School of Pharmacy, Lanzhou University, Lanzhou, China
- 3School of Biology and Food Engineering, Fuyang Normal University, Fuyang, China
- 4State Key Laboratory of Natural Medicines, China Pharmaceutical University, Nanjing, China
- 5Collaborative Innovation Center for Northwestern Chinese Medicine, Lanzhou University, Lanzhou, China
Non-small cell lung cancer (NSCLC) is a worldwide disease with a high morbidity and mortality rate, which is most derived from its metastasis. Some studies show that the epithelial–mesenchymal transition (EMT) process promotes lung cancer cell migration and invasion, leading to NSCLC metastasis. Total flavonoid aglycones extract (TFAE) isolated from Scutellaria baicalensis was reported to inhibit tumor growth and induce apoptosis. In this study, we found that baicalein, wogonin, and oroxylin-A were the active compounds of TFAE. After reconstructing with these three compounds [baicalein (65.8%), wogonin (21.2%), and oroxylin-A (13.0%)], the reconstructed TFAE (reTFAE) inhibited the EMT process of A549 cells. Then, bioinformatic technology was employed to elucidate the potential pharmacodynamic mechanism network of reTFAE. We identified the relationship between reTFAE and PI3K/Akt signaling pathways, with TWIST1 as the key protein. LY294002, the inhibitor of the PI3K/Akt signaling pathway, and knock-down TWIST1 could significantly enhance the efficacy of reTFAE, with increasing expression of epithelial markers and decreasing expression of mesenchymal markers in A549 cells at the same time. Furthermore, stable isotope dimethyl-labeled proteomics technology was conducted to complement the follow-up mechanism that the EMT-inhibition process may be realized through the glycolysis pathway. In conclusion, we claim that TWIST1-targeted flavonoids could provide a new strategy to inhibit EMT progress for the treatment of NSCLC.
Introduction
Lung cancer is a worldwide disease with a high morbidity and mortality rate. According to the World Health Organization (2020), there were 2.207 million new lung cancer patients worldwide in 2020, only less than breast cancer, and lung cancer ranked first in the global causes of cancer-related death (World Health Statistics, 2020; Thai et al., 2021). Among the histological subtypes of lung cancer collectively, non-small cell lung cancer (NSCLC) accounts for 80%–85%, and most lung cancer patients are diagnosed at advanced stages, so traditional chemotherapy and radiotherapy have limited efficacy (Siegel et al., 2017). NSCLC microenvironment has limited nutrition, so cancer cells often obtain the energy and substances required for proliferation and cell growth through the “Warburg Effect” glucose metabolism changes (Gatenby and Gillies, 2004; Lopez-Lazaro, 2008), that is, cancer cells tend to produce energy through glycolysis rather than oxidative phosphorylation (Heiden et al., 2009; Muz et al., 2015; Eales et al., 2016). However, the metabolic change leads to further tumor cell growth, aggravation of epithelial–mesenchymal transition (EMT) process, and resistance to treatment. During the treatment of NSCLC, 40% of patients have metastases (Dwyer-Nield et al., 2017), and the main biological process is EMT (Tulchinsky et al., 2018). In the EMT process, cells gradually become spindle-shaped, losing the unique epithelial characteristics, intercellular adhesion, and cytoskeletal structure of epithelial tissue, and possessing polarity, individual migration ability, and invasion ability (Nieto et al., 2016). Cellular EMT process is related to TGFR, Wnt/β-Catenin, PI3K/Akt, Ras, MAPK, NF-κB pathways, and the transcription factors involved are TWIST (TWIST1 and TWIST2), Slug, ZEB (ZEB1 and ZEB2), and Snail (Shaw and Solomon, 2011; Otsuki et al., 2018). EMT transcription factors promote the development of drug resistance (Staalduinen et al., 2018; Liu et al., 2020). In the process of EMT, the activated PI3K/Akt pathway promotes glucose uptake and glycolysis (Semenza, 2010; Tennant et al., 2010), and phosphorylated Akt can phosphorylate TWIST1, then up-regulates the TGFβ-Smad pathway, which further promotes the EMT process (Xue et al., 2012). TWIST1 was known to bind to the E-box sequence of E-cadherin and inhibited the expression of E-cadherin to induce EMT in a variety of tumors (Ren et al., 2016; Zhu et al., 2016; Ding et al., 2019). Inhibition of TWIST1 was known to induce growth inhibition and apoptosis of EGFR-mutant NSCLC cells (Yochum et al., 2019).
Some herbals and their extracts have been reported to have therapeutic effects on cancer, including Scutellaria baicalensis (Sato et al., 2013; Bie et al., 2017; Zhang et al., 2018; Hsiao et al., 2019), whose main components are flavonoids. We previously reported that Total Flavonoid Aglycones Extract (TFAE) isolated from S. baicalensis inhibited tumor growth and induced apoptosis, in which the main components, baicalein, wogonin, and oroxylin-A accounted for 41.4%, 13.3%, and 8.2%, respectively (Wang et al., 2016). Several studies have investigated the treatment of lung cancer with S. baicalensis components, including wogonin, baicalin, and oroxylin-A (Wei et al., 2016; Yu et al., 2017; Wang et al., 2021). This inspired us to reconstruct a compound combination with baicalein, wogonin, and oroxylin-A.
In this study, we reconstructed TFAE with baicalein (65.8%), wogonin (21.2%), and oroxylin-A (13.0%) and named it as reTFAE (Figure 1A), and found that the reTFAE could significantly inhibit EMT of A549 cells. We constructed the reTFAE-PI3K/AKT-TWIST1 correlation in inhibiting EMT progress, and revealed that this phenomenon may be related to glycolysis pathway through stable isotope dimethyl-labeled proteomics technology. Moreover, we proved that TWIST1 may be a potential target to inhibit the EMT process of lung cancer, and hope to explore a new combined medication strategy through the compatibility study of active compounds in herbal extracts.
FIGURE 1. reTFAE inhibited EMT processes of A549 cells. (A) The structure and proportion of the components of reTFAE. (B) The invasion-inhibiting ability of reTFAE on A549 cells (n = 3). (C,D) Effects of different concentrations of reTFAE on the expression of EMT-related proteins (n = 3). (E,F) reTFAE inhibited A549 cells’ migration (n = 3). Statistical differences were determined by a two-sided Student’s t-test; Compared with 0 μg/μl TFAE. ***p < 0.001; **p < 0.01; *p < 0.05.
Materials and Methods
Cell Culture
The A549 cells were obtained from the National Collection of Authenticated Cell Culture and were cultured at 37°C under 5% CO2 in Dulbecco’s modified Eagle’s medium (DMEM; Sigma, D6429) supplemented with 10% (vol./vol.) fetal bovine serum (FBS; Gibco, 10270-106) and 1% (vol./vol.) penicillin-streptomycin (PS; Gibco, 15140-122). Culture medium was refreshed and the cells were passaged when the confluency was 80%.
Preparation of Cisplatin, LY294002, and reTFAE
As shown in Table1, 4 mg/ml cisplatin (Solarbio, D8810) stock solution, 20 mM LY294002 (AbMole, M1925) stock solution, and 128 mg/ml reTFAE solution were prepared by DMSO and were used for the later experiments. For reTFAE, the three main components and their mass percentages were baicalein (65.8%, Sichuan Weikeqi Biological, wkq-00289), wogonin (21.2%, Sichuan Weikeqi Biological, wkq-00245), and oroxylin-A (13.0%, Sichuan Weikeqi Biological, wkq-00180). We weighed the amount of each component according to their mass percentage, then mixed the three compounds, and dissolved them to the desired volume with DMSO.
TABLE 1. Compounds.
Wound Healing Assay and Cell Invasion Assay
For wound healing assay: 5.0 × 106 A549 cells were plated in six-well plates. After 7 h, the cells formed a single layer with nearly 100% confluency. Confluent cell layers were scratched using 200-μl tips to generate wounds with approximately 800 µm width. Cells were washed three times with PBS to remove detached cells and debris. Cells were incubated with complete medium with or without 8 μg/ml reTFAE, 4 μg/ml Cisplatin, or 20 μM LY294002 for 24 h. At the 0th and 24th hour after scratching, three fields of view for each well were selected randomly to take pictures. ImageJ software was used to scan the scratch areas and to calculate the wound healing area. Wound healing rate = (wound area at 0th hour − wound area at 24th hour)/wound area at 0th hour × 100%.
For cell invasion assay (Bauer et al., 2005; Li et al., 2020; Nishida et al., 2020): the transwell containing polycarbonate membranes with 8 μm pore size was precoated with 50 µl of Matrigel (1.0 mg/ml) and incubated at 37°C overnight. A total of 1.0 × 106 A549 cells and reTFAE (2, 4, and 8 μg/ml) in serum-free DMEM medium were plated in the upper chamber, and 150 µl of DMEM supplemented with 10% FBS in the lower chamber. After 24-h incubation, cells inside the upper chamber might pass through the matrix layer and transfer to the lower membrane surface. The upper chamber was taken out and washed two times by wash buffer. Cells in the upper surface of the membrane were completely removed with cotton swabs, and cells migrated to the lower membrane surface were fixed with 4% paraformaldehyde for 5 min. Then, the upper chamber was transferred to test board containing 100 µl Cell Dissociation Solution/Calcein-AM. After 1-h incubation, fluorescence value was detected at 485-nm excitation wavelength and 520-nm emission wavelength. Cell invasion rate = (number of invaded cells in experimental wells/number of invaded cells in control wells) ×100%.
Western Blot
The A549 cells under different conditions were collected after washing three times by PBS and centrifugation at 8,600 rcf for 3 min. The cells were lysed in RIPA Lysis Buffer (CWBIO, CW2333S) containing EDTA-free cOmplete (Roche, 4693132001) with sonication on ice. The cell lysates were collected by centrifugation (2.15 × 104 rcf, 20 min) at 4°C. The protein concentration was detected by Pierce™ BCA Protein Assay Kit (Thermo Fisher Scientific, 23225) and adjusted to 2 mg/ml. Twenty micrograms of the lysates was separated by 10% SDS-PAGE and then transferred onto a PVDF membrane in an ice bath at 300 mA for 1 h. The membrane was blotted with 5% nonfat milk in TBST at RT for 1 h and then incubated with the following antibodies at 4°C overnight: TWIST1 (1:1,000, Cell Signaling Technology, 46702S), E-cadherin (1:1,000, Cell Signaling Technology, 3,195), Vimentin (1:1,000, Cell Signaling Technology, 5,741), N-cadherin (1:1,000, Santa, Sc-59987), AKT (1:1,000, Cell Signaling Technology, 4,685), p-AKT (Ser473) (1:2,000, Cell Signaling Technology, 4,060), β-Actin (1:5,000, Proteintech, 20536-1-AP), and GAPDH (1:5,000, Proteintech, 10494-1-AP). Then, the membrane was incubated with secondary antibodies at RT for 1 h. The secondary antibodies are as follows: Goat anti-Rabbit IgG (1:5,000, Proteintech, SA00001-2) and Goat anti-mouse IgG (1:5,000, Proteintech, SA00001-1). Finally, the immunoreactive bands were visualized by the Clarity™ Western ECL Substrate (BIO-RAD, 1705061) using the Tanon4600 Automatic chemiluminescence image analysis system (Tanon, Shanghai, China). ImageJ software was used to analyze the thickness of protein bands. Relative expression of target protein = target protein band/β-Actin band. AKT phosphorylation level = p-AKT band/AKT band.
Transfection With siRNA
A total of 6.0 × 105 A549 cells were cultured in medium without PS and were plated in six-well plates. After 7 h, the cells with nearly 70% confluency were transfected with TWIST1 siRNA. The siRNA sequences were as follows: TWIST1 (sense: 5′-CAAGAUUCAGACCCUCAAGTT-3′, antisense: 5′-CUUGAGGGUCUGAAUCUUGTT-3′); GAPDH (sense: 5′-AATGGGCAGCCGTTAGGAAA-3′, antisense: 5′-TGAAGGGGTCATTGATGGCA-3′). Negative control siRNA, which were confirmed not to interact with any mRNA sequence else, was used to balance siRNA where necessary. siRNA sequences were designed and synthesized by GenePharma (Shanghai, China). The interference was performed based on the protocol of Lipofectamine® 3,000 (Thermo Fisher Scientific, L3000150). In detail, for each well, 5 μl of 15-μm siRNA primers dissolved in 125 μl of opti-MEM (Gibco) was mixed with 7.5 μl of Lipofectamine® 3,000 (Thermo Fisher Scientific, L3000150) dissolved in 125 μl of opti-MEM. After 20-min mixture, replace with 1 ml of new medium and add 250 μl of mixture dropwise to cells. After 8-h incubation, add 1 ml of medium and culture for another 48 h. The cells were collected separately and the knockdown efficiency was examined by qRT-PCR.
Quantitative Real-Time PCR
Total RNA was extracted and purified with Eastep® Super Total RNA Extraction Kit (Promega, LS1040) and 1 μg of total RNA was used as the cDNA synthesis template for reverse transcription with the RevertAid First Strand cDNA Synthesis Kit (Thermo Scientific, k1622). The reverse transcription was performed in 20 μl volume and the conditions were as follows: 42°C for 30 min; 85°C for 10 min. qRT-PCR was performed by MX3000P (Agilent) using UltraSYBR Mixture (Low ROX) (CWBIO, CW2601), and reaction was performed in 20 μl volume and the conditions were as follows: 95°C for 3 min; 95°C for 15 s; 55°C for 30 s, 72°C for 30 s; 40 cycles. The primers were as follows: TWIST1 (forward: 5′-TCGGACAAGCTGAGAGCAAGATTCA-3′, reverse: 5′-TCCATCCTCCAGACCGAGAAGG-3′), GAPDH (forward: 5′-CATGAGAAGTATGACAACAGCCT-3′, reverse: 5′-AGTCCTTCCACGATACCAAAGT-3′). CT values were used to evaluate the relative mRNA expression by 2−ΔΔCt method and β-Actin served as an internal control. The primers were designed and synthesized by GenePharma (Shanghai, China).
Protein In-Solution Digestion and Dimethyl Labeling
A549 cells were grown to 80% confluence in a 15-cm dish. In our previous experiments, 64 μg/ml reTFAE had a more stable and stronger EMT-inhibition effect than 8 μg/ml reTFAE. To obtain a more accurate and reliable proteomic result, A549 cells were treated with 64 μg/ml reTFAE for 24 h. The cells were collected after washing three times by PBS and centrifugation at 8,600 rcf for 3 min. The cells were lysed in 0.1% Triton X-100 (Sigma-Aldrich, T8787)-100 mM TEAB (Sigma-Aldrich, T7408) containing EDTA-free cOmplete with sonication on ice. The cell lysates were collected and the protein concentration was detected by Pierce™ BCA Protein Assay Kit. Ten-microliter cell lysates (3 mg/ml) were reacted with 30 µl of 8 M urea (Sigma, BCBZ1744) and 2 µl of 200 mM DTT (Sigma-Aldrich, 43815) at 65°C for 15 min in the dark. Then, 2 µl of 400 mM iodoacetamide (Sigma, I1149) was added to react in the dark for 30 min at 35°C. Two microliters of 200 mM DTT was added to react with the remaining iodoacetamide in the dark for 15 min at 65°C. One hundred microliters of 100 mM TEAB, 2 µl of 0.2 μg/μl trypsin (Promega, V528A), and 1.5 µl of 100 mM CaCl2 were added and the trypsin digestion was performed at 37°C overnight. For dimethyl labeling: for the control (“light”) group, the peptides were reacted with 6 µl of 4% CH2O (Sigma-Aldrich, F1635) and 6 µl of 0.6 M NaBH3CN (Sigma-Aldrich, 42077); for the reTFAE (“heavy”) group, the peptides were reacted with 6 µl of 4% 13CD2O (Sigma-Aldrich, 596,388) and 6 µl of 0.6 M NaBD3CN (Sigma-Aldrich,190020). The reaction was incubated at 22°C for 1 h. The reaction was quenched by 24 µl of 1% ammonia (Sigma-Aldrich, 221228) and 12 µl of formic acid (Sigma-Aldrich, F0507). The “light” and “heavy” samples were combined and the desalination was performed by Pierce™ C18 Tips (Thermo Scientific,87782). LC-MS/MS analysis.
Samples were analyzed by LC-MS/MS on Q Exactive series Orbitrap mass spectrometers (Thermo Fisher Scientific). In positive-ion mode, full-scan mass spectra were acquired over the m/z ratio from 350 to 1,800 using the Orbitrap mass analyzer with a mass resolution of 7,000. MS/MS fragmentation was performed in a data-dependent mode, of which the TOP 20 most intense ions were selected for MS/MS analysis with a resolution of 17,500 under HCD’s collision mode. The isolation window was 2.0 m/z units, the default charge was 2+, the normalized collision energy was 28%, the maximum IT was 50 ms, and dynamic exclusion was 20.0 s. Under +57.0215 Da’s cysteine modification, the LC-MS/MS data were analyzed by ProLuCID. The isotopic modifications were set as static modifications on the N-terminal of a peptide and lysins, and 28.0313 and 34.0631 Da were for light and heavy labeling, respectively. The CIMAGE software was applied for quantitation, proteins with an average ratio (light/heavy) above 1.5 were selected for further KEGG analysis.
Network Pharmacology
In this study, the common mechanism of reTFAE in the treatment of lung cancer was studied based on network pharmacology. The active ingredients and related targets of reTFAE were integrated from TCMSP, BATMAN-TAM, STP, and Pubchem databases. The standard names of these targets were united by UniProt database. Targets of lung cancer were enriched through GeneCards, NCBI (Gene), Therapeutic Target Database, and DisGeNET (v7.0) databases. Then, the intersection targets of reTFAE and disease were obtained. The STRING network and the Cytoscape 3.7.2 were used to construct a protein–protein interaction (PPI) network, and the DAVID database was used to perform KEGG analysis. Then, Cytoscape 3.6.1 was used to build the “Ingredient-Target-Signal Pathway” network.
Statistical Analysis
Student’s t-test was used to compare experimental data. We analyzed the data in GraphPad Prism (GraphPad Software), using the unpaired, two-tailed t-test module. Statistical significance was considered when a p value was below 0.05. *p < 0.05; **p < 0.01; ***p < 0.001. N.S., not significant.
Results
reTFAE Inhibited EMT Processes of A549 Cells
To explore reTFAE’s effects in the EMT process of lung cancer, A549 cells were treated with reTFAE (2, 4, and 8 μg/ml) for 24 h. Compared to 0 μg/ml reTFAE, the invasion rate of cells treated by reTFAE concentration-dependently decreased (Figure 1B). In addition, reTFAE dose-dependently inhibited the expression of TWIST1 and N-cadherin proteins, and reduced AKT phosphorylation level, while increasing the expression of E-cadherin (Figures 1C,D). Compared with the wound area at the 0th hour, the wound was dose-dependently healed by reTFAE at the 48th hour (Figures 1E,F), which means that reTFAE inhibited the migration of A549 cells.
reTFAE Inhibited the EMT Process via the PI3K/Akt Pathway
The above results indicated that reTFAE exerted a concentration-dependent inhibitory effect on the EMT process of A549 cells. To explore the connection between reTFAE and lung cancer, network pharmacology analysis was performed. Three active ingredients and 65 related targets of reTFAE were integrated from TCMSP, BATMAN-TAM, STP, and Pubchem databases. The targets of lung cancer were enriched through GeneCards, NCBI (Gene), Therapeutic Target Database, and DisGeNET (v7.0) databases. Venn diagram showed the intersection of 62 reTFAE and lung cancer targets (Figure 2A). STRING network and Cytoscape 3.7.2 were used to construct the PPI network (Figure 2B); AKT1 was the central protein. DAVID database was used to perform KEGG analysis, reTFAE regulated pathways with p < 0.01 were as follows: TNF signaling pathway, p53 signaling pathway, PI3K/Akt signaling pathway, VEGF signaling pathway, HIF-1 signaling pathway, and MAPK signaling pathway (Figure 2C). After literature research, the PI3K/Akt signaling pathway attracted our attention, which could significantly activate tumor metastasis and directly regulate the expression of TWIST1, one of the key markers in EMT. To investigate whether reTFAE affected the EMT process of A549 cells by acting on the PI3K/Akt signaling pathway, we incubated A549 cells with the PI3K/Akt pathway inhibitor, LY294002. In this part, cisplatin was used as a positive control.
FIGURE 2. reTFAE inhibited the PI3K/Akt pathway to inhibit EMT processes of A549 cells. (A) Intersection protein targets of reTFAE and lung cancer. reTFAE targets and lung cancer targets shared 62 proteins. (B) Protein–protein interaction (PPI) network of 62 proteins. AKT1, TP53, VEGFA, and JUN showed more connections with other proteins. (C) KEGG pathway analysis of 62 proteins. p < 0.01 was the cutoff for displaying the pathways. (D,E) Both reTFAE and PI3K/AKT pathway inhibitors inhibited the migration of A549 cells (n = 3). (F) Both reTFAE and PI3K/AKT pathway inhibitors inhibited the invasion of A549 cells (n = 3). (G,H) The influence of PI3K/Akt-suppression and reTFAE on the expression of EMT-related proteins (n = 3). Statistical differences were determined by a two-sided Student’s t-test. Compared with control, N.S. not significant, ***p < 0.001; **p < 0.01; *p < 0.05.
Compared with normal cells, inhibition of the PI3K/Akt pathway showed poor migration and invasion capabilities (Figures 2D–F), and increased expression of E-cadherin and decreased expression of Vimentin, N-cadherin and TWIST1, and reduced AKT phosphorylation level (Figures 2G,H), which were related with EMT progress. Therefore, we believed that reTFAE inhibited the EMT process of A549 cells by inhibiting the PI3K/Akt signaling pathway.
reTFAE Inhibited TWIST1 in the EMT Process of A549 Cells
To identify the key target of reTFAE, Cytoscape 3.6.1 was used to build the “Ingredient-Target- Pathway” network. Among proteins involved in the network, TWIST1 attracted our attention, which was reported to be directly regulated by PI3K-Akt and connected with both baicalein and wogonin (Xue et al., 2012; Yao et al., 2017; Yuan et al., 2020; Zeng et al., 2020) (Figure 3A). To investigate whether TWIST1 played an important role in reTFAE’s effect on the EMT process, we applied the small interfering RNA (siRNA) method to knock down the TWIST1 gene in A549 cells.
FIGURE 3. reTFAE inhibited the EMT process of A549 cells via inhibiting TWIST1. (A) Relationships of reTFAE components–KEGG pathways–Protein targets. (B) Relative TWIST1 mRNA level after siRNA interference (n = 3). (C). The inhibitions of siTWIST1 and reTFAE on A549 cell invasion rate (n = 3). (D,E) The inhibitions of siTWIST1 and reTFAE on wound healing (n = 3). (F,G) The influence of siTWIST1 and reTFAE on the expression of EMT-related proteins (n = 3). Negative control siRNA was confirmed not to interact with any mRNA sequence else, and was used to balance siRNA where necessary. Statistical differences were determined by a two-sided Student’s t-test; Compared with control, N.S., not significant, ***p < 0.001; **p < 0.01; *p < 0.05.
Compared with normal cells, TWIST1-knockdown A549 cells showed significant lower mRNA level (Figure 3B), and reTFAE-treated/siTWIST1/reTFAE and siTWIST1-treated cells showed poor migration and invasion capabilities, increased expression of epithelial marker (E-cadherin), and decreased expression of mesenchymal markers (Vimentin and N-cadherin), EMT-related transcription factor (TWIST1), and protein in the PI3K/Akt pathway (p-AKT) (Figures 3C–E). siTWIST1 showed similar effects to reTFAE, and the synergistic application of reTFAE and siTWIST1 showed a more obvious effect on EMT-related protein expression (Figures 3F,G). Therefore, we claimed that reTFAE inhibited the EMT process of A549 cells via inhibiting TWIST1.
The Effect of reTFAE on Protein Expression of A549 Cells
As above, we found that reTFAE inhibited the PI3K/Akt pathway and TWIST, yet how reTFAE led to the decline of EMT is not clear. To explain a further mechanism, we collected A549 cells treated with reTFAE for 24 h and performed a stable-isotope dimethyl labeling proteomics experiment. The proteins of A549 cells treated with DMSO or reTFAE were marked with “light” or “heavy”, respectively, and then mixed, desalted, and analyzed by LC-MS/MS (Figure 4A). Proteins with Average Ratio Light/Heavy (L/H) > 1.5 were filtered out, and 158 proteins downregulated by reTFAE appeared at least twice in repeated experiments (Figure 4B). KEGG pathway analysis was performed on the 158 proteins, and reTFAE downregulated pathways with p < 0.01 were as follows: Carbon metabolism, Glycolysis/Gluconeogenesis, Biosynthesis of amino acid, and Pentose phosphate pathway (Figure 4C). These four pathways contain 11, 9, 7, and four proteins, respectively (Table 1). In short, reTFAE interferes with the sugar metabolism of A549 cells.
FIGURE 4. reTFAE interfered with the glucose metabolism of A549 cells. (A) Overall scheme of In-solution dimethyl labeling experiment, the control or reTFAE group were labeled with “light” or “heavy”, respectively. (B) Venn diagram showing the 158 proteins with an averaged ratio >1.5 appearing twice in triple experiments. (C) KEGG pathway analysis of the 158 proteins. p < 0.01 was the cutoff for displaying the pathways.
Discussion
As reported, TFAE of S. baicalensis has an inhibitory effect on NSCLC (Wang et al., 2016). In this study, we configured reTFAE and found that it inhibited the EMT process of A549 cells by affecting the PI3K/Akt pathway and TWIST1. As “knockdown TWIST1” and “reTFAE incubation” showed similar EMT-inhibition tendency, we believe that reTFAE inhibits the PI3K/Akt pathway and TWIST1, thereby inhibiting the EMT process of A549 cells. Actually, some evidence showed that the TWIST1-related axis may participate in the EMT process by activating the Wnt/β-catenin signaling pathway, thereby accelerating the process of lung cancer (Pan et al., 2020). TWIST1 promotes EMT and metastasis through serine phosphorylation of p38, c-Jun N-terminal kinase (JNK), and Erk1/2 (Jun et al., 2011). As our research provides direct evidence, TWIST1-targeted flavonoid may provide a new strategy to inhibit EMT progress.
Furthermore, stable-isotope dimethyl labeling proteomics was employed to detail the pharmacodynamic network of lung cancer cells treated with reTFAE. In this research, reTFAE downregulated the carbon metabolism, glycolysis/gluconeogenesis, biosynthesis of amino acid, pentose phosphate pathway, etc. Actually, tumor cells obtain energy mainly through the process of glycolysis, which promotes the EMT process of tumor cells (Liberti and Locasale, 2016); when nutrients were depleted, cancer cells tend to obtain metabolic materials through the pentose phosphate pathway (Georgakopoulos-Soares et al., 2020). Studies had shown that TWIST1 can promote the glycolysis process (Yang et al., 2015; Lovisa et al., 2020; Wang et al., 2020). In this study, ALDOA and GPI were enriched in both glycolysis and pentose phosphate pathways, while PKM and LDHA were enriched in glycolysis (Table2). The relationships between the above genes and the EMT process had been reported: PKM expresses pyruvate kinase and catalyzes the transfer of phosphorylation groups from phosphoenolpyruvate to ADP to generate ATP and pyruvate (Dombrauckas et al., 2005). PKM is expressed in fetal tissues and cancers, and participates in the EMT process of human colon cancer cells (Wong et al., 2015). PKM interacts with TGFβ-induced factor homeobox 2 (TGIF2) to inhibit the transcription of E-Cadherin (Hamabe et al., 2014). Glucose-6-phosphate isomerase edited by GPI catalyzes the conversion of glucose-6-phosphate to fructose-6-phosphate. GPI can also be produced by cancer cells to promote the EMT process (Huang et al., 2019). GPI promotes the EMT process of breast cancer cells by inhibiting miR-200 and inducing ZEB1/2 (Teng et al., 2014). Silencing GPI promotes the transition of human lung fibroblasts from the mesenchymal to the epithelial state (Funasaka et al., 2007). ALDOA expresses fructose-bisphosphate aldolase A, which is involved in glycolysis and gluconeogenesis, and is also overexpressed in cancer (Chang et al., 2018). ALDOA is highly expressed in lung squamous cell carcinoma (LSCC) and depletion of ALDOA in lung squamous carcinoma cells reduces cell motility capabilities (Du et al., 2014). Overexpression of ALDOA in colon cancer cells leads to the EMT progress (Ye et al., 2018). In pancreatic cancer and bladder cancer cells, ALDOA-silencing increases E-Cadherin and decreases N-Cadherin expression (Ji et al., 2016; Li et al., 2019). Lactate dehydrogenase-A (LDHA) enzyme converts pyruvate into lactic acid. LDHA promotes the EMT process, but the specific mechanism is not clear (Hou et al., 2019; Tirpe et al., 2019). Our hypothesis is that reTFAE had a relationship with glycolysis or PI3K/Akt-related proteins to exert an EMT-inhibiting effect, yet it needs further verification.
TABLE 2. Pathways interfered by reTFAE, and proteins enriched in each pathway.
However, the details of the synergy between the three molecules need to be further studied. Therefore, next, we will explore the detailed pharmacodynamic mechanism of the three molecules.
To conclude, in this work, we first determined that reTFAE composed of three compounds in a definite proportion had the activity of inhibiting the invasion and metastasis of lung cancer A549 cells. Then, combining with network pharmacology and molecular biology technology, the PI3K-AKT-TWIST1 axis was inhibited by reTFAE to inhibit EMT progress. Furthermore, chemoproteomics was used to elucidate the changes of downstream protein networks after TFEA inhibited the PI3K-AKT-TWIST1 axis, and it was found that the glycolysis pathway may play an important role. For clinical drug development, this study shows that multiple components can be used to treat lung cancer. Herbal medicine containing these three components can be used in the treatment of lung cancer. At the same time, the expression of EMT-related indicators in the body can be used as quality monitoring during the development of lung cancer treatment drugs. In summary, reTFAE can inhibit the PI3K/Akt pathway with the core of TWIST1, thereby inhibiting the glycolytic pathway to suppress EMT in A549 cells (Figure 5). TWIST1-targeted flavonoid provided a new strategy to inhibit EMT progress for the treatment of NSCLC.
FIGURE 5. EMT process induced by reTFAE. Baicalein and wogonin inhibit the PI3K/Akt pathway, thereby inhibiting TWIST1, which inhibits the glycolysis pathway and inhibits the EMT process of A549 cells.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: ProteomeXchange Consortium via the PRIDE partner repository with the dataset identifier PXD030097.
Author Contributions
J-YD and H-JC designed all the experiments. H-JC, WZ, and X-LX conducted most of the experiments. WZ, S-JS, and C-YT designed siRNA primers. WZ, P-JD, and F-LT participated in the analysis of network pharmacology. WZ, C-HJ, T-TF, and SZ participated in the Western blot experiment. J-YD and WZ performed LC-MS/MS analysis. All the authors participated in manuscript writing.
Funding
This work is financially supported by the Science and Technology Project of Hebei Education Department (QN2018060), the National Key Research and Development Program of China (2018YFC1706300 and 2018YFC17063005), the Science and Technology Planning Project of Gansu Province (20JR10RA586), the Fundamental Research Funds for the Central Universities (lzujbky-2021-kb40), and the Project for Longyuan Youth Innovation and Entrepreneurship Talent.
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.
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.
Acknowledgments
We thank North China University of Science and Technology and Lanzhou University for experimental platform. We thank the BioRxiv platform for allowing the article to be published in preprinted version.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2021.821485/full#supplementary-material
References
Bauer, T. W., Fan, F., Liu, W., Johnson, M., Parikh, N. U., Parry, G. C., et al. (2005). Insulinlike Growth Factor-I-Mediated Migration and Invasion of Human colon Carcinoma Cells Requires Activation of C-Met and Urokinase Plasminogen Activator Receptor. Ann. Surg. 241, 748–8. doi:10.1097/01.sla.0000160699.59061.92
Bie, B., Sun, J., Guo, Y., Li, J., Jiang, W., Yang, J., et al. (2017). Baicalein: A Review of its Anti-cancer Effects and Mechanisms in Hepatocellular Carcinoma. Biomed. Pharmacother. 93, 1285–1291. doi:10.1016/j.biopha.2017.07.068
Chang, Y. C., Yang, Y. C., Tien, C. P., Yang, C. J., and Hsiao, M. (2018). Roles of Aldolase Family Genes in Human Cancers and Diseases. Trends Endocrinol. Metab. 29, 549–559. doi:10.1016/j.tem.2018.05.003
Ding, X., Li, F., and Zhang, L. (2019). Knockdown of Delta-like 3 Restricts Lipopolysaccharide-Induced Inflammation, Migration and Invasion of A2058 Melanoma Cells via Blocking Twist1-Mediated Epithelial-Mesenchymal Transition. Life Sci. 226, 149–155. doi:10.1016/j.lfs.2019.04.024
Dombrauckas, J. D., Santarsiero, B. D., and Mesecar, A. D. (2005). Structural Basis for Tumor Pyruvate Kinase M2 Allosteric Regulation and Catalysis. Biochemistry 44, 9417–9429. doi:10.1021/bi0474923
Du, S., Guan, Z., Hao, L., Song, Y., Wang, L., Gong, L., et al. (2014). Fructose-Bisphosphate Aldolase A Is a Potential Metastasis-Associated Marker of Lung Squamous Cell Carcinoma and Promotes Lung Cell Tumorigenesis and Migration. PLoS One 9, e85804. doi:10.1371/journal.pone.0085804
Dwyer-Nield, L., Hickey, G. A., Friedman, M., Choo, K., Mcarthur, D. G., Tennis, M. A., et al. (2017). The Second-Generation PGI2 Analogue Treprostinil Fails to Chemoprevent Tumors in a Murine Lung Adenocarcinoma Model. Cancer Prev Res (Phila) 10, 671-679. doi:10.1158/1940-6207
Eales, K. L., Hollinshead, K. E., and Tennant, D. A. (2016). Hypoxia and Metabolic Adaptation of Cancer Cells. Oncogenesis 5, e190. doi:10.1038/oncsis.2015.50
Funasaka, T., Hu, H., Yanagawa, T., Hogan, V., and Raz, A. (2007). Down-regulation of Phosphoglucose Isomerase/autocrine Motility Factor Results in Mesenchymal-To-Epithelial Transition of Human Lung Fibrosarcoma Cells. Cancer Res. 67, 4236–4243. doi:10.1158/0008-5472.CAN-06-3935
Gatenby, R. A., and Gillies, R. J. (2004). Why Do Cancers Have High Aerobic Glycolysis? Nat. Rev. Cancer 4, 891–899. doi:10.1038/nrc1478
Georgakopoulos-Soares, I., Chartoumpekis, D. V., Kyriazopoulou, V., and Zaravinos, A. (2020). EMT Factors and Metabolic Pathways in Cancer. Front. Oncol. 10, 499. doi:10.3389/fonc.2020.00499
Hamabe, A., Konno, M., Tanuma, N., Shima, H., Tsunekuni, K., Kawamoto, K., et al. (2014). Role of Pyruvate Kinase M2 in Transcriptional Regulation Leading to Epithelial-Mesenchymal Transition. Proc. Natl. Acad. Sci. U S A. 111, 15526–15531. doi:10.1073/pnas.1407717111
Hou, X. M., Yuan, S. Q., Zhao, D., Liu, X. J., and Wu, X. A. (2019). LDH-A Promotes Malignant Behavior via Activation of Epithelial-To-Mesenchymal Transition in Lung Adenocarcinoma. Biosci. Rep. 39, BSR20181476. doi:10.1042/BSR20181476
Hsiao, Y. H., Lin, C. W., Wang, P. H., Hsin, M. C., and Yang, S. F. (2019). The Potential of Chinese Herbal Medicines in the Treatment of Cervical Cancer. Integr. Cancer Ther. 18, 1534735419861693. doi:10.1177/1534735419861693
Huang, H. C., Wen, X. Z., Xue, H., Chen, R. S., Ji, J. F., and Xu, L. (2019). Phosphoglucose Isomerase Gene Expression as a Prognostic Biomarker of Gastric Cancer. Chin. J. Cancer Res. 31, 771–784. doi:10.21147/j.issn.1000-9604.2019.05.07
Ji, S., Zhang, B., Liu, J., Qin, Y., Liang, C., Shi, S., et al. (2016). ALDOA Functions as an Oncogene in the Highly Metastatic Pancreatic Cancer. Cancer Lett. 374, 127–135. doi:10.1016/j.canlet.2016.01.054
Jun, H., Jian, Z., Junjiang, F., Tao, H., Jun, Q., Li, W., et al. (2011). Phosphorylation of Serine 68 of Twist1 by MAPKs Stabilizes Twist1 Protein and Promotes Breast Cancer Cell Invasiveness. Cancer Res. 71, 3980–3990. doi:10.1158/0008-5472.CAN-10-2914
Li, J., Wang, F., Gao, H., Huang, S., Cai, F., and Sun, J. (2019). ALDOLASE A Regulates Invasion of Bladder Cancer Cells via E-Cadherin-EGFR Signaling. J. Cell. Biochem. 120, 13694–13705. doi:10.1002/jcb.28642
Li, M., Wei, L., Zhou, W., He, Z., Ran, S., and Liang, J. (2020). miR-200a Contributes to the Migration of BMSCs Induced by the Secretions of E. faecalis via FOXJ1/NFκB/MMPs axis. Stem Cell Res Ther 11, 317. doi:10.1186/s13287-020-01833-1
Liberti, M. V., and Locasale, J. W. (2016). Correction to: 'The Warburg Effect: How Does it Benefit Cancer Cells? Trends Biochem. Sci. 41, 287–218. doi:10.1016/j.tibs.2016.01.004
Liu, T., Zhao, X., Zheng, X., Zheng, Y., Dong, X., Zhao, N., et al. (2020). The EMT Transcription Factor, Twist1, as a Novel Therapeutic Target for Pulmonary Sarcomatoid Carcinomas. Int. J. Oncol. 56, 750–760. doi:10.3892/ijo.2020.4972
Lopez-Lazaro, M. (2008). The Warburg Effect: Why and How Do Cancer Cells Activate Glycolysis in the Presence of Oxygen? Acamc 8, 305–312. doi:10.2174/187152008783961932
Lovisa, S., Fletcher-Sananikone, E., Sugimoto, H., Hensel, J., Lahiri, S., Hertig, A., et al. (2020). Endothelial-to-mesenchymal Transition Compromises Vascular Integrity to Induce Myc-Mediated Metabolic Reprogramming in Kidney Fibrosis. Sci. Signal. 13, eaaz2597. doi:10.1126/scisignal.aaz2597
Muz, B., de la Puente, P., Azab, F., and Azab, A. K. (20152015). The Role of Hypoxia in Cancer Progression, Angiogenesis, Metastasis, and Resistance to Therapy. Hypoxia (Auckl) 3, 83–92. doi:10.2147/HP.S93413
Nieto, M. A., Huang, R. Y.-J., Jackson, R. A., and Thiery, J. P. (2016). EMT: 2016. Cell 166, 21–45. doi:10.1016/j.cell.2016.06.028
Nishida, K., Kuwano, Y., and Rokutan, K. (2020). The MicroRNA-23b/27b/24 Cluster Facilitates Colon Cancer Cell Migration by Targeting FOXP2. Cancers (Basel) 12. doi:10.3390/cancers12010174
Otsuki, Y., Saya, H., and Arima, Y. (2018). Prospects for New Lung Cancer Treatments that Target EMT Signaling. Dev. Dyn. 247, 462–472. doi:10.1002/dvdy.24596
Pan, J., Fang, S., Tian, H., Zhou, C., Zhao, X., Tian, H., et al. (2020). lncRNA JPX/miR-33a-5p/Twist1 axis Regulates Tumorigenesis and Metastasis of Lung Cancer by Activating Wnt/β-Catenin Signaling. Mol. Cancer 19, 9. doi:10.1186/s12943-020-1133-9
Ren, H., Du, P., Ge, Z., Jin, Y., Ding, D., Liu, X., et al. (2016). TWIST1 and BMI1 in Cancer Metastasis and Chemoresistance. J. Cancer 7, 1074–1080. doi:10.7150/jca.14031
Sato, D., Kondo, S., Yazawa, K., Mukudai, Y., Li, C., Kamatani, T., et al. (2013). The Potential Anticancer Activity of Extracts Derived from the Roots of Scutellaria Baicalensis on Human Oral Squamous Cell Carcinoma Cells. Mol. Clin. Oncol. 1, 105–111. doi:10.3892/mco.2012.14
Semenza, G. L. (2010). HIF-1: Upstream and Downstream of Cancer Metabolism. Curr. Opin. Genet. Dev. 20, 51–56. doi:10.1016/j.gde.2009.10.009
Shaw, A. T., and Solomon, B. (2011). Targeting Anaplastic Lymphoma Kinase in Lung Cancer. Clin. Cancer Res. 17, 2081–2086. doi:10.1158/1078-0432.CCR-10-1591
Siegel, R. L., Miller, K. D., and Jemal, A. (2017). Cancer Statistics, 2017. CA Cancer J. Clin. 67, 7–30. doi:10.3322/caac.21387
Teng, Y., Mei, Y., Hawthorn, L., and Cowell, J. K. (2014). WASF3 Regulates miR-200 Inactivation by ZEB1 through Suppression of KISS1 Leading to Increased Invasiveness in Breast Cancer Cells. Oncogene 33, 203–211. doi:10.1038/onc.2012.565
Tennant, D. A., Durán, R. V., and Gottlieb, E. (2010). Targeting Metabolic Transformation for Cancer Therapy. Nat. Rev. Cancer 10, 267–277. doi:10.1038/nrc2817
Thai, A. A., Solomon, B. J., Sequist, L. V., Gainor, J. F., and Heist, R. S. (2021). Lung Cancer. Lancet 398, 535–554. doi:10.1016/S0140-6736(21)00312-3
Tirpe, A. A., Gulei, D., Ciortea, S. M., Crivii, C., and Berindan-Neagoe, I. (2019). Hypoxia: Overview on Hypoxia-Mediated Mechanisms with a Focus on the Role of HIF Genes. Int. J. Mol. Sci. 20, 6140. doi:10.3390/ijms20246140
Tulchinsky, E., Demidov, O., Kriajevska, M., Barlev, N. A., and Imyanitov, E. (2018). EMT: A Mechanism for Escape from EGFR-Targeted Therapy in Lung Cancer. Biochim. Biophys. Acta Rev. Cancer 1871, 29–39. doi:10.1016/j.bbcan.2018.10.003
van Staalduinen, J., Baker, D., Ten Dijke, P., and van Dam, H. (2018). Epithelial-mesenchymal-transition-inducing Transcription Factors: New Targets for Tackling Chemoresistance in Cancer? Oncogene 37, 6195–6211. doi:10.1038/s41388-018-0378-x
Vander Heiden, M. G., Cantley, L. C., and Thompson, C. B. (2009). Understanding the Warburg Effect: The Metabolic Requirements of Cell Proliferation. Science 324, 1029–1033. doi:10.1126/science.1160809
Wang, L., Zhang, J., Shan, G., Liang, J., Jin, W., Li, Y., et al. (2021). Establishment of a Lung Cancer Discriminative Model Based on an Optimized Support Vector Machine Algorithm and Study of Key Targets of Wogonin in Lung Cancer. Front. Pharmacol. 12, 728937. doi:10.3389/fphar.2021.728937
Wang, X. X., Yin, G. Q., Zhang, Z. H., Rong, Z. H., Wang, Z. Y., Du, D. D., et al. (2020). TWIST1 Transcriptionally Regulates Glycolytic Genes to Promote the Warburg Metabolism in Pancreatic Cancer. Exp. Cell Res. 386, 111713. doi:10.1016/j.yexcr.2019.111713
Wang, Y., Cao, H. J., Sun, S. J., Dai, J. Y., Fang, J. W., Li, Q. H., et al. (2016). Total Flavonoid Aglycones Extract in Radix Scutellariae Inhibits Lung Carcinoma and Lung Metastasis by Affecting Cell Cycle and DNA Synthesis. J. Ethnopharmacol. 194, 269–279. doi:10.1016/j.jep.2016.07.052
Wei, L., Yao, Y., Zhao, K., Huang, Y., Zhou, Y., Zhao, L., et al. (2016). Oroxylin A Inhibits Invasion and Migration through Suppressing ERK/GSK-3β Signaling in Snail-Expressing Non-small-cell Lung Cancer Cells. Mol. Carcinog 55, 2121–2134. doi:10.1002/mc.22456
Wong, N., Ojo, D., Yan, J., and Tang, D. (2015). PKM2 Contributes to Cancer Metabolism. Cancer Lett. 356, 184–191. doi:10.1016/j.canlet.2014.01.031
World Health Statistics (2020). International Agency for Research on Cancer. Available at: https://www.who.int/data/collections.
Xue, G., Restuccia, D. F., Lan, Q., Hynx, D., Dirnhofer, S., Hess, D., et al. (2012). Akt/PKB-mediated Phosphorylation of Twist1 Promotes Tumor Metastasis via Mediating Cross-Talk between PI3K/Akt and TGF-β Signaling Axes. Cancer Discov. 2, 248–259. doi:10.1158/2159-8290.CD-11-0270
Yang, L., Hou, Y., Yuan, J., Tang, S., Zhang, H., Zhu, Q., et al. (2015). Twist Promotes Reprogramming of Glucose Metabolism in Breast Cancer Cells through PI3K/AKT and P53 Signaling Pathways. Oncotarget 6, 25755–25769. doi:10.18632/oncotarget.4697
Yao, Y., Zhao, K., Yu, Z., Ren, H., Zhao, L., Li, Z., et al. (2017). Wogonoside Inhibits Invasion and Migration through Suppressing TRAF2/4 Expression in Breast Cancer. J. Exp. Clin. Cancer Res. 36, 103. doi:10.1186/s13046-017-0574-5
Ye, F., Chen, Y., Xia, L., Lian, J., and Yang, S. (2018). Aldolase A Overexpression Is Associated with Poor Prognosis and Promotes Tumor Progression by the Epithelial-Mesenchymal Transition in colon Cancer. Biochem. Biophys. Res. Commun. 497, 639–645. doi:10.1016/j.bbrc.2018.02.123
Yochum, Z. A., Cades, J., Wang, H., Chatterjee, S., Simons, B. W., O'brien, J. P., et al. (2019). Targeting the EMT Transcription Factor TWIST1 Overcomes Resistance to EGFR Inhibitors in EGFR-Mutant Non-small-cell Lung Cancer. Oncogene 38, 656–670. doi:10.1038/s41388-018-0482-y
Yu, M., Qi, B., Xiaoxiang, W., Xu, J., and Liu, X. (2017). Baicalein Increases Cisplatin Sensitivity of A549 Lung Adenocarcinoma Cells via PI3K/Akt/NF-Κb Pathway. Biomed. Pharmacother. 90, 677–685. doi:10.1016/j.biopha.2017.04.001
Yuan, R., Chang, J., and He, J. (2020). Effects of Twist1 on Drug Resistance of Chronic Myeloid Leukemia Cells through the PI3K/AKT Signaling Pathway. Cell Mol Biol (Noisy-le-grand) 66, 81–85. doi:10.14715/cmb/2020.66.6.15
Zeng, Q., Zhang, Y., Zhang, W., and Guo, Q. (2020). Baicalein Suppresses the Proliferation and Invasiveness of Colorectal Cancer Cells by Inhibiting Snailinduced Epithelialmesenchymal Transition. Mol. Med. Rep. 21 (6), 2544–2552. doi:10.3892/mmr.2020.11051
Zhang, X. W., Liu, W., Jiang, H. L., and Mao, B. (2018). Chinese Herbal Medicine for Advanced Non-small-cell Lung Cancer: A Systematic Review and Meta-Analysis. Am. J. Chin. Med. 46, 923–952. doi:10.1142/S0192415X18500490
Keywords: epithelial–mesenchymal transition, reconstructed TFAE, stable isotope dimethyl-labeled proteomics, Twist1, PI3K/Akt signaling pathway
Citation: Cao H-J, Zhou W, Xian X-L, Sun S-J, Ding P-J, Tian C-Y, Tian F-L, Jiang C-H, Fu T-T, Zhao S and Dai J-Y (2022) A Mixture of Baicalein, Wogonin, and Oroxylin-A Inhibits EMT in the A549 Cell Line via the PI3K/AKT-TWIST1-Glycolysis Pathway. Front. Pharmacol. 12:821485. doi: 10.3389/fphar.2021.821485
Received: 24 November 2021; Accepted: 31 December 2021;
Published: 09 February 2022.
Edited by:
Sirajudheen Anwar, University of Hail, Saudi ArabiaReviewed by:
Onur Bender, Ankara University, TurkeyAntonio Speciale, University of Messina, Italy
Germano Aguiar Ferreira, Hemocentro Foundation of Ribeirão Preto, Brazil
Copyright © 2022 Cao, Zhou, Xian, Sun, Ding, Tian, Tian, Jiang, Fu, Zhao and Dai. 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: Jian-Ye Dai, daijy@lzu.edu.cn
†These authors have contributed equally to this work and share first authorship