Fangyu Xiao
Jihua Qiu
Ying Zhao1,3,4*- 1Guangzhou University of Chinese Medicine, Guangzhou, China
- 2South China Agricultural University, Guangzhou, China
- 3Department of Gynecology, First Affiliated Hospital of Guangzhou University of Chinese Medicine, Guangzhou, China
- 4First Clinical Medical College of Guangzhou University of Chinese Medicine, Guangzhou, China
Objective: This study focuses on whether vine tea contains potentially toxic components that trigger hepatotoxicity as a mechanism of action, which further provides some reference for the consumption and guides future product development of vine tea.
Methods: The chemical components of vine tea were collected from the reported literature and the toxicological information matched with the CTD database was collected, and the dataset of potential toxic components was established. The toxic components were submitted to the PharmMapper server to obtain potential targets. At the same time, the relevant targets were searched in the CTD database and GeneCards database with keywords such as “Hepatic Toxicity,” “Liver Damage,” and “Drug-induced liver injury.” After intersection, the potential hepatotoxic targets of vine tea were obtained. The protein interactions of potential hepatotoxic targets of vine tea were analyzed by the STRING database. Protein–protein interaction (PPI) networks were constructed by Cytoscape3.6.1 software. The GO molecular function and KEGG pathway of hepatotoxic targets were enriched by the R package to screen the key targets. The role of the components and key targets was analyzed by the LEDOCK program. The data from GEO database were mined for the functional correlation characterized by cell transcriptional expression caused by vine tea as a disturbance factor.
Results: This study has searched 34 potential toxic components and 57 potential hepatotoxic targets of vine tea, and the result showed that these targets were mainly involved in oxidative stress, cell metabolism, and apoptosis to affect the liver.
Conclusion: Vine tea has the interrelationship of multi-components, multi-targets, and multi-pathways. At the cellular level, the toxic components of vine tea, mainly flavonoids, may promote oxidative stress, promote oxidation to produce free radicals, guide apoptosis, and affect cell metabolism and other cytotoxic mechanisms. However, this hepatotoxicity is related to the dose, duration of vine tea, and individual differences. This study revealed the potential hepatotoxic components of vine tea and provides a reference for further research and development of related functional products.
Introduction
Vine tea is made from the tender leaves of Ampelopsis grossedentata (Hand. -Mazz.) W.T. Wang, a genus of Ampelopsis in the family Vitaceae, and the chemical composition of vine tea is mainly flavonoids and polyphenolic substances (Jia et al., 2021), with a total flavonoid content of up to 40% or more (He et al., 2000), represented by dihydromyricetin (30%) (Wang et al., 2014). Vine tea has been classified as a new food resource by the Chinese regulatory authority since 2013 and is generally taken orally as a raw herb for medicine. Vine tea is a typical tea-like plant in which the tender leaves are simply stir-fried during production to inhibit the action of intracellular enzymes, thus maximizing the retention of endogenous components of vine tea. Vine tea has been consumed in China for thousands of years and is considered a therapeutic tea or herbal tea in the Tujia, Yao, and Hakka regions of China for anti-inflammatory purposes (Tong et al., 2020), treatment of sore throat (Gao et al., 2009), and prevention of high fat (Fan et al., 2020).
There has been a misconception that Chinese herbal medicines are natural and less toxic, whereas chemically purified drugs are more toxic, when in fact, several have found that Chinese herbal medicines cause drug-induced liver injury (Sun et al., 2018; Zhao et al., 2018); that is, the chemical components of Chinese herbal medicines or their metabolites have toxic effects on the liver, and clinical data show that 42.9% of drug-induced liver injuries are caused by Chinese herbal medicines (Fengqin et al., 2011).
Presently, vine tea is becoming more popular among consumers worldwide because of its therapeutic functions and taste, but the assessment of its potential toxicity and drug interactions has not been definitional, and only few components are considered safe for consumption, such as dihydromyricetin (Zhang et al., 2021), so it is debatable whether vine tea has potential toxicological effects on humans after consumption, for it contains rich flavonoids. In the past two years, a case of a patient with multi-organ functional impairment due to vine tea poisoning admitted to the First Affiliated Hospital of Guangxi Medical University was reported (Faxin et al., 2020), which shows that the research on the mechanism of hepatotoxic effects of vine tea should be paid attention.
However, the present research on vine tea mainly focuses on chemical composition, nutritional composition, and pharmacology, and only some studies on toxicology are available. Moreover, traditional toxicological studies generally use animal models as the basis and pathology and immunology as the techniques to study the toxic targets of drugs, which not only have a long experimental period and high cost but also the correlation of the studied toxic targets is poor and it is more difficult to summarize the core targets. Until the concept of network toxicology (Fan et al., 2011) was proposed, the potential core targets, GO molecular functions, and KEGG signaling were analyzed by constructing herbal–component–target networks and combined with molecular docking techniques to verify the interactions between vine tea components and important targets, as well as by mining gene chips or RNA-seq to analyze the functions characterized by the transcriptional expression of vine tea components as perturbation factors. This study mainly used network toxicology and transcriptomics to analyze the potential toxic components of vine tea and reveal the potential mechanism of action for the occurrence of hepatotoxicity in vine tea, which further provides some reference for the consumption and guides future product development of vine tea.
In this study, first, a dataset of potentially toxic components of vine tea was established and these components were used to obtain binding targets through the PharmMapper server. On the other hand, the CTD database and the GeneCards database were searched to obtain targets with high correlation with keywords, such as “Hepatictoxicity,” “Liver Damage,” and “Drug-induced liver injury.” Next, potential hepatotoxic targets were obtained after the intersection of these two kinds of target collections. Then, the protein–protein interaction network was constructed for the aforementioned obtained hepatotoxic targets and acquiring the key targets with the toxic components of vine tea by combining GO and KEGG analysis. Finally, the relationships between the relevant toxic components of vine tea and the key targets were verified by a molecular docking technique and dynamics simulation. Subsequently, the RNA-seq expression profiles acquired from the experiment with hepatotoxic components as perturbation factors were further subjected to an enrichment analysis of the biological process (BP) module to observe the differences between the enrichment analysis results and the biological processes involved in the toxic targets acted by vine tea.
Materials and Methods
Establishment of the Dataset of Toxic Components of Vine Tea
The domestic and foreign literature (Wang et al., 1999; Du et al., 2004; Ming et al., 2015; Vieira Carneiro, 2016; Gao et al., 2017) was looked up and the chemical composition of vine tea was collected, and the chemical composition information was matched with the PubChem database (Kim et al., 2021) to obtain the 3D stereo structure and SMILES numbers of the compounds. The SMILES numbers were entered into the SWISS database (Daina et al., 2017) for predicting the ADMET properties of the compounds, and the compounds with intestinal absorption and compliance with Lipinski’s rule were screened as potential active components. The toxicological information of the active components was further queried through the CTD database (Davis et al., 2021) to establish a dataset of potentially toxic components of vine tea.
Prediction of Potential Targets of Vine Tea and Acquisition of Liver-Related Targets
The 3D stereo structure of vine tea toxicity components dataset was submitted to the PharmMapper server (Wang et al., 2017) for the prediction of relevant action targets. The prediction results were matched with the UniProt database (UniProt, 2021) to filter human-derived targets and annotate the target information as official symbols to obtain the potential action target set A of vine tea toxicity components. These keywords, such as “Hepatictoxicity,” “Liver Damage,” and “Drug-induced liver injury,” were imported into the GeneCards database (Stelzer et al., 2016) and CTD database to obtain the relevant targets of hepatotoxicity. The targets with a score value greater than 200 in the returned results of the GeneCards database and those with reported direct evidence in the returned results of the CTD database were used as filtering conditions, and the hepatotoxicity target set B was further established by deleting duplicate and false-positive genes. The predicted potential targets of vine tea toxicity components set A and the hepatotoxicity target set B were intersected to obtain the potential hepatotoxicity targets of vine tea toxicity components.
Construction of Protein–Protein Interaction Network for Hepatotoxic Targets
The potential hepatotoxic targets of vine tea were analyzed by the STRING online server (Szklarczyk et al., 2020) for protein–protein interactions (PPIs) and importing the results into Cytoscape 3.6.1 software for visualization analytics. CytoHubba and NetworkAnalyzer plug-ins were used to calculate the nodes and edges of the entire PPI network. First, the degree calculated by the CytoHubba plug-in was demonstrated by node color. Second, the degree of closeness centrality calculated by the NetworkAnalyzer plug-in was represented by node size. Third, the degree of the interaction between the nodes was represented by a circular-ranked approach which represents a further analysis of the importance of nodes in the network and key target proteins.
GO Enrichment Analysis and KEGG Pathway Enrichment Analysis
An enrichment analysis of GO molecule function and KEGG pathway was performed by the Bioconductor (www.bioconductor.org/) bioinformatics analysis package of R software (https://www.r-project.org/). Running the R 4.0.2 software, the R package (org.Hs.eg.db) was referenced, transforming potential hepatotoxic target gene names of vine tea toxic components for gene ID numbers, and the GO biomolecular functions and KEGG pathways of the targets were enriched by the clusterProfiler (Yu, 2018).
Molecular Docking and Molecular Dynamics Simulation
The top 25 important nodes of the PPI network and the genes with corrected p < 0.01 for KEGG pathway hits were used as filtering criteria to screen key target proteins for potential hepatotoxicity of vine tea. Protein crystal structures with more comprehensive information were searched in the PDB database (Burley et al., 2021), with human, inclusion of natural ligand molecules, and resolution ≤3.5 Å as the priority criteria for screening, and the stereo structures of protein crystals were downloaded. The LEDOCK program (Zhang and Zhao, 2016) was used to perform a semi-flexible molecular docking of vine tea toxic components with key target proteins, which the docking program set the receptor protein as a rigid molecule and the ligand as a flexible molecule, and the lowest energy conformation was further searched using simulated annealing and genetic algorithm crossover operations. Compared with the natural ligand of the protein crystal and the ligand with a significant effect reported in the literature, the RMSD value of the natural ligand after docking should be less than 2 Å to verify the accuracy of the docking program parameter settings. The output file of the docking task was sorted by the energy space conformation of ligands with a score system, and selecting the top-ranked group of receptors and ligands for molecular dynamics simulations to verify the reliability of the docking procedure. Using the lowest energy conformation after docking as the input file for molecular dynamics simulations, GROMACS software was used to perform fully flexible simulations of the protein and ligand molecules to further analyze the interaction between the protein and ligand molecules in the aqueous solvent system. The topology file of the protein and ligand molecules was constructed using the CHARMM36 force field (Huang et al., 2017), a dodecahedral unit box was defined for periodic simulations. Afterward, the water model of the TIP4P force field was used as the solvent to fill the unit box, while chloride or sodium ions were added to make the whole system neutral. Before running the kinetic simulation, the whole system was energy minimized, the system temperature was controlled to 300 K and pressure to 1bar, and after the system passed the NVT and NPT system equilibrium, the molecular dynamics simulation would be run for ten ns.
Transcriptome Analysis
The aforementioned analysis would obtain the important hepatotoxic components of vine tea. The expression profiles of the relevant microarray and RNA-seq were searched and downloaded from the GEO database using the components as perturbation factor, when expression profiles of compounds that have been reported to be hepatotoxic as perturbation factors were also obtained as reference. The expression profiles of different perturbers were grouped and analyzed for gene expression. Based on the highly conserved gene expression profiles of different species (Alam et al., 2020), if the expression profiles were for nonhuman species, the homologous transformation of genes was performed through the Ensembl database. The expression profiles were further subjected to a gene set enrichment analysis (GSEA), using the biological process (BP) of Gene Ontology (GO) as a functional module for enrichment, and the process was implemented through the gseGO of clusterProfiler. The enrichment results were obtained GO functions with significant levels by screening conditions of p < 0.05 and enrichment score >50. The GO functions enriched with vine tea components as perturbation factors were compared with the GO functions of reference expression profiles for similarity, respectively, and the correlation between GO functional modules was analyzed by the GOSemSim package. The similarity comparison used the GO function similarity algorithm explored by Wang et al. (2007). Generally speaking, the correlation of 0.6–0.8 indicates that both are similar, and the correlation >0.8 is considered as homology. The GO functions of vine tea components as perturbation factors and the GO function correlation >0.6 of the reference expression profile were screened as functional modules with similarity. The GO functions related to liver injury were selected for analysis. Hit gene sets corresponding to the similar modules were compared for gene function similarity to verify the functions characterized by transcriptional expression of different components of vine tea after perturbation of cells.
Results
Dataset of Toxic Components of Vine Tea
Searching the reported literature relating to vine tea, a total of 95 chemical components were collected, mainly including flavonoid, flavonoid glycoside, isoflavone, limonoid, fatty acid, and volatile component. The ADMET properties of the chemical components were predicted by the Swiss database, and 60 active components that conformed to human intestinal absorption and Lipinski’s rule were filtered out. Moreover, matching the toxic information of compounds in the CTD database, and 34 toxic components with direct evidence were screened out. The obtained information of potential toxic components of vine tea is exhibited in Table 1.
TABLE 1. Potentially toxic components of vine tea.
Potential Hepatotoxic Targets of the Toxic Components of Vine Tea
The 3D stereo structures of potential toxic components of vine tea were submitted to the PharmMapper server, and the system performed target prediction based on the principle of molecular reverse docking. The results returned 1,020 hit targets. After deleting nonhuman targets and duplicate values, a total of 551 potential targets were obtained, and annotating them as official names. A total of 699 targets related to hepatotoxicity, liver injury, and drug-induced liver injury were retrieved from the GeneCards database and CTD database. Fifty seven potential hepatotoxic targets of vine tea toxic components were obtained after intersecting with the action targets of vine tea toxic components.
Protein–Protein Interaction Network
Information on potential hepatotoxic targets of vine tea was submitted to the STRING database for analyzing protein–protein interactions (PPI). The topological properties of the PPI network allowed the analysis of the interaction mechanism of the potential hepatotoxic targets of vine tea. The results were imported into Cytoscape 3.6.1 software for visualization, as shown in Figure 1. The entire PPI network was constructed by 55 nodes and 308 edges. The node size was represented by the degree of closeness centrality, the node color was represented by degree, and the nodes of the entire network were ranked by the degree. The larger the node, the darker the color, the higher the ranking, the more important the node was in the network. The highest degree was ALB, which interacted with 38 targets (69%) in the network. Moreover there were also six interactive targets with more than 19 edges, namely TP53, VEGFA, CAT, MMP9, MPO, and PPARG. Those main functions are to regulate cell microenvironment, growth, proliferation, apoptosis, immunity, angiogenesis, and other physiological conditions, which are also important targets for the development of hepatotoxicity.
FIGURE 1. Hepatic toxicity targets protein–protein interaction network. Nodes are hepatotoxic targets; edges are target protein interactions.
Gene Ontology Enrichment Analysis and KEGG Pathway Enrichment Analysis
Gene Ontology Enrichment Analysis
The GO molecular functions of the potential hepatotoxic targets of vine tea were enriched by the ClusterProfiler, setting the molecular function p-value cutoff to 0.05 and Q-value cutoff to 0.05. The top 20 GO molecular functions were plotted in a bar graph after enrichment. From Figure 2, it can be seen that there were 19 GO molecular functions with p < 0.001, and all of these molecular functions were related to drug-induced liver injury, among which there were nine molecular functions with ≥6 hits. The coenzyme binding with the highest number of targets (19%) may lead to drug-induced liver injury or cirrhosis by binding to pyridoxal phosphate, flavin adenine dinucleotide, and nicotinamide adenine dinucleotide. In addition, molecular functions, such as heme binding, organic acid binding, fatty acid binding, carboxylic acid binding, oxygen binding, and heparin binding, may result in toxicity and damage to the liver by affecting the metabolism of drugs or alcohol by the liver, inducing the incidence of negative blood-related effects (Dadarkar et al., 2011), lipid peroxidation (Das et al., 2010), and activating inflammation. Moreover, direct effects on antioxidant activity and nuclear receptor activity may cause an imbalance in hepatocyte microenvironmental homeostasis or damage by free radicals in vivo by decreasing the superoxide dismutase activity and mineralocorticoid receptor activity.
FIGURE 2. GO enrichment analysis of hepatotoxic target.
KEGG Pathway Enrichment Analysis
The Kyoto encyclopedia of genes and genomes (KEGG) pathway enrichment analysis was calculated by a Bioconductor package for potential hepatotoxic targets of vine tea with a p-value cutoff of 0.05 and Q-value cutoff of 0.05. The enrichment results are shown in Table 2. The enrichment results obtained 27 signaling pathways, including 13 cancer pathways, four cell proliferation or differentiation or apoptosis- related pathways, one oxidative stress pathway, one drug resistance pathway, four immunoloregulation related pathways, and three cell metabolism- related pathways. Among these signaling pathways, endocrine resistance, fluid shear stress and atherosclerosis, HIF-1 signaling pathway, longevity regulating pathway, thyroid hormone signaling pathway, hepatitis C, glutathione metabolism, carbon metabolism, and arginine biosynthesis were proved that these maybe the main pathways of hepatotoxicity in vine tea.
TABLE 2. KEGG enrichment pathway of hepatotoxic target genes.
Molecular Docking and Molecular Dynamics Simulation Analysis Results
Sixteen key target proteins of potential hepatotoxicity of vine tea were screened by the total conditions, and the corresponding protein crystal structures were obtained by the PDB database. Among them, nine protein crystal structures with natural ligands were extracted and pretreated for the docking task, and the specific docking parameter information is shown in Table 3. The RMSD of natural ligands and corresponding proteins after docking were all within 0.2 nm, indicating that the pretreatment and parameter settings before docking had a certain accuracy and could be used for the virtual docking task between the toxic components of vine tea and protein crystals. The 34 toxic components of vine tea and 16 protein crystals were docking calculated. The docking of 16 positive ligands with active effects to the corresponding proteins as a control, the docking results were represented by the scoring (kcal/mol), and the scoring equation included the sum of several effects such as electrostatic interaction, Van der Waals’ force, hydrogen bonding contributions, intermolecular conflicts, and intramolecular conflicts of the ligands. Generally speaking, a docking system score of ≥4.25 indicates that the ligand interacts with the protein, a score of ≥5.0 indicates that the ligand has some binding activity with the active region of the protein, and a score of ≥7.0 indicates that the ligand has strong binding activity with the active region of the protein (Hsin et al., 2013). The results in Table 4 lists the docking scores of key target proteins with positive ligands and vine tea toxic components (scored with an absolute value ≥6.0), while the docking scores not listed had 231 groups with an absolute value ≥4.25, 151 groups with an absolute value ≥5.0, and only eight groups with an absolute value ≥7.0. The entire docking task had 560 groups. Nearly half of the docking tasks indicated that the ligands interacted with the target proteins, mainly focusing on 15 vine tea toxicity components and 16 key target proteins, which further indicated that vine tea may manifest hepatotoxicity through multi-component-multi-target effects. Further analyzed in Table 4, among the vine tea toxic components whose absolute values of docking scores were greater than the positive ligands were mainly dihydromyricetin, hesperitin, kaempferol, myricetin, phloretin, phloridzin, quercetin, and taxifolin. Among them, the absolute values of the docking scores of phloridzin and 16 target proteins were all greater than 5.0, and they had strong binding effects with the active regions of seven target proteins. The energy minimum conformation of its docking with 1DGF was used as the input file of GROMACS software, and the results are shown in Figure 3 after molecular dynamics simulations of 10 ns. The RMSD variation of protein residues or protein residues-ligands within 10 ns of a 5 nm expansion centered on the ligand is shown in Figure 3A. It was obvious that the ligand was stable with the active binding region of the protein, its RMSD values were within 0.2 nm and were stable above or below 0.10 nm. Figure 3B shows the variation of the hydrogen bonding value between the ligand and the protein residue in 10 ns. Its average value was 4, the bond length was less than 0.35 nm, and the bond angle was in the range of 30°, indicating that the preprocessing and parameter setting of the simulation task was more reasonable. Moreover, it also showed that the results of the ligand–protein interaction and the docking task were in agreement within 10 ns, which verified the reliability of the LEDOCK program for this study.
FIGURE 3. Results of 10 ns simulation of 1DGF and Phloridzin. Graph (A) indicates the RMSD variation of protein residue-ligand within 10 ns with the ligand as the center expanding 5 nm; graph (B) indicates the change in the value of hydrogen bonding generated by this ligand with protein residues within 10 ns.
TABLE 3. Molecular docking parameters.
TABLE 4. Molecular docking of potentially toxic components of vine tea.
Transcriptome Analysis Results
The results of the aforementioned analysis obtained eight major toxic components of vine tea. When these components were used as perturbation factors to search the GEO database for relevant expression profiles, only a few components’ expression profiles such as quercetin, taxifolin, phloridzin, and kaempferol were found, whose database platform numbers were GSE62805, GSE59704, GSE38138, and GSE145665, respectively. At the same time, referring to the expression profile no. GSE5789, the experimental protocol for this data was obtained with the transcriptional data using toxic drugs (2,3,7,8-tetrachlorodibenzo-p-dioxin, 3,3′,4,4′,5-pentachlorobiphenyl, 2,2′,4,4′,5,5′-hexachlorobiphenyl, and 2,3,4,7,8-pentachlorodibenzofuran) as perturbers that caused liver damage. After the expression profiles were enriched for the GO functional biological process (BP) modules separately, the BP modules corresponding to the enrichment of vine tea components were compared with the BP modules enriched for the reference expression profiles for similarity. Eight similar common BP modules were found, namely gland growth, angiogenesis, morphogenesis of blood vessels, epithelial cell differentiation, development of circulatory system, sprouting angingenesis, and muscle organ growth, which were more consistent with the gene functions of the important targets analyzed previously, and the results are shown in Figure 4A. The BP modules associated with liver injury were further analyzed, and the gene sets hit by the BP modules of expression profiles and the corresponding similar BP modules of various components were compared for similarity. The gene sets of GO functional similar modules with correlation coefficients greater than 0.6 were found to have a functional similarity in the range of 0.6–0.9, indicating high similarity between the BP modules enriched by the transcriptional expression of vine tea components for the treatment of perturbation factors and the characteristics of transcriptional expression of liver injury, and the results are represented in Table 5. This revealed that these four vine tea components may promote biological processes at the cellular level such as the response to reactive oxygen species, response to toxic components, response to oxidative stress, lipid metabolism, and apoptosis. Subsequently, the genes hit by these highly similar modules were intersected with the hepatotoxic targets obtained previously, which obtained seven common target genes, such as HMOX1, PPARG, PTPRC, G6PD MMP9, MMP2, and GSTM1, as shown in Figure 4B. Further validating the reliability of the screening of the aforementioned key targets, and suggesting that vine tea components played an important role in angiogenesis, cell proliferation, and cell differentiation.
FIGURE 4. Venn diagram of transcriptome analysis. Graph (A) indicates the distribution of GO functions enriched by vine tea components as perturbation factors similar to the GO functions enriched by the reference expression profile; graph (B) indicates the intersection of transcriptional targets hit by hepatotoxic targets and liver injury-related GO functions.
TABLE 5. Gene set correlation and hit count of GO similar functional modules.
Discussion
In this study, we analyzed the GO molecular functions and KEGG pathway of the same targets by integrating the action targets of the toxic components of vine tea and the hepatotoxicity-related targets and validated them by molecular docking techniques. By mining the expression profile of GEO, we analyzed the toxicological effects of the GO functional modules enriched by the transcriptional expression of vine tea components as perturbation factors at the cellular level, to predict the potential hepatotoxic action mechanism of vine tea. The liver is involved in the physiological processes of body metabolism, hormone balance, and toxin removal, and is an essential organ of the human body. Most drugs are absorbed orally through the gastrointestinal tract and then metabolized by the liver and excreted from the kidneys; if a few drugs are metabolized by the liver to form end-toxicants or transient active substances, then the drugs are exposed to the liver will lead to pharmacologic liver damage. Therefore, the active components of drugs should be actively concerned whether they can cause hepatotoxicity while exerting their medicinal effects, and some active components may cause hepatotoxicity through direct cellular toxicity, oxidative stress, immune response, and cellular metabolism (Bleibel et al., 2007). Although some drugs are relatively safe, they may also cause liver injury (Galan et al., 2005). There are reports in the literature that Chinese herb medicines can cause liver injury (S et al., 2014). Earlier literature (Yang et al., 2006) demonstrated that the consumption of biological doses of vine tea was safe through safety toxicology experiments in rats, however, the available literature (Sahu and Gray, 1996; Dzoyem et al., 2013; Skibola and Smith, 2000; Galati and O'Brien, 2004) showed that at certain doses, durations, and individual differences in the effects of vine tea, kaempferol, quercetin, myricetin, taxifolin, dihydromyricetin, physcion, dihydrokaempferol, phloretin, and other components can induce mutagenic cytotoxicity.
The role of chemical components such as flavonoid substances and polyphenols in vine tea in vivo is complex. Although there is a large body of literature supporting the antioxidant function of dietary polyphenol or flavonoid components in vivo, the antioxidant activity of these substances in vivo can be expressed by the substances themselves either directly as reactive oxygen radical scavengers or indirectly as regulators of intracellular antioxidant enzymes and pro-oxidant enzymes, when involved in peroxidases containing systems with redox metals, these substances have been shown to act as pro-oxidants, directing the formation of reactive oxygen radicals (ROS) and phenoxyl radicals, further promoting the destruction of lipids, DNA, or their biomolecules (Decker, 1997), and increasing the concentrations of substances, which can mediate the breakage of cellular DNA and activate the cellular caspase apoptotic pathway (Hadi et al., 2000), leading to non-specific necrosis of cells and cause cytotoxicity. For example, quercetin, a toxic component of vine tea with a catechol B ring, is absorbed orally by the body and can be distributed to liver organs (de Boer et al., 2005), where it is oxidized and metabolized to o-quinone by tyrosine oxidase and catalase in hepatocytes (Oyama et al., 1999) and bound to macromolecules as an electrophilic body, and may also react with O2 redox to generate O2—or OH radicals, and these oxidized quinones or quinone methyl intermediates may also react with glutathione to form adducts (Galati et al., 2001) when the resulting cytotoxicity is related to the quercetin-dependent concentration. It has also been shown that the opposite occurs, as ROS in hepatocytes metabolizes flavonoids such as quercetin into products that covalently bind to cellular DNA and proteins (Walle et al., 2003), and similarly, kaempferol and myricetin has been shown to cause ROS formation in systems containing peroxidases involved in the presence of redox metals (Galati and O'Brien, 2004). In contrast, the transcriptional analysis indicated that four components of vine tea characterize the cellular response to reactive oxygen species, response to oxidative stress, response to toxic substances, and BP function in response to injury at the cellular level, which is consistent with the aforementioned analysis.
Another mechanism of cytotoxicity that causes apoptosis is that quercetin and myricetin directly affect the mitochondrial membrane potential in a pro-oxidant effect with the involvement of transition metals (Said Ahmad et al., 1992), and this effect leads to the rupture of the mitochondrial membrane and the release of cytochrome C and triggers the release of mitochondrial toxicity, which is an irreversible marker of cellular mortality, and that the toxic components of vine tea, represented by palmitic acid, represented by saturated fatty acids, induces lipid apoptosis (Jinju et al., 2016). The present work demonstrated that the interaction of vine tea toxic components with superoxide dismutase (SOD2) and catalase (CAT) may affect the antioxidant activity of hepatocytes and disrupt the intracellular balance of scavenging toxic substances such as superoxide anion radicals, indirectly leading to the accumulation of toxins and attacking intracellular biomolecules and inducing apoptosis pathway. In contrast, 11 (32.35%) of the vine tea toxic components act on superoxide dismutase and 15 (44.12%) on catalase, so vine tea may affect the activity of superoxide dismutase and catalase downstream of the longevity regulating pathway, reducing the detoxification capacity of hepatocytes and increasing the potential of toxins to damage cells. The transcriptional analysis revealed that the enriched BP module was highly similar in the regulation of the apoptotic process, cytokine stimulation response, and regulation of protein serine/threonine kinase activity, which verified that the characteristics of this component as a perturber- induced cell expression were consistent with the behavior of the apoptotic mechanism.
In addition, this study demonstrated that vine tea toxic components can act on heme oxygenase 1 (HMOX1), glutathione S-transferase Mu 1 (GSTM1), whose involvement in phase II metabolism is considered as a pathway for cellular clearance of exogenous toxic substances, and it is important to note that inhibition of these enzymes may lead to accumulation of exogenous toxic substances. Among them affecting the glutathione metabolic pathway can lead to depletion of glutathione, thus affecting the content of components such as bilirubin and glutathione, which is consistent with the literature reporting that flavonoids cause glutathione oxidation in rat hepatocytes (Galati et al., 2002). On the other hand, glutathione S-transferase P is closely linked to CYP450 in the signaling pathway for the metabolic action of exogenous compounds, and the CYP enzyme (CYP2E1) acting on toxic components of vine tea is involved in molecular functions such as heme binding, tetrapyrrole binding, sulfide binding, oxygen binding, NADP binding, and oxidoreductase activity, and the functions of the transcriptome- enriched BP module in cellular detoxification, reaction to inorganic and organic substances, fatty acid and lipid metabolism indicate that the toxic components of vine tea can perturb cellular bio metabolism and may affect the metabolism of exogenous compounds in hepatocytes and prolong the duration of pharmacological action of toxic components and bioactive substances of vine tea, increasing the chance of Drug-drug interaction (DDI), which in turn causes some toxic side effects in the liver (Jana and Paliwal, 2007). Here, it should be noted that the physical differences of different people make the specific reactions triggered by the differences in the expression and catalytic activity of metabolic enzymes such as CYP enzymes in sensitive individuals, and most of them can tolerate the dose and duration of drugs well, only for a very small number of physically sensitive healthy individuals hide dangerous or even fatal signals. Both network toxicology and transcriptomic analyses support the significant relevance of vine tea components in angiogenesis and cell differentiation through the interactions of vascular endothelial growth factor A (VEGFA), neurogenic locus notch homolog protein 1 (NOTCH1), 72 kDa type IV collagenase (MMP2), catenin beta-1 (CTNNB1), receptor tyrosine-protein kinase erbB-2 (ERBB2), and other proteins interact with each other, whether this effect predicts that the toxic components of vine tea can block bile entry into the bile ducts and cause bile stasis by affecting vascular and peripheral capillary vessels is to be verified by further experiments. In conclusion, the types of liver damage induced by drugs vary, depending on the nature, dosage, and duration of the drug itself, and the toxic components of vine tea have both significant pharmacological activity and potential for cytotoxicity, so it is of some significance to explore the mechanism of hepatotoxic effects of vine tea in depth.
In summary, this work demonstrates that vine tea has a multi-component, multi-target, and multi-pathway interrelationships and that the flavonoid-based toxic components of vine tea may promote oxidative stress and pro-oxidative generation of free radicals and attack the mitochondria at the cellular level, causing cytotoxicity, triggering apoptotic mechanisms, or affecting the metabolic activities of cells. These adverse effects are dose-dependent toxic manifestations, and although the general population has a good degree of tolerance at certain doses and durations, the potential pitfalls hidden in a few healthy individuals with sensitive constitutions should not be overlooked. Such herb tea beverages do not require SFDA approval for market circulation, and their therapeutic functions and taste are increasingly favored by consumers worldwide, and the assessment of potential toxicity and drug interactions has not been definitional, so mixing and abuse with other herbs should be avoided, and it is recommended that a safe dosage standard for consumption of herbal-based tea beverages be established. This work is based on database mining and analysis to investigate whether vine tea is potentially hepatotoxic using a network toxicology and molecular dynamics approach. Although the method is becoming increasingly popular, its accuracy and reliability are related to the subjective choice of the user, together with drawbacks such as the limitations of the database, which require the user to have a deep basic knowledge of network toxicology and molecular dynamics theoretical approach. However, it also illustrates the important contribution of the method to the development of biology, which will become a core research technique in modern biology and further generate new interdisciplinary disciplines. Finally, the metabolic indexes and action targets predicted in this paper are informative but do not fully reflect the overall profile of vine tea. Because of the potential applications in material and pharmaceutical sciences for constituents containing therapeutic functions such as dihydromyricetin, vine tea is of increasing interest to researchers. Further follow-up studies and validation of the toxicological effects, serum pharmacology, and serum medicinal chemistry of various constituents of vine tea are essential to guide future product development and the possibility of reducing the cost of new formulations.
Data Availability Statement
The datasets supporting the conclusions of this article are available in a public database from PubChem, SMILSS, CTD, PharmMapper, UniProt, GeneCards, STRING, GEO. The accession numbers of gene chip are GSE62805, GSE59704, GSE38138, GSE145665, GSE5789.
Author Contributions
All authors made a significant contribution to the work reported, whether that is in the conception, study design, software, acquisition of data, analysis, and interpretation, or in all these areas; took part in drafting, revising or critically reviewing the article; gave final approval of the version to be published; have agreed on the journal to which the article has been submitted; and agree to be accountable for all aspects of the work.
Funding
This work was financially supported by the National Academic School of Chinese Medicine Inheritance Studio Construction Project of China (National Chinese Medicine Human Education Letter (2019) No. 62), Science and Technology Program Project of Guangdong Province, China (No. 2016A020215180), and Luo Yuankai Zhi Kidney Fertility Pill Middle Youth Research Fund Project of China (No. 20190803).
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.
References
Alam, T., Agrawal, S., Severin, J., Young, R. S., Andersson, R., Arner, E., et al. (2020). Comparative Transcriptomics of Primary Cells in Vertebrates. Genome Res. 30 (7), 951–961. doi:10.1101/gr.255679.119
Bleibel, W., Kim, S., D'Silva, K., Lemmer, E. R., D Silva, K., D Silva, K., et al. (2007). Drug-Induced Liver Injury: Review Article. Dig. Dis. Sci. 52 (10), 2463–2471. doi:10.1007/s10620-006-9472-y
Burley, S. K., Bhikadiya, C., Bi, C., Bittrich, S., Chen, L., Crichlow, G. V., et al. (2021). RCSB Protein Data Bank: Powerful New Tools for Exploring 3D Structures of Biological Macromolecules for Basic and Applied Research and Education in Fundamental Biology, Biomedicine, Biotechnology, Bioengineering and Energy Sciences. Nucleic Acids Res. 49 (D1), D437–D451. doi:10.1093/nar/gkaa1038
Dadarkar, S. S., Fonseca, L. C., Mishra, P. B., Lobo, A. S., Doshi, L. S., Dagia, N. M., et al. (2011). Phenotypic and Genotypic Assessment of Concomitant Drug-Induced Toxic Effects in Liver, Kidney and Blood. J. Appl. Toxicol. 31 (2), 117–130. doi:10.1002/jat.1562
Daina, A., Michielin, O., and Zoete, V. (2017). SwissADME: a Free Web Tool to Evaluate Pharmacokinetics, Drug-Likeness and Medicinal Chemistry Friendliness of Small Molecules. Sci. Rep. 7 (1), 42717. doi:10.1038/srep42717
Das, J., Ghosh, J., Manna, P., and Sil, P. C. (2010). Acetaminophen Induced Acute Liver Failure via Oxidative Stress and JNK Activation: Protective Role of Taurine by the Suppression of Cytochrome P450 2E1. Free Radic. Res. 44 (3), 340–355. doi:10.3109/10715760903513017
Davis, A. P., Grondin, C. J., Johnson, R. J., Sciaky, D., Wiegers, J., Wiegers, T. C., et al. (2021). Comparative Toxicogenomics Database (CTD): Update 2021. Nucleic Acids Res. 49 (D1), D1138–D1143. doi:10.1093/nar/gkaa891
de Boer, V. C., Dihal, A. A., van der Woude, H., Arts, I. C., Wolffram, S., Alink, G. M., et al. (2005). Tissue Distribution of Quercetin in Rats and Pigs. J. Nutr. 135 (7), 1718–1725. doi:10.1093/jn/135.7.1718
Decker, E. A. (1997). Phenolics: Prooxidants or Antioxidants? Nutr. Rev. 55 (11), 396–398. doi:10.1111/j.1753-4887.1997.tb01580.x
Du, Q., Chen, P., Jerz, G., and Winterhalter, P. (2004). Preparative Separation of Flavonoid Glycosides in Leaves Extract of Ampelopsis Grossedentata Using High-Speed Counter-current Chromatography. J. Chromatogr. A. 1040 (1), 147–149. doi:10.1016/j.chroma.2004.03.062
Dzoyem, J. P., Hamamoto, H., Ngameni, B., Ngadjui, B. T., and Sekimizu, K. (2013). Antimicrobial Action Mechanism of Flavonoids from Dorstenia Species. Drug Discov. Ther. 7 (2), 66–72. doi:10.5582/ddt.2013.v7.2.66
Fan, L., Qu, X., Yi, T., Peng, Y., Jiang, M., Miao, J., et al. (2020). Metabolomics of the Protective Effect of Ampelopsis Grossedentata and its Major Active Compound Dihydromyricetin on the Liver of High-Fat Diet Hamster. Evid. Based Complement. Alternat Med. 2020, 3472578. doi:10.1155/2020/3472578
Fan, X., Zhao, X., Jin, Y., Shen, X., and Liu, C. (2011). Network Toxicology and its Application to Traditional Chinese Medicine. Zhongguo Zhong Yao Za Zhi 36 (21), 2920–2922. doi:10.4268/cjcmm20112104
Faxin, L., Changyi, X., Zong, N., and Jijin, Z. (2020). A Case of Multi-Organ Dysfunction Caused by Vine tea tea in Guangxi. Chin. J. Gen. Pract. 18 (02), 335–337. doi:10.16766/j.cnki.issn.1674-4152.001240
Fengqin, H., Na, H., Xia, L., Tailing, W., and Guiqiang, W. (2011). Clinical Characteristics of 33 Patients with Drug Induced Chronic Liver Injury. Chin. J. Pract. Intern. 31 (12), 941–943. doi:CNKI:SUN:SYNK.0.2011-12-019
Galan, M. V., Potts, J. A., Silverman, A. L., and Gordon, S. C. (2005). The burden of Acute Nonfulminant Drug-Induced Hepatitis in a United States Tertiary Referral center [corrected]. J. Clin. Gastroenterol. 39 (1), 64–67. doi:10.1097/mcg.0000148380.81683.74
Galati, G., Moridani, M. Y., Chan, T. S., and O'Brien, P. J. (2001). Peroxidative Metabolism of Apigenin and Naringenin versus Luteolin and Quercetin: Glutathione Oxidation and Conjugation. Free Radic. Biol. Med. 30 (4), 370–382. doi:10.1016/s0891-5849(00)00481-0
Galati, G., and O'Brien, P. J. (2004). Potential Toxicity of Flavonoids and Other Dietary Phenolics: Significance for Their Chemopreventive and Anticancer Properties. Free Radic. Biol. Med. 37 (3), 287–303. doi:10.1016/j.freeradbiomed.2004.04.034
Galati, G., Sabzevari, O., Wilson, J. X., and O'Brien, P. J. (2002). Prooxidant Activity and Cellular Effects of the Phenoxyl Radicals of Dietary Flavonoids and Other Polyphenolics. Toxicology 177 (1), 91–104. doi:10.1016/S0300-483X(02)00198-1
Gao, J., Liu, B., Ning, Z., Zhao, R., Zhang, A., and Wu, Q. (2009). CHARACTERIZATION AND ANTIOXIDANT ACTIVITY OF FLAVONOID-RICH EXTRACTS FROM LEAVES OF AMPELOPSIS GROSSEDENTATA. J. Food Biochem. 33 (6), 808–820. doi:10.1111/j.1745-4514.2009.00253.x
Gao, Q., Ma, R., Chen, L., Shi, S., Cai, P., Zhang, S., et al. (2017). Antioxidant Profiling of Vine tea (Ampelopsis Grossedentata): Off-Line Coupling Heart-Cutting HSCCC with HPLC-DAD-QTOF-MS/MS. Food Chem. 225, 55–61. doi:10.1016/j.foodchem.2016.11.122
Hadi, S. M., Asad, S. F., Singh, S., and Ahmad, A. (2000). Putative Mechanism for Anticancer and Apoptosis-Inducing Properties of Plant-Derived Polyphenolic Compounds. Iubmb Life 50 (3), 167–171. doi:10.1080/152165400300001471
He, G. X., Pei, G., Zhou, T. D., and Zhou, X. X. (2000). Determination of Total Flavonoids and Dihydromyricetin in Ampelopsis Grossedentala (Hand-Mazz) W. T. Wang. Zhongguo Zhong Yao Za Zhi 25 (7), 423–425.
Hsin, K. Y., Ghosh, S., and Kitano, H. (2013). Combining Machine Learning Systems and Multiple Docking Simulation Packages to Improve Docking Prediction Reliability for Network Pharmacology. PLoS One 8 (12), e83922. doi:10.1371/journal.pone.0083922
Huang, J., Rauscher, S., Nawrocki, G., Ran, T., Feig, M., de Groot, B. L., et al. (2017). CHARMM36m: an Improved Force Field for Folded and Intrinsically Disordered Proteins. Nat. Methods 14 (1), 71–73. doi:10.1038/nmeth.4067
Jana, S., and Paliwal, J. (2007). Molecular Mechanisms of Cytochrome P450 Induction: Potential for Drug-Drug Interactions. Curr. Protein Pept. Sci. 8 (6), 619–628. doi:10.2174/138920307783018668
Jia, C., Li, J., Zhang, M., Ma, W., Zhao, S., Liu, R., et al. (2021). Antioxidant Properties of the Extracts of Vine tea (Ampelopsis Grossedentata) with the Different Color Characteristics and Inhibition of Rapeseed and sunflower Oil Oxidation. LWT 136, 110292. doi:10.1016/j.lwt.2020.110292
Jinju, W., Chunyan, Z., Bin, X., Hong, L., Jianxiong, J., Chunhong, F., et al. (2016). The Roles of GRP78 in Inhibiting PA-mediated Hepatocellular Carcinoma Cells Apoptosis. Chin. J. Gerontol. 36 (9), 2056–2059. doi:10.3969/j.issn.1005-9202.2016.09.003
Kim, S., Chen, J., Cheng, T., Gindulyte, A., He, J., He, S., et al. (2021). PubChem in 2021: New Data Content and Improved Web Interfaces. Nucleic Acids Res. 49 (D1), D1388–D1395. doi:10.1093/nar/gkaa971
Ming, F. U., Xiao-Ying, L. I., Deng-Yu, W., and Meng-Qi, G. (2015). Flavonoid Constituents of Leaves of Ampelopsis Grossedentata (Hand-Mazz)W.T.Wang. Chin. Pharm. J. 50 (7), 574–578.
Oyama, Y., Noguchi, S., Nakata, M., Okada, Y., Yamazaki, Y., Funai, M., et al. (1999). Exposure of Rat Thymocytes to Hydrogen Peroxide Increases Annexin V Binding to Membranes: Inhibitory Actions of Deferoxamine and Quercetin. Eur. J. Pharmacol. 384 (1), 47–52. doi:10.1016/S0014-2999(99)00649-4
S, G., Xm, S., Y, X., Bs, Z., and Wl, L. (2014). Research on Medicine - Induced Liver Toxicity Mechanisms. J. Harbin Univ. Commerce(Natural Sci. Edition) 30 (03), 257–259+270. doi:10.19492/j.cnki.1672-0946.2014.03.001
Sahu, S. C., and Gray, G. C. (1996). Pro-oxidant Activity of Flavonoids: Effects on Glutathione and Glutathione S-Transferase in Isolated Rat Liver Nuclei. Cancer Lett. 104 (2), 193–196. doi:10.1016/0304-3835(96)04251-6
Said Ahmad, M., Fazal, F., Rahman, A., Hadi, S. M., and Parish, J. H. (1992). Activities of Flavonoids for the Cleavage of DNA in the Presence of Cu(II): Correlation with Generation of Active Oxygen Species. Carcinogenesis 13 (4), 605–608. doi:10.1093/carcin/13.4.605
Skibola, C. F., and Smith, M. T. (2000). Potential Health Impacts of Excessive Flavonoid Intake. Free Radic. Biol. Med. 29 (3-4), 375–383. doi:10.1016/s0891-5849(00)00304-x
Stelzer, G., Rosen, N., Plaschkes, I., Zimmerman, S., Twik, M., Fishilevich, S., et al. (2016). The GeneCards Suite: from Gene Data Mining to Disease Genome Sequence Analyses. Curr. Protoc. Bioinformatics 54 (1), 1–33. doi:10.1002/cpbi.5
Sun, W., Zeng, C., Liu, S., Fu, J., Hu, L., Shi, Z., et al. (2018). Ageratina Adenophora Induces Mice Hepatotoxicity via ROS-NLRP3-Mediated Pyroptosis. Sci. Rep. 8 (1), 16032. doi:10.1038/s41598-018-34492-7
Szklarczyk, D., Gable, A. L., Nastou, K. C., Lyon, D., Kirsch, R., Pyysalo, S., et al. (2020). The STRING Database in 2021: Customizable Protein-Protein Networks, and Functional Characterization of User-Uploaded Gene/measurement Sets. Nucleic Acids Res. 49 (D1), D605–D612. doi:10.1093/nar/gkaa1074
Tong, H., Zhang, X., Tan, L., Jin, R., Huang, S., and Li, X. (2020). Multitarget and Promising Role of Dihydromyricetin in the Treatment of Metabolic Diseases. Eur. J. Pharmacol. 870, 172888. doi:10.1016/j.ejphar.2019.172888
UniProt (2021). UniProt: the Universal Protein Knowledgebase in 2021. Nucleic Acids Res. 49 (D1), D480–D489. doi:10.1093/nar/gkaa1100
Vieira Carneiro, R. C. (2016). Volatile Compounds in Vine Tea (Ampelopsis grossedentata). [dissertation/master’s thesis]. Virginia: Virginia Tech.
Walle, T., Vincent, T. S., and Walle, U. K. (2003). Evidence of Covalent Binding of the Dietary Flavonoid Quercetin to DNA and Protein in Human Intestinal and Hepatic Cells. Biochem. Pharmacol. 65 (10), 1603–1610. doi:10.1016/S0006-2952(03)00151-5
Wang, D., Liu, J., Lu, J., and Zheng, S. (1999). ChemInform Abstract: Two New Limonoids of Ampelopsis grossedentata Hand.—Mazz. ChemInform. 30 (27), 240–242. doi:10.1002/chin.199927203
Wang, J., He, L., Zheng, N., Kong, Q., Yu, Q., and Yu, Z. (2014). Dihydromyricetin in Ampelosis Grossedentata Collected from Different Habitats. Chin. traditional patent Med. 36 (1), 145–147. doi:10.3969/j.issn.1001-1528.2014.01.035
Wang, J. Z., Du, Z., Payattakool, R., Yu, P. S., and Chen, C. F. (2007). A New Method to Measure the Semantic Similarity of GO Terms. Bioinformatics 23 (10), 1274–1281. doi:10.1093/bioinformatics/btm087
Wang, X., Shen, Y., Wang, S., Li, S., Zhang, W., Liu, X., et al. (2017). PharmMapper 2017 Update: a Web Server for Potential Drug Target Identification with a Comprehensive Target Pharmacophore Database. Nucleic Acids Res. 45 (W1), W356–W360. doi:10.1093/nar/gkx374
Yang, W., Ping, L., and Ming-Hui, D. (2006). Evaluation on Immuno Toxicology of Feng Ming Ampelopsis Grossedentata Exstract. Mod. Prev. Med. 33 (07), 1100–1101. doi:10.3969/j.issn.1003-8507.2006.07.017
Yu, G. (2018). Using Meshes for MeSH Term Enrichment and Semantic Analyses. Bioinformatics 34 (21), 3766–3767. doi:10.1093/bioinformatics/bty410
Zhang, H., Caprioli, G., Hussain, H., Khoi Le, N. P., Farag, M. A., and Xiao, J. (2021). A Multifaceted Review on Dihydromyricetin Resources, Extraction, Bioavailability, Biotransformation, Bioactivities, and Food Applications with Future Perspectives to Maximize its Value. eFood 2 (4), 164–184. doi:10.53365/efood.k/143518
Zhang, N., and Zhao, H. (2016). Enriching Screening Libraries with Bioactive Fragment Space. Bioorg. Med. Chem. Lett. 26 (15), 3594–3597. doi:10.1016/j.bmcl.2016.06.013
Keywords: vine tea, hepatotoxicity, network toxicology, molecular docking, molecular dynamics, GEO database
Citation: Xiao F, Qiu J and Zhao Y (2022) Exploring the Potential Toxicological Mechanisms of Vine Tea on the Liver Based on Network Toxicology and Transcriptomics. Front. Pharmacol. 13:855926. doi: 10.3389/fphar.2022.855926
Received: 16 January 2022; Accepted: 01 March 2022;
Published: 22 March 2022.
Edited by:
Sajjad Gharaghani, University of Tehran, IranReviewed by:
Karim Abbasi, Sharif University of Technology, IranJen-Tsung Chen, National University of Kaohsiung, Taiwan
Copyright © 2022 Xiao, Qiu and Zhao. 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: Ying Zhao, ZHJ6aGFveWluZ0AxMjYuY29t
†These authors have contributed equally to this work and share first authorship