• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Furan Derivatives and Polyketides from the Fungus Irpex lacteus

    2021-04-10 01:43:08MengWangZhengHuiLiMasahikoIsakaJiKaiLiuTaoFeng
    Natural Products and Bioprospecting 2021年2期

    Meng Wang·Zheng-Hui Li·Masahiko Isaka·Ji-Kai Liu·Tao Feng

    Abstract Eight new furan derivatives, irpexins A-H (1- 8), two new polyketides, irpexins I and J (9 and 10), together with nine known compounds were isolated from the fermentation of Irpex lacteus.The structures and absolute confi gurations were elucidated on the basis of extensive spectroscopic methods and Mosher ester reaction.All compounds shows no cytotoxicity to human MCF-7 and Hela cancer cell lines at the concentration of 10 μM.

    Keywords Irpex lacteus· Furan derivatives · Polyketides

    1 Introduction

    Fungi are important sources of natural products for the discovery of drugs and agricultural chemicals [1- 5], which have attracted great attention from pharmaceutical chemists and pharmacologists.Irpex lacteusis a pathogenic wooddecaying fungus belonging to the family Polyporaceae [6].Its crude polysaccharide fraction has long been used as a traditional Chinese medicine for the treatment of chronic glomerulonephritis in clinic [7].Chemists have researched the fermentation products of this fungus and the chemical constituents of fruiting bodies.In 1981, Hayashi et al.isolated three new nematicidal metabolites including two furan derivatives from their lactose liquid medium [8].Our previous chemical investigations on the liquid fermentation of this fungus reported new carbon skeletons irlactins A—D that featured a 6/6 backbone [9], as well as a series of tremulane sesquiterpenoids [10,11].Later, studies on the fruiting bodies reported a series of triterpenes, irpeksins A—E, which inhibit the production of NO by RAW264.7 macrophages induced by LPS, and a new skeletal triterpenoid irpexolidal with a 6/5/6/5/6/5-fused polycyclic skeletal system [12,13].Recently, Duan et al.isolated four new bioactive metabolites from the endophytic fungusIrpex lacteusDR10-1, including one tremulane sesquiterpene irpexlacte A and three new furan derivatives irpexlactes B-D [14].In the current study, eight new furan derivatives, irpexins A-H (1- 8), two new polyketides, irpexins I and J (9 and 10), together with nine known compounds were obtained from cultures ofI.lacteusin rice medium.Their structures have been established by extensive spectroscopic methods, as well as modified Mosher’s method.All new compounds were evaluated for their cytotoxic activities against two human MCF-7 and Hela cell lines.

    2 Results and Discussion

    Compound 1 was isolated as a colorless oil.Its molecular formula C10H14O3was determined on the basis of positive high resolution (HR) ESI—MS atm/z183.10158 [M + H] + (calcd for C10H15O3+, 183.10157), corresponding to four degrees of unsaturation.In the 1 H NMR spectrum (Table 1), one methyl group atδH0.96 (3H, t,J= 7.5 Hz, H-10), one aldehyde proton atδH9.50 (1H, s, H-1), and two olefi nic protons atδH7.19 (1H, d,J= 3.5 Hz, H-3) and 6.28 (1H, d,J= 3.5 Hz, H-4) were readily identified.The 13 C NMR and DEPT data (Table 2) revealed ten carbon resonances ascribable for one CH3, three CH2, four CH, and two C.Of them, signals atδC151.8 (s, C-2), 124.1 (d, C-3), 108.9 (d, C-4), 163.8 (s, C-5) were suggested to establish a furan moiety conjugated to a aldehyde group atδC177.0 (d, C-1).These data were very similar to these of irpexlacte C [14].The main difference was that the position of the hydroxyl group at C-7 in irpexlacte C was moved to C-8 in 1, which was verified by the 1 H- 1 H COSY correlation ofδH1.51 (2H, m, H-9) withδH3.56 (2H, m, H-8) and of H-8 with H-10, together with the HMBC correlations from H-10 toδC30.3 (t, C-9) and 72.2 (d, C-8).Therefore, a planar structure of 1 was fi gured out as shown in Fig.1.The absolute confi guration of 1 was determined by the Mosher reaction.Compound 1 reacted withR-Mosher reagent andS-Mosher reagent to obtain compounds 1a and 1b (for their 1 H NMR spectra, see Figs.S8 and S9 in the Supporting Information), respectively.In the 1 H NMR spectrum, the chemical shifts ofδH2.76—2.71 (2H, m, H-6) andδH 2.07—2.03 (2H, m, H-7) in 1a were higher thanδH2.63—2.56 (2H, m, H-6) andδH2.02—1.96 (2H, m, H-7) in 1b, while the chemical shifts ofδH1.67—1.65 (2H, m, H-9) andδH0.85—0.83 (3H, m, Me-10) were lower thanδH1.74—1.73 (2H, m, H-9) andδH0.95—0.92 (3H, m, Me-10) in 1b.Due to the shielding effect of phenyl, this result suggested that —C9-C10— in 1a was on the same side as phenyl inR-Mosher, and —C9-C10— in 1b is not on the same side as phenyl inS-Mosher corresponding to A-type and B-type in Fig.2, respectively.At this point, the absolute confi guration of C-8 was determined asR.Compound 1 was, therefore, established as irpexin A, as depicted.

    Compound 2 was isolated as a colorless oil.Its molecular formula C10H12O3was determined on the basis of HRESI-MS, corresponding to fi ve degrees of unsaturation.The13C NMR and DEPT data of 2 were very similar to these of irpexlacte C [14], except for one double bond in 2.In1H-1H COSY spectrum, the correlation ofδH6.84 (1H, m, H-9) withδH2.21 (1H, m, H-8) andδH5.06 (2H, m, H-10) revealed a double bond was formed between C-9 and C-10 in 2.The optical rotation of 2 ([α] = - 38.0), contrary to that of 1 ([α] = + 44.8) and those analogues reported in the literature [15], indicated that the absolute confi guration of the C-6 in the side chain wasS.Compound 2 was, therefore, established as irpexin B, as depicted.

    Compound 3 was isolated as a colorless oil.Its molecular formula C10H16O3was determined on the basis of HR-ESIMS, corresponding to three degrees of unsaturation.The13C NMR and DEPT data (Table 2) of 3 were very similar to these of irpexlacte C [14], except that the aldehyde unit was replaced by a hydroxymethyl (δC57.5) in 3, as supported by the HMBC correlations fromδH4.59 (1H, s, H-1) toδC153.3 (s, C-2) and 108.4 (d, C-3).The optical rotation of 3 ([α] = - 16.1) suggested C-6 to beSform.Compound 3 was, therefore, established as irpexin C, as depicted.

    Compound 4 was isolated as a colorless oil.Its molecular formula C10H20O4was determined on the basis of HR-ESI—MS atm/z251.12544 [M + Na] + (calcd for C12H20O4Na + , 251.12538), corresponding to three degrees of unsaturation.The 13 C NMR and DEPT data of 4 were very similar to these of 1 (Table 2), except that the aldehyde unit in 1 was replaced by a dimethoxymethine group as supported by the HMBC correlations fromδH5.38 (1H, s, H-1) toδC149.0 (s, C-2), 109.1 (d, C-3) and two OCH3(δC52.8).The optical rotation of 4 ([α] = + 37.5), related to that of 1, revealed the absolute confi guration of C-8 to beR.Compound 4 was, therefore, established as irpexin D, as depicted.

    Compound 5 was isolated as a colorless oil.Its molecular formula C12H18O3was determined on the basis of the HR-ESI—MS data atm/z211.13280 [M + H] + (calcd for C12H19O3+, 211.13287), corresponding to four degrees of unsaturation.The 13 C NMR and DEPT data of 5 were very similar to these of 1 (Table 2), except for two additional carbon resonances linked to a carbonyl carbon forming a propionyl group, replacing the formyl group in 1.It was supported by the 1 H- 1 H COSY correlation ofδH1.20 (3H, t,J= 7.4 Hz, Me-1) withδH2.80 (2H, m, H-2), together with the HMBC correlations fromδH7.11 (1H, d, H-5), H-1, and H-2 toδC189.7 (s, C-3) and 151.2 (s, C-4).The optical rotation data ([α] = + 14.6) indicated C-10 to beRform.Compound 5 was, therefore, established as irpexin E, as depicted.

    Compound 6 was isolated as a colorless oil.Its molecular formula C12H18O3was determined on the basis of HR-ESIMS data.The 1 H and 13 C NMR spectra were very close to those of 5 (Tables 1 and 2), except for the hydroxyl group at C-10 (δC72.3) in 5 moved to the C-11 (δC67.7) in 6, as verified by the 1 H- 1 H COSY correlation ofδH3.84 (1H, m, H-11) withδH1.51 (2H, m, H-10) andδH1.21 (3H, d,J= 6.2 Hz, Me-12).The calculated optical rotation data ([α] = - 10.7) matched well with the experimental data([α] = - 14.3) for 6 suggested C-10 to beRform (see Supporting Information).Compound 6 was, therefore, established as irpexin I, as depicted.

    Table 1 1 H NMR (600 MHz) Data of 1 - 10 (δin ppm and Jin Hz)

    Table 2 13C NMR (150 MHz) Data of 1- 10 (δ in ppm and J in Hz)

    Fig.1 Chemical structures of compounds 1- 19

    Compound 7 was isolated as a colorless oil.Its molecular formula C12H16O4was determined on the basis of the positive high resolution (HR) ESI-MS at m/z 225.11211 [M + H] + (calcd for C12H17O4+, 225.11214), corresponding to fi ve degrees of unsaturation.The13C and DEPT spectra displayed 12 carbon resonances which were closely related to that of 5 (Table 2), except that the signals atδC31.3 (t, C-2) andδC72.3 (d, C-10) in 5 were replaced by signals ofδC69.3 (d, C-2) andδC209.0 (s, C-10) in 7, respectively.The 1 H- 1 H COSY correlation ofδH 1.47 (3H, d,J= 7.0 Hz, Me-1) withδH 4.83 (1H, q,J= 7.0 Hz, H-2) suggested that C-2 in 5 was substituted by a hydroxyl group.In addition, the1H-1H COSY correlation ofδH3.03 (2H, t,J= 7.3 Hz, H-8) withδH2.84 (2H, t,J= 7.3 Hz, H-9), andδH2.46 (2H, q,J= 7.3 Hz, H-11) withδH1.08 (3H, t,J= 7.3 Hz, Me-12), as well as the HMBC correlation from H-12 to C-10 andδC35.9 (t, C-11) suggested that C-10 in 5 was substituted by a ketone group.The optical rotation ([α] = - 14.0), comparing with those analogues reported in the literature [16], indicated C-2 to beSform.Therefore, the structure of 7 was established as irpexin G, as shown.

    Fig.2 Four types of products of compound 1 reacting with R- and S-type Mosher reagents

    Compound 8 was isolated as a colorless oil.Its molecular formula C13H18O4was determined on the basis of HR-ESIMS.The 13 C NMR and DEPT data of 8 were very similar to these of 7 (Table 2), except that the hydroxy at C-2 was replaced by a methoxy in 8, as supported by a HMBC correlation fromδH3.23 (3H, s,OCH3) toδC81.4 (d, C-2).Detailed analysis of 2D NMR data suggested that other parts of 8 were the same to those of 7.Therefore, compound 8 was established as irpexin K, as depicted.

    Compound 9 was isolated as a colorless oil.Its molecular formula C10H14O4was determined on the basis of HR-ESIMS atm/z221.07843 [M + Na] + (calcd for C10H14O4Na+, 221.07843), corresponding to four degrees of unsaturation.In the 1 H NMR spectrum (Table 1), one methyl group atδH1.06 (t,J= 7.3 Hz, Me-9), one methoxy group atδH3.20 (s,OCH3), two olefi nic protons atδH6.22 (d,J= 5.7 Hz, H-2) and 7.12 (d,J= 5.7 Hz, H-3) were identified.The 13 C and DEPT spectra (Table 2) displayed 10 carbon resonances corresponding to one CH3, one OCH3, three CH2, two CH (two olefi nic carbon atδC153.6 and 124.7), and three C (one ketone carbon atδC210.0, one ester carbonyl atδC169.7, one oxygenated sp3quaternary carbon atδC110.4) (Table 2).Considering one double bond, one ketone carbon and ester carbonyl, compound 9 should possess a single ring framework.Primary analysis of 2D NMR data indicated a fi ve-membered unsaturated lactone as built by C-1 to C-4 [17], as well as a 3-pentone moiety link to C-4 (Fig.3).In addition, one HMBC correlation from OMe toδC110.4 (s, C-4) indicated one methoxy at C-4.The optical rotation ([α] = 0) indicated that 9 should be racemic.Compound 9 was, therefore, identifed as irpexin I, as depicted.

    Fig.3 Key HMBC and 1H- 1 H COSY correlations of compounds 1,9, and 10

    Compound 10 was isolated as a colorless oil.Its molecular formula C12H20O2was determined on the HR-ESIMS atm/z197.15361 [M + H] + (calcd for C12H21O2+, 197.15361), corresponding to three degrees of unsaturation.In the 1 H NMR spectrum (Table 1), two methyl group atδH0.98 (t,J= 7.6 Hz, Me-1) and 0.97 (t,J= 7.6 Hz, Me-12), one oxygenated proton atδH3.56 (1H, m, H-10) were identified.The13C and DEPT spectra (Table 2) displayed 12 carbon resonances corresponding to two CH3, six CH2, one CH (one oxygenated methine atδC72.9), and three C (one ketone carbon atδC210.0, two olefi nic carbon atδC173.2 and 142.0) (Table 2).Considering one double bond and one ketone carbon, compound 10 should possess a single ring framework.The gross structure of 10 was deduced by comprehensive analysis of its 2D NMR spectra.In the1H-1H COSY spectrum, three fragments were revealed as shown with bold lines in Fig.3.In the HMBC spectrum, the correlation fromδH2.21 (2H, q,J= 7.6 Hz, H-2) toδC142.0 (s, C-3), 210.0 (s, C-4), and 173.2 (s, C-7), and fromδH2.60 (1H, m, H-8a) and 2.52 (1H, m, H-8b), as well as fromδH2.51 (2H, m, H-6) to C-4 established a fi ve-membered ring by C-3, C-4, C-5, C-6, and C-7.At the same time, the ethyl branch was connected to C-3, and the other branch composed by C8, C9, C10, C11, and C12 was connected to C-7.According to 1 H- 1 H COSY spectrum, a hydroxy attached to C-10 (δC72.9) was established.The optical rotation data ([α] = - 9.7) suggestedSform of C-10.Compound 10 was, therefore, identified as irpexin J, as depicted.

    Other known compounds were identified by comparison of their NMR spectroscopic data with the literature, namely irpexlacte B (11) [14], irpexlacte C (12) [14], irpexlacte D (13) [14], 1-(5-(hydroxymethyl)-2-furanyl)-pentanone (14) [18], 5-(3-oxopentyl)-2-furancarboxaldehyde (15) [19], (1R,2R)-1-(5-(pent-4-en-1-yl)furan-2-yl)propane-1,2-diol (16) [20], (1R,2R)-1-(5-pentylfuran-2-yl)propane-1,2-diol (17) [20], 2,5,8-decanetrione (21) [19], 2-butyl-3hydroxy-2-cyclopenten-1-one (19) [22].

    All new compounds were evaluated for their cytotoxicity against human MCF-7 and Hela cancer cell lines.However, no compounds showed cytotoxicity against two cell lines at the concentration of 10 μM.

    3 Experimental Section

    3.1 General Experimental Procedures

    Optical rotations were measured with a Horiba SEPA-300 polarimeter.IR spectra were obtained with a Tenor 27 spectrophotometer using KBr pellets.1D and 2D spectra were run on a Bruker Avance III 600 MHz spectrometer with TMS as an internal standard.Chemical shifts (δ) were expressed in ppm with reference to the solvent signals.Mass spectra were recorded on an Agilent 6200 Q-TOF MS system.Column chromatography (CC) was performed on silica gel (200—300 mesh, Qingdao Marine Chemical Ltd., Qingdao, People’s Republic of China), RP-18 gel (20—45 μm, Fuji Silysia Chemical Ltd., Japan), and Sephadex LH-20 (Pharmacia Fine Chemical Co., Ltd., Sweden).Medium Pressure Liquid Chromatography (MPLC) was performed on a Biotage SP1 equipment, and columns packed with RP-18 gel.Preparative High Performance Liquid Chromatography (prep-HPLC) was performed on an Agilent 1260 liquid chromatography system equipped with Zorbax SB-C18 columns (5 μm, 9.4 mm × 150 mm or 21.2 mm × 150 mm) and a DAD detector.Fractions were monitored by TLC (GF 254, Qingdao Haiyang Chemical Co., Ltd.Qingdao), and spots were visualized by heating silica gel plates sprayed with 10% H2SO4in EtOH.

    3.2 Fungal Material

    Irpex lacteuswas collected from the Wangtianshu Scenic Area, Xishuangbanna, Yunnan Province in July 2014, and authenticated by Prof.Yu-Cheng Dai of Beijing Forestry University.A voucher specimen ofI.lacteuswas deposited at the Higher Fungi Chemistry Group of School of Pharmaceutical Sciences, South-Central University for Nationalities (No.HFG 201407-SCUN201804.2).

    3.3 Fermentation Condition

    This strain was cultured on PDA medium for 8 days, and then was cut into small pieces to incubate on solid rice medium.The rice medium was inoculated in 500-mL Erlenmeyer flasks, each containing 100 g of rice medium and 100 mL of water.Flask cultures were put into a high pressure sterilizing pot and sterilized at 121 ℃ for 15 min.Then, the strain ofI.lacteuswas cultured in the rice medium, and two hundred 500-mL Erlenmeyer flasks were dark incubated fi xedly at 25 °C for 30 days.

    3.4 Extraction and Isolation

    The culture ofI.lacteusin rice medium (20 kg) was extracted fi ve times with Me2CO to give a crude extract, which was partitioned into water and EtOAc layers.The EtOAc layer (72 g) was subjected to CC over silica gel (80-100 mesh) eluted with a solvent system of CHCl3/MeOH (from 1:0 to 0:1) to obtain nine fractions (A-I).Fraction C (6 g) was separated by MPLC over RP-18 silica gel column eluted with MeOH/H2O (from 5:95 to 100:0, v/v) to give fi ve sub-fractions (C1-C5).Fraction C1was repeatly separated by CC over silica gel eluted with petroleum ether/ethyl acetate (15:1) to gain 18 (4.6 mg), and then purified by prep-HPLC (CH3CN/H2O from 5:95 to 25:75 in 25 min, v/v) to give 13 (0.8 mg,tR= 21.6 min).Fraction C2was separated by CC over Sephadex LH-20 (methanol), and then purified by prep-HPLC (CH3CN/H2O = 16:84, v/v) to give 2 (2.2 mg,tR= 20.1 min),12 (23.7 mg,tR= 18.0 min) and 19 (2.0 mg,tR= 24.1 min).Fraction C3was repeatly separated by CC over silica gel eluted with petroleum ether/ ethyl acetate (15:1), and then purified by prep-HPLC (CH3CN/H2O from 8:92 to 20:80 in 40 min, v/v) to give 7 (4.5 mg,tR= 20.8 min) and 15 (0.6 mg,tR= 17.3 min).Fraction C4was separated by CC over Sephadex LH-20 (methanol), and then purified by prep-HPLC (MeOH/H2O from 25:75 to 55:45 in 30 min, v/v) to give 1 (24.8 mg,tR= 22.6 min) and 4 (2.5 mg,tR= 17.4 min).

    Fraction D (8 g) was separated by MPLC over RP-18 silica gel column eluted with MeOH/H2O (from 5:95 to 100:0, v/v) to give eleven sub-fractions (D1-D11).Fraction D2was purified by prep-HPLC (CH3OH/H2O = 26:74, v/v) to give 3 (2.5 mg,tR= 18.2 min) and 9 (1.0 mg,tR= 14.6 min).Fraction D3was separated by CC over Sephadex LH-20 (methanol), and then purified by prep-HPLC (CH3CN/H2O = 21:79, v/v) to give 6 (2.4 mg,tR= 20.6 min) and 10 (5.0 mg,tR= 25.0 min), and by prep-HPLC (CH3CN/H2O = 26:73, v/v) to give 11 (14.1 mg,tR= 16.2 min) and 14 (2.4 mg,tR= 18.9 min).Fraction D4was repeatly separated by CC over silica gel eluted with petroleum ether/acetone (10:1), and then purified by prep-HPLC (CH3CN/H2O = 32:68, v/v) to give 7 (3.3 mg,tR= 15.4 min).Fraction D6was separated by CC over silica gel eluted with petroleum ether/acetone (10:1), and then purified by prep-HPLC (CH3CN/H2O = 32:68, v/v) to give 16 (2.2 mg,tR= 17.5 min) and 17 (12.8 mg,tR= 19.8 min).Fraction E (7.2 g) was separated by MPLC over RP-18 eluted with MeOH/H2O (from 5/95 to 100/0, v/v) to give fi ve subfractions (E1-E5).Fraction E3was subjected to CC over silica gel (200—300 mesh) and Sephadex LH-20 (methanol), and then purified by prep-HPLC (CH3CN/H2O from 20:80 to 70:30 in 20 min, v/v) to give 8 (3.2 mg,tR= 8.9 min).

    3.5 Spectroscopic Data of Compounds

    3.5.1 Irpexin A (1)

    3.5.2 Irpexin B (2)

    3.5.3 Irpexin C (3)

    3.5.4 Irpexin D (4)

    3.5.5 Irpexin E (5)

    3.5.6 Irpexin F (6)

    3.5.7 Irpexin G (7)

    3.5.8 Irpexin H (8)

    3.5.9 Irpexin I (9)

    3.5.9 .1 Irpexin J (10)Colorless oil;- 9.7 (c0.12, MeOH); UV (MeOH) λmax(log ε) 235 (3.18) nm.IR (KBr)νmax: 3338, 2943, 2831, 1647, 1452, 1114, 1031 cm-1;1H NMR (600 MHz,CDCl3) and13C NMR (150 MHz,CDCl3) spectroscopic data, see Tables 1 and 2; HR-ESI—MS:m/z197.15361 [M + H]+(calcd for C12H21O2+, 197.15361).

    3.6 Mosher Reaction

    Compound 1 (25 μmol) was added to a 10 mL round bottom flask, then 2.5 mL of dichloromethane and 50 μmol of pyridine were added.Finally, 50 μmol of (R) or (S) -MTPA-Cl was added dropwise, which was detected by the thin-layer chromatography.After 3 h of stirring at room temperature, the reaction was terminated, and the reaction solution was washed twice with 1 M diluted hydrochloric acid.The organic phase was concentrated and separated by HPLC to obtain compounds 1a and 1b respectively.

    3.7 Cytotoxicity Assay

    Two human cancer cell lines, breast cancer MCF-7 and Hela, were used in the cytotoxic assay.All the cells were cultured in RPMI-1640 or DMEM medium (Hyclone, USA), supplemented with 10% fetal bovine serum (Hyclone, USA) in 5% CO2at 37 ℃.The cytotoxicity assay was performed according to the MTT (3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyl tetrazolium bromide) method in 96-well microplates [23].Briefl y, 100 μL adherent cells were seeded into each well of 96-well cell culture plates and allowed to adhere for 12 h before drug addition, while suspended cells were seeded just before drug addition with initial density of 1 × 10 5 cells/mL.Each tumor cell line was exposed to the test compound at concentrations of 0.0625, 0.32, 1.6, 8, and 40 μM in triplicates for 48 h, with taxol as a positive control.After compound treatment, cell viability was detected and cell growth curve was graphed.IC50value was calculated by Reed and Muench’s method [24].

    AcknowledgementsThis work was financially supported by the joint research project from the National Natural Science Foundation of China (Grant No.21961142008) and Thailand Research Fund (Grant No.DBG6280008).The authors thank Analytical & Measuring Centre, South-Central University for Nationalities for the spectra measurements.

    Compliance with Ethical Standards

    Conflict of interestThe authors declare no confl ict of interest.

    Open AccessThis article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material.If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.To view a copy of this licence, visit http://creat iveco mmons.org/licen ses/by/4.0/.

    日韩国内少妇激情av| 桃色一区二区三区在线观看| 卡戴珊不雅视频在线播放| 99九九线精品视频在线观看视频| 五月伊人婷婷丁香| 精品一区二区三区人妻视频| 久久久久久久久中文| 国产日本99.免费观看| 99在线视频只有这里精品首页| 校园人妻丝袜中文字幕| 久久综合国产亚洲精品| 亚洲国产欧美人成| 色哟哟·www| 国产男人的电影天堂91| 国产精品人妻久久久影院| 日本与韩国留学比较| 国产高潮美女av| 国产精品亚洲美女久久久| 亚洲中文日韩欧美视频| 欧美最黄视频在线播放免费| 亚洲欧美日韩高清在线视频| 欧美中文日本在线观看视频| 桃色一区二区三区在线观看| 在线观看午夜福利视频| 国产精品爽爽va在线观看网站| 一个人免费在线观看电影| 久久天躁狠狠躁夜夜2o2o| 色综合站精品国产| 国产一区二区在线av高清观看| 色综合亚洲欧美另类图片| 日本与韩国留学比较| 亚洲精品日韩在线中文字幕 | 久久99热6这里只有精品| 国产淫片久久久久久久久| 欧美日韩乱码在线| 国产精品爽爽va在线观看网站| 丰满的人妻完整版| 我的女老师完整版在线观看| 人人妻人人澡欧美一区二区| 97碰自拍视频| 国内精品美女久久久久久| 亚洲一区二区三区色噜噜| 亚洲国产精品合色在线| 久久久久国产网址| 国产av麻豆久久久久久久| 91在线观看av| 啦啦啦观看免费观看视频高清| 黄色欧美视频在线观看| 一卡2卡三卡四卡精品乱码亚洲| 国产亚洲欧美98| 禁无遮挡网站| 男女视频在线观看网站免费| 免费大片18禁| 露出奶头的视频| 麻豆国产97在线/欧美| 女生性感内裤真人,穿戴方法视频| 九九久久精品国产亚洲av麻豆| 欧美bdsm另类| 神马国产精品三级电影在线观看| 国产爱豆传媒在线观看| 国产老妇女一区| 网址你懂的国产日韩在线| 亚洲性久久影院| 国产亚洲精品综合一区在线观看| 国产aⅴ精品一区二区三区波| 中文在线观看免费www的网站| 免费看日本二区| 久久久国产成人精品二区| 国产精品一区www在线观看| 人妻少妇偷人精品九色| 欧美激情在线99| 69人妻影院| 日韩成人av中文字幕在线观看 | 99热这里只有精品一区| 97超级碰碰碰精品色视频在线观看| 亚洲国产精品久久男人天堂| 在线国产一区二区在线| 级片在线观看| 精品久久久久久久久久久久久| 欧美成人免费av一区二区三区| 成人亚洲精品av一区二区| 亚洲成a人片在线一区二区| 色av中文字幕| 51国产日韩欧美| 亚洲av不卡在线观看| 国产视频内射| 夜夜夜夜夜久久久久| 亚洲人与动物交配视频| 三级毛片av免费| 真实男女啪啪啪动态图| 免费看日本二区| 麻豆av噜噜一区二区三区| 国产精品久久久久久久久免| 精品无人区乱码1区二区| 女人被狂操c到高潮| 欧美成人免费av一区二区三区| 亚洲久久久久久中文字幕| 又爽又黄a免费视频| 国产在视频线在精品| 日本三级黄在线观看| 内地一区二区视频在线| 高清毛片免费看| 国产在线精品亚洲第一网站| 色尼玛亚洲综合影院| 夜夜夜夜夜久久久久| 淫妇啪啪啪对白视频| 国产精品99久久久久久久久| 免费一级毛片在线播放高清视频| 搡老熟女国产l中国老女人| 亚洲精品一卡2卡三卡4卡5卡| 久久精品91蜜桃| av天堂在线播放| 久久99热这里只有精品18| 国产精品一区二区三区四区久久| 亚洲最大成人av| or卡值多少钱| 欧美在线一区亚洲| 日韩亚洲欧美综合| 成人鲁丝片一二三区免费| 午夜视频国产福利| 久久久久国产网址| 国产乱人视频| av在线观看视频网站免费| 国产蜜桃级精品一区二区三区| 变态另类丝袜制服| 国产黄a三级三级三级人| 亚洲欧美日韩卡通动漫| 国产精品久久视频播放| 精品欧美国产一区二区三| 别揉我奶头~嗯~啊~动态视频| 极品教师在线视频| 国产大屁股一区二区在线视频| 国产伦精品一区二区三区四那| 欧美绝顶高潮抽搐喷水| 你懂的网址亚洲精品在线观看 | 丝袜美腿在线中文| 亚洲成a人片在线一区二区| 中文字幕久久专区| 亚洲成人av在线免费| 男女下面进入的视频免费午夜| 免费不卡的大黄色大毛片视频在线观看 | 欧美中文日本在线观看视频| 免费观看精品视频网站| 日韩一区二区视频免费看| 身体一侧抽搐| 99久久中文字幕三级久久日本| 亚洲内射少妇av| 亚洲国产精品成人综合色| 别揉我奶头~嗯~啊~动态视频| 欧美3d第一页| 联通29元200g的流量卡| 久久精品国产亚洲av涩爱 | 国产精品国产高清国产av| 国产精品久久久久久久电影| 男人舔奶头视频| 最新中文字幕久久久久| 国产精品免费一区二区三区在线| 激情 狠狠 欧美| 亚洲精品久久国产高清桃花| 国产日本99.免费观看| 真实男女啪啪啪动态图| 亚洲人成网站在线观看播放| 色综合色国产| 人妻夜夜爽99麻豆av| 成人三级黄色视频| 成年女人永久免费观看视频| 别揉我奶头~嗯~啊~动态视频| 一本精品99久久精品77| 日日摸夜夜添夜夜添av毛片| 插逼视频在线观看| 中文字幕av成人在线电影| 高清日韩中文字幕在线| 国产精品日韩av在线免费观看| 午夜亚洲福利在线播放| 在现免费观看毛片| 亚洲不卡免费看| 国产精品不卡视频一区二区| 在线观看免费视频日本深夜| 波多野结衣高清无吗| 香蕉av资源在线| 少妇裸体淫交视频免费看高清| 一级毛片我不卡| 国产v大片淫在线免费观看| 精品久久久久久久人妻蜜臀av| 三级经典国产精品| 在线观看av片永久免费下载| 3wmmmm亚洲av在线观看| 亚洲18禁久久av| 成人特级av手机在线观看| 国产亚洲精品久久久com| 欧美丝袜亚洲另类| 内地一区二区视频在线| 欧美高清成人免费视频www| 亚洲中文字幕日韩| 卡戴珊不雅视频在线播放| 韩国av在线不卡| 久久人妻av系列| 国产精品一区二区三区四区免费观看 | 在线观看66精品国产| 久久欧美精品欧美久久欧美| 日韩中字成人| 亚洲精品色激情综合| 女人十人毛片免费观看3o分钟| 99久久精品国产国产毛片| 啦啦啦啦在线视频资源| 久久久久久久亚洲中文字幕| 看免费成人av毛片| 久久欧美精品欧美久久欧美| 亚洲av免费在线观看| 午夜亚洲福利在线播放| 国产中年淑女户外野战色| 亚洲人成网站高清观看| 国产伦一二天堂av在线观看| .国产精品久久| 免费电影在线观看免费观看| 午夜老司机福利剧场| 在线免费观看不下载黄p国产| 人妻丰满熟妇av一区二区三区| 精品少妇黑人巨大在线播放 | 久久6这里有精品| 亚洲国产精品成人久久小说 | 久久精品国产亚洲av天美| 欧美+日韩+精品| 中文字幕av成人在线电影| 国产成人91sexporn| 亚洲av熟女| 国产麻豆成人av免费视频| 啦啦啦观看免费观看视频高清| 一级黄色大片毛片| 此物有八面人人有两片| 国产精品国产三级国产av玫瑰| 一级毛片久久久久久久久女| 好男人在线观看高清免费视频| 免费黄网站久久成人精品| 国产在线男女| 午夜福利高清视频| 色吧在线观看| 男女视频在线观看网站免费| 国产精品美女特级片免费视频播放器| 蜜桃久久精品国产亚洲av| 日韩在线高清观看一区二区三区| 国产极品精品免费视频能看的| 午夜a级毛片| 我要搜黄色片| av女优亚洲男人天堂| 成年av动漫网址| 国内精品一区二区在线观看| 日本免费一区二区三区高清不卡| 久久精品夜夜夜夜夜久久蜜豆| 欧美性猛交黑人性爽| 99久国产av精品国产电影| 免费黄网站久久成人精品| 日韩一区二区视频免费看| 亚洲精品一卡2卡三卡4卡5卡| 人人妻,人人澡人人爽秒播| 国产一区亚洲一区在线观看| 日日摸夜夜添夜夜爱| 非洲黑人性xxxx精品又粗又长| 精品久久久久久久久久免费视频| 欧美最新免费一区二区三区| 亚洲国产欧洲综合997久久,| 老司机影院成人| 内地一区二区视频在线| 内射极品少妇av片p| 中文亚洲av片在线观看爽| 久久久久久伊人网av| 日本免费a在线| 99热这里只有是精品在线观看| 99热这里只有是精品50| 小说图片视频综合网站| 91麻豆精品激情在线观看国产| 在线观看美女被高潮喷水网站| 非洲黑人性xxxx精品又粗又长| 天堂影院成人在线观看| 免费av不卡在线播放| 婷婷色综合大香蕉| а√天堂www在线а√下载| 成年女人毛片免费观看观看9| 国产成人一区二区在线| 寂寞人妻少妇视频99o| 色吧在线观看| 成人无遮挡网站| 亚洲av不卡在线观看| 变态另类成人亚洲欧美熟女| 亚洲国产精品sss在线观看| 超碰av人人做人人爽久久| 非洲黑人性xxxx精品又粗又长| 国产色婷婷99| 日本-黄色视频高清免费观看| 色吧在线观看| 国产成人a∨麻豆精品| 男人狂女人下面高潮的视频| 免费av观看视频| 中文字幕av成人在线电影| 国产精品三级大全| 亚洲欧美日韩高清专用| 乱码一卡2卡4卡精品| 国产精品国产高清国产av| 国产午夜精品久久久久久一区二区三区 | 97人妻精品一区二区三区麻豆| 变态另类丝袜制服| 18禁裸乳无遮挡免费网站照片| 亚洲人成网站在线播| 中文资源天堂在线| 欧美+亚洲+日韩+国产| 久久人人爽人人爽人人片va| 久久久久国内视频| 蜜桃久久精品国产亚洲av| 成人高潮视频无遮挡免费网站| 日本 av在线| 国模一区二区三区四区视频| 国产成人a区在线观看| 亚洲专区国产一区二区| 亚洲一级一片aⅴ在线观看| 国产精品一区www在线观看| 美女被艹到高潮喷水动态| 色在线成人网| av在线天堂中文字幕| 麻豆乱淫一区二区| 欧美日韩国产亚洲二区| 亚洲激情五月婷婷啪啪| 一个人免费在线观看电影| 成人二区视频| 国语自产精品视频在线第100页| av免费在线看不卡| 国产淫片久久久久久久久| 成人精品一区二区免费| 国产午夜精品久久久久久一区二区三区 | 十八禁网站免费在线| 国产69精品久久久久777片| 女同久久另类99精品国产91| 亚洲内射少妇av| 18禁裸乳无遮挡免费网站照片| 中文在线观看免费www的网站| 国产午夜福利久久久久久| 少妇的逼好多水| 国产日本99.免费观看| 国产高清不卡午夜福利| 有码 亚洲区| 俺也久久电影网| av在线老鸭窝| 不卡一级毛片| 精品久久久久久久久av| 国产黄色视频一区二区在线观看 | 国产精品一及| 男插女下体视频免费在线播放| 国产综合懂色| 亚洲电影在线观看av| 国产成年人精品一区二区| 久久国内精品自在自线图片| 国产在线精品亚洲第一网站| 永久网站在线| 99在线人妻在线中文字幕| 国产精品嫩草影院av在线观看| 久久久久久久亚洲中文字幕| 亚洲电影在线观看av| 久久人人爽人人爽人人片va| 欧美色视频一区免费| 亚洲丝袜综合中文字幕| 嫩草影视91久久| 俄罗斯特黄特色一大片| 女的被弄到高潮叫床怎么办| 国产精品国产高清国产av| 丰满的人妻完整版| 最近在线观看免费完整版| 国产亚洲精品av在线| avwww免费| 久久国内精品自在自线图片| 午夜视频国产福利| 搡老岳熟女国产| 国产精品久久视频播放| 99热这里只有是精品50| 久久热精品热| 综合色丁香网| 国产精品永久免费网站| 高清午夜精品一区二区三区 | 九九久久精品国产亚洲av麻豆| 亚洲欧美清纯卡通| 97超碰精品成人国产| 黄片wwwwww| 午夜免费激情av| 久久久久免费精品人妻一区二区| 少妇被粗大猛烈的视频| 久久久久国内视频| 欧美潮喷喷水| 久久国产乱子免费精品| 插阴视频在线观看视频| 成人欧美大片| 国产午夜精品久久久久久一区二区三区 | 亚洲av成人av| 久久精品国产亚洲网站| 伊人久久精品亚洲午夜| av免费在线看不卡| 99热6这里只有精品| 欧美最新免费一区二区三区| 中文字幕熟女人妻在线| 黄色欧美视频在线观看| 久久这里只有精品中国| 99热全是精品| 日本黄大片高清| 嫩草影院新地址| av在线老鸭窝| 欧美日韩在线观看h| 高清毛片免费观看视频网站| 国产午夜精品论理片| 国产精品免费一区二区三区在线| 国产视频内射| 国产视频一区二区在线看| 中文字幕熟女人妻在线| 直男gayav资源| 国产三级中文精品| 亚洲高清免费不卡视频| 国产精品不卡视频一区二区| 日日撸夜夜添| 久久99热6这里只有精品| 国产视频一区二区在线看| 听说在线观看完整版免费高清| av黄色大香蕉| 日产精品乱码卡一卡2卡三| 亚洲精品久久国产高清桃花| 国产伦在线观看视频一区| 最新中文字幕久久久久| 国产中年淑女户外野战色| 在线观看免费视频日本深夜| 舔av片在线| 又黄又爽又刺激的免费视频.| 最好的美女福利视频网| 国产亚洲精品久久久com| 国产探花极品一区二区| 久久精品影院6| 色av中文字幕| 又粗又爽又猛毛片免费看| 性色avwww在线观看| 久久精品国产亚洲av涩爱 | 一级a爱片免费观看的视频| 在线观看午夜福利视频| 赤兔流量卡办理| 午夜a级毛片| 美女 人体艺术 gogo| .国产精品久久| av在线天堂中文字幕| 老女人水多毛片| 日韩欧美精品v在线| 国产片特级美女逼逼视频| 久久久久国产网址| 国产精品日韩av在线免费观看| 中文在线观看免费www的网站| 无遮挡黄片免费观看| 好男人在线观看高清免费视频| 国语自产精品视频在线第100页| 国产黄a三级三级三级人| 亚洲激情五月婷婷啪啪| 国产精品一区二区免费欧美| 日本爱情动作片www.在线观看 | 午夜影院日韩av| 久久久色成人| 三级经典国产精品| 亚洲自偷自拍三级| 特大巨黑吊av在线直播| 男女啪啪激烈高潮av片| 最近手机中文字幕大全| 久久精品国产亚洲av涩爱 | 国产一区二区亚洲精品在线观看| 国产在视频线在精品| 成人三级黄色视频| 久99久视频精品免费| 午夜福利在线观看免费完整高清在 | 蜜臀久久99精品久久宅男| 日韩亚洲欧美综合| 人人妻人人澡人人爽人人夜夜 | 国产蜜桃级精品一区二区三区| 欧美中文日本在线观看视频| 嫩草影视91久久| 禁无遮挡网站| 国产精品国产高清国产av| a级毛片免费高清观看在线播放| 久久久久久久亚洲中文字幕| 国产不卡一卡二| 99热这里只有是精品在线观看| 免费在线观看影片大全网站| 卡戴珊不雅视频在线播放| 午夜日韩欧美国产| 亚洲熟妇熟女久久| 久久久a久久爽久久v久久| 精品国产三级普通话版| 久久草成人影院| 欧美日本亚洲视频在线播放| 麻豆成人午夜福利视频| 欧美最新免费一区二区三区| 欧美在线一区亚洲| 国产精品乱码一区二三区的特点| 深爱激情五月婷婷| 人人妻人人看人人澡| АⅤ资源中文在线天堂| 国产成人a区在线观看| 久久精品久久久久久噜噜老黄 | 国产精品av视频在线免费观看| 欧美色欧美亚洲另类二区| 精品人妻偷拍中文字幕| 国产亚洲av嫩草精品影院| 精品福利观看| 免费人成视频x8x8入口观看| 久久久久免费精品人妻一区二区| 久久久久久国产a免费观看| 久久久a久久爽久久v久久| 国产伦精品一区二区三区视频9| 91精品国产九色| 自拍偷自拍亚洲精品老妇| 亚洲最大成人手机在线| 亚洲av五月六月丁香网| a级毛片免费高清观看在线播放| 国产高清三级在线| 黄片wwwwww| 麻豆av噜噜一区二区三区| 波多野结衣高清无吗| 乱系列少妇在线播放| 日韩欧美 国产精品| 亚洲av免费高清在线观看| 精品不卡国产一区二区三区| 69av精品久久久久久| 欧美绝顶高潮抽搐喷水| 欧美性猛交╳xxx乱大交人| 午夜视频国产福利| 天堂√8在线中文| 老女人水多毛片| 久久这里只有精品中国| 91精品国产九色| 无遮挡黄片免费观看| 三级毛片av免费| 国产高潮美女av| 亚洲av免费高清在线观看| 国产熟女欧美一区二区| 国产精品久久视频播放| 色综合色国产| 女人十人毛片免费观看3o分钟| 久久精品国产亚洲av涩爱 | 久久久久久国产a免费观看| 欧美最黄视频在线播放免费| 日韩欧美三级三区| 十八禁国产超污无遮挡网站| 菩萨蛮人人尽说江南好唐韦庄 | 国产极品精品免费视频能看的| 国产成人影院久久av| 欧美日韩综合久久久久久| 欧美成人a在线观看| 免费av不卡在线播放| 日韩av在线大香蕉| 色综合色国产| 国产精品久久久久久av不卡| 国产精品1区2区在线观看.| 男女视频在线观看网站免费| 日本熟妇午夜| 午夜亚洲福利在线播放| 国产免费男女视频| 天美传媒精品一区二区| 欧美人与善性xxx| 久久久精品欧美日韩精品| 国产伦精品一区二区三区四那| 床上黄色一级片| 亚洲三级黄色毛片| 日本精品一区二区三区蜜桃| 天堂影院成人在线观看| 我的女老师完整版在线观看| av黄色大香蕉| 日韩欧美三级三区| 成人美女网站在线观看视频| 丰满人妻一区二区三区视频av| 亚洲不卡免费看| av国产免费在线观看| 欧美中文日本在线观看视频| 蜜臀久久99精品久久宅男| 秋霞在线观看毛片| 免费看a级黄色片| 久久99热6这里只有精品| 亚洲精品日韩av片在线观看| 久久久国产成人精品二区| 97人妻精品一区二区三区麻豆| 国产单亲对白刺激| 18禁黄网站禁片免费观看直播| 国产精品久久久久久久久免| 看非洲黑人一级黄片| 蜜桃久久精品国产亚洲av| 九九久久精品国产亚洲av麻豆| 午夜久久久久精精品| 日韩三级伦理在线观看| 亚洲在线观看片| 国产一区二区在线av高清观看| av专区在线播放| 日韩强制内射视频| 国产午夜精品论理片| 最近视频中文字幕2019在线8| 人人妻人人看人人澡| 综合色丁香网| 欧美+日韩+精品| 免费观看人在逋| 成人高潮视频无遮挡免费网站| 狠狠狠狠99中文字幕| 国产精品国产高清国产av| 午夜福利在线观看吧| 变态另类成人亚洲欧美熟女| 最后的刺客免费高清国语| 亚洲成人av在线免费| 99riav亚洲国产免费| 国产精品久久久久久久电影| 美女cb高潮喷水在线观看| 欧美三级亚洲精品| 男人和女人高潮做爰伦理| 日本-黄色视频高清免费观看| 亚洲欧美日韩高清专用| ponron亚洲| 狂野欧美白嫩少妇大欣赏| 亚洲av中文av极速乱| 又爽又黄无遮挡网站| 日韩精品有码人妻一区| 久久久a久久爽久久v久久|