Integration of transcriptomics and proteomics to elucidate inhibitory effect and mechanism of rosmarinic acid from Perilla frutescens (L.) Britt. in treating Trichophyton mentagrophytes
Chinese Medicine volume 18, Article number: 67 (2023)
Dermatophyte caused by Trichophyton mentagrophytes is a global disease with a growing prevalence that is difficult to cure. Perilla frutescens (L.) Britt. is an edible and medicinal plant. Ancient books of Traditional Chinese Medicine and modern pharmacological studies have shown that it has potential anti-fungi activity. This is the first study to explore the inhibitory effects of compounds from P. frutescens on Trichophyton mentagrophytes and its mechanism of action coupled with the antifungal activity in vitro from network pharmacology, transcriptomics and proteomics.
Five most potential inhibitory compounds against fungi in P. frutescens was screened with network pharmacology. The antifungal activity of the candidates was detected by a broth microdilution method. Through in vitro antifungal assays screening the compound with efficacy, transcriptomics and proteomics were performed to investigate the pharmacological mechanisms of the effective compound against Trichophyton mentagrophytes. Furthermore, the real-time polymerase chain reaction (PCR) was applied to verify the expression of genes.
The top five potential antifungal compounds in P. frutescens screened by network pharmacology are: progesterone, luteolin, apigenin, ursolic acid and rosmarinic acid. In vitro antifungal assays showed that rosmarinic acid had a favorable inhibitory effect on fungi. The transcriptomic findings exhibited that the differentially expressed genes of fungus after rosmarinic acid intervention were mainly enriched in the carbon metabolism pathway, while the proteomic findings suggested that rosmarinic acid could inhibit the average growth of Trichophyton mentagrophytes by interfering with the expression of enolase in the glycolysis pathway. Comparison of real-time PCR and transcriptomics results showed that the trends of gene expression in glycolytic, carbon metabolism and glutathione metabolic pathways were identical. The binding modes and interactions between rosmarinic acid and enolase were preliminary explored by molecular docking analysis.
The key findings of the present study manifested that rosmarinic acid, a medicinal compound extracted from P. frutescens, had pharmacological activity in inhibiting the growth of Trichophyton mentagrophytes by affecting its enolase expression to reduce metabolism. Rosmarinic acid is expected to be an efficacious product for prevention and treatment of dermatophytes.
In recent years, an increasing number of cases of dermatophytes have been reported, in which Trichophyton mentagrophytes responsible for human dermatophytosis such as tinea corpora, tinea cruris and tinea capitis [1,2,3]. Although many drugs like terbinafine and azoles have shown desirable curative effects in treating intractable recurrent skin tinea disease, abuse of these drugs is liable to produce drug-resistant fungi [4,5,6] like Trichophyton mentagrophytes, thus significantly reducing their therapeutic effects in clinical practice [7,8,9]. Another study manifested that Trichophyton mentagrophytes could be isolated from patients with beriberi . In addition, treatment with fluconazole and other drugs contributes to an increase in 4β-hydroxycholesterol , which in turn negatively impacts the body.
With the development of Chinese pharmacology, an increasing number of natural bioactivity components used in combination with the commonly-used antifungal drugs can improve the susceptibility of dermatophytes to antimicrobials and even reverse drug resistance. Pharmacological studies have exhibited that natural compounds like curcumin, eugenol, magnoflorine, and geraniol can inhibit drug-resistant dermatophytes with different mechanisms from the existing antifungals [12,13,14,15,16]. Therefore, the use of newly-discovered natural antifungal compounds will probably mitigate the evolution of drug-resistant fungi, reduce the drug dose and improve the therapeutic efficacy in the treatment of dermatophytes. Perilla frutescens (L.) Britt. is an annual herb of the Labiatae family [17, 18]. The medicinal part is mainly the leaves, which have analgesic, sedative and detoxifying effects. In the Northern Song Dynasty, the Taiping Shenghui Fang written by Wang Huaiyin recorded the use of P. frutescens as the main medicine, Zisu Powder, for the treatment of women's beriberi. In addition, the Compendium of Materia Medica compiled by Li Shizhen, a world-renowned medical scientist in the Ming dynasty has recorded that P. frutescens can effectively treat diarrhea, abdominal distension, dermatophytosis, etc. In Chinese medicinals, the stems, leaves, and seeds of P. frutescens are commonly used as effective medicines for treating pain, colds, cough, nausea, poisoning, asthma and constipation . Modern studies have shown that P. frutescens has anti-bacterial , anti-inflammatory [21,22,23], anti-cancer , anti-viral [25, 26], and anti-aging effects . Most of these studies have focused on the inhibitory effects of P. frutescens extracts on the common bacteria and plant fungi in daily foods [20, 28]. However, studies about the effects of P. frutescens extracts on dermatophytes, a type of fungus present on the body surface of animals, are rather limited . Investigators determined the minimum inhibitory doses and minimum fungicidal dose of Trichophyton mentagrophytes by using airtight boxes, in which P. frutescens extracts had a good inhibitory effect [30, 31]. Nonetheless, the potential mechanism of how the active ingredients of P. frutescens extracts inhibit the growth of dermatophytes remains unclear. However, the network pharmacology through the disease-target-drug interaction network, contributes to the discovery of new functions of known compounds .
Given all the setbacks above, by utilizing network pharmacology, transcriptomics and proteomics, the present study aimed to screen the potential compounds from P. frutescens in inhibiting the growth of Trichophyton dermatophytes and elucidating the potential molecular mechanisms. In this study, firstly, network pharmacology was utilized to screen the potential compounds and then performed in vitro fungal inhibition assays to verify the therapeutic efficacy of the potential compounds. Besides, transcriptomics and proteomics were adopted to explore the potential compounds’ targets of action, and real-time PCR and molecular docking analysis were used to verify the gene expression. This study would probably provide a potentially more effective treatment for Trichophyton mentagrophytes. The whole study was illustrated in the flowchart of Fig. 1.
Materials and methods
Network pharmacology analysis
Active compounds collection and targets prediction
P. frutescens was used as a keyword to search the terms in the Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform (TCMSP) . Compounds that met the following criteria were selected and considered as P. frutescens active ones: drug-likeness (DL) ≥ 0.18 . The targets corresponding to the active compounds were extracted as P. frutescens targets. P. frutescens and subsequent tinea targets were normalized in the UniProt protein database to prevent bias caused by different nomenclatural methods .
Tinea targets collection
Tinea was used as a keyword to search in the Comparative Toxicogenomics Database (CTD)  to obtain tinea-related targets, and the intersection of P. frutescens and tinea targets was obtained by EXCEL. The Venn diagram was plotted through the online website www. bioinformatics.com.cn.
Drug-target-disease network construction
Cytoscape (version 3.9.0)  was used to construct the P. frutescens-target-tinea network to visualize the complex relationships between compounds, targets and diseases. Analysis Network of Cytoscape was used to perform the network characteristics analysis and to screen the major active compounds.
Chemicals and reagents
Progesterone, ursolic acid, apigenin, luteolin and miconazole (≥ 95%) were purchased from Aladdin. Rosmarinic acid (≥ 95%) was purchased from Shanghai Puyu Kemao Co., Ltd. (Shanghai, China) and Dimethyl sulfoxide (DMSO) (≥ 99.8%) from Shandong Yousuo Chemical Technology Co., Ltd. (Linyi, Shandong, China).SDS-PAGE gel preparation kits, Prestained Protein Marker II (10–200 kDa), RPMI-1640, Radio Immunoprecipitation Assay (RIPA) lysate, Cocktail protease inhibitors, phenylmethylsulfonyl fluoride (PMSF), phosphorylation protease inhibitors, bicinchoninic acid (BCA) kit, Coomassie brilliant blue, loading buffer, Trizol, SweScript RT II First Strand cDNA Synthesis Kit and 2 × SYBR Green qPCR Master Mix were all purchased from Servicebio (Wuhan, Hubei, China). Fetal Bovine Serum (FBS) was purchased from Zhejiang Tianhang Biotechnology Co., Ltd. (Huzhou, Zhejiang, China). Potato Dextrose Agar (PDA) was purchased from Qingdao Hope Bio-Technology Co., Ltd. (Qingdao, Shandong, China) Clinical strains 1078, 1079, 130, and 1413 were shared by Wuhan No. 1 Hospital (Wuhan, Hubei, China). Trichophyton mentagrophytes ATCC9533 was purchased from Shanghai Bioresource Collection Center (Shanghai, China).
Preparation of drug working solution and spore suspension
Progesterone, ursolic acid, apigenin and luteolin were dissolved using DMSO to make a final concentration of 12.8 mg/mL mother solution and then filtered for sterilization. An appropriate amount of rosmarinic acid powder was weighed and dissolved in RPMI-1640 with 10% FBS to make the final concentration of 8 mg/mL and then filtered for sterilization. Trichophyton dermatophytes were cultivated on PDA at 28 ºC and were periodically transferred at 60-day intervals for preservation, filtered by 8 layers of sterile gauze with an appropriate amount of distilled water, and counted by blood cell counting plates. The concentration of spore suspension was finally diluted to 2 × 104 CFU/mL by adding RPMI-1640 with 10% FBS.
Determination of minimum inhibitory concentration and minimum bactericidal concentration
The antifungal activity of the drug was tested by the broth microdilution method according to Clinical and Laboratory Standards Institute (CLSI) M38-A2. A sterile 96-well plate was taken, and 100 μL RPMI-1640 containing 10% FBS was added to the well 1 of each row as a blank control. 200 μL of rosmarinic acid working solution containing spore suspension was added to well 2. 100 μL of spore suspension was added to wells 3–11 respectively. Wells 3–11 were two-fold diluted. Well 12 in each row did not contain drugs and was used as the growth control. The drug-sensitive plates were incubated at 28 ºC for 6 days under constant temperature and light. Compared with the growth control, the minimum inhibitory concentration (MIC) was determined to be the lowest drug concentration where higher than 80% growth was inhibited. The maximum concentration of progesterone, ursolic acid, apigenin and luteolin were set as 128 μg/mL two-fold dilution, and miconazole concentration in the positive group was determined to be 8 μg/mL, and the interference was eliminated by adding DMSO of equal concentration in the control group, with DMSO concentrations ≤ 1% in all systems. The minimum bactericidal concentration (MBC) was determined by the plate recovery method, and 100 μL liquid absorbed from the integral multiples MIC wells measured on the 6th day was coated onto PDA plates and incubated at 28 ºC for 5 days. The drug concentration corresponding to the colony-free growth plate was used as MBC. All experiments were performed in three biological replicates.
Transcriptomic analysis and data processing
After 4 h treatment with 1MIC rosmarinic acid, the mycelia of well-grown Trichophyton mentagrophytes were collected for transcriptomic sequencing, which was performed by Beijing Biomarker Technologies Co.,Ltd. Briefly, total RNA was extracted from the control and rosmarinic acid groups and its quality was detected, then mRNA was purified and fragmented, and finally double-stranded cDNA was synthesized using mRNA as the template. RNA sequencing library was constructed by PCR amplification. Upon completion, library quality was assessed on the Agilent Bioanalyzer 2100. After qualified library screening, TruSeq PE Cluster Kit v3-cBot-HS (Illumia) was used to cluster index-coded samples on cBot Cluster Generation System. Library preparations were then sequenced using the Illumina Hiseq 2000 and paired-end reads were generated. Trinity software was used to assemble all reads , and the gene expression level of each sample was estimated by RNA-Seq by Expectation–Maximization (RSEM) . Specific steps are available in Additional file 1: Method 1. Sequencing results have been uploaded to the National Center for Biotechnology Information, and sequence number PRJNA847061 provides access to all RNA-seq data.
In this experiment, the DESeq R package was used to analyze differentially expressed genes (DEGs) in the two samples. In the process of DEGs analysis, the recognized and effective Benjamini–Hochberg method was used to correct the P-value of significance obtained from the original hypothesis test, and the corrected P-value was finally adopted. False Discovery Rate (FDR) was used as a key indicator for DEGs screening to reduce the false positives caused by independent statistical hypothesis testing of the expression values of numerous genes. In the screening process, FDR less than 0.05 and Fold Change (FC) greater than or equal to 2 were considered as the screening criteria. Here, FC represents the ratio of expression between two samples or groups. Functional annotation of DEGs was performed using eggNOG , Gene Ontology (GO, http://www.geneontology.org/) and Kyoko Encyclopedia of Genes and Genomes (KEGG, http://www.genome.jp/kegg/) databases, enrichment analysis used the clusterProfiler R package , the heat map and volcano plot used gplots and ggplot2 R package respectively.
Proteomic analysis and data processing
Rosmarinic acid working solution was added to the 5-day cultured Trichophyton mentagrophytes solution to adjust the concentration to 1 MIC, and RPMI-1640 containing 10% FBS was used as the control. The samples were incubated at 28 ºC for 4 h with constant temperature shock, then centrifuged at 5000 × g for 15 min. After the supernatant was discarded, the samples were washed with ultrapure water 3 times. RIPA lysate containing, coupled with its instructions specified Cocktail protease inhibitors, PMSF and phosphorylation protease inhibitors were then added, and the samples were centrifuged at 12,000 × g for 15 min after shaking for 10 min on ice. After BCA kit was used for protein quantification, the lowest protein concentration sample was used as the reference to adjust the protein concentration of each sample, and bovine serum protein was used as the standard. After adding loading buffer, the samples were boiled in water bath for 5 min, followed by SDS-PAGE of 8% separation glue and 5% concentration glue. Finally, Coomassie brilliant blue was used for staining after electrophoresis.
Strip sequencing was performed by Biotree Biotech Co.,Ltd. (Shanghai, China). In brief, the bands on SDS-PAGE were cut off, decolorized and dehydrated, and the polypeptide chains were extracted after trypsin digestion, dried and desalted by C18 column, and then analyzed by nanoLC-MS/MS. The off-machine data were searched and analyzed by Proteome Discoverer (version 188.8.131.525, Thermo Fisher Scientific) and the built-in Sequest HT search engine. After the abundance was normalized, differential expression was calculated, and clusterProfiler and ggplot2 R package were used for enrichment analysis and mapping. Details are available in Additional file 1: Method 2.
Quantitative real-time PCR assays
The well-grown Trichophyton mentagrophytes were treated with 1 MIC rosmarinic acid for 4 h to collect mycelia. After grinding with liquid nitrogen, the RNA samples were extracted using Trizol. SweScript RT II First Strand cDNA Synthesis Kit was used to reverse transcribe RNA to obtain cDNA. Amplification was performed using 2 × SYBR Green qPCR Master Mix on qTower 3G Real-time PCR Instrument (Analytik Jena, Germany). PCR reaction conditions were as follows: 40 cycles of predenaturation at 95 ºC for 30 s, denaturation at 95 ºC for 15 s, annealing at 65 ºC for 30 s, and extension at 72 ºC for 30 s. 18S RNA gene was used as reference, and the differential expression levels of target genes and reference genes were calculated by 2(−ΔΔCt) method. All primer sequences used in this experiment are available in Additional file 1: Table S1.
Molecular Docking Studies
Virtual molecular docking was employed to analyze the potential binding modes of enolase and rosmarinic acid (Fig. 2C). Four recognized enolase inhibitors were used for comparison [42,43,44,45]. All compound structures were obtained from the PubChem database (https://pubchem.ncbi.nlm.nih.gov/). AP-III-a4 is nonsubstrate analog that directly binds to enolase and inhibits its activity. Hex and D-(-)-3-Phosphoglyceric acid disodium are substrate-competitive enolase inhibitors and POMHEX increases the permeability of HEX into cells and tissue. The crystal structure of enolase was downloaded from the Uniport, which was selected and saved as pdb format. The ligand and receptor were split by Pymol. Autodock Tools 1.5.7 was used to transform pdb to pdbqt format files with gird boxes adjusted to cover the entire pocket for the preparation of virtual docking. The structures of chemicals were collected from Pubchem as sdf format and transformed into pdbqt format . Autodock vina 1.2.0 was used to simulate the potential interactions among the selected compounds and the targets .
The data were analyzed using GraphPad Prism statistical software. All results are presented as mean ± SEM deviation. One-way analysis of variance was performed between multiple groups when the homogeneity of variance and normality were met. P < 0.05 was considered statistically significant.
Screening of P. frutescens active ingredients and targets and establishment of the drug-target-disease network
328 P. frutescens compounds were screened from TCMSP. Twenty-nine active compounds were screened by DL parameters (Additional file 1: Table S2). The targets corresponding to the active compounds were extracted and entered into EXCEL, with removed duplicate entries, 253 P. frutescens targets were obtained after merging.
A total of 1,262 Tinea targets were collected through the CTD database, and several different Tinea diseases were obtained by Tiena search (Additional file 1: Table S3). A total of 99 related targets were obtained through the mapping of P. frutescens and Tinea targets. Intersection of P. frutescens and Tinea targets was then plotted into the Venn diagram (Fig. 2A).
In the P. frutescens-target-Tinea network (Fig. 2B), the Degree of a single target is the number of the connected nodes. The network analysis was made by the Cytoscape, and the active compounds were sorted by Degree. The top 5 active compounds, namely, progesterone, luteolin, apigenin, ursolic acid and rosmarinic acid were listed by Betweenness Centrality, Closeness Centrality and Degree (Additional file 1: Table S4).
Antifungi activity analysis
In this experiment, MIC refers to the drug concentration at which the growth rate after drug treatment was less than 20% of the blank group. The inhibitory effect of luteolin on Trichophyton dermatophytes was shown in Fig. 3A, which revealed no significant inhibitory effect on all five experimental fungi. Ursolic acid and apigenin demonstrated similar inhibitory effects on Trichophyton dermatophytes (Fig. 3B, C). With increasing concentrations, the growth rate of 130 and 1413 slowed down but had no significant inhibitory effect. The inhibitory effect of progesterone on the experimental fungi was shown in Fig. 4, with MIC values of 32 μg/mL for 1078, 1079, 130 and 1413 and 64 μg/mL for ATCC 9533 in a dose-dependent manner (Fig. 4A, B). It was considered that the drug concentration for no colony growthte incubation, it was found that colonies still grew on the plate, indicating that progesterone could only inhibit fungal growth, but had no bactericidal activity. The inhibitory effect of rosmarinic acid on the experimental fungi was shown in Fig. 5 in a dose-dependent manner, with MIC values of 1000 μg/mL for 1078 and 1079, 125 μg/mL for 130 μg/mL and 250 μg/mL for 1413 (Fig. 5A, B). MBC values of 2000 μg/mL for rosmarinic acid against 1078 and 1079, 250 μg/mL for 130 and 250 μg/mL for 1413 were observed after plate incubation (Fig. 5C).
To explore the inhibitory mechanism of Trichophyton mentagrophytes by rosmarinic acid, RNA-seq was used to identify differential mRNA expression. The base quality values measured by Illumina sequencing platform were expressed as Phred scores (Additional file 1: Fig. S1 and Table S5). The quality score of reads was Q30 (base recognition accuracy ≥ 99.9%) and the error rate of single base sequencing was < 0.1% (Additional file 1: Table S6). Since the first few bases at the 5' end of the reads were random primer sequences with some bias, there was a significant fluctuation in the front end of the base distribution map (Additional file 1: Fig. S2)0.15,946 Unigene and 38,006 transcripts were obtained using Trinity assembly (Additional file 1: Fig. S3 and Table S7). A total of 32,798 annotated transcripts were obtained by comparing transcripts with eggNOG database using the eggNOG-Mapper tool (Additional file 1: Table S8). The reads obtained by sequencing were compared with Unigene library with Bowtie. Based on the comparison results, RSEM was used to estimate the expression level, and the expression abundance of the corresponding Unigene was represented by the trimmed mean of M values value (Additional file 1: Fig. S4, S5).
There were significant differences in gene expression between the control group and the rosmarinic acid treatment group. The genes were clustered according to the same or similar expression patterns (Fig. 6A), where the color-key abscissa bar indicates the multiple difference of the normalized genes and the larger absolute value indicates the more obvious differential expression. The ordinate represents the density of differential gene enrichment in the corresponding value and the trend line in the heat map represents the degree of gene deviation from the median value. The volcano plot showed the differential gene expression versus priority, where each point represents a gene, and the negative value of the logarithm of the difference Fold (horizontal coordinate: log2 Fold Change) indicates a decrease of gene expression, while the positive number indicates an increase (Fig. 6B and Additional file 1: Table S9). The larger logarithm of error incidence (vertical coordinate: log10 FDR) indicates higher confidence in the results.
DEGs were mapped to the GO database. The enrichment results of DEGs were obtained from Biological Process (BP), Cellular Component (CC), and Molecular Function (MF) of GO terms (Fig. 7A). There were 110 statistically significant GO terms, including 83 BP terms, 7 CC terms and 20 MF terms. DEGs in BP were mainly enriched in detoxification, cell modified amino acid metabolism and glutathione metabolism; DEGs in CC were mainly enriched in plasma membrane, mitochondrial nuclear mass and integral components of mitochondrial tricarboxylic acid cycle enzyme complex; DEGs in MF were mainly enriched in oxidoreductase activity, cation transmembrane transporter activity and proton transmembrane transporter activity.
The topGO R package was used to rank the pathways to visualize the GO nodes and their hierarchical relationships for differential gene enrichment. The most significant nodes were the hexose catabolic process, siderophore metabolic process, and iron assimilation under the homeostatic process. Based on the directed acyclic graph of CC, the most significant DEGs enrichment nodes were found in mitochondrial nucleocapsid and plasma membrane components. The directed acyclic graph of MF showed that the main enriched nodes of DEGs were oxidoreductase activity and transmembrane transporter activity (Fig. 7B–D).
Figure 8A showed the top 10 pathways of KEGG pathway terms, most of which were related to metabolisms, such as the metabolisms of Glutathione, cysteine, methionine, amino sugar, nucleotide sugar, starch and sucrose. The significantly enriched metabolic pathways in KEGG were selected and analyzed by Gene Set Enrichment Analysis. The results exhibited that the expression of genes enriched in glutathione metabolic pathway in the blank group was higher than that in the rosmarinic acid treatment group, while the expression levels of genes enriched in other metabolic pathways, such as amino acid sugar and nucleotide sugar metabolism, biosynthesis of secondary metabolites and carbohydrate metabolism, were lower in the blank group than in the rosmarinic acid treatment group (Fig. 8B).
A large number of protein bands were found to disappear or appear around 75 kDa, 55 kDa, 42 kDa, and 30 kDa after treatment with rosmarinic acid compared to the control group by SDS-PAGE electrophoresis (Fig. 9A), which was possibly related to the inhibition of fungal growth caused by the increased or decreased expression of these proteins. SDS-PAGE subjected to qualitative proteomic assays further clarifies the mechanism of rosmarinic acid action. 17 of all 35 proteins were identified as differentially expressed, including 4 down-regulated proteins and 13 up-regulated ones (Fig. 9B).
The GO and KEGG pathway analysis were then used to further investigate the mechanism of rosmarinic acid against Trichophyton mentagrophytes (fold change ≥ 2 and P < 0.05). There were 221 statistically significant GO terms in Fig. 9C, including 190 BP terms, 20 CC terms and 11 MF terms. Differential proteins in BP were mainly enriched in glucose catabolism, NADH regeneration, fructose-6-phosphate glycolysis and glucose-6-phosphate glycolysis; differential proteins in CC were mainly enriched in fungal cell wall, cell wall, outer envelope structure and extracellular region; the differential proteins in MF were mainly enriched in glyceraldehyde-3-phosphate dehydrogenase (NAD +) (phosphorylation) activity, glyceraldehyde-3-phosphate dehydrogenase (NAD (P) +) (phosphorylation) activity and promoter specific chromatin binding and oxidoreductase activity.
Among the top 10 pathways of KEGG pathway terms (Fig. 9D), the glycolysis/glucose production pathway was highlighted. The significantly down-regulated differentially expressed protein X5CH36 in Additional file 1: Table S10 was shown on Uniport as triglyceride 2-phosphodehydrase, another name for enolase.
Real-time PCR verification results
Comparison of real-time PCR and transcriptomics results showed that the expression trends of genes in glycolytic, carbon metabolism and glutathione metabolic pathways were identical (Fig. 10). Real-time PCR results of genes in the rosmarinic acid treatment group showed significant low or high fold changes compared with the control group, although no significant differences between these groups were observed in some transcriptomics results (Fig. 10). Therefore, the transcriptomics results can be considered reliable.
Molecular docking studies
Molecular docking was used to explore the possible binding sites and modes of action between enolase and the active compounds including AP-III-a4, Hex, D-(-)-3-Phosphoglyceric acid disodium, POMHEX and rosmarinic acid. As shown in Table 1, the binding energies of Hex, POMHEX and Rosmarinic acid with enolase were all less than − 5.0 kcal/mol, indicating their good binding affinity with enolase. Moreover, rosmarinic acid showed better binding affinity than AP-III-a4, Hex, D-(-)-3-Phosphoglyceric acid disodium and POMHEX. The 3D and 2D interaction diagrams of the compounds with enolase were shown in Fig. 11.
Although many drugs are clinically available to treat dermatophytes, the emergence of drug-resistant fungi has reduced efficacy of existing drugs and even caused recurrent skin tinea disease. This may be caused by the massive abuse of existing drugs, or related to the fact that mutations in genes in the ergosterol pathway cause their resistance . Therefore, it is necessary to discover more effective compounds from Traditional Chinese medicine to provide new treatment protocols. However, the complex composition of Traditional Chinese medicine makes it challenging to discover such effective compounds. With the rapid development of genomic technologies such as genomics, transcriptomics, proteomics and metabolomics, it has become possible to elucidate complex biological phenomena and obtain large amounts of genetic data, which can help explore natural compounds from Traditional Chinese medicine in treating dermatophytes . In this study, the TCMSP database was used to screen potentially active compounds and transcriptomics and proteomics to make predictions about the inhibitory mechanisms of Trichophyton mentagrophytes growth.
P. frutescens is a medicinal and edible plant widely distributed in subtropical regions. P. frutescens is used as a medicine in China, Japan and Southeast Asian countries due to its significant medicinal functions and minor adverse effects. Rosmarinic acid is the most abundant bioactive polyphenolic component of P. frutescens leaves . The important effects of rosmarinic acid and its derivatives on human health have been demonstrated in quite several in vitro/in vivo pharmaceutical and clinical studies. In recent years, the biological activity of rosmarinic acid has drawn widespread attention from the academic community, as it doesn’t merely inhibit the growth of methicillin-resistant Staphylococcus aureus, but also the growth of other drug-resistant fungi through its combined use with antibiotics. Besides, studies have also exhibited that rosmarinic acid has certain inhibitory effects on some fungi . This study found that rosmarinic acid from P. frutescens could inhibit the average growth of Trichophyton mentagrophytes by affecting their metabolism.
Inhibition assays in vitro are commonly used to assess the ability of drugs of inhibiting microorganisms. In line with other studies, MIC and MBC are the main parameters for evaluating the antimicrobial ability of drugs . In this study, among five compounds screened by network pharmacology: progesterone, luteolin, apigenin, ursolic acid and rosmarinic acid, only progesterone and rosmarinic acid exhibited desired antifungal ability in vitro. The results of MIC and spot test indicated that progesterone has a lower MIC but no MBC (Fig. 4), which is consistent with the previous study. However, the rosmarinic acid could not only inhibit the growth of Trichophyton mentagrophytes, which is one of filamentous fungi, in a dose-dependent manner, but also killed it at a concentration of 2MIC (Fig. 5). The previous studies indicated the rosmarinic acid could play a significant role in inhibiting the growth of filamentous fungi, and filamentous fungi treated with high concentrations of rosmarinic acid cannot grow on nutrient-rich plates . It is consistent with the results of the study.
At present, commercially available antifungal drugs usually work by interfering with ergosterol synthesis, formation of water channels by ergosterol binding on cell membranes to interfere with cell wall synthesis and interference with nucleic acid replication . However, the recent emergence of drug-resistant fungi in the world is mainly due to mutations in the squalene epoxidase gene and squalene epoxidase is a key link in the ergosterol synthesis pathway [7,8,9]. Therefore, there are significant implications in the antifungal field by developing a new pathway. As we all known, there was no report of antifungal drugs through affecting enolase expression. However, the proteomic results in this study indicated that the expression of enolase significantly decreased while Trichophyton mentagrophytes treated with rosmarinic acid (Fig. 9). There were studies shown that inhibition of enolase could have significant effects on cell growth . Enolase is a key glycolytic metalloenzyme involved in carbon metabolism and is mainly responsible for the key step of catalyzing the conversion of 2-phosphoglycerate (2-PG) to phosphoenolpyruvate (PEP) in glycolysis, which provides substrate support for subsequent oxidative phosphorylation and is an important protein in glycolysis . Most of the substrates they react to are carbohydrates or their derivatives, such as glycolic acid and amino acid . In addition, the activation of the mitogen-activated protein kinase MAPK pathway could increase the expression of enolase . The results of this study showed that there were 52 transcripts enriched in the MAPK pathway, and they were highly expressed in the rosmarinic acid treatment group (Fig. 8B). These results suggested rosmarinic acid could decrease the expression of enolase by blocking the upstream pathway of enolase synthesis. Related results were shown in Table S10 in Additional file 1.
Enolase serves as a prototype for metalloenzymes with labile metal ions that play a key role in catalytic turnover . In the present study, the molecular docking results suggested that rosmarinic acid and the substrate-competitive inhibitor of enolase Hex and POMHEX could bind with the Magnesium ion in enolase via metal-acceptor interaction (Fig. 11). The reported studies have showed that the amino acid residues GLU211, LYS345, HIS159, GLU168 and ASP246 could be potential binding sites of the substrate 2-PG and enolase [57, 58]. Meanwhile, in the present study, the hydrogen bond interaction formed the rosmarinic acid and HIS159 residue. The bond of the GLU168, ASP246 and ASP323 residue with the benzene ring of rosmarinic acid was affected by Pi-anion and Pi-cation interaction. The hydrogen bond interaction formed the Hex and the GLU211 residue. And the bond of the HIS159 residue with the Hex and POMHEX were all affected by Pi-Alkyl interaction (Table 1). These results suggested that rosmarinic acid could be a substrate-competitive inhibitor of enolase.
The findings of this study demonstrated that rosmarinic acid has a certain degree of inhibition and lethal effect on the growth of Trichophyton mentagrophytes. And the preliminary mechanisms exploration indicated that the rosmarinic acid might affect the normal metabolism and glycolysis of Trichophyton mentagrophytes by inhibiting the expression of enolase, resulting in an insufficient intracellular supply of energy to fungi, thus inhibiting the growth of Trichophyton mentagrophytes due to insufficient energy supply for normal physiological activities (Fig. 12). This provided more extensive research ideas for the development of antifungal drugs.
Availability of data and materials
The datasets during and/or analyzed during the current study available from the corresponding author on reasonable request.
Clinical and Laboratory Standards Institute
Comparative toxicogenomics database
Differentially expressed genes
Fetal bovine serum
False discovery rate
Kyoko Encyclopedia of Genes and Genomes
Mitogen-activated protein kinase
Minimum bactericidal concentration
Minimum inhibitory concentration
Polymerase chain reaction
Potato dextrose agar
Radio immunoprecipitation assay
RNA-Seq by Expectation–Maximization
Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform
Rao AG, Datta N. Tinea corporis due to Trichophyton mentagrophytes and Trichophyton tonsurans mimicking tinea imbricata. Indian J Dermatol Venereol Leprol. 2013;79(4):554. https://doi.org/10.4103/0378-6323.113109.
Khosravi AR, Mansouri P, Naraghi Z, Shokri H, Ziglari T. Unusual presentation of tinea cruris due to Trichophyton mentagrophytes var. mentagrophytes. J Dermatol. 2008;35(8):541–5. https://doi.org/10.1111/j.1346-8138.2008.00517.x.
Muller VL, Kappa-Markovi K, Hyun J, Georgas D, Silberfarb G, Paasch U, Uhrlass S, Nenoff P, Schaller J. Tinea capitis et barbae caused by Trichophyton tonsurans: a retrospective cohort study of an infection chain after shavings in barber shops. Mycoses. 2021;64(4):428–36. https://doi.org/10.1111/myc.13231.
Aharaz A, Jemec GBE, Hay RJ, Saunte DML. Tinea capitis asymptomatic carriers: what is the evidence behind treatment? J Eur Acad Dermatol Venereol. 2021;35(11):2199–207. https://doi.org/10.1111/jdv.17462.
Feng J, Liu F, Wu F, De Deng Q, Zeng HM, Kong TQ, Chen J, Sang H. Tinea infection with scutula-like lesions caused by Microsporum gypseum in a SLE patient: case report and literature review. Mycopathologia. 2013;176(3–4):255–8. https://doi.org/10.1007/s11046-013-9652-1.
Kong QT, Du X, Yang R, Huang SY, Sang H, Liu WD. Chronically recurrent and widespread tinea corporis due to Trichophyton rubrum in an immunocompetent patient. Mycopathologia. 2015;179(3–4):293–7. https://doi.org/10.1007/s11046-014-9834-5.
Hsieh A, Quenan S, Riat A, Toutous-Trellu L, Fontao L. A new mutation in the SQLE gene of Trichophyton mentagrophytes associated to terbinafine resistance in a couple with disseminated tinea corporis. J Mycol Med. 2019;29(4):352–5. https://doi.org/10.1016/j.mycmed.2019.100903.
Khurana A, Masih A, Chowdhary A, Sardana K, Borker S, Gupta A, Gautam RK, Sharma PK, Jain D. Correlation of in vitro susceptibility based on MICs and squalene epoxidase mutations with clinical response to terbinafine in patients with tinea corporis/cruris. Antimicrob Agents Chemother. 2018. https://doi.org/10.1128/AAC.01038-18.
Sacheli R, Harag S, Dehavay F, Evrard S, Rousseaux D, Adjetey A, Seidel L, Laffineur K, Lagrou K, Hayette MP. Belgian National Survey on Tinea Capitis: epidemiological considerations and highlight of terbinafine-resistant T. mentagrophytes with a mutation on SQLE gene. J Fungi (Basel). 2020. https://doi.org/10.3390/jof6040195.
Parekh M, Ramaiah G, Pashilkar P, Ramanujam R, Johnston P, Ilag LL. A pilot single centre, double blind, placebo controlled, randomized, parallel study of Calmagen(R) dermaceutical cream and lotion for the topical treatment of tinea and onychomycosis. BMC Complement Altern Med. 2017;17(1):464. https://doi.org/10.1186/s12906-017-1970-2.
Lutjohann D, Stellaard F, Kerksiek A, Lotsch J, Oertel BG. Serum 4beta-hydroxycholesterol increases during fluconazole treatment. Eur J Clin Pharmacol. 2021;77(5):659–69. https://doi.org/10.1007/s00228-020-03041-5.
Apisariyakul A, Vanittanakom N, Buddhasukh D. Antifungal activity of turmeric oil extracted from Curcuma longa (Zingiberaceae). J Ethnopharmacol. 1995;49(3):163–9. https://doi.org/10.1016/0378-8741(95)01320-2.
de Oliveira PF, Mendes JM, de Oliveira LE. Investigation on mechanism of antifungal activity of eugenol against Trichophyton rubrum. Med Mycol. 2013;51(5):507–13. https://doi.org/10.3109/13693786.2012.742966.
Luo N, Jin L, Yang C, Zhu Y, Ye X, Li X, Zhang B. Antifungal activity and potential mechanism of magnoflorine against Trichophyton rubrum. J Antibiot (Tokyo). 2021;74(3):206–14. https://doi.org/10.1038/s41429-020-00380-4.
Pereira Fde O, Mendes JM, Lima IO, Mota KS, Oliveira WA, Lima Ede O. Antifungal activity of geraniol and citronellol, two monoterpenes alcohols, against Trichophyton rubrum involves inhibition of ergosterol biosynthesis. Pharm Biol. 2015;53(2):228–34. https://doi.org/10.3109/13880209.2014.913299.
Salehi Z, Fatahi N, Taran M, Izadi A, Badali H, Hashemi SJ, Rezaie S, DaieGhazvini R, Ghaffari M, Aala F, et al. Comparison of in vitro antifungal activity of novel triazoles with available antifungal agents against dermatophyte species caused tinea pedis. J Mycol Med. 2020;30(2):100935. https://doi.org/10.1016/j.mycmed.2020.100935.
Fan Y, Cao X, Zhang M, Wei S, Zhu Y, Ouyang H, He J. Quantitative comparison and chemical profile analysis of different medicinal parts of Perilla frutescens (L.) Britt. from different varieties and harvest periods. J Agric Food Chem. 2022;70(28):8838–53. https://doi.org/10.1021/acs.jafc.2c03104.
Wang P, Jin B, Lian C, Guo K, Ma C. Comparative analysis of polycyclic aromatic hydrocarbons and halogenated polycyclic aromatic hydrocarbons in different parts of Perilla frutescens (L.) Britt. Molecules. 2022. https://doi.org/10.3390/molecules27103133.
Committee SP. China Pharmacopoeia Committee (CPC) PotPsRoC. Beijing: China Medical Science and Technology Press; 2020.
Zhao Y, Li H, Zhang Z, Ren Z, Yang F. Extraction, preparative monomer separation and antibacterial activity of total polyphenols from Perilla frutescens. Food Funct. 2022;13(2):880–90. https://doi.org/10.1039/d1fo02282b.
Jeon IH, Kim HS, Kang HJ, Lee HS, Jeong SI, Kim SJ, Jang SI. Anti-inflammatory and antipruritic effects of luteolin from Perilla (P. frutescens L.) leaves. Molecules. 2014;19(6):6941–51. https://doi.org/10.3390/molecules19066941.
Kangwan N, Pintha K, Lekawanvijit S, Suttajit M. Rosmarinic acid enriched fraction from perilla frutescens leaves strongly protects indomethacin-induced gastric ulcer in rats. Biomed Res Int. 2019;2019:9514703. https://doi.org/10.1155/2019/9514703.
Kitamura N, Nishino M, Fujii A, Hashizume K, Nakamura J, Kondo H, Ohuchi A, Hase T, Murase T. Perilla extract improves frequent urination in spontaneously hypertensive rats with enhancement of the urothelial presence and anti-inflammatory effects. Int J Urol. 2018;25(3):298–304. https://doi.org/10.1111/iju.13516.
Pintha K, Chaiwangyen W, Yodkeeree S, Suttajit M, Tantipaiboonwong P. Suppressive effects of rosmarinic acid rich fraction from perilla on oxidative stress, inflammation and metastasis ability in A549 cells exposed to PM via C-Jun, P-65-Nf-Kappab and Akt signaling pathways. Biomolecules. 2021. https://doi.org/10.3390/biom11081090.
Kuo YT, Liu CH, Corona A, Fanunza E, Tramontano E, Lin LT. The methanolic extract of perilla frutescens robustly restricts ebola virus glycoprotein-mediated entry. Viruses. 2021. https://doi.org/10.3390/v13091793.
Tang WF, Tsai HP, Chang YH, Chang TY, Hsieh CF, Lin CY, Lin GH, Chen YL, Jheng JR, Liu PC, et al. Perilla (Perilla frutescens) leaf extract inhibits SARS-CoV-2 via direct virus inactivation. Biomed J. 2021;44(3):293–303. https://doi.org/10.1016/j.bj.2021.01.005.
Bae JS, Han M, Shin HS, Kim MK, Shin CY, Lee DH, Chung JH. Perilla frutescens leaves extract ameliorates ultraviolet radiation-induced extracellular matrix damage in human dermal fibroblasts and hairless mice skin. J Ethnopharmacol. 2017;195:334–42. https://doi.org/10.1016/j.jep.2016.11.039.
Bais HP, Walker TS, Schweizer HP. Root specific elicitation and antimicrobial activity of rosmarinic acid in hairy root cultures of Ocimum basilicum. Plant Physiol Biochem. 2002;40(11):983–95. https://doi.org/10.1016/S0981-9428(02)01460-2.
Xiaoqin L, Fuzu W, Shiling Z. Effects of anti-bacteria and fungus of Perilla frutescens (L.) britton in vitro. Nat Prod Res Dev. 2000;12(1):42–5. https://doi.org/10.16333/j.1001-6880.2000.01.010.
Inouye S, Nishiyama Y, Uchida K, Hasumi Y, Yamaguchi H, Abe S. The vapor activity of oregano, perilla, tea tree, lavender, clove, and geranium oils against a Trichophyton mentagrophytes in a closed box. J Infect Chemother. 2006;12(6):349–54. https://doi.org/10.1007/s10156-006-0474-7.
Inouye S, Uchida K, Yamaguchi H. In-vitro and in-vivo anti-Trichophyton activity of essential oils by vapour contact. Mycoses. 2001;44(3–4):99–107. https://doi.org/10.1046/j.1439-0507.2001.00618.x.
Zhang R, Zhu X, Bai H, Ning K. Network pharmacology databases for traditional chinese medicine: review and assessment. Front Pharmacol. 2019;10:123. https://doi.org/10.3389/fphar.2019.00123.
Ru J, Li P, Wang J, Zhou W, Li B, Huang C, Li P, Guo Z, Tao W, Yang Y, et al. TCMSP: a database of systems pharmacology for drug discovery from herbal medicines. J Cheminform. 2014;6:13. https://doi.org/10.1186/1758-2946-6-13.
Tao W, Xu X, Wang X, Li B, Wang Y, Li Y, Yang L. Network pharmacology-based prediction of the active ingredients and potential targets of Chinese herbal Radix Curcumae formula for application to cardiovascular disease. J Ethnopharmacol. 2013;145(1):1–10. https://doi.org/10.1016/j.jep.2012.09.051.
UniProt C. UniProt: the universal protein knowledgebase in 2021. Nucleic Acids Res. 2021;49(D1):D480–9. https://doi.org/10.1093/nar/gkaa1100.
Davis AP, Grondin CJ, Johnson RJ, Sciaky D, Wiegers J, Wiegers TC, Mattingly CJ. Comparative toxicogenomics database (CTD): update 2021. Nucleic Acids Res. 2021;49(D1):D1138–43. https://doi.org/10.1093/nar/gkaa891.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504. https://doi.org/10.1101/gr.1239303.
Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng Q, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29(7):644–52. https://doi.org/10.1038/nbt.1883.
Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics. 2011;12:323. https://doi.org/10.1186/1471-2105-12-323.
Huerta-Cepas J, Forslund K, Coelho LP, Szklarczyk D, Jensen LJ, von Mering C, Bork P. Fast genome-wide functional annotation through orthology assignment by eggNOG-mapper. Mol Biol Evol. 2017;34(8):2115–22. https://doi.org/10.1093/molbev/msx148.
Wu T, Hu E, Xu S, Chen M, Guo P, Dai Z, Feng T, Zhou L, Tang W, Zhan L, et al. clusterProfiler 4.0: a universal enrichment tool for interpreting omics data. Innovation (Camb). 2021;2(3):100141. https://doi.org/10.1016/j.xinn.2021.100141.
Jung DW, Kim WH, Park SH, Lee J, Kim J, Su D, Ha HH, Chang YT, Williams DR. A unique small molecule inhibitor of enolase clarifies its role in fundamental biological processes. ACS Chem Biol. 2013;8(6):1271–82. https://doi.org/10.1021/cb300687k.
Lin YH, Satani N, Hammoudi N, Yan VC, Barekatain Y, Khadka S, Ackroyd JJ, Georgiou DK, Pham CD, Arthur K, et al. An enolase inhibitor for the targeted treatment of ENO1-deleted cancers. Nat Metab. 2020;2(12):1413–26. https://doi.org/10.1038/s42255-020-00313-3.
Ramos ML, Justino LL, Gil VM, Burrows HD. NMR and DFT studies of the complexation of W(VI) and Mo(VI) with 3-phospho-D-glyceric and 2-phospho-D-glyceric acids. Dalton Trans. 2009;43:9616–24. https://doi.org/10.1039/b905933d.
Zheng-Hong MFMDBWLY-HPBABP: Enolase inhibitors and methods of treatment there with. In: vol. WO2016145113A1, A61K31/662 edn. WIPO (PCT): UNIV TEXAS (US); 2016.
Kim S, Chen J, Cheng T, Gindulyte A, He J, He S, Li Q, Shoemaker BA, Thiessen PA, Yu B, et al. PubChem in 2021: new data content and improved web interfaces. Nucleic Acids Res. 2021;49(D1):D1388–95. https://doi.org/10.1093/nar/gkaa971.
Eberhardt J, Santos-Martins D, Tillack AF, Forli S. AutoDock vina 1.2.0: new docking methods, expanded force field, and python bindings. J Chem Inf Model. 2021;61(8):3891–8. https://doi.org/10.1021/acs.jcim.1c00203.
Zhang W, Huai Y, Miao Z, Qian A, Wang Y. Systems pharmacology for investigation of the mechanisms of action of traditional chinese medicine in drug discovery. Front Pharmacol. 2019;10:743. https://doi.org/10.3389/fphar.2019.00743.
Wei MC, Wang CS, Liou RM, Yang YC. Development and validation of a green and sustainable procedure for the preparation of Perilla frutescens extracts. Food Chem. 2022;369:130929. https://doi.org/10.1016/j.foodchem.2021.130929.
Zhu C, Liao B, Ye X, Zhou Y, Chen X, Liao M, Cheng L, Zhou X, Ren B. Artemisinin elevates ergosterol levels of Candida albicans to synergise with amphotericin B against oral candidiasis. Int J Antimicrob Agents. 2021;58(3):106394. https://doi.org/10.1016/j.ijantimicag.2021.106394.
Clemons KV, Stover EP, Schar G, Stathis PA, Chan K, Tokes L, Stevens DA, Feldman D. Steroid metabolism as a mechanism of escape from progesterone-mediated growth inhibition in Trichophyton mentagrophytes. J Biol Chem. 1989;264(19):11186–92.
Zhang AY, Camp WL, Elewski BE. Advances in topical and systemic antifungals. Dermatol Clin. 2007. https://doi.org/10.1016/j.det.2007.01.002.
Krucinska J, Falcone E, Erlandsen H, Hazeen A, Lombardo MN, Estrada A, Robinson VL, Anderson AC, Wright DL. Structural and functional studies of bacterial enolase, a potential target against gram-negative pathogens. Biochemistry. 2019;58(9):1188–97. https://doi.org/10.1021/acs.biochem.8b01298.
Piast M, Kustrzeba-Wojcicka I, Matusiewicz M, Banas T. Molecular evolution of enolase. Acta Biochim Pol. 2005;52(2):507–13.
Coleman JP, Smith CJ. Microbial metabolism: reference module in biomedical sciences; 2014. https://doi.org/10.1016/B978-0-12-801238-3.05146-1
Chen S, Duan G, Zhang R, Fan Q. Helicobacter pylori cytotoxin-associated gene A protein upregulates alpha-enolase expression via Src/MEK/ERK pathway: implication for progression of gastric cancer. Int J Oncol. 2014;45(2):764–70. https://doi.org/10.3892/ijo.2014.2444.
Li P, Hammes-Schiffer S. Substrate-to-product conversion facilitates active site loop opening in yeast enolase: a molecular dynamics study. ACS Catal. 2019;9(10):8985–90. https://doi.org/10.1021/acscatal.9b03249.
Thompson EJ, Paul A, Iavarone AT, Klinman JP. Identification of thermal conduits that link the protein-water interface to the active site loop and catalytic base in enolase. J Am Chem Soc. 2021;143(2):785–97. https://doi.org/10.1021/jacs.0c09423.
We extend our thanks to Dr. Bi-lin Dong for sharing the fungi and to Shanghai Biotree Biotech Co.,Ltd for proteomic data analysis.
This study was supported by the Hubei University of Chinese Medicine “Construction of graduate student innovation base (Graduate Workstation of Guangzhou Wellhealth Bio-Pharmaceutical Co., Ltd.)” Project 125 and Science and Technology Project of Hubei Provincial Department of Education for Young and Middle-aged scientists (Q20222009).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Library preparation for Transcriptome sequencing and transcriptomics data quality control and differential expression analysis. Method 2. Sample preparation of SDS-PAGE and data analysis. Fig. S1. Distribution of base sequencing error rates. Fig. S2. Distribution of ATCG content. Fig. S3. Muscari Unigene length Distribution Fig. S4. Comparative plot of the trimmed mean of M (TMM) density distribution. Fig. S5. Box line plot of TMM. Table S1. List of primer sequences for real-time PCR. Table S2. Active compounds of P. frutescens. Table S3. Diseases associated with Tiena and CTD ID. Table S4. The top five active compounds of P. frutescens. Table S5. Relationship between base score and error rate. Table S6. Sample sequencing data. Table S7. Assembly results. Table S8. Unigene annotated. Table S9. Differential genes. Table S10. Differential proteins.
About this article
Cite this article
Xu, Yd., Guo, Yj., Mao, Hr. et al. Integration of transcriptomics and proteomics to elucidate inhibitory effect and mechanism of rosmarinic acid from Perilla frutescens (L.) Britt. in treating Trichophyton mentagrophytes. Chin Med 18, 67 (2023). https://doi.org/10.1186/s13020-023-00772-2
- Perilla frutescens
- Rosmarinic acid
- Trichophyton mentagrophytes