Skip to main content

ORIGINAL RESEARCH article

Front. Pharmacol., 21 September 2018
Sec. Cardiovascular and Smooth Muscle Pharmacology

Luteolin Inhibits Vascular Smooth Muscle Cell Proliferation and Migration by Inhibiting TGFBR1 Signaling

\r\nYu-Ting Wu,Yu-Ting Wu1,2Ling Chen,Ling Chen1,2Zhang-Bin Tan,Zhang-Bin Tan1,2Hui-Jie Fan,Hui-Jie Fan1,2Ling-Peng Xie,Ling-Peng Xie1,2Wen-Tong Zhang,Wen-Tong Zhang1,2Hong-Mei Chen,Hong-Mei Chen1,2Jun Li,Jun Li1,2Bin Liu,*Bin Liu2,3*Ying-Chun Zhou,*Ying-Chun Zhou1,2*
  • 1School of Traditional Chinese Medicine, Southern Medical University, Guangzhou, China
  • 2Department of Traditional Chinese Medicine, Nanfang Hospital, Southern Medical University, Guangzhou, China
  • 3Guangzhou Institute of Cardiovascular Disease, The Second Affiliated Hospital of Guangzhou Medical University, Guangzhou, China

Vascular smooth muscle cell (VSMC) proliferation and migration play a critical role in the development of arterial remodeling during various vascular diseases including atherosclerosis, hypertension, and related diseases. Luteolin is a food-derived flavonoid that exerts protective effects on cardiovascular diseases. Here, we investigated whether transforming growth factor-β receptor 1 (TGFBR1) signaling underlies the inhibitory effects of luteolin on VSMC proliferation and migration. We found that luteolin reduced the proliferation and migration of VSMCs, specifically A7r5 and HASMC cells, in a dose-dependent manner, based on MTS and EdU, and Transwell and wound healing assays, respectively. We also demonstrated that it inhibited the expression of proliferation-related proteins including PCNA and Cyclin D1, as well as the migration-related proteins MMP2 and MMP9, in a dose-dependent manner by western blotting. In addition, luteolin dose-dependently inhibited the phosphorylation of TGFBR1, Smad2, and Smad3. Notably, adenovirus-mediated overexpression of TGFBR1 enhanced TGFBR1, Smad2, and Smad3 activation in VSMCs and partially blocked the inhibitory effect of luteolin on TGFBR1, Smad2, and Smad3. Moreover, overexpression of TGFBR1 rescued the inhibitory effects of luteolin on the proliferation and migration of VSMCs. Additionally, molecular docking showed that this compound could dock onto an agonist binding site of TGFBR1, and that the binding energy between luteolin and TGFBR1 was -10.194 kcal/mol. Simulations of molecular dynamics showed that TGFBR1-luteolin binding was stable. Collectively, these data demonstrated that luteolin might inhibit VSMC proliferation and migration by suppressing TGFBR1 signaling.

Introduction

Atherosclerosis is a chronic inflammatory disease of the arteries and is the leading cause of morbidity and mortality worldwide (Weber and von Hundelshausen, 2017). VSMCs are one of the principal components of the arterial walls. Recent evidence indicates that abnormal VSMC proliferation and migration play a vital role in the pathologies and development of AS (Lacolley et al., 2012). Therefore, suppressing the proliferation and migration of these cells could play an important role in preventing the pathological process of AS and might become a novel therapeutic strategy.

Luteolin is a common flavonoid that is part of a group of potentially chemoprotective compounds that widely exist in many Huoxue Huatan Chinese herbal medicines (Radix Salviae, Carthami Flos, Spatholobus suberectus Dunn, Platycodon grandiflorus, Inulae Flos, Asteris Radix Et Rhizoma, etc.), vegetables (carrot, cabbage, artichoke, and celery), fruits (apple), and plant-based beverages (tea) (Wang et al., 2009; Luo et al., 2017). Preclinical in vitro and in vivo studies have provided reasonable evidence of the cardioprotective effects of luteolin (Lopez-Lazaro, 2009). For example, experimental studies have found that it has multiple pharmacological effects including anti-inflammatory, anti-oxidant, anti-proliferative, anti-migrative, anti-apoptotic, and other pharmacological activities, most of which participate in the pathological process of AS and mediate CVD protection (Jia et al., 2015; Nile et al., 2018). Moreover, previous studies have demonstrated that luteolin treatment inhibits the proliferation of VSMCs (Lang et al., 2012; Jiang et al., 2013; Xu et al., 2015). However, the underlying mechanisms associated with its inhibitory effects on VSMC proliferation and migration are still unclear.

TGFBR1, a specific TGF-β receptor, mediates the physiological and pathological function of TGF-β, which was demonstrated to be a major component of the atherosclerotic process (Serralheiro et al., 2017). TGFBR1 is essential for mediating TGF-β-induced effects on nearly all cell types of the cardiovascular system, including endothelial cells, VSMCs, myofibroblasts, and macrophages (Goumans and Ten Dijke, 2018). It was reported that a TGFBR1 inhibitor suppressed intimal hyperplasia in a mouse model of transplant arteriosclerosis and reduced VSMC proliferation in vitro (Sun et al., 2014). Thus, targeting TGFBR1 might represent a strategy for AS treatment. However, it remains to be determined whether luteolin prevents VSMC proliferation and migration through TGFBR1 inhibition. In this study, an in vitro model of VSMC proliferation and migration was utilized to investigate if TGFBR1 signaling is involved in the inhibitory effects of luteolin on VSMC proliferation and migration.

Materials and Methods

Materials

Antibodies targeting TGFBR1, proliferative cell nuclear antigen (PCNA), Cyclin D1, MMP2, and MMP9 were purchased from Abcam (Cambridge, MA, United States). The p-TGFBR1 (Ser 165) antibody was purchased from MyBioSource (San Diego, CA, United States). Antibodies specific for Smad2, Smad3, p-Smad2 (Ser 465/Ser 467), p-Smad3 (Ser 423/Ser 425), BAX, BCL-2, and GAPDH were purchased from Cell Signaling Technology (Beverly, MA, United States). Crystal Violet Staining Solution was purchased from KeyGEN BioTECH (Nanjing, China). Luteolin was purchased from Chengdu Mansite Biotechnology (Chengdu, China). Other reagents used in this study were obtained from reagent companies.

A7r5 and HASMC Cell Culture

A7r5 and HASMC cell lines were purchased from the Institute of Biochemistry and Cell Biology Cell Bank (Shanghai Institutes for Biological Sciences, Chinese Academy of Sciences, Shanghai, China). A7r5 and HASMC cells were incubated in DMEM (Gibco, Grand Island, NY, United States) supplemented with 10% fetal bovine serum (Gibco) at 37°C in an atmosphere of 5% CO2 and 95% air.

Adenovirus Preparation

To verify the effect of luteolin, we used recombinant adenovirus, including those used for overexpression of TGFBR1 and control vectors expressing a GFP marker (Ad-GFP), which were purchased from Vigene Biosciences (Shandong, China). A7r5 and HASMC cells were infected with adenovirus for 4 h, at a multiplicity of infection (MOI) of 100, with serum-free DMEM. Medium was replaced with DMEM containing 10% fetal bovine serum for 48 h before performing different treatments.

Cell Proliferation Assay

The Cell Titer 96® Aqueous One Solution Cell Proliferation Assay (MTS, Promega, Madison, WI, United States) was used to detect the effect of luteolin on the proliferation of A7r5 and HASMC cells according to the manufacturer’s protocol. Briefly, A7r5 and HASMC cells in logarithmic growth phase were harvested and seeded in 96-well plates at a density of 5 × 103 cells/well. Upon adherence to the plates, cells were starved with serum-free DMEM overnight. Next, they were treated with different concentrations of luteolin (10, 20, and 40 μM) and incubated for 48 h. MTS solution was directly added to each well (20 μl/well), and cells were cultured for 1 h; the optical density (OD) at 490 nm was measured for each well using an automatic microplate reader (Gene Company, Hong Kong, China).

5-Ethynyl-2′-Deoxyuridine Proliferation Assay

EdU was used to detect the effect of luteolin on the proliferation of A7r5 and HASMC cells according to the Cell-LightTM EdU Apollo®488 In Vitro Imaging Kit (Guangzhou RiboBio, China) instructions. Briefly, cells were treated with MTS as described in Section “Cell Proliferation Assay.” Then, EdU (100 mM, 100 μl/well) was added to the 96-well plates, which were incubated for 12 h. Next, A7r5 and HASMC cells were fixed with 4% paraformaldehyde. After 30 min, 1× Apollo® staining reaction liquid was added (100 μl/well), and cells were incubated for 30 min at room temperature; after 10 min of permeabilization with 0.5% Triton X-100, the cells were stained with 1× Hoechst 33342 (100 μl/well) for 30 min. The ratio of EdU-positive cells (EdU-stained cells/Hoechst-stained cells × 100%) was determined using a fluorescence microscope (Nikon Eclipse Ti-S, Tokyo, Japan).

Transwell Migration Assay

A Transwell chamber was used to detect the migrative ability of A7r5 and HASMC cells. Briefly, A7r5 and HASMC cells were seeded into the upper chamber at a density of 1 × 105 cells/well and treated with different concentrations of luteolin (10, 20, and 40 μM) with serum-free DMEM. The lower chambers, which were devoid of cells, were filled with DMEM with 10% fetal bovine serum and placed in 24-well plates. After incubation for 12 h, cells on the upper surface had migrated through the micropores (8-μM pores), and Transwell chambers were washed with phosphate-buffered saline three times; the lower side was fixed with 4% paraformaldehyde. After 20 min, the migrated cells were stained with crystal violet staining solution for 15 min. Five visual fields were randomly selected from each Transwell filter and captured at 400× magnification using an inverted fluorescence microscope (Nikon Eclipse Ti-S). The average number of cells that migrated through the Transwell filter was assayed using ImageJ software.

Wound Healing Migration Assay

Wound healing assays were used to detect the migrative ability of A7r5 and HASMC cells. Briefly, A7r5 and HASMC cells were treated with different concentrations of luteolin (10, 20, and 40 μM) for 24 h in DMEM containing 1% fetal bovine serum. Then, monolayers were wounded with a sterile 10-μl pipette tip, and suspended cells were washed away with phosphate-buffered saline twice. Images were captured of each well at 200× magnification at 0 and 48 h. Migration distance was estimated based on the widths of the wounds at 0 and 48 h, which were assayed using ImageJ software.

Western Blot Analysis

In brief, protein extracts from cultured cells were prepared in RIPA buffer containing protease and phosphatase inhibitors (all purchased from Beyotime, Shanghai, China). Lysates were centrifuged at 12,000 rpm at 4°C for 15 min to obtain the supernatant. Approximately 15 μg of protein was separated on a 10% SDS-PAGE gel and transferred to a polyvinylidene difluoride membrane (Millipore, Bedford, MA, United States), which was blocked with Tris-buffered saline containing 5% non-fat milk. Target proteins were incubated with the corresponding primary antibodies overnight at 4°C as follows: p-TGFBR1, 1:1,000; total TGFBR1, 1:1,000; p-Smad2, 1:1,000; p-Smad3, 1:1,000; total Smad2, 1:1,000; total Smad3, 1:1,000; Cyclin D1, 1:10,000; PCNA, 1:1,000; MMP2, 1:1,000; MMP9, 1:1,000; BAX, 1:1,000; BCL-2, 1:1,000; GAPDH, 1:1,000. Membranes were then washed three times with Tris-buffered saline containing 0.1% Tween 20 (TBST) and incubated with horseradish peroxidase-goat anti-rabbit IgG (Beverly, MA, United States) for 2 h at 4°C. GAPDH was used as an internal control. Band intensity was quantified using ImageJ software.

Molecular Docking and Simulation of Molecular Dynamics

Molecular docking, which was analyzed by Autodock Vina (Scripps Research Institute, United States) (Dos Santos et al., 2018), was used to analyze the binding mechanism between TGFBR1 (PDB ID: 1PY5) and luteolin (ZINC ID: 5280445). The molecular docking simulation protein of TGFBR1 was prepared by removing water molecules and bound ligands. YASARA was used to perform the energy minimization of ligands. In this study, the agonist-binding domain of TGFBR1 corresponded to the protein-ligand binding site of TGFBR1. The star conformation for the MD simulation was considered the best conformation.

YASARA was used to perform MD simulation (Land and Humble, 2018). AMBER 03 forcefield was used to run all simulations. Briefly, 0.9% NaCl, as the solvation of the receptor-ligand complex, was placed in a dodecahedron box; the distance between the box and the solute as 5 Å. The initiation of simulated annealing minimizations was set at 298 K, with velocities scaling down 0.9 every 10 steps lasting for 5 ps. After the energy was minimized, the temperature of the system was adjusted using a Berendsen thermostat to minimize the influence of temperature control. In addition, velocities were rescaled only every 100 simulation steps, whenever the mean of the last 100 detected temperatures converged. Finally, 100-ns MD simulations were conducted at a rate of 2 fs, and the coordinates of the complexes were saved every 10 ps.

Data Analysis

All data are expressed as the mean ± SD and were analyzed using the SPSS statistical package (version 13.0, SPSS Inc., Chicago, IL, United States). One-way analysis of variance factorial analysis was used to compare measurement data among the groups. If the variances were homogenous, the mean values were compared by performing an F-test. The differences between two groups were analyzed using a Bonferroni test. If the variances were not homogenous, Welch’s test and Dunnett’s T3 analysis were adopted. P < 0.05 indicated statistically significant differences. Each experiment was repeated three or more times.

Results

Luteolin Inhibits VSMC Proliferation Without Inducing Apoptosis

To detect the effect of luteolin on VSMC proliferation, A7r5 and HASMC cells were treated as described in the experimental section. Then, MTS assays and EdU staining were used to detect the proliferation of VSMCs. As shown in Figures 1A,B, treatment with different concentrations of luteolin (10, 20, and 40 μM) significantly inhibited A7r5 and HASMC cell proliferation in a dose-dependent manner (P < 0.05). Consistently, the ratios of EdU-positive cell numbers for A7r5 and HASMC cells were also dose-dependently reduced by luteolin (P < 0.05, Figures 1C–E). In addition, we further detected the inhibitory effect of luteolin through the expression of proliferation-related proteins including Cyclin D1 and PCNA. CyclinD1 is a regulatory protein of the cell cycle and plays a key role in progression from the G1 to S phase in the mammalian cell. PCNA is a well-defined regulator of DNA replication and cell cycle control and is used as a proliferation index of a broader spectrum of cells (Zhao et al., 2013; Kim et al., 2014; Xu et al., 2015). As shown in Figures 1F–H, this compound significantly inhibited the expression of Cyclin D1 and PCNA in a dose-dependent manner (P < 0.05). Furthermore, to determine if luteolin inhibits cell proliferation by inducing apoptosis, we detected the expression of apoptosis-related proteins by western blotting. As shown in Figures 1I–K, we found that luteolin could not induce VSMC apoptosis. These results indicated that luteolin inhibited VSMC proliferation without inducing apoptosis.

FIGURE 1
www.frontiersin.org

FIGURE 1. Luteolin suppresses vascular smooth muscle cell (VSMC) proliferation without inducing apoptosis. (A,B) A7r5 and HASMC cells were incubated with different concentrations of luteolin (10, 20, and 40 μM) for 48 h, and then VSMC proliferation was tested by MTS assays (n = 6). (C) EdU proliferation assay. The red regions (EdU) in the images represent proliferating cells, and the blue regions (Hoechst 33324) represent the nuclei of all cells. (D,E) Percentage of EdU-positive VSMCs (the proliferation ratio was tested by assessing Hoechst 33324/EdU staining; n = 3). (F) A7r5 and HASMC cells were incubated with luteolin (10, 20, and 40 μM) for 24 h, and the expression of Cyclin D1 and PCNA was tested by western blotting. (GH) Relative expression levels of Cyclin D1 and PCNA (n = 3). (I) Expression levels of BAX and BCL-2 were tested by western blotting. (J) Relative expression levels of BAX (n = 3). (K) Relative expression levels of BCL-2 (n = 3). Data are presented as the mean ± SD. P < 0.05, P < 0.01, compared to the control group.

Luteolin Inhibits VSMC Migration

To detect the effect of luteolin on VSMC migration, wound healing and Transwell assays were performed using A7r5 and HASMC cells. As shown in Figures 2A–C, luteolin inhibited A7r5 and HASMC cell migration in a dose-dependent manner (P < 0.05). Consistently, as shown in Figures 2D–F, the migrative ability of A7r5 and HASMC cells was also dose-dependently reduced by luteolin (P < 0.05). In addition, we further detected the expression levels of migration-related proteins (MMP2 and MMP9) to verify these effects. MMPs are a large family of calcium-dependent zinc-containing endopeptidases, which are responsible for tissue remodeling and degradation of the extracellular matrix (ECM). Previous studies have demonstrated that MMP2 and MMP9 are most closely involved in VSMC migration in both in vitro and in vivo models (Kenagy et al., 1997; Ding et al., 2011; Guo et al., 2013). Selective gene silencing of MMP2 and MMP9, or interruption of MMP2 and MMP9 gene expression by chemical inhibitors, inhibit the migration of VSMCs (Pinney et al., 2003; Hlawaty et al., 2007; Turner et al., 2007). As shown in Figures 2G–K, this compound inhibited the expression of MMP2 and MMP9 in a dose-dependent manner (P < 0.05). These findings indicated that luteolin inhibits VSMC migration.

FIGURE 2
www.frontiersin.org

FIGURE 2. Luteolin suppresses the migration of vascular smooth muscle cells (VSMCs). (A) A7r5 and HASMC cells were incubated with different concentrations of luteolin (10, 20, and 40 μM) for 24 h, and the migration of VSMCs was tested by wound healing assays; the migration distance was estimated based on the width of the scrape at 0 and 48 h. (B,C) Percentage of relative VSMC migration (n = 3). (D) A7r5 and HASMC cells were harvested in the logarithmic growth phase, treated with different concentrations of luteolin (10, 20, and 40 μM), and tested by performing Transwell assays for 12 h. (E,F) The number of cells in each field of view (n = 5). (G) A7r5 and HASMC cells, respectively, were incubated with luteolin (10, 20, and 40 μM) for 24 h, and the expression levels of MMP2 and MMP9 were tested by western blotting. (HK) Relative expression levels of MMP2 and MMP9 (n = 3). Data are presented as the mean ± SD. P < 0.05, ∗∗P < 0.01, compared to the control group.

Luteolin Inhibits Activation of the TGFBR1 Signaling Pathway

Here, we investigated whether luteolin inhibits TGFBR1 signaling. As shown in Figures 3A–C, administration of luteolin significantly inhibited the phosphorylation of TGFBR1 in a dose-dependent manner without altering total levels of TGFBR1 (P < 0.05). Smad2 and Smad3 are the downstream signaling molecules of TGFBR1, and as shown in Figures 3A,D–G, the phosphorylation levels of Smad2 and Smad3 were also reduced with treatment in a dose-dependent manner without altering total levels of Smad2 and Smad3 (P < 0.05). In addition, we further evaluated the effects of luteolin on TGF-β-induced TGFBR1 and Smad activation in VSMCs. As shown in Supplementary Figures S1A,B, the phosphorylation levels of TGFBR1 were significantly enhanced under TGF-β treatment (P < 0.05), suggesting that the binding of TGF-β activates TGFBR1. Consistently, as shown in Supplementary Figures S1A,C,D, the phosphorylation levels of Smad2 and Smad3 were also significantly enhanced in TGF-β-stimulated VSMCs (P < 0.05). Furthermore, luteolin (40 μM) significantly reduced the phosphorylation levels of TGFBR1, Smad2, and Smad3 (P < 0.05). These results strongly indicate that luteolin inhibited the activation of TGFBR1/Smad signaling.

FIGURE 3
www.frontiersin.org

FIGURE 3. Luteolin suppresses activation of the TGFBR1 signaling pathway. (A) A7r5 and HASMC cells, respectively, were incubated with luteolin (10, 20, and 40 μM) for 1 h, and the expression levels of p-TGFBR1, TGFBR1, p-Smad2, Smad2, p-Smad3, and Smad3 were tested by western blotting. (BG) Relative phosphorylation levels of TGFBR1, Smad2, and Smad3 (n = 3). Data are presented as the mean ± SD. P < 0.05, ∗∗P < 0.01, compared to the control group.

Overexpression of TGFBR1 Partially Abolishes the Inhibitory Effect of Luteolin on Smads

To determine if the inhibitory effect of luteolin on Smad signaling was due to suppression of TGFBR1, we overexpressed human TGFBR1 in A7r5 and HASMC cells using a recombinant adenovirus (Ad-TGFBR1). After Ad-TGFBR1 and negative control adenovirus (Ad-NC) treatment was administered for 48 h, A7r5 and HASMC cells were treated with luteolin (40 μM). As shown in Figures 4A–C, the phosphorylation levels of TGFBR1 were significantly enhanced in TGFBR1-overexpressing VSMCs (P < 0.05). The phosphorylation levels of Smad2 and Smad3, which are downstream molecules of TGFBR1, were also measured by western blotting. As presented in Figures 4A,D–G, overexpression of TGFBR1 reduced the inhibitory effects of luteolin on the phosphorylation of Smad2 and Smad3 in A7r5 and HASMC cells (P < 0.05). Previous studies have demonstrated that luteolin exerts an inhibitory action on VSMC proliferation through suppressing Akt signaling (Lang et al., 2012; Xu et al., 2015). Thus, in this study, we investigated whether luteolin inhibits Akt in a TGFBR1-dependent manner. Interestingly, we found that the phosphorylation level of Akt was reduced in TGFBR1-overexpressing VSMCs (Supplementary Figures S2A,B; P < 0.05). Treatment with luteolin (40 μM) inhibited the phosphorylation level of Akt; overexpression of TGFBR1 enhanced the inhibitory effects of luteolin on Akt (Supplementary Figures S2A,B; P < 0.05). These results indicated that luteolin inhibited Smads trough a TGFBR1-dependent manner.

FIGURE 4
www.frontiersin.org

FIGURE 4. Overexpression of TGFBR1 partially rescues TGFBR1 activation after luteolin treatment. (A) A7r5 and HASMC cells, respectively, were incubated with a TGFBR1 overexpression adenovirus vector and GFP control adenovirus for 48 h. Cells were then treated with luteolin (40 μM) for 1 h. The expression levels of p-TGFBR1, TGFBR1, p-Smad2, Smad2, p-Smad3, and Smad3 were tested by western blotting. (BG) Relative phosphorylation levels of TGFBR1, Smad2, and Smad3 (n = 3). Data are presented as the mean ± SD. P < 0.05, ∗∗P < 0.01, compared to the control group.

Overexpression of TGFBR1 Partially Blocks the Inhibitory Effect of Luteolin on VSMC Proliferation

To further verify that the inhibitory effects of luteolin on VSMC proliferation are associated with inhibiting the activation of TGFBR1, we overexpressed TGFBR1 in A7r5 and HASMC cells and then detected the anti-proliferative effects of luteolin. As shown in Figures 5A,B, overexpression of TGFBR1 significantly enhanced the proliferation of luteolin-treated A7r5 and HASMC cells (P < 0.05). As shown in Figures 5C–E, consistently, overexpression of TGFBR1 also significantly reduced the inhibitory effect of luteolin on the EdU-positive cell ratio in both cell lines (P < 0.05). Next, we further detected the expression levels of cell cycle-related proteins including Cyclin D1 and PCNA. As shown in Figures 5F–J, the expression of Cyclin D1 and PCNA was partially rescued by overexpression of TGFBR1 (P < 0.05). These findings indicated that the inhibitory effects of luteolin on VSMC proliferation are at least partly due to inhibition of TGFBR1 activation.

FIGURE 5
www.frontiersin.org

FIGURE 5. Overexpression of TGFBR1 partially blocks the inhibitory effect of luteolin on vascular smooth muscle cell (VSMC) proliferation. (A,B) A7r5 and HASMC cells, respectively, were incubated with TGFBR1 overexpression adenovirus vector and GFP control adenovirus for 48 h and then treated with luteolin (40 μM) for 48 h. Then, the proliferation of VSMCs was tested by MTS assays (n = 6). (C) EdU proliferation assay. The red regions (EdU) in the images represent proliferating cells, and the blue regions (Hoechst 33324) represent the nuclei of all cells. (D,E) Percentage of EdU-positive VSMCs (n = 3). (F) A7r5 and HASMC cells, respectively, were incubated with TGFBR1 overexpression adenovirus vector and GFP control adenovirus for 48 h and then treated with luteolin (40 μM) for 24 h. the expression levels of Cyclin D1 and PCNA were tested by western blotting. (GJ) Relative expression levels of Cyclin D1 and PCNA (n = 3). Data are shown as the mean ± SD. P < 0.05, ∗∗P < 0.01.

Overexpression of TGFBR1 Partially Blocks the Inhibitory Effect of Luteolin on VSMC Migration

Next, we further determined whether the inhibitory effect of luteolin on VSMC migration was due to inhibited activation of TGFBR1. TGFBR1 was overexpressed in A7r5 and HASMC cells using Ad-TGFBR1, and the anti-migration effects of luteolin were measured by wound healing and Transwell assays. As shown in Figures 6A–C, overexpression of TGFBR1 significantly enhanced the migration of luteolin-treated A7r5 and HASMC cells (P < 0.05). As shown in Figures 6D–F, consistently, overexpression of TGFBR1 also partially blocked the inhibitory effect of luteolin on the number of migrating cells, and this difference was statistically significant (P < 0.05). Our findings also suggested that the inhibitory effects of luteolin on VSMC migration were due to, at least in part, inhibition of TGFBR1 activation.

FIGURE 6
www.frontiersin.org

FIGURE 6. Overexpression of TGFBR1 partially blocks the inhibitory effect of luteolin on vascular smooth muscle cell (VSMC) migration. (A) A7r5 and HASMC cells were incubated with TGFBR1 overexpression adenovirus vector and GFP control adenovirus for 48 h and then treated with luteolin (40 μM) for 24 h. Then, the migration of VSMCs was tested by wound healing assays. (B,C) Percentage of relative VSMC migration (n = 3). (D) A7r5 and HASMC cells were incubated with TGFBR1 overexpression adenovirus vector and GFP control adenovirus for 48 h. Then, Transwell assays were performed after treatment with luteolin (40 μM) for 12 h. (E,F) The number of cells in each field of view (n = 5). Data are shown as the mean ± SD. P < 0.05, ∗∗P < 0.01.

Luteolin Directly Targets TGFBR1

To explore the interaction between luteolin and TGFBR1, we performed molecular docking analysis, using the docking program Autodock Vina. The binding energy of the TGFBR1-luteolin complex was found to be -10.194 kcal/mol, which indicated good binding ability. The three-dimensional binding conformation of the TGFBR1-luteolin complex is presented in Figure 7A. We found six hydrogen bonds that formed between luteolin and GLU-245, LYS-232, ASP-351, SER-280, SER-287, and ALA-230 of TGFBR1. Next, to further verify the results of molecular docking, the best conformation of TGFBR1-luteolin was taken as the start conformation for MD simulation by YASARA. The surface visualization models of the TGFBR1-luteolin complex are shown in Figure 7B. Luteolin steadily presented at the center of the TGFBR1 binding site until the end of the MD simulation. In addition, as shown in Figure 7C, both the heavy atoms RMSD track of the TGFBR1-luteolin complex (Figure 7C, red line) and unbound TGFBR1 (TGFBR1-free) (Figure 7C, blue line) mildly fluctuated around 2 Å during 0–100 ns (Figure 7C). The track of potential energy of the TGFBR1-luteolin complex was higher than that of unbound TGFBR1 (Figure 7D). These results suggested that binding between TGFBR1 and luteolin is stable, and that luteolin might directly target TGFBR1.

FIGURE 7
www.frontiersin.org

FIGURE 7. Molecular docking and simulation of the molecular dynamics. (A) Three-dimensional crystal structure of luteolin (ZINC18185774) in complex with TGFBR1 (PDB ID: 1PY5). Luteolin is shown in green, and the hydrogen bonds are indicated by sky-blue lines. (B) Surface presentation of the TGFBR1-luteolin complex crystal structure at 0 and 100 ns. (C) Plots of root mean square deviation (RMSD) of heavy atoms of TGFBR1 unbound (blue) and TGFBR1-luteolin complex (red). (D) Potential energy profiles of TGFBR1-free (blue) and TGFBR1-luteolin complex (red) during the 100-ns molecular dynamics simulation.

Discussion

Under physiological conditions, VSMCs regulate vascular tone through contraction and relaxation and maintain normal vascular function by secreting and releasing vascular regulatory factors. However, VSMC migration from the media to the subendothelial space and abnormal proliferation represent common pathological bases of AS, hypertension, and vascular restenosis. Luteolin is a common flavonoid that widely exists in many Chinese herbal medicines, vegetables, and fruits. Several studies have demonstrated that luteolin has multiple cardio-protective effects via complex signal transduction pathways and target effectors. In addition, epidemiological evidence has shown that intake of these compounds helps prevent CVD (Luo et al., 2017). These results suggest that generous dietary intake of flavonoids may be an effective primary prevention strategy against CVD. Relevant studies have indicated that luteolin prevents VSMC remodeling, using in vivo models. For example, Zhang et al. (2009) reported that flavonoids inhibit neointima formation. Su et al. (2015) reported that luteolin effects decline in media thickness of vascular wall and attenuates hypertensive vascular remodeling through inhibiting the proliferation and migration of VSMCs. In this study, we have demonstrated the inhibitory effect of luteolin on VSMC proliferation and migration. Our study will enrich the current knowledge on the role of luteolin on VSMC proliferation and migration, similar to previous studies (Lang et al., 2012; Jiang et al., 2013; Xu et al., 2015).

TGF-β is a multifunctional cytokine that regulates various cell functions (Agrotis et al., 2005; Weiss and Attisano, 2013) including smooth muscle cell proliferation and migration (Goumans et al., 2009). It has been demonstrated that TGF-β binds TGFBR2 and then induces the phosphorylation of TGFBR1 at Ser165 (Souchelnytskyi et al., 1996). In the present study, we found that luteolin significantly inhibits phosphorylation of TGFBR1 at Ser165 in a dose-dependent manner. Activated TGFBR1 specifically recognizes and phosphorylates Smad2/3, following which phosphorylated Smad2/3 translocates to the nucleus with Smad4 to activate various transcriptional programs through binding to a Smad-binding element (Weiss and Attisano, 2013). Therefore, we further investigated the inhibitory effects of luteolin on the phosphorylation of Smad2/3. We found that treatment with luteolin reduced the phosphorylation levels of Smad2 and Smad3 in a dose-dependent manner. Our findings are in agreement with a previous study showing that luteolin inhibits Smad2 activation in hepatic stellate cells (Li et al., 2015). Then, we further determined whether the inhibitory effects of luteolin on the activation of Smad2/3 were due to TGFBR1 inhibition by overexpressing this receptor. We found that the inhibitory effects of this compound were partly counteracted by TGFBR1 overexpression. These data suggested that luteolin inhibits activation of Smad2/3 via inhibition of TGFBR1.

Next, we further determined whether the inhibitory action of luteolin on VSMC proliferation and migration is associated with its ability to inhibit the activation of TGFBR1 using a TGFBR1-overexpression adenovirus vector. We found that such inhibitory effects were partially blocked with exogenous TGFBR1 expression. These results demonstrated that TGFBR1 is the pharmacological target of luteolin with respect to VSMC proliferation and migration. In addition, the essential involvement of TGFBR1 has been detected in a series of pathophysiological processes associated with CVDs, such as cardiac fibrosis, vascular restenosis, and inflammation (Kajdaniuk et al., 2013; Goumans and Ten Dijke, 2018). Our findings suggest that luteolin could be used to regulate these processes.

The present findings suggested that luteolin might target TGFBR1. Thus, we tested its ability to target TGFBR1 directly by performing molecular docking and MD simulation to assess the mechanism through which luteolin and TGFBR1 interact. The binding energy of the luteolin-TGFBR1 complex indicated strong binding. Moreover, luteolin interacted with GLU-245, LYS-232, ASP-351, SER-280, SER-287, and ALA-230 of TGFBR1 via hydrogen bonds. These residues are the key residues involved in inhibition of the binding domain of TGFBR1 (Sawyer et al., 2004). Furthermore, MD simulation results demonstrated that the luteolin-TGFBR1 binding conformation was stable. These results strongly indicated that luteolin directly targets TGFBR1, acting as a potential TGFBR1 inhibitor.

Multi-targets and complex signaling is an important characteristic of natural products. Previous studies have demonstrated that luteolin inhibits the proliferation and migration of VSMCs by regulating Akt (Lang et al., 2012; Xu et al., 2015). To determine whether luteolin inhibits Akt through suppressing TGFBR1 signaling, we investigated the effect of luteolin on the phosphorylation level of Akt both in TGFBR1-overexpressing and negative control VSMCs. Interestingly, we found the phosphorylation level of Akt to be reduced in TGFBR1-overexpressing VSMCs. Previous studies have shown that treatment with TGF-β increased Akt activation in type 2 diabetic nephropathy mice model and cultured human podocytes (Madne and Dockrell, 2018; Zhang et al., 2018). Different roles of TGFBR1 in Akt signaling might be cell-type-dependent. We also found that luteolin inhibited the phosphorylation level of Akt; overexpression of TGFBR1 significantly enhanced the inhibitory effect of luteolin on Akt, while the activation of Smad2 and Smad3 was increased and the proliferation of VSMC was further enhanced. These results strongly indicate that luteolin performs its inhibitory function on Akt in a TGFBR1-independent manner. Inhibition of TGFBR1 signaling plays a more important role in the inhibitory actions of luteolin than suppression of Akt.

In this study, we demonstrated that luteolin inhibits VSMC proliferation and migration by suppressing TGFBR1 signaling without inducing cell injury. Our findings provide novel insights into the anti-atherosclerotic mechanisms of luteolin. However, there are some limitations to our study. Firstly, we investigated this new potential mechanism of luteolin only in vitro, and did not further validate it in animal studies. Secondly, we did not investigate the effective therapeutic dosage and the safety of long-term administration of luteolin. Our future studies will further investigate the molecular mechanisms of luteolin both in vivo and in vitro.

In summary, a series of experiments was designed to clarify whether TGFBR1 signaling underlies the inhibitory effects of luteolin on VSMC proliferation and migration. The main findings in the present study are as follows. (i) Luteolin inhibited the proliferation and migration of A7r5 and HASMC cells, without affecting apoptosis. (ii) Luteolin decreased TGFBR1 and Smad2/3 phosphorylation. (iii) Overexpression of TGFBR1 relieved the inhibitory effects of luteolin on VSMC proliferation and migration, and phosphorylation of TGFBR1 pathway-associated proteins. (iv) Molecular simulation demonstrated strong and stable binding between luteolin and TGFBR1.

Author Contributions

Y-TW supervised the entire work and performed the cell cultures and proliferation and migration assays. LC and Z-BT performed the cell cultures and adenovirus transfection. L-PX and W-TZ performed the molecular docking and simulation of molecular dynamics. H-MC and JL analyzed the data. Y-CZ and BL conceived and designed the experiments and critically revised the manuscript. All authors discussed the results and contributed to manuscript writing.

Funding

This work was supported by the National Natural Science Foundation of China (81673805, 81704058, and 81774100), Science and Technology Planning Project of Guangdong Province (2014A020221013, 2014A020221059, 2015A020211040, and 2016A020226026), Traditional Chinese Medicine Bureau of Guangdong Province (20161167), Climbing Program Special Funds of Guangdong Province (pdjh2017a0094), and National Training Program of Innovation and Entrepreneurship for Undergraduates (201712121006).

Conflict of Interest Statement

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.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2018.01059/full#supplementary-material

Abbreviations

AS, atherosclerosis; CVD, cardiovascular diseases; DMEM, Dulbecco’s modified Eagle’s medium; Edu, 5-ethynyl-2′-deoxyuridine; GFP, green fluorescent protein; TGFBR1, transforming growth factor beta receptor 1; VSMC, vascular smooth muscle cell.

References

Agrotis, A., Kalinina, N., and Bobik, A. (2005). Transforming growth factor-beta, cell signaling and cardiovascular disorders. Curr. Vasc. Pharmacol. 3, 55–61. doi: 10.2174/1570161052773951

CrossRef Full Text | Google Scholar

Ding, Q., Chai, H., Mahmood, N., Tsao, J., Mochly-Rosen, D., and Zhou, W. (2011). Matrix metalloproteinases modulated by protein kinase Cepsilon mediate resistin-induced migration of human coronary artery smooth muscle cells. J. Vasc. Surg. 53, 1044–1051. doi: 10.1016/j.jvs.2010.10.117

PubMed Abstract | CrossRef Full Text | Google Scholar

Dos Santos, R. N., Ferreira, L. G., and Andricopulo, A. D. (2018). Practices in molecular docking and structure-based virtual screening. Methods Mol. Biol. 1762, 31–50. doi: 10.1007/978-1-4939-7756-7_3

PubMed Abstract | CrossRef Full Text | Google Scholar

Goumans, M. J., Liu, Z., and Ten Dijke, P. (2009). TGF-beta signaling in vascular biology and dysfunction. Cell Res. 19, 116–127. doi: 10.1038/cr.2008.326

PubMed Abstract | CrossRef Full Text | Google Scholar

Goumans, M. J., and Ten Dijke, P. (2018). TGF-beta signaling in control of cardiovascular function. Cold Spring Harb. Perspect. Biol. 10:a022210 doi: 10.1101/cshperspect.a022210

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo, J., Dhaliwall, J. K., Chan, K. K., Ghanim, H., Al Koudsi, N., Lam, L., et al. (2013). In vivo effect of insulin to decrease matrix metalloproteinase-2 and -9 activity after arterial injury. J. Vasc. Res. 50, 279–288. doi: 10.1159/000351611

PubMed Abstract | CrossRef Full Text | Google Scholar

Hlawaty, H., San Juan, A., Jacob, M. P., Vranckx, R., Letourneur, D., and Feldman, L. J. (2007). Inhibition of MMP-2 gene expression with small interfering RNA in rabbit vascular smooth muscle cells. Am. J. Physiol. Heart Circ. Physiol. 293, H3593–H3601. doi: 10.1152/ajpheart.00517.2007

PubMed Abstract | CrossRef Full Text | Google Scholar

Jia, Z., Nallasamy, P., Liu, D., Shah, H., Li, J. Z., Chitrakar, R., et al. (2015). Luteolin protects against vascular inflammation in mice and TNF-alpha-induced monocyte adhesion to endothelial cells via suppressing IKappaBalpha/NF-kappaB signaling pathway. J. Nutr. Biochem. 26, 293–302. doi: 10.1016/j.jnutbio.2014.11.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, D., Li, D., and Wu, W. (2013). Inhibitory effects and mechanisms of luteolin on proliferation and migration of vascular smooth muscle cells. Nutrients 5, 1648–1659. doi: 10.3390/nu5051648

PubMed Abstract | CrossRef Full Text | Google Scholar

Kajdaniuk, D., Marek, B., Borgiel-Marek, H., and Kos-Kudla, B. (2013). Transforming growth factor beta1 (TGFbeta1) in physiology and pathology. Endokrynol. Pol. 64, 384–396. doi: 10.5603/EP.2013.0022

PubMed Abstract | CrossRef Full Text | Google Scholar

Kenagy, R. D., Hart, C. E., Stetler-Stevenson, W. G., and Clowes, A. W. (1997). Primate smooth muscle cell migration from aortic explants is mediated by endogenous platelet-derived growth factor and basic fibroblast growth factor acting through matrix metalloproteinases 2 and 9. Circulation 96, 3555–3560. doi: 10.1161/01.CIR.96.10.3555

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, M. H., Ham, O., Lee, S. Y., Choi, E., Lee, C. Y., Park, J. H., et al. (2014). MicroRNA-365 inhibits the proliferation of vascular smooth muscle cells by targeting cyclin D1. J. Cell. Biochem. 115, 1752–1761. doi: 10.1002/jcb.24841

PubMed Abstract | CrossRef Full Text | Google Scholar

Lacolley, P., Regnault, V., Nicoletti, A., Li, Z., and Michel, J. B. (2012). The vascular smooth muscle cell in arterial pathology: a cell that can take on multiple roles. Cardiovasc. Res. 95, 194–204. doi: 10.1093/cvr/cvs135

PubMed Abstract | CrossRef Full Text | Google Scholar

Land, H., and Humble, M. S. (2018). YASARA: a tool to obtain structural guidance in biocatalytic investigations. Methods Mol. Biol. 1685, 43–67. doi: 10.1007/978-1-4939-7366-8_4

PubMed Abstract | CrossRef Full Text | Google Scholar

Lang, Y., Chen, D., Li, D., Zhu, M., Xu, T., Zhang, T., et al. (2012). Luteolin inhibited hydrogen peroxide-induced vascular smooth muscle cells proliferation and migration by suppressing the Src and Akt signalling pathways. J. Pharm. Pharmacol. 64, 597–603. doi: 10.1111/j.2042-7158.2011.01438.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, J., Li, X., Xu, W., Wang, S., Hu, Z., Zhang, Q., et al. (2015). Antifibrotic effects of luteolin on hepatic stellate cells and liver fibrosis by targeting AKT/mTOR/p70S6K and TGFbeta/Smad signalling pathways. Liver Int. 35, 1222–1233. doi: 10.1111/liv.12638

PubMed Abstract | CrossRef Full Text | Google Scholar

Lopez-Lazaro, M. (2009). Distribution and biological activities of the flavonoid luteolin. Mini Rev. Med. Chem. 9, 31–59. doi: 10.2174/138955709787001712

CrossRef Full Text | Google Scholar

Luo, Y., Shang, P., and Li, D. (2017). Luteolin: a flavonoid that has multiple cardio-protective effects and its molecular mechanisms. Front. Pharmacol. 8:692. doi: 10.3389/fphar.2017.00692

PubMed Abstract | CrossRef Full Text | Google Scholar

Madne, T. H., and Dockrell, M. E. C. (2018). TGFbeta1-mediated PI3K/Akt and p38 MAP kinase dependent alternative splicing of fibronectin extra domain A in human podocyte culture. Cell. Mol. Biol. 64, 127–135. doi: 10.14715/cmb/2018.64.5.22

PubMed Abstract | CrossRef Full Text | Google Scholar

Nile, S. H., Keum, Y. S., Nile, A. S., Jalde, S. S., and Patel, R. V. (2018). Antioxidant, anti-inflammatory, and enzyme inhibitory activity of natural plant flavonoids and their synthesized derivatives. J. Biochem. Mol. Toxicol. 32:e22002. doi: 10.1002/jbt.22002

PubMed Abstract | CrossRef Full Text | Google Scholar

Pinney, S. P., Chen, H. J., Liang, D., Wang, X., Schwartz, A., and Rabbani, L. E. (2003). Minocycline inhibits smooth muscle cell proliferation, migration and neointima formation after arterial injury. J. Cardiovasc. Pharmacol. 42, 469–476. doi: 10.1097/00005344-200310000-00003

PubMed Abstract | CrossRef Full Text | Google Scholar

Sawyer, J. S., Beight, D. W., Britt, K. S., Anderson, B. D., Campbell, R. M., Goodson, T., et al. (2004). Synthesis and activity of new aryl- and heteroaryl-substituted 5,6-dihydro-4H-pyrrolo[1,2-b]pyrazole inhibitors of the transforming growth factor-beta type I receptor kinase domain. Bioorg. Med. Chem. Lett. 14, 3581–3584. doi: 10.1016/j.bmcl.2004.04.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Serralheiro, P., Soares, A., Costa Almeida, C. M., and Verde, I. (2017). TGF-beta1 in Vascular Wall Pathology: unraveling Chronic Venous Insufficiency Pathophysiology. Int. J. Mol. Sci. 18:E2534 doi: 10.3390/ijms18122534

PubMed Abstract | CrossRef Full Text | Google Scholar

Souchelnytskyi, S., Ten Dijke, P., Miyazono, K., and Heldin, C. H. (1996). Phosphorylation of Ser165 in TGF-beta type I receptor modulates TGF-beta1-induced cellular responses. EMBO J. 15, 6231–6240. doi: 10.1002/j.1460-2075.1996.tb01013.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Su, J., Xu, H. T., Yu, J. J., Gao, J. L., Lei, J., Yin, Q. S., et al. (2015). Luteolin Ameliorates Hypertensive Vascular Remodeling through Inhibiting the Proliferation and Migration of Vascular Smooth Muscle Cells. Evid. Based Complement. Alternat. Med. 2015:364876. doi: 10.1155/2015/364876

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, Y., Ye, P., Wu, J., Liu, Z., Zhang, A., Ren, L., et al. (2014). Inhibition of intimal hyperplasia in murine aortic allografts by the oral administration of the transforming growth factor-beta receptor I kinase inhibitor SD-208. J. Heart Lung Transplant. 33, 654–661. doi: 10.1016/j.healun.2014.02.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Turner, N. A., Hall, K. T., Ball, S. G., and Porter, K. E. (2007). Selective gene silencing of either MMP-2 or MMP-9 inhibits invasion of human saphenous vein smooth muscle cells. Atherosclerosis 193, 36–43. doi: 10.1016/j.atherosclerosis.2006.08.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, L., Lee, I. M., Zhang, S. M., Blumberg, J. B., Buring, J. E., and Sesso, H. D. (2009). Dietary intake of selected flavonols, flavones, and flavonoid-rich foods and risk of cancer in middle-aged and older women. Am. J. Clin. Nutr. 89, 905–912. doi: 10.3945/ajcn.2008.26913

PubMed Abstract | CrossRef Full Text | Google Scholar

Weber, C., and von Hundelshausen, P. (2017). CANTOS trial validates the inflammatory pathogenesis of atherosclerosis: setting the stage for a new chapter in therapeutic targeting. Circ. Res. 121, 1119–1121. doi: 10.1161/CIRCRESAHA.117.311984

PubMed Abstract | CrossRef Full Text | Google Scholar

Weiss, A., and Attisano, L. (2013). The TGFbeta superfamily signaling pathway. Wiley Interdiscip. Rev. Dev. Biol. 2, 47–63. doi: 10.1002/wdev.86

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, T., Zhu, H., Li, D., Lang, Y., Cao, L., Liu, Y., et al. (2015). Luteolin inhibits angiotensin II-stimulated VSMC proliferation and migration through downregulation of Akt phosphorylation. Evid. Based Complement. Alternat. Med. 2015:931782. doi: 10.1155/2015/931782

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, H., Han, M., and Wen, J. (2009). [Flavonoids extracts from Inula britannica inhibit oxidative stress in vessel after balloon injury]. Zhongguo Zhong Yao Za Zhi 34, 615–619.

Google Scholar

Zhang, Y., Wang, B., Guo, F., Li, Z., and Qin, G. (2018). Involvement of the TGFbeta1- ILK-Akt signaling pathway in the effects of hesperidin in type 2 diabetic nephropathy. Biomed. Pharmacother. 105, 766–772. doi: 10.1016/j.biopha.2018.06.036

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, Y., Lv, M., Lin, H., Cui, Y., Wei, X., Qin, Y., et al. (2013). Rho-associated protein kinase isoforms stimulate proliferation of vascular smooth muscle cells through ERK and induction of cyclin D1 and PCNA. Biochem. Biophys. Res. Commun. 432, 488–493. doi: 10.1016/j.bbrc.2013.02.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: luteolin, vascular smooth muscle cell, proliferation, migration, TGFBR1, Smad2/3

Citation: Wu Y-T, Chen L, Tan Z-B, Fan H-J, Xie L-P, Zhang W-T, Chen H-M, Li J, Liu B and Zhou Y-C (2018) Luteolin Inhibits Vascular Smooth Muscle Cell Proliferation and Migration by Inhibiting TGFBR1 Signaling. Front. Pharmacol. 9:1059. doi: 10.3389/fphar.2018.01059

Received: 13 May 2018; Accepted: 03 September 2018;
Published: 21 September 2018.

Edited by:

Changhua Wang, Wuhan University, China

Reviewed by:

Javier Angulo, Hospital Universitario Ramón y Cajal, Spain
Nazareno Paolocci, Johns Hopkins University, United States

Copyright © 2018 Wu, Chen, Tan, Fan, Xie, Zhang, Chen, Li, Liu and Zhou. 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: Bin Liu, xmhoolv@163.com Ying-Chun Zhou, zhychun@126.com, orcid.org/0000-0003-2623-0472

Disclaimer: 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.