ORIGINAL RESEARCH article

Front. Chem., 17 August 2021
Sec. Medicinal and Pharmaceutical Chemistry
https://doi.org/10.3389/fchem.2021.736070

Aspulvinones Suppress Postprandial Hyperglycemia as Potent α-Glucosidase Inhibitors From Aspergillus terreus ASM-1

www.frontiersin.orgChangjing Wu1,2,3, www.frontiersin.orgXiang Cui1, www.frontiersin.orgLuzhen Sun1, www.frontiersin.orgJiajia Lu1, www.frontiersin.orgFeng Li1, www.frontiersin.orgMinghui Song1, www.frontiersin.orgYunxia Zhang1, www.frontiersin.orgXinqi Hao2, www.frontiersin.orgCongkui Tian3*, www.frontiersin.orgMaoping Song2* and www.frontiersin.orgXiaomeng Liu1,4*
  • 1College of Life Sciences and Agronomy, Zhoukou Normal University, Zhoukou, China
  • 2College of Chemistry and Molecular Engineering, Zhengzhou University, Zhengzhou, China
  • 3Wuling Mountain Institute of Natural Medicine, Hubei Minzu University, Enshi, China
  • 4College of Public Health, Xinxiang Medical University, Xinxiang, China

Chemical investigation of Aspergillus terreus ASM-1 fermentation resulted in the isolation of three new prenylated aspulvinones V–X (13), together with the previously reported analogs, aspulvinone H (4), J-CR (5), and R (6). Their structures were elucidated by various spectroscopic methods including HRESIMS and NMR, and the absolute configurations of 2 and 3 were determined by ECD comparison. Compounds 16 were evaluated for α-glucosidase inhibitory effects with acarbose as positive control. As a result, compounds 1 and 4 exhibited potent α-glucosidase inhibitory activities with IC50 values of 2.2 and 4.6 µM in mixed-type manners. The thermodynamic constants recognized the interaction between inhibitors and α-glucosidase was hydrophobic force-driven spontaneous exothermic reaction. The CD spectra also indicate that the compounds 1 and 4 changed the enzyme conformation. Furthermore, compound 4 significantly suppressed the increases in postprandial blood glucose levels in the C57BL/6J mice.

Introduction

Diabetes mellitus is chronic metabolic disease with worldwide concerns, which causes a major challenge for the health system (Kharroubi and Darwish, 2015). The high prevalence of diabetes has focalized much efforts for novel therapeutic alternatives (Ghosh et al., 2016). Nowadays, alleviating postprandial hyperglycemia is one of the first-line therapeutical strategies for the treatment of diabetes and its complications (Taylor et al., 2021). α-Glucosidase inhibitors (AGIs), such as acarbose, miglitol, and voglibose, are usually employed for controlling postprandial blood glucose levels by delaying the intestinal digestion of carbohydrates (Hossain et al., 2020). However, utilization of clinical AGIs often have some shortcomings such as side-effects including abdominal discomfort and flatulence, limited efficacy, failure in metabolism adjustment (Calcutt et al., 2009). Therefore, much effort has been focused on searching for natural AGIs with better safety and efficacy from natural sources in the past decade (Deng et al., 2015; Zhang et al., 2020).

Aspergillus terreus ML-44 is marine-derived fungi previously isolated from the fresh gut of pacific oyster. Our former study reported five terretonins isolated from ML-44 fermentation, including a new one, which showed weak anti-inflammatory activity (Wu et al., 2019). In order to exploit the potential of strain ML-44 in the medical field, the diethyl sulfate (DES) mutagenesis strategy (Fang et al., 2014) was applied to strain ML-44 in this study, mutant strain ASM-1 was screened out with different phenotypic morphology of colonies. HPLC-DAD-UV analysis of the mutant fermentations comparing to parent strain exhibited a series of metabolites with unique ultra-violet absorption were observed in the mutant ASM-1. Subsequent HPLC-guided chemical investigation of ASM-1 fermentation resulted in the isolation of six aspulvinone derivatives (16), including three new ones (13) (Figure 1). Aspulvinones and the analogs have been reported with various biological activities, such as inhibiting antibacterial (Machado et al., 2021), luciferase (Cruz et al., 2011), anti-influenza A viral (Gao et al., 2013), anticancer (Sun et al., 2019), anti-DPPH radicals (Zhang et al., 2015), as well as α-glucosidase inhibitory activity (Dewi et al., 2014; Wang et al., 2016; Zhang et al., 2016). However, there has been no systematic report on the mechanism for inhibition of α-glucosidase, structure–activity relationships, and hpyerglycemic effect in vivo by natural aspulvinones. In this study, the α-glucosidase inhibitory activities of compounds 16 were evaluated in vitro, in silico, and in vivo. Herein, we report the isolation, structure elucidation, and the α-glucosidase inhibitory activities of the isolated aspulvinones.

FIGURE 1
www.frontiersin.org

FIGURE 1. Chemical structures of Compounds 16.

Material and Methods

General Experimental Procedures

Sephadex™ LH-20 (GE Healthcare, Uppsala, Sweden), and YMC*GEL® ODS-A-HG (12 nm S-50 μm, YMC Co., Ltd., Kyoto, Japan) were used for column chromatography. The MPLC was performed on a QuikSep chromatographic system (H&E, Beijing, China), and a Gemini C18 column (21.2 × 250 mm, column temperature: 26°C) was used for separation and purification. Optical rotations were measured on a JASCO P-2000 digital polarimeter (JASCO, Tokyo, Japan). UV spectra were recorded on a PerkinElmer Lambda 25 spectrophotometer. Electronic Circular Dichroism (ECD) data were taken on a Chirascan circular dichroism spectrometer (Applied Photophysics, Surrey, United Kingdom). HR-ESI-MS was measured on Agilent 6520 Q-TOP mass spectrometer (Agilent, CA, United States), and all 1D and 2D NMR spectra were obtained on a Bruker-500 (500 MHz 1H and 125 MHz 13C-NMR) NMR spectrometer. A SynergyHTX micro plate reader (BioTek, VT, United States) was used to read optical density (OD). The intrinsic fluorescence spectra (280–500 nm) were measured using Perkin Elmer LS55 fluorescence spectrophotometer (United Kingdom).

Chemical Mutagenesis of A. terreus ML-44 and Mutant Selection

The DES mutagenesis procedure was referred to the method that we previously reported (Fang et al., 2014), and with proper modifications: DES was dissolved in DMSO to obtain a 20% (v/v) solution, which was further mixed with spore suspension of A. terreus ML-44 in a ratio of 1:9 (v/v). The mixture was treated with assistance of ultrasonic wave (40 KHz) at room temperature. Each 80 μl portion of the treated spore suspensions was sampled and spread on PDA plates at 1 and 2 h of treatment followed by incubation at 28°C for 5–7 days. Mutants from the test groups were obtained by selection of colonies with different colonial morphology, and the genetic stability were verified by passing three generations.

The initial ML-44 strain and mutants were activated by incubation at 28°C for 3–5 days, and further inoculated into 100 ml of liquid medium (glucose 2%, maltose 1%, mannitol 2%, glutamic acid 1%, peptone 0.5%, and yeast extract 0.3% in distilled water) in an Erlenmeyer (250 ml) and fermented at 28°C on a rotary shaker at 200 rpm for 12 days. Each 100 ml of the fermentation broth was extracted with equal volumes of EtOAc with assistance of ultrasonic wave (40 KHz) for 30 min. The EtOAc extractions were concentrated in vacuo at 37°C, followed by re-dissolved in 1.0 ml methanol which were used for further chemical analysis. HPLC-PDAD-UV analysis was performed using an analytical Kromasil C18 column (5 μm, 100 Å, 4.6 × 250 mm; Akzo Nobel) on an Agilent 1100 HPLC system equipped with photo-diode array detector (G1316A). After filtered with 0.45 μm membrane, the extraction solution in methanol (1.0 ml) was injected (10 μl) into the column and eluted with a MeOH-H2O linear gradient (20%→100% MeOH in 30 min followed by 5 min with isocratic 100% MeOH) mobile phase (flow rate 1 ml/min). The acquired PDAD data were processed with Agilent OpenLAB software.

Chemical Investigation to Mutant ASM-1

The mutant ASM-1 was inoculated into ten Erlenmeyer (500 ml) each containing 200 ml of sterile liquid medium and cultured at 28°C for 48 h on a rotary shaker at 200 rpm providing a seed culture (2 L). The seed culture was inoculated into a fermentation cylinder containing the same sterile liquid medium (70 L), and was cultured at 28°C for 12 days, with sterile compressed air passes from the bottom of the cylinder keeping a positive pressure of 0.15 MPa. The whole broth (65 L) was filtered to separate into the filtrate and the mycelial cake. The filtrate (60 L) was subjected to an AB-8 macroporous resin column (column volume, CV 2.4 L), eluted by water and 95% ethanol successively. The water elute (3 CVs) was discarded, and the 95% ethanol elute (3 CVs) was gathered. The mycelial cake was extracted two times with 95% ethanol (5 L) assisted by ultra-sonication for 2 h, followed by filtration giving the ethanol extract. All the ethanol solutions were combined and concentrated to a water suspension, which was further extracted with EtOAc to afford a total of 60.5 g of EtOAc extract.

The EtOAc extract (60.5 g) of the mutant ASM-1 was subjected to silica gel column chromatography by stepwise elution with b. p. 60–90°C petroleum ether (P)- dichloro-methane (D)-methanol (M) to obtain 9 fractions. HPLC analysis showed that the newly produced metabolites were contained in Fr-4 (1.9 g, eluted by D), Fr-5 (5 g, eluted by DM 98:2), and Fr-7 (2.3 g, eluted by DM 90:10) (Supplementary Figure S3). Subsequent repeated preparative reverse phase HPLC separation led to the purification of compounds. Fr-4 was subjected to reduced pressure ODS column chromatography (cc) to give subfraction Fr-4–8 (0.3 g, eluted by 90%M), which was further separated by HPLC (Methanol-H2O (0.1%HCl) 90:10, 10 ml/min) to afford 1 (25 mg, tR = 36.0 min) and 6 (2.1 mg, tR = 24.5 min); Fr-5–4 (1.6 g, eluted by 80%M), the subfraction of Fr-5 by ODS cc, afforded 3 (12 mg, tR = 64.0 min) and 4 (218 mg, tR = 67.5 min) with preparative HPLC separation (Methanol-H2O (0.1%HCl) 65:35, 10 ml/min); Fr-7 was separated by ODS cc to give subfraction Fr-7–9 (0.6 g, eluted by 80%M), which provided 2 (28 mg, tR = 28.4 min) and 5 (16 mg, tR = 18.6 min) after preparative HPLC separation (Methanol-H2O (0.1%HCl) 77:23, 10 ml/min).

Aspulvinone V (1): yellow solid (MeOH), UV (MeOH) λmax (log ε): 203 (4.36), 239 (4.11), 380 (4.23). Positive HR-ESI-MS: m/z measured 499.2480 [M + H]+, calcd for C32H35O5 [M + H]+ 499.2484. 1H NMR and 13C NMR spectroscopic data, see Table 1.

TABLE 1
www.frontiersin.org

TABLE 1. 500 MHz 1H and 125 MHz13C NMR data of compounds 13, 5 in CD3OD.a

Aspulvinone W (2): yellow solid (MeOH) [α]D24‒32.8 (c 0.23, MeOH). UV (MeOH) λmax (log ε): 203 (4.35), 243 (3.98), 374 (4.15). Positive HR-ESI-MS: m/z measured 449.1963 [M + H]+, calcd for C27H29O6 [M + H]+ 449.1964. 1H NMR and 13C NMR spectroscopic data, see Table 1.

Aspulvinone X (3): yellow solid (MeOH) [α]D24‒27.9 (c 0.26, MeOH). UV (MeOH) λmax (log ε): 203 (4.35), 242 (3.98), 376 (4.16). Positive HR-ESI-MS: m/z measured 449.1958 [M + H]+, calcd for C27H29O6 [M + H]+ 449.1964. 1H NMR and 13C NMR spectroscopic data, see Table 1.

α-Glucosidase Inhibitory Assay

α-Glucosidase (EC:3.2.1.20, MAL12) from Saccharomyces cerevisiae was dissolved in 0.1 mol/L PBS solutions with a pH of 6.8, and diluted to be a 1.0 U/ml solution. The substrate p-nitrophenyl-β-D-glucopyranoside (pNPG) was dissolved in PBS to be a 1 mM solution. Acarbose and the compounds were dissolved in mehanol and further diluted to a series of concentrations from 0.1 μmol/L to 10 mmol/L. In vitro α-glucosidase inhibitory assay was performed according to a method described previously with some modification (Dan et al., 2019). Briefly, 20 μl of 1.0 U/ml enzyme solution and 10 μl of acarbose or compound solution, was mixed with 50 μl PBS solution in 96-well plate, and the mixed solution was incubated at 37°C for 10 min 20 μl of 1 mmol/L pNPG was subsequently added and further incubated at 37°C for 15 min, after which 100 μl of 1 M Na2CO3 solution was added to terminate the reaction. The absorbance of p-nitrophenol was monitored at 405 nm. All samples were analysed in triplicate, and acarbose was used as positive control. The negative control was performed by adding PBS instead of α-glucosidase, the blank was prepared by adding solvent without tested compounds. The inhibition rate was calculated as Eq. 1:

IR% =[(Ac-As)/Ac] × 100%(1)

where Ac represents the absorbance of control without sample solution, and As denotes the absorbance of sample.

Enzymatic Kinetics of α-Glucosidase

pNPG with a concentration range of 100–4,000 μM and α-glucosidase were incubated with different concentrations of inhibitor for 10 min, respectively. 20 μl of 1.0 U/ml enzyme solutions were first mixed with 10 μl of different concentrations of inhibitors, then 50 μl PBS solutions were added, and the mixed solutions were incubated at 37°C for 10 min. Subsequently, 20 μl of pNPG solutions (1.25, 2.5, 5, 10 and 20 mM) were added, and the mixed solutions were further incubated at 37°C for 25 min, the absorbance of reaction solution was measured at 405 nm every 3 min. The kinetic parameters, Michaelis–Menten (Km) and maximum velocity (Vmax), were found using Lineweaver–Burk plots to check the mode of α-glucosidase inhibition for compounds 1 and 4. The dissociation constants between inhibitor and enzyme (Ki) were calculated from Dixon plots. Two inhibition constants, KI or KIS, for inhibitor binding with either free or enzyme-substrate complex, were calculated from secondary plots of the slopes of the straight lines (Vmax/Km) or vertical intercept (1/Vmax) verse the concentration of inhibitors, respectively (Jenis et al., 2019).

Fluorescence Quenching Analysis

α-Glucosidase (1 μM) was mixed with different concentrations of inhibitors (0–15 μM) at 20, 31, and 37°C, respectively, and the fluorescence spectra of mixed solutions were determined after equilibration for 5 min at an excitation wavelength of 250 nm. Both the excitation and emission slits were set at 10 nm. The quenching rate constant (Kq), binding constant (Ka), the number of the binding sites (n), and thermodynamic parameters enthalpy change (ΔH) and entropy change (ΔS) were calculated according to the Stern-Volmer Eq. 2 and the van’t Hoff Eqs 35, which were listed as follows (Xu et al., 2019):

F0/F = 1 + Ksv[Q] = 1 + Kqτ0[Q](2)
log((F0-F)/F) = logKa+ nlogQ(3)
lnKa= (-1/T(ΔH/R)) + ΔS/R(4)
Δ= ΔH-TΔS(5)

Where F0 and F represent the fluorescence intensities in the absence or presence of inhibitor, [Q] denote the concentration of inhibitor, τ0 is the constant of the lifetime of fluorophore (10–8 s) and R is the gas constant of 8.31 J/(mol × K).

Circular Dichroism Spectroscopy

The CD measurements were performed in a wavelength range of 190–250 nm at a speed of 60 nm/min. All measurements were carried out at 20°C using 1.0 mm path length quartz cuvette and sodium phosphate buffer (pH 6.8) was considered as a blank. The concentration of α-glucosidase was 1.25 μM, whereas the molar ratios of inhibitors (25 and 50 μM) to α-glucosidase were 20:1 and 40:1. All the results were expressed as ellipticity in mill degrees.

Molecular Docking

Autodock Vina software was used in docking calculations to investigate the modes of glucosidase inhibition for individual aspulvinones. The 3D structures of aspulvinones were generated and then energetically minimized with MM2 force field to a minimum Root Mean Square (RMS) gradient of 0.005 using Chem3D Ultra 2017 (Version 17.0.0.206). The crystal structure of Saccharomyces cerevisiae isomaltase (PDB ID: 3A4A; Resolution 1.6 Å) was retrieved from protein Data Bank (www.rcsb.org/pdb/), and was prepared by removing the water molecules and original inhibitors. The ligand and protein pdb files were prepared and grid box formation was accomplished using AutoDock Tools. AutoGrid was used in order to prepare the grid map using a grid box. The dimension of the grid size was set to 50 × 50 × 50, and the grid box center was designated in coordinates x = 21.285, y = −0.64, and z = 18.475. Nine output poses were generated and evaluated by their calculated free energy of binding. The best pose of each ligand was determined by its Affinity score (kcal/mol), which was visualized by Discovery Studio Visualizer v21.1.0.20298 (Accelrys, San Diego, United States) to analyze the interactions between the target enzyme and the inhibitor.

Oral Disaccharide Tolerance Test

Female C57BL/6J mice, 6 weeks old, weighting 16–20 g, were obtained from the Skbex Biotechenology Company Limited (Henan, China). The animals were housed in Experimental Animal Center of Zhoukou Normal University under 12 h light-dark cycle at controlled temperature (22 ± 1°C), and provided with a standard pellet diet and water ad libitum. The mice were adapted to diet and general conditions of vivarium for 1 week before the experiment. After an 16 h fasting, the mice were divided into three groups randomly (eight mice each group). Sucrose or maltose, as well as the inhibitors (compound 4 and acarbose), was dissolved in 0.5% sodium carboxymethyl cellulose (CMC–Na) solution. Compound 4 was tested at dose of 25 mg/kg, whereas acarbose was evaluated at dose of 50 mg/kg of body weight (BW). Mice were fasted 16 h and then intragastrically (i.g.) administrated the inhibitors, 15 min later, 2 g/kg BW of sucrose or 2 g/kg BW of maltose solution was given i.g. to the animals. Blood samples were taken from the tail vein at 0, 30, 60, and 120 min after sucrose or maltose loading, and blood glucose was measured with Accu-Chek Active glucometer (Roche, Germany). Area under the curve (AUC) over 0–120 min was calculated by the trapezoidal method. All animals were cared under the frame of the China Council on Animal Care and all procedures were approved by the Health Sciences Animal Welfare Committee of Zhoukou Normal University.

Results and Discussion

Isolation and Structure Determination of Aspulvinones

The mutant ASM-1 was obtained with different colonial morphologies with the parent strain ML-44 (Supplementary Figure S1), through treatment of ML-44 spores with 2% (v/v) DES under ultrasonic assistance for 1 h, and its EtOAc extract HPLC profile showed a series of chromatographic peaks of newly produced secondary metabolites with unique UV absorption spectra comparing to the wild strain ML-44 (Supplementary Figure S2). This mutant was deposited at the China General Microbiological Culture Collection Center under the accession number CGMCC No. 22417. In order to clarify the newly produced secondary metabolites, large scale fermentation of mutant ASM-1 and HPLC-guided separation were performed (Supplementary Figure S3). The EtOAc extract of the mutant ASM-1 was separated by silica gel column chromatography and repeated preparative reverse phase HPLC separation under HPLC-PDAD-UV monitoring, resulted in the isolation of aspulvinones 16. Compound 4 was determined to be aspulvinone H by comparison its MS and NMR data with that of literature (Nagia et al., 2012). Compound 6 was identified as aspulvinone R, which was recently isolated from a marine sponge-associated fungus A. flavipes KUFA1152 as the first example of triprenylated aspulvinones (Machado et al., 2021).

Compound 1 was obtained as yellow solid. The molecular formula was determined to be C32H34O5 on the basis of a HRESIMS peak at m/z 499.2480 [M + H]+ (calcd. 499.2484), indicating 16 degrees of unsaturation. The illustration of 13C NMR, DEPT, and HSQC spectra came up with 32 resonances, which were indicative of one ketone carbonyl, twelve sp2 and one sp3 quaternary carbons, ten sp2 methine, two sp3 methylene, six methyl groups in 1. The remaining unsaturation was thus attributed to four rings. The 1H NMR spectroscopic data of 1 (Table 1) showed a series of protons signals, and were affiliated to relevant carbons via HSQC analysis. The detailed 1D (Table 1) and 2D (Figure 2) NMR analyses of 1 indicated the same pulvinone nucleus as 4 and 6. The 1H NMR signals δ6.36 (d, J = 9.7 Hz, H-18) and δ5.68 (d, J = 9.7 Hz, H-19) indicated that one of the prenyl occurred cyclization, which was verified by the HMBC correlations between the above two protons with related carbons (Figure 2). In addition, there were two linear prenyl groups (C-1′ to C-5′, and C-1″ to C-5″) in the molecular of 1 according the NMR data analysis. Thus, compound 1 should also be a triprenylated pulvinone. Subsequently, the HMBC correlations between H-18 to C-13 and C-15, H-1″ to C-17 and H-2″ to C-16, H-19 to C-14, H-13 and H-17 to C-2, indicated that a 1,3,4,5-tetrasubstituted benzene ring bound to the γ-butenolide core directly at C-2. Thus, the another one prenylated benzene ring linked to the core via C-5, which was also confirmed by relative HMBC signals. At this stage, the planar structure has been constructed as 1. The relatively small chemical shifts of C-2 (δC 102.2) and C-5 (δC 109.5) established the Z geometry of the ∆4,5-double bond (Campbell et al., 1985), which was the same as that of compound 6. Since 1 has never been reported previously, it was named as aspulvinone V in the order of the names for this series of prenylated pulvinones from A. terreus.

FIGURE 2
www.frontiersin.org

FIGURE 2. Key COSY and HMBC correlations of compounds 13, and CD spectra of compound 2, 3, and 5.

Compound 2 was obtained as yellow solid. Its molecular formula was established as C27H28O6 on the basis of a HRESIMS peak at m/z 449.1963 [M + H]+ (calcd. 449.1964), which indicated 14 degrees of unsaturation and 16 amu more than compound 4 (C27H28O5). The 1H and 13C NMR data (Table 1) indicated that the structure of 2 is very similar to compound 4. The most significant differences in the NMR data exist in the high-field shift effect of the signal at C-1′ (δH 3.20, 3.15, δC 31.2), and the presence of a O-substituted sp3 methine instead of a sp2 methine at C-2′ (δH 4.59, δC 91.1), a O-substituted sp3 quaternary carbon instead of a sp2 one at C-3′ (δC 72.4). These data indicated that compound 2 bears a dihydrofuran ring fused to the benzene ring, as opposed to the linear prenyl present in compound 4. This structure was further supported by COSY, HSQC, and HMBC spectra. COSY correlations between H2-1′ and H-2′, as well supported the presence of one dihydrofuran ring systems. HMBC correlation from H2-1′ to C-8 and C-9, H-10 to C-8, H-7, and H-11 to C-9 confirmed the presence of a dihydrofuran fused to one of the benzene rings. Furthermore, HMBC signals of H-5 to C-7 and C-11 indicated the dihydrofuran fused benzene ring was linked to C-5.

Compound 3 was isolated as yellow solid with the same molecular formula as 2 base on HRESIMS peak at m/z 449.1958 [M + H]+ (calcd. 449.1964). The 1H and 13C NMR data of 3 (Table 1) were similar to those of 2 except for the high-field shift effect of the signal at C-18 (δH 2.81, δC 23.4), the presence of a sp3 methylene instead of a sp2 methine at C-19 (δH 1.82, δC 33.8), and a O-substituted sp3 quaternary carbon instead of a sp2 one at C-20 (δC 75.5). These data indicated that compound 3 bears a tetrahydropyran ring fused to the benzene ring, instead of the linear prenyl present in compound 2. The location of furan ring and pyran ring were confirmed through HMBC analysis (Figure 2). In addition, both 2 and 3 have the Z geometry for the ∆4,5-double bond according to 13C NMR data. The transisomer of compound 3 (trans-3) has been reported previously, which was determined the absolute configuration as R by the ECD calculations (Sun et al., 2018). The CD spectra of 2 and 3 was tested, and the negative cotton effect at 350 nm and general spectrum shapes were both consistent with trans-3 (Figure 2), indicating the R configuration at C-2′ for compounds 2 and 3.

By illustration of 1D and 2D NMR data, compound 5 was deduced to have the same structure with aspulvinone J-CR, whereas there exists considerable discrepancy between our 13C NMR data and that of the literature (Cruz et al., 2011), especially for C2–C5 and C12. Taking consideration of the relatively strong acidity of 4-OH (predicted pKa value of 4.50 ± 1.00 calculated by Advanced Chemistry Development (ACD/Labs) Software V11.02), it is presumed that pulvinone derivatives may be dissociative in relative high pH solutions, which would cause changes in NMR and UV spectral characteristics. Based on this assumption, 5 sodium was prepared by adding NaHCO3 aqueous solution to 5 methanol solution in molar ratio 1:1. The 13C NMR data of 5 sodium exhibited significant difference with 5 (Figure 3A), while it is identical with the literature data. Therefore, the NMR data for aspulvinone J-CR and other analogs in the literature should be for their sodium. In addition, the UV spectra of the 5 and its sodium represented different maximum absorption peaks both in number and intensity (Figure 3B), 5 has one absorption peak at 376 nm, while 5 sodium exhibited two absorption peaks at 327 and 376 nm, and the latter with a relatively low absorbance exists as a shoulder peak of the former.

FIGURE 3
www.frontiersin.org

FIGURE 3. Spectral comparisons between 5 and 5 sodium. (A)13C-NMR data differences; (B) UV spectra of 5 and 5 sodium.

α-Glucosidase Inhibitory Activities

All compounds 16 showed potent inhibitions towards α-glucosidase with IC50s ranging from 2.2 to 44.3 μM (Table 2). It was reported that the transisomer of compound 3 inhibited α-glucosidase with IC50 of 24.8 μM (Sun et al., 2018). The inhibitory potencies varied with the modification of benzene rings. In term of diprenylated aspulvinones, compound 4 (IC50 4.6 μM) is appropriate ten times more potent than 2, 3 and its transisomer, and 5, indicating that the linear prenyl is significant to the α-glucosidase inhibitory activity, while the configuration of ∆4,5-double bond has little influence. However, Aspulvinone E (IC50 2.70 μM) was reported with higher inhibitory activity than its transisomer, isoaspulvinone E (IC50 8.92 μM), and the ∆4,5-double bond stereochemistry significantly affected the inhibition activity to α-glucosidase for non-prenylated pulvinones (Dewi et al., 2014). Therefore, we presumed that non-prenylated and prenylated pulvinones had different binding modes with α-glucosidase. On the other hand, both of compound 1 (IC50 2.2 μM) and 4 were more potent than 6 (IC50 10.8 μM). It was speculated that the fork-like structure of two linear prenyl groups has steric hindrance effect hampering the binding of 6 with enzyme.

TABLE 2
www.frontiersin.org

TABLE 2. α-Glucosidase inhibitory activity of aspulvinones 16.

Compounds 1 and 4, the most potent inhibitors, were selected for enzyme kinetic studies to elucidate the inhibition mode. In the Lineweaver–Burk double-reciprocal plots, as shown in Figure 4, the plots of 1/V versus 1/[S] give a group of straight lines with different slopes that intersect at the second quadrant for both of 1 and 4, suggesting that both of them are mixed-type inhibitors (Dan et al., 2019). Therefore, both compounds could bind to free enzyme (EI), and interfere with the formation of the α-glucosidase-pNPG (ES) intermediate through forming an α-glucosidase-pNPG-inhibitor (ESI) complex (Wikul et al., 2012; Wu et al., 2014). The inhibition constant for the inhibitor binding with free enzyme (KI) was determined by a plot of the slope (Km/Vm) versus the inhibitor concentration, and the inhibition constant for the inhibitor binding with enzyme–substrate complex (KIS) was obtained from the vertical intercept (1/Vm) versus the inhibitor concentration (Supplementary Figure S5) (Sheng et al., 2018). The results are shown in Table 2: the KI values of both 1 and 4 are smaller than their KIS values, which suggest that them have higher affinity with the free enzyme than with the enzyme-substrate complex.

FIGURE 4
www.frontiersin.org

FIGURE 4. Inhibition mechanism of 1 and 4 against α-glucosidase. Lineweaver–Burk plot of 1 (A) and 4 (D). Fluorescence spectra of α-glucosidase with 1 (B) and 4 (E). CD spectra of α-glucosidase with 1 (C) and 4 (F).

Fluorescence Quenching Mechanism and Binding Characterizations

Subsequently, the interaction between the inhibitors and α-glucosidase was investigated by fluorescence spectroscopy and circular dichroism (CD) spectroscopy. As shown in Figure 4, with the increased concentrations of 1 and 4, the intrinsic fluorescence intensity of α-glucosidase decreased gradually, indicating that the inhibitors interacted with α-glucosidase and then quenched its intrinsic fluorescence. Comparing to the maximum scattering collision quenching constant of the biomacromolecule (2 × 1010 L/mol/s), the quenching rate constants (Kq) were much larger, demonstrating that the fluorescence quenching process was static quenching predominantly (Table 3). The number of binding sites (n) were all close to one at the three incubation temperatures, indicating that both 1 and 4 interact with α-glucosidase at only one binding site. The binding constants (Ka) at the three temperatures were in the order of 105 and 104 L/mol for 1 and 4, respectively, indicating that there were high binding affinities existed in the complex of α-glucosidase with the both compounds, especially as for 1. In addition, the thermodynamic parameters (∆S, ∆H and ∆G) were calculated, showing that ∆H and ∆S were positive, while ∆G was negative (Table 3). The binding process could be defined to be thermodynamically favorable and spontaneous, which was driven mainly by a hydrophobic force (Ross and Subramanian, 1981). In the CD spectroscopy analysis, after the addition of 25 and 50 mM of 1 or 4, the absorption of the two negative peaks at 209 and 222 nm decreased, which demonstrated a loss of the α-helix structure (Figure 4) (Xu et al., 2019). In addition, with an increase in molar ratios of inhibitors to α-glucosidase (from 20:1 to 40:1), the loss of the α-helix structure increased, associating with the decreased α-glucosidase activity. The transformation from the α-helix to other conformations in the presence of inhibitors indicated a partial unfolding of the α-glucosidase structure, causing alterations of the secondary structure of α-glucosidase, and thereby some hydrogen bonding networks might be destroyed. These alterations may prevent the binding of the substrate to α-glucosidase or hamper the formation of an active center, eventually resulted in the dysfunction of the enzyme (Xu et al., 2019).

TABLE 3
www.frontiersin.org

TABLE 3. Quenching constants (Ksv), binding constants (Ka), number of binding sites (n), and thermodynamic parameters for the α-glucosidase-inhibitor system.

Molecular Docking

Since there is no crystal structure for the commercially available Saccharomyces cerevisiae α-glucosidase for preparing the protein for docking, the constructed homology models based on the isomaltases or itself were often used to perform the molecular docking (Xu et al., 2019; Khosravi et al., 2020). In this study, the crystal structure of isomaltase from Saccharomyces cerevisiae (PDB ID: 3A4A; Resolution 1.6 Å) was adopted for silico docking of to confirm the interaction. Compounds 1 and 4 exhibits a strong binding affinity with the protein by the low binding energy of −10.9 and −9.6 kcal/mol, respectively. As shown in Figure 5, both compounds could bond at the gate of the hydrophobic pocket, and partially inserted into the binding pocket. In this bonding mode, the ligands could hamper the substrate loading into the catalytic pocket in EI complex formation, or cause structural modification of α-glucosidase leading to the dysfunction in ESI complex formation (Kim et al., 2005). According to the molecular docking results, the binding pocket involves the amino acid residues Asp307, Pro312, Tyr158, Thr310, Arg315, and Lys156 for the both inhibitors, and additional Phe303, Phe314 and Val319 for 1, whereas additional Val308, Ser311, and His280 for 4. There were the hydrogen bond interactions between the carbonyl group of 1 with Thr310 and Arg315 (the distance: 2.26 and 2.62 Å), and the rest interactions were all hydrophobic effect including alkyl, Pi-alkyl and Pi-anion (Figure 5C). In comparison, there are relatively more hydrogen bonds and less hydrophobic interaction for 4 (Figure 5D), indicating the different binding force compositions between the two inhibitors.

FIGURE 5
www.frontiersin.org

FIGURE 5. Docking binding model of inhibitors with Saccharomyces cerevisiae isomaltase (3A4A). Predicted dock conformation of the isomaltase to inhibitors 1 (A) and 4 (C); 2D interaction diagrams between 3A4A and inhibitors 1 (B) and 4 (D).

Effect of Compound 4 on Postprandial Hyperglycemia in vivo

The intestinal α-glucosidase inhibitory activity in vivo was evaluated by oral sucrose and oral maltose tolerance tests in female C57BL/6J mice. Acarbose (50 mg/kg BW) was used as a positive control, and compound 4 was chosen for its potent inhibitory activity and high yield. In the oral sucrose tolerance test (Figure 6A), after oral administration of sucrose (2 g/kg of BW), the blood glucose level rapidly increased from 5.00 ± 0.07 mM to a maximum of 14.24 ± 0.45 mM in 30 min, and then recovered to the pretreatment level at 120 min. In the treatment group, 4 significantly suppressed the blood glucose rise at 30 and 60 min comparing to that of the negative control group, and led to 13.2% decrease of the AUC at a dose of 25 mg/kg BW comparable to that of acarbose (11.8% decrease) at dose of 50 mg/kg BW (Figure 6B). Similarly, in the sucrose tolerance test, compound 4 treatment resulted in a significant decrease in the postprandial blood glucose peak versus the negative control group (Figure 6C) and he AUC for postprandial plasma glucose was reduced by 19.7% in 2 h after 4 administration, which was more potent than acarbose (16.2%) (Figure 6D). These results strongly confirmed that 4 could alleviate the postprandial hyperglycemia through inhibiting intestinal α-glucosidase. Therefore, natural aspulvinones can be regarded as potential candidate for hpyerglycemic agents.

FIGURE 6
www.frontiersin.org

FIGURE 6. Effects of acarbose and 4 on postprandial blood glucose levels in female C57BL/6J mice. (A) Blood glucose concentrations after oral administration of 4 and sucrose. (B) Area under the curve (AUC) after oral administration of 4 and sucrose for 2 h. (C) Blood glucose concentrations after an oral administration of 4 and maltose. (D) Area under the curve (AUC) after oral administration of 4 and maltose for 2 h. Each value represents the mean standard deviation (n = 8). Asterisks indicate a significant difference (*p < 0.05, **p < 0.01, and ***p < 0.001) compared with the control group.

Conclusion

In summary, we conducted the DES mutagenesis on the marine-derived A. terreus ML-44, and a mutant strain ASM-1 was obtained by morphological and HPLC analyses. Six aspulvinone secondary metabolites were isolated from the ASM-1 culture, including three new ones. Their structures including the absolute configurations were elucidated by various spectroscopic methods and ECD comparison. All compounds were evaluated for α-glucosidase inhibitory activity with acarbose as positive control. Among them, compounds 1 and 4 exhibited potent α-glucosidase inhibitory activities with IC50 values of 2.2 and 4.6 µM in mixed-type manners. The thermodynamic and molecular docking studies recognized the interaction between inhibitors and α-glucosidase was spontaneous exothermic reaction driven mainly by hydrophobic forces. Furthermore, 4 significantly suppressed the increases in postprandial blood glucose levels in the C57BL/6J mice more potently than acarbose at a smaller dosage. The results suggested that aspulvinones could be promising candidates for further pharmacologic research. In addition, the mechanism of the mutagenesis of the strain ASM-1 from strain ML-44 deserve further investigation through genome and transcriptome analyses, which may make contribution to understanding the metabolic regulation of aspulvinones biosynthesis.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.

Ethics Statement

The animal study was reviewed and approved by the Health Sciences Animal Welfare Committee of Zhoukou Normal University.

Author Contributions

CW and XC performed isolation, structure determination and bioassays of the compounds and wrote the manuscript. LS and JL performed the mutation and fermentation of the fungus, extraction of the culture broths and isolation of the compounds. FL, MS, and YZ. performed the bioassays. XH carried out spectroscopic tests. CT illustrated the inhibitory activity data and revised the manuscript. MS and XL designed the study and revised the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (No. 81803425), China Postdoctoral Science Foundation (No. 2019M662552).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

We thank Jinwei Ren from Institute of Microbiology, Chinese Academy of Sciences, for his professional technological support on NMR and HRESIMS tests.

Supplementary Material

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

References

Calcutt, N. A., Cooper, M. E., Kern, T. S., and Schmidt, A. M. (2009). Therapies for Hyperglycaemia-Induced Diabetic Complications: from Animal Models to Clinical Trials. Nat. Rev. Drug Discov. 8, 417–430. doi:10.1038/nrd2476

PubMed Abstract | CrossRef Full Text | Google Scholar

Campbell, A. C., Maidment, M. S., Pick, J. H., and Stevenson, D. F. M. (1985). Synthesis of (E)- and (Z)-pulvinones. J. Chem. Soc. Perkin Trans. 1 1, 1567–1576. doi:10.1039/P19850001567

CrossRef Full Text | Google Scholar

Cruz, P. G., Auld, D. S., Schultz, P. J., Lovell, S., Battaile, K. P., MacArthur, R., et al. (2011). Titration-based Screening for Evaluation of Natural Product Extracts: Identification of an Aspulvinone Family of Luciferase Inhibitors. Chem. Biol. 18, 1442–1452. doi:10.1016/j.chembiol.2011.08.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Dan, W.-J., Zhang, Q., Zhang, F., Wang, W.-W., and Gao, J.-M. (2019). Benzonate Derivatives of Acetophenone as Potent α-glucosidase Inhibitors: Synthesis, Structure-Activity Relationship and Mechanism. J. Enzyme Inhib. Med. Chem. 34, 937–945. doi:10.1080/14756366.2019.1604519

PubMed Abstract | CrossRef Full Text | Google Scholar

Deng, Y.-T., Lin-Shiau, S.-Y., Shyur, L.-F., and Lin, J.-K. (2015). Pu-erh tea Polysaccharides Decrease Blood Sugar by Inhibition of α-glucosidase Activity In Vitro and in Mice. Food Funct. 6, 1539–1546. doi:10.1039/c4fo01025f

PubMed Abstract | CrossRef Full Text | Google Scholar

Dewi, R. T., Tachibana, S., Fajriah, S., and Hanafi, M. (2014). α-Glucosidase Inhibitor Compounds from Aspergillus terreus RCC1 and Their Antioxidant Activity. Med. Chem. Res. 24, 737–743. doi:10.1007/s00044-014-1164-0

CrossRef Full Text | Google Scholar

Fang, S.-M., Wu, C.-J., Li, C.-W., and Cui, C.-B. (2014). A Practical Strategy to Discover New Antitumor Compounds by Activating Silent Metabolite Production in Fungi by Diethyl Sulphate Mutagenesis. Mar. Drugs 12, 1788–1814. doi:10.3390/md12041788

PubMed Abstract | CrossRef Full Text | Google Scholar

Gao, H., Guo, W., Wang, Q., Zhang, L., Zhu, M., Zhu, T., et al. (2013). Aspulvinones from a Mangrove Rhizosphere Soil-Derived Fungus Aspergillus terreus Gwq-48 with Anti-influenza A Viral (H1N1) Activity. Bioorg. Med. Chem. Lett. 23, 1776–1778. doi:10.1016/j.bmcl.2013.01.051

PubMed Abstract | CrossRef Full Text | Google Scholar

Ghosh, A. K., Reddy, B. S., Yen, Y.-C., Cárdenas, E. L., Rao, K. V., Downs, D., et al. (2016). Design of Potent and Highly Selective Inhibitors for Human β-secretase 2 (Memapsin 1), a Target for Type 2 Diabetes. Chem. Sci. 7, 3117–3122. doi:10.1039/c5sc03718b

PubMed Abstract | CrossRef Full Text | Google Scholar

Hossain, U., Das, A. K., Ghosh, S., and Sil, P. C. (2020). An Overview on the Role of Bioactive α-glucosidase Inhibitors in Ameliorating Diabetic Complications. Food Chem. Toxicol. 145, 111738. doi:10.1016/j.fct.2020.111738

PubMed Abstract | CrossRef Full Text | Google Scholar

Jenis, J., Baiseitova, A., Yoon, S. H., Park, C., Kim, J. Y., Li, Z. P., et al. (2019). Competitive α-glucosidase Inhibitors, Dihydrobenzoxanthones, from the Barks of Artocarpus Elasticus. J. Enzyme Inhib. Med. Chem. 34, 1623–1632. doi:10.1080/14756366.2019.1660653

PubMed Abstract | CrossRef Full Text | Google Scholar

Kharroubi, A. T., and Darwish, H. M. (2015). Diabetes Mellitus: the Epidemic of the century. Wjd 6, 850–867. doi:10.4239/wjd.v6.i6.850

PubMed Abstract | CrossRef Full Text | Google Scholar

Khosravi, A., Vaezi, G., Hojati, V., and Abdi, K. (2020). Study on the Interaction of Triaryl-Dihydro-1,2,4-Oxadiazoles with α-glucosidase. DARU J. Pharm. Sci. 28, 109–117. doi:10.1007/s40199-019-00322-y

CrossRef Full Text | Google Scholar

Kim, Y.-M., Jeong, Y.-K., Wang, M.-H., Lee, W.-Y., and Rhee, H.-I. (2005). Inhibitory Effect of pine Extract on α-glucosidase Activity and Postprandial Hyperglycemia. Nutrition 21, 756–761. doi:10.1016/j.nut.2004.10.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Machado, F. P., Kumla, D., PereiraSousa, J. A. E., Sousa, E., Dethoup, T., Freitas-Silva, J., et al. (2021). Prenylated Phenylbutyrolactones from Cultures of a marine Sponge-Associated Fungus Aspergillus flavipes KUFA1152. Phytochemistry 185, 112709. doi:10.1016/j.phytochem.2021.112709

PubMed Abstract | CrossRef Full Text | Google Scholar

Nagia, M. M., El-Metwally, M., Shaaban, M., El-Zalabani, S. M., and Hanna, A. G. (2012). Four Butyrolactones and Diverse Bioactive Secondary Metabolites from Terrestrial Aspergillus flavipes MM2: Isolation and Structure Determination. Org. Med. Chem. Lett. 2, 9. doi:10.1186/2191-2858-2-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Ross, P. D., and Subramanian, S. (1981). Thermodynamics of Protein Association Reactions: Forces Contributing to Stability. Biochemistry 20, 3096–3102. doi:10.1021/bi00514a017

PubMed Abstract | CrossRef Full Text | Google Scholar

Sheng, Z., Ge, S., Xu, X., Zhang, Y., Wu, P., Zhang, K., et al. (2018). Correction: Design, Synthesis and Evaluation of Cinnamic Acid Ester Derivatives as Mushroom Tyrosinase Inhibitors. Med. Chem. Commun. 9, 897. doi:10.1039/C8MD00099A10.1039/c8md90024h

CrossRef Full Text | Google Scholar

Sun, K., Zhu, G., Hao, J., Wang, Y., and Zhu, W. (2018). Chemical-epigenetic Method to Enhance the Chemodiversity of the marine Algicolous Fungus, Aspergillus terreus OUCMDZ-2739. Tetrahedron 74, 83–87. doi:10.1016/j.tet.2017.11.039

CrossRef Full Text | Google Scholar

Sun, W., Luan, S., Qi, C., Tong, Q., Yan, S., Li, H., et al. (2019). Aspulvinone O, a Natural Inhibitor of GOT1 Suppresses Pancreatic Ductal Adenocarcinoma Cells Growth by Interfering Glutamine Metabolism. Cell Commun Signal 17, 111. doi:10.1186/s12964-019-0425-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Taylor, S. I., Yazdi, Z. S., and Beitelshees, A. L. (2021). Pharmacological Treatment of Hyperglycemia in Type 2 Diabetes. J. Clin. Invest. 131, e142243. doi:10.1172/JCI142243

CrossRef Full Text | Google Scholar

Wang, C., Guo, L., Hao, J., Wang, L., and Zhu, W. (2016). α-Glucosidase Inhibitors from the Marine-Derived Fungus Aspergillus flavipes HN4-13. J. Nat. Prod. 79, 2977–2981. doi:10.1021/acs.jnatprod.6b00766

PubMed Abstract | CrossRef Full Text | Google Scholar

Wikul, A., Damsud, T., Kataoka, K., and Phuwapraisirisan, P. (2012). (+)-Pinoresinol Is a Putative Hypoglycemic Agent in Defatted Sesame (Sesamum indicum) Seeds Though Inhibiting α-glucosidase. Bioorg. Med. Chem. Lett. 22, 5215–5217. doi:10.1016/j.bmcl.2012.06.068

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, C.-J., Cui, X., Xiong, B., Yang, M.-S., Zhang, Y.-X., and Liu, X.-M. (2019). Terretonin D1, a New Meroterpenoid from marine-derived Aspergillus terreus ML-44. Nat. Product. Res. 33, 2262–2265. doi:10.1080/14786419.2018.1493583

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, P.-P., Zhang, K., Lu, Y.-J., He, P., and Zhao, S.-Q. (2014). In Vitro and In Vivo Evaluation of the Antidiabetic Activity of Ursolic Acid Derivatives. Eur. J. Med. Chem. 80, 502–508. doi:10.1016/j.ejmech.2014.04.073

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, Y., Xie, L., Xie, J., Liu, Y., and Chen, W. (2019). Pelargonidin-3-O-rutinoside as a Novel α-glucosidase Inhibitor for Improving Postprandial Hyperglycemia. Chem. Commun. 55, 39–42. doi:10.1039/c8cc07985d

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, L.-H., Feng, B.-M., Zhao, Y.-Q., Sun, Y., Liu, B., Liu, F., et al. (2016). Polyketide Butenolide, Diphenyl Ether, and Benzophenone Derivatives from the Fungus Aspergillus flavipes PJ03-11. Bioorg. Med. Chem. Lett. 26, 346–350. doi:10.1016/j.bmcl.2015.12.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, P., Li, X.-M., Wang, J.-N., Li, X., and Wang, B.-G. (2015). New Butenolide Derivatives from the marine-derived Fungus Paecilomyces variotii with DPPH Radical Scavenging Activity. Phytochemistry Lett. 11, 85–88. doi:10.1016/j.phytol.2014.11.014

CrossRef Full Text | Google Scholar

Zhang, X., Li, G., Wu, D., Yu, Y., Hu, N., Wang, H., et al. (2020). Emerging Strategies for the Activity Assay and Inhibitor Screening of Alpha-Glucosidase. Food Funct. 11, 66–82. doi:10.1039/c9fo01590f

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: Aspergillus terreus, secondary metabolites, aspulvinone, structure elucidation, α-glucosidase inhibitory effect

Citation: Wu C, Cui X, Sun L, Lu J, Li F, Song M, Zhang Y, Hao X, Tian C, Song M and Liu X (2021) Aspulvinones Suppress Postprandial Hyperglycemia as Potent α-Glucosidase Inhibitors From Aspergillus terreus ASM-1. Front. Chem. 9:736070. doi: 10.3389/fchem.2021.736070

Received: 04 July 2021; Accepted: 04 August 2021;
Published: 17 August 2021.

Edited by:

Yuanyuan Lu, China Pharmaceutical University, China

Reviewed by:

Yi Wang, Zhejiang University, China
Reshma Rani, Amity University, India
Jianhua Shao, Yangzhou University, China

Copyright © 2021 Wu, Cui, Sun, Lu, Li, Song, Zhang, Hao, Tian, Song and Liu. 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: Xiaomeng Liu, lxmxm_99@126.com; Maoping Song, mpsong@zzu.edu.cn; Congkui Tian, pkutianck@163.com

These authors have contributed equally to this work

Download