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

    Simultaneous determination of 35 constituents and elucidation of effective constituents in a multi-herb Chinese medicine formula Xiaoer-Feire-Kechuan

    2022-01-19 06:39:50ZhnpengShngLuluXuYiKungYnLinShungLiuLongSunToBoMinYeXueQio
    Journal of Pharmaceutical Analysis 2021年6期

    Zhnpeng Shng,Lulu Xu,Yi Kung,Yn Lin,Shung Liu,Long Sun,To Bo,Min Ye,**,Xue Qio,*

    a State Key Laboratory of Natural and Biomimetic Drugs,School of Pharmaceutical Sciences,Peking University,Beijing 100191,China

    bSchool of Pharmaceutical Sciences,Guizhou Medical University,Guizhou,550025,China

    c Thermo Fisher Scientific,Beijing,100102,China

    Keywords:

    Multi-herb formulas

    Effective constituents

    Xiaoer-Feire-Kechuan formula(XFK)Parallel reaction monitoring(PRM)Cyclooxygenase-2(COX-2)

    Peer review under responsibility of Xi’an Jiaotong University.

    A B S T R A C T

    Xiaoer-Feire-Kechuan(XFK)is an 11-herb Chinese medicine formula to treat cough and pulmonary inflammation.The complicated composition rendered its chemical analysis and effective-component elucidation.In this study,we combined quantitative analysis and bioactivity test to reveal the antiinflammatory constituents of XFK.First,UPLC-DAD and UHPLC/Q-Orbitrap-MS methods were established and validated to quantify 35 analytes(covering 9 out of 11 herbs)in different XFK formulations.Parallel reaction monitoring mode built in Q-Orbitrap-MS was used to improve the sensitivity and selectivity.Then,anti-inflammatory activities of the 35 analytes were analyzed using in vitro COX-2 inhibition assay.Finally,major analytes forsythosides H,I,A(8-10),and baicalin(15)(total contents varied from 21.79 to 91.20 mg/dose in different formulations)with significant activities(inhibitory rate≥80%)were proposed as the anti-inflammatory constituents of XFK.The present study provided an effective strategy to discover effective constituents of multi-herb formulas.

    1.Introduction

    Traditional Chinese medicines(TCMs)have attracted increasing attention due to their curative effects on complex and chronic diseases[1,2].In clinic,TCMs are mainly used in the form of formulas.The component herbs in a formula may target different symptoms of complex diseases,which might play synergistic roles in the drug efficacy[3,4].In total,1,933 TCM formulas are recorded in the Chinese Pharmacopoeia(2015 edition),and many of them are composed of more than 10 herbs[5],which are called “multi-herb formulas”or“Dafufang”[6].Up to now,only a few of reports have focused on the quality control of multi-herb formulas[6-9].For instance,12(6 out of 12 component herbs involved),16(5 out of 10 herbs involved),and 41(15 out of 19 herbs involved)analytes were determined in Kangjing formula,Yougui pill,and Niuhuangshangqing pill,respectively[6,7,9].However,due to their complex components,it is still challenging to comprehensively analyze the constituents and to discover the effective ones for multi-herb formulas.

    Liquid chromatography with diode array detector(LC-DAD)and liquid chromatography tandem mass spectrometry(LC/MS)are popularly used methods to determine the contents of constituents in complex samples[9-14].Especially,parallel reaction monitoring(PRM)scan mode built in quadrupole(Q)-Orbitrap-MS combines the mass isolation capability of the quadrupole and the high resolution of Orbitrap detector,which could avoid false compliant and non-compliant results in complex samples[15-17].It emerges as a promising method to monitor characteristic analytes in multi-herb formulas.To further elucidate the effective constituents of multiherb formulas,feasible pharmacological models that are related to the therapeutic effect should be developed[18].For example,cyclooxygenase-2(COX-2)is one of the main isozymes responsible for inflammation[19].It has been proved to play roles in lung inflammation[20],such as pneumonia,bronchitis,and asthma[21,22].Therefore,COX-2 inhibitory assay can be used to reveal anti-inflammatory constituents from formulas.

    Xiaoer-Feire-Kechuan(XFK)is a multi-herb formula composed of 10 herbal and 1 mineral herbs(Table 1).It is a patent TCM formula to treat bronchitis,pneumonia,and cough in children[5].Both oral solution and granule formulations have been developed.XFK oral solution is recorded in the Chinese Pharmacopoeia(2015 edition),while only ephedrine and pseudoephedrine from Mahuang(MH)are used as quality control markers[5].Several studies determined the contents of chemical constituents in XFK[23,24].For example,9 analytes from MH,Huangqin(HQ),Jinyinhua(JYH),Gancao(GC),and Lianqiao(LQ)were quantified using an 83-min HPLC-DAD method[24].To fully evaluate the quality of XFK formula,it is important to monitor characteristic analytes of each herb and to identify its major effective constituents.

    Table 1The herbs of Xiaoer-Feire-Kechuan(XFK)formula.

    In the present work,we reported an integrated method to reveal the anti-inflammatory constituents of XFK formula.The contents of 35 characteristic analytes in 18 batches of XFK formulations were determined using ultra performance liquid chromatography(UPLC)-DAD and ultra-high performance liquid chromatography(UHPLC)/Q-Orbitrap-MS.Furthermore,the COX-2 inhibitory activities of the 35 analytes,separate herbs,and different XFK formulations were investigated to discover the main anti-inflammatory constituents of XFK.

    2.Experimental

    2.1.Chemicals and reagents

    Methanol,acetonitrile(Fisher Scientific,Fair Lawn,NJ,USA),and formic acid(Sigma-Aldrich,St.Louis,MO,USA)were of LC/MS grade.De-ionized water was prepared using a Milli-Q water purification system(Millipore,Burlington,MA,USA).COX-2 inhibitor screening kit was purchased from Beyotime Biotechnology(Shanghai,China).Reference standards 8-10,13,16,23,and 26 were isolated from LQ;15,17,18,and 34 were from HQ[11];20,21,and 22 were purchased from the National Institute for the Control of Biological and Pharmaceutical Products of China(Beijing,China);1,2,4,5,11,12,14,19,24,27,28,31-33,and 35 were purchased from Chengdu DeSiTe Biological Technology Co.,Ltd.(Chengdu,China);3,6,7,25,29,30,and internal standards(IS1 and IS2)were purchased from Chengdu Must Bio-technology Co.,Ltd.(Chengdu,China).Their structures are shown in Fig.1.Purity of all these standards was above 98% by HPLC analysis.

    Fig.1.Chemical structures of analytes 1-35 and the internal standards(IS1 and IS2).

    Herbs including Mahuang(MH),HQ,JYH,LQ,Kuxingren(KXR),GC,Zhimu(ZM),Banlangen(BLG),Maidong(MD),Yuxingcao(YXC),XFK oral solutions XFK1-XFK10(OS[a],10 mL/dose),and granules XFK13-XFK16(GR[a],3 g/dose)were supplied by company a.Granules XFK11 and XFK12(GR[b],3 g/dose)were from company b and granules XFK17 and XFK18(GR[c],4 g/dose)from company c.

    2.2.UPLC-DAD method for major components

    2.2.1.Preparation of calibration standard and sample solutions

    Reference standards(1-14,16,17,and 19)were dissolved in 50% methanol to prepare a mixed standard solution 1.Reference standards 15 and 18 were dissolved in 50% methanol to prepare a mixed standard solution 2.Their concentrations ranged from 29.5-217.0μg/mL of each analyte.The mixed reference solutions were respectively diluted by 2,4,8,16,32,64,128,256,and 512-fold using 50% methanol to prepare a series of calibration samples.XFK oral solution(0.2 mL)was accurately diluted by 25-fold with 50% methanol.The fine powder of XFK granules(200 mg)were extracted using 10 mL of 50% methanol in an ultrasonic bath for 5 min.All samples were filtered through 0.22μm membranes before use.

    2.2.2.Chromatographic conditions

    A Waters UPLC H-Class system(Waters Technologies,Corp.,Milford,MA,USA)was employed.Samples were separated on an Acquity HSS T3 column(2.1 mm × 100 mm,1.8μm,Waters Corporation,Milford,MA,USA)and eluted using mobile phase A(water containing 0.1% formic acid)and B(acetonitrile).The gradient program was as follows:0 min,5% B;2 min,5% B;5 min,10% B;5.5 min,12% B;13 min,18% B;20 min,40% B;25 min,100% B.The flow rate was 400μL/min and the column temperature was set at 55°C.An aliquot of 2 μL was injected for analysis.The detector wavelengths were selected according to the UV absorption of each analyte(Table 2).

    2.3.UHPLC/Q-Orbitrap-MS method for minor components

    2.3.1.Preparation of calibration standard,IS,and sample solutions

    An appropriate amount of each reference standard(20-35)was dissolved in 50% methanol to prepare a mixed standard solution containing 6-30μg/mL of each analyte.The mixed standard solution was then serially diluted by 2,4,8,16,32,64,128,256,512,and 1024-fold using 50% methanol.The series of calibration solutions were then diluted by 2-fold using the mixed internal standard solution(containing 400 ng/mL of IS1 and 400 ng/mL of IS2),respectively.The XFK oral solution(0.5 mL)was accurately diluted by 100-fold with 50% methanol.The fine powder of XFK granule(40 mg)was extracted using 10 mL of 50% methanol in an ultrasonic bath for 5 min.The sample solutions were then diluted with the mixed internal standard solution by 2-fold,respectively.All samples were filtered through 0.22μm membranes before use.

    2.3.2.Chromatographic and mass spectrometry conditions

    A Thermo Vanquish UHPLC system(Thermo Fisher Scientific,San Jose,CA,USA)was employed.Samples were separated on a Phenyl-Hexyl(2.1 mm × 100 mm,1.8μm,Agilent Technologies,Waldbronn,Germany)and eluted using mobile phase A(water containing 0.2% formic acid)and B(acetonitrile).The gradient program was as follows:0 min,2% B;2 min,2% B;7 min,13% B;13 min,28% B;16 min,80% B.An aliquot of 2μL was injected for analysis.The flow rate was 400μL/min and column temperature was set at 50°C.

    Mass spectral data acquisition was performed on a Q-Exactive Focus hybrid Q-Orbitrap mass spectrometer equipped with a heated electrospray ionization source(HESI)(Thermo Scientific,San Jose,USA).The parameters were set as follows:spray voltage,-3.5 kV;sheath gas,45 arb;auxiliary gas,10 arb;capillary temperature,350°C;auxiliary temperature,300°C;S-lens RF level,55 V.Polarity switch negative(-)/positive(+)and PRM mode were employed.The scan windows of each analyte were set from 1.5 to 2.0 min based on their retention times to ensure sufficient data points. MS/MS resolution was set at 17,500 FWHM.Quantitative product ion and collision energy of each analyte were optimized using the MS Tune software(Thermo Scientific,Wilmington,DE,USA)and are provided in Table 2 and Fig.S1.Data were processed using XcaliburTM4.1 software(Thermo Scientific,Wilmington,DE,USA).

    Table 2The information for UV and PRM parameters for quantitation of 35 analytes.

    2.4.In vitro COX-2 inhibition assay

    The inhibitory activities of different XFK formulations(10μg/mL),10 herbs(MH,HQ,GC,JYH,LQ,ZM,KXR,BLG,MD,and YXC,10 μg/mL),and 35 quantified analytes(10μM)were tested using a COX-2 inhibitor screening kit(Beyotime Biotechnology,Shanghai,China)according to the manufacturer's instructions.Celecoxib(50 nM)was used as the positive control.All experiments were carried out in triplicate.Preparation of XFK extracts is shown in Supplementary data.

    3.Results and discussion

    3.1.Selection of analytes

    Only two analytes from MH were determined according to theChinese Pharmacopoeia (2015 edition) [5]. To comprehensively evaluate the quality of XFK, characteristic constituents from each herbwere selected for analysis. As a result, 35 analytes from 9 herbs in XFK formula were determined. (R,S)-goitrin, a characteristic component inBLG, was excluded due to its low concentration in XFK. In the presentmethod, 19 major analytes were quantified using UPLC-DAD(Figs. 2A-C), and 16 minor ones were determined by UHPLC/QOrbitrap-MS (Fig. 2D).

    Fig.2.Typical UPLC-UV chromatograms and UHPLC-PRM-MS/MS ion chromatograms.(A)UPLC-UV chromatograms of mixed reference standards(1-19).(B)UPLC-UV chromatograms of oral solution XFK-8.(C)UPLC-UV chromatograms of granule XFK-13.(D)PRM ion chromatograms of XFK-13,showing minor analytes 20-35 and internal standards IS1 and IS2.

    3.2.Optimization of separation method for UHPLC/Q-Orbitrap-MS analysis

    Analytes 20-35 were determined using UHPLC/Q-Orbitrap-MS.They contain hydrophilic,hydrophobic,basic and acidic compounds,which are challenging for chromatographic separation.For example,the three ephedra alkaloids(20,21,and 22)are difficult to be fully separated using LC/MS-compatible mobile phases[5].Thus,the separation method was optimized.Different stationary phases were tested,including Acquity CSH C18(2.1 mm ×100 mm,1.7μm,Waters Corporation,Milford,MA,USA),Acquity HSS T3 C18(2.1 mm ×100 mm,1.8μm,Waters Corporation,Milford,MA,USA),Acquity Cortecs C18(2.1 mm×100 mm,1.6μm,Waters Corporation,Milford,MA,USA),Xterra MS C18(2.1 mm ×150 mm,3.5μm,Waters),and Phenyl-Hexyl(2.1 mm × 100 mm,1.8μm,Agilent).All columns exhibited good peak shape for analytes 24,26,27 and 34,but only the Phenyl-Hexyl column could effectively separate analytes 20,21,and 22(Fig.S2).By increasing the concentration of formic acid in the mobile phase(0.2%),the three analytes(20,21,and 22)showed satisfactory peak shape and resolution(Fig.S2F).The optimized stationary and mobile phases were used for the follow-up experiments.

    3.3.Optimization of MS conditions for UHPLC/Q-Orbitrap-MS analysis

    To improve the sensitivity,three different scan modes were compared,including full scan/data-dependent MS2(FS/ddMS2),target-selected ion monitoring(t-SIM)/ddMS2,and PRM(Fig.3).In FS/ddMS2,the precursor ions were detected by full scan mode and then delivered to the high energy collision-induced dissociation(HCD)cell via C-trap.The top N abundant ions in each scan were fragmented successively,and the fragment ions were detected by Orbitrap-MS to confirm the analyte.In t-SIM/ddMS2,precursor ions were selectively detected and delivered to the HCD cell,and their product ions were monitored by Orbitrap-MS for analyte confirmation.In PRM mode,predefined precursor ions were selected by the quadrupole and delivered directly to the HCD cell without detection.The product ions were then detected by Orbitrap-MS.

    Data point(DP),sensitivity,and selectivity of the three scan modes were compared using representative analytes(Fig.3).The PRM mode allowed shorter scan time and duty cycle.For alkaloids 20 and 21,DP and signal-to-noise(S/N)values increased significantly when using PRM mode.The duty cycles(calculated as shown in Fig.S3)were 1.80,0.60,and 0.12 s for FS/ddMS2,t-SIM/ddMS2and PRM,respectively.Similarly,the S/N value for the glycoside 24 and the phenylethanoid 26 in PRM was much higher than that in t-SIM/ddMS2and FS/ddMS2.Our results also indicated the highest specificity of the PRM mode(Fig.3E).For analytes 28 and 31(precursor ion m/z 549.1602),false positive ions could be observed in FS/ddMS2and t-SIM/ddMS2modes,but not observed in PRM mode(Fig.S4).

    Fig.3.Comparison of FS/ddMS2,t-SIM/ddMS2,and PRM scan modes by analyzing a mixed reference standard sample.(A-C)Working principles and sensitivities for the three scan modes.The flow of precursor ions(bigger dots)and product ions(smaller dots)are indicated using green and red arrows,respectively.(D)MS parameters and duty cycles for analytes 20,21,24,and 26.(E)Extracted ion chromatograms(FS/ddMS2and t-SIM/ddMS2)and PRM chromatogram for analytes 28 and 31.Concentrations of the four analytes were 100 ng/mL.Images for C-trap and Orbitrap detector were obtained from the producer’s website(http://www.thermo.com)and modified by the authors.DP:data point;S/N:signal to noise;CE:collision energy.

    3.4.Method validation

    3.4.1.Linearity,dynamic ranges,and limits of detection

    Calibration curves of analytes 1-19 detected by UPLC-UV were constructed by plotting the peak areas(Y)against the concentrations(X).Calibration curves of analytes 20-35 detected by UPLC/QOrbitrap-MS were constructed by plotting the analyte/IS peak area ratio(Y)against the correspondent concentration(X).ISs were used to ensure precision of the analyses.Phenylpropanolamine(IS1)and daidzin(IS2)corresponded to the analytes 20-22 and 23-35,respectively.Standard calibration curves of the 19 analytes in UPLCDAD analysis showed good linearity within the range of 0.08-246.00μg/mL(r2> 0.9995).The analyte/IS peak area ratio of the 16 analytes in UHPLC/Q-Orbitrap-MS analysis showed good correlation with concentrations(r2>0.99)within the range of 2.00-27840.00 ng/mL(Table 2).The limits of detection(LOD,S/N=3)for UPLC-DAD and UHPLC/Q-Orbitrap-MS methods were 40.00 to 3760.00 ng/mL and 0.02 to 5.90 ng/mL,respectively.

    3.4.2.Precision,repeatability,and stability

    Intra-and inter-day precisions were assessed by testing a sample solution in the same day for six times and on three consecutive days.The relative standard deviation(RSD)values for intra-and inter-day precisions ranged from 0.2% to 5.6% and 0.8% to 7.2%,respectively,indicating acceptable precision of the method.The repeatability was described by analyzing six samples(XFK-6)prepared using the same method.The results indicated that the sample preparation method was repeatable with RSD values ranging from 0.1% to 8.6%.The stability was evaluated by analyzing the same sample solution at 10°C after 0,2,4,8,12,and 24 h.RSD values of the analytes within 24 h ranged from 0.1% to 5.8%,indicating the analytes were stable.The detailed data are listed in Table S1.

    3.4.3.Accuracy

    The accuracy was validated by spiking the reference solutions to a real XFK sample.For analytes 1-19,accuracy analysis was conducted at 100% concentration level.For analytes 20-35,recovery was conducted at 80%,100%,and 120% concentratio levels.Samples XFK-6 and XFK-5 were used for accuracy study of 1-14,16,17,19-35 and 15,18,respectively.Recoveries were calculated by the formula:recovery(%)= (detected amount-original amount)/spiked amount×100%.Recoveries of the 19 analytes detected by UPLC-DAD ranged from 84.5%-113.8% with RSD values ranging from 0.2% to 3.1%.Recoveries of the 16 analytes detected by UHPLC/Q-Orbitrap-MS ranged from 72.4%-118.5% with RSD values ranging from 0.5% to 10.7%.The data are shown in Table S2.

    3.4.4.Sample analysis

    The validated method was applied to analyze 18 batches of XFK formula,including 10 batches of XFK oral solution and 8 batches of granules from three pharmaceutical companies(Figs.4A and S5).Due to different formulation methods and packages,we converted the concentrations(mg/g or mg/mL)into the contents in a single dose(mg/3 g for GR[a]and GR[b],mg/4 g for GR[c],and mg/10 mL for OS[a])to facilitate the comparison among different formulations.All of the oral solution samples met the requirements of the Chinese Pharmacopoeia(20 and 21≥1.8 mg/dose)[5].The contents of analytes 20 and 21 varied from 1.09 mg/dose to 2.21 mg/dose in granules from different pharmaceutical companies.Baicalin(15)from HQ was the most abundant component in all samples(38.92±4.35 mg/dose for oral solution and 17.03±7.87 mg/dose for granules).The total contents of the 35 analytes in oral solutions showed slight variations ranging from 145.92-206.69 mg/dose,while significant variations (49.59-133.22 mg/dose) were observed in different granules.The results are shown in Table S3.

    The quantitation results were then analyzed by principal component analysis(PCA)using SIMCA-P software(version 13.0).The first and second principal components accounted for 68.8% and 18.8% of the variation,respectively.Different formula samples were grouped in different clusters in Fig.4B.GR[a]and GR[c]were closer due to their similar chemical contents.Partial least squares discrimination analysis(PLS-DA)was then used to explore the variables that contributed to the grouping of the samples(Fig.S6).As shown in Fig.4C,contents of 5 and 33 from ZM,9 from LQ,15 and 34 from HQ,28 and 31 from GC,and 24 from KXR showed the highest intra-group variance,as suggested by the largest variable importance in projection(VIP)values(>1.20).For example,the contents of 5 and 33 from ZM were much higher in GR[b](16.90±1.66 mg/dose)than in other formulas(5.72±0.95 mg/dose for OS[a],2.87±0.17 mg/dose for GR[a],and 3.06±0.16 mg/dose for GR[c]).The data indicated the different qualities for the crude drugs used to prepare the XFK formula.

    Fig.4.Contents of 35 analytes in XFK and their principal component analysis.(A)Contents of 35 analytes in four different XFK formulas;(B)PCA scatter plots for 18 batches of formulas;(C)variable importance in projection(VIP)values for 35 analytes in different XFK formulations.OS[a],oral solution from company a,GR[a],granules from company a,GR[b],granules from company b,GR[c],granules from company c.Red asterisk represented the analytes with the highest intra-group variance.

    3.5.COX-2 inhibitory activities of chemical markers

    In present study,in vitro COX-2 inhibitory activities of the 35 analytes,10 herbs,and 4 different formulations were investigated to discover the main anti-inflammatory constituents of XFK.XFK formula,along with JYH,LQ,HQ,YXC,and ZM,exhibited potent COX-2 inhibitory activities(inhibitory rate≥60%)(Fig.5A).These five herbs might be responsible for the anti-inflammatory activity of XFK.For single analytes,phenolic acids(1-4,11,12,14)from JYH,phenylethanoid glycosides(8,9,10,23,26)from LQ,xanthone(5)from ZM,and flavonoid glycosides(15 from HQ,25 and 29 from JYH,30 from YXC)exhibited significant inhibition activities at 10 μM level(inhibitory rate≥ 80%,Fig.5B).

    Fig.5.COX-2 inhibitory activities XFK formula,herbs,and 35 analytes.(A)COX-2 inhibitory activities of XFK formula and herbs at 10 μg/mL,*P<0.05;(B)COX-2 inhibitory activities of 35 analytes at 10μM.

    Among the COX-2 inhibitors,8,9,10,and 15 were the major components.Their contents reached 81.70±7.47 mg/dose for OS[a],50.41±0.18 mg/dose for GR[b],23.28±2.00 mg/dose for GR[a],and 28.33±6.13 mg/dose for GR[c].Combining their contents and bioactivities,8,9,10,and 15 could be the major anti-inflammatory constituents for XFK.This was further supported by the activities of different formulations.For example,XFK oral solution exhibited higher inhibitory rate(78%±4%)than granules(average inhibitory rates of three granules at 50%±5%).Accordingly,the four effective constituents are higher in XFK oral solution.

    4.Conclusions

    In this study, quantitative analyses and bioactivity test werecombined to elucidate the anti-inflammatory constituents in XFK.Firstly, UPLC-DAD and UHPLC/Q-Orbitrap-MS methods were established and validated to quantify 35 analytes in different XFK formulations. The total contents of the 35 analytes in different XFKformulations showed significant variations ranging from75.69-269.46 mg/dose. Further COX-2 inhibitory assay revealed thatJYH, LQ, HQ, YXC, and ZM might be responsible for the antiinflammatory activity of XFK formula. Four major analytes 8, 9, 10,and 15 exhibited high abundance (total contents varied from 21.79 to91.20 mg/dose in different formulations) and potent COX-2 inhibitionactivities (inhibitory rate ≥ 80%) were proposed as the major effectivecomponents of XFK. The work also provided an effective strategy fordiscovery of effective constituents in multi-herb formulas.Declaration of competing interest

    The authors declare that there are no conflicts of interest.

    Acknowledgements

    The work was supported by the National Key Research and Development Program of China(Grant No.:2018YFC1707304,2018YFC1707301),Beijing Natural Science Foundation(Grant No.:JQ18027),and National Natural Science Foundation of China(Grant No.:81725023).We thank Sunflower Pharmaceutical Co.,Ltd.for providing the samples.

    Appendix A.Supplementary data

    Supplementary data to this article can be found online at https://doi.org/10.1016/j.jpha.2021.01.003.

    蜜臀久久99精品久久宅男| 精品久久久久久久人妻蜜臀av| 亚洲人成网站高清观看| 亚洲色图av天堂| 亚洲精品乱码久久久久久按摩| 菩萨蛮人人尽说江南好唐韦庄 | 99久久中文字幕三级久久日本| 噜噜噜噜噜久久久久久91| 国产成人影院久久av| 免费黄网站久久成人精品| 欧美人与善性xxx| 久久精品人妻少妇| 亚洲最大成人av| 在线免费十八禁| 老女人水多毛片| 网址你懂的国产日韩在线| 中文欧美无线码| 嫩草影院精品99| 成人av在线播放网站| 丝袜喷水一区| 不卡一级毛片| 亚洲内射少妇av| 久久国内精品自在自线图片| 日韩大尺度精品在线看网址| 亚洲成人久久性| 亚洲精品色激情综合| 亚洲婷婷狠狠爱综合网| 乱码一卡2卡4卡精品| 亚洲色图av天堂| 欧美bdsm另类| 亚洲五月天丁香| 欧美极品一区二区三区四区| 人妻制服诱惑在线中文字幕| 26uuu在线亚洲综合色| 超碰av人人做人人爽久久| 美女 人体艺术 gogo| 久久这里有精品视频免费| 色哟哟哟哟哟哟| 亚洲一级一片aⅴ在线观看| 精品久久久久久久末码| 日本一本二区三区精品| 国产免费一级a男人的天堂| 国产成年人精品一区二区| 18禁裸乳无遮挡免费网站照片| 国产91av在线免费观看| 99视频精品全部免费 在线| 九九久久精品国产亚洲av麻豆| 99国产精品一区二区蜜桃av| 国产av一区在线观看免费| 99久久九九国产精品国产免费| 黄片wwwwww| 日本爱情动作片www.在线观看| 日本黄色片子视频| 黄色配什么色好看| 久久久成人免费电影| 亚洲精品影视一区二区三区av| 午夜精品国产一区二区电影 | 国产一区二区亚洲精品在线观看| 一级黄片播放器| 欧美变态另类bdsm刘玥| 亚洲成人精品中文字幕电影| 欧美在线一区亚洲| 国产大屁股一区二区在线视频| 女同久久另类99精品国产91| 五月玫瑰六月丁香| 中文字幕精品亚洲无线码一区| 亚洲人成网站在线观看播放| 亚洲国产精品sss在线观看| 久久综合国产亚洲精品| 免费黄网站久久成人精品| 国产精品一区二区在线观看99 | 国语自产精品视频在线第100页| 看非洲黑人一级黄片| 岛国在线免费视频观看| 亚洲av电影不卡..在线观看| 人体艺术视频欧美日本| 久久中文看片网| 插阴视频在线观看视频| 亚洲三级黄色毛片| av天堂在线播放| 青青草视频在线视频观看| 免费看av在线观看网站| 亚洲色图av天堂| 久久精品国产亚洲av香蕉五月| 六月丁香七月| 亚洲欧美日韩高清在线视频| 亚洲av中文av极速乱| 国产视频首页在线观看| 久久精品国产99精品国产亚洲性色| 国产精品99久久久久久久久| av国产免费在线观看| 岛国在线免费视频观看| 国产在线精品亚洲第一网站| 日日啪夜夜撸| 亚洲欧美日韩卡通动漫| 国产在线精品亚洲第一网站| 两个人的视频大全免费| 在线观看美女被高潮喷水网站| 蜜桃亚洲精品一区二区三区| 夜夜看夜夜爽夜夜摸| 又粗又硬又长又爽又黄的视频 | 国产午夜福利久久久久久| 五月玫瑰六月丁香| 热99re8久久精品国产| 欧美日本视频| 精品国产三级普通话版| 久久久久九九精品影院| 小蜜桃在线观看免费完整版高清| 男女下面进入的视频免费午夜| 久久韩国三级中文字幕| 久久久精品大字幕| 亚洲欧美精品专区久久| 日韩欧美精品免费久久| 欧美变态另类bdsm刘玥| 亚洲欧美成人综合另类久久久 | 国产高清三级在线| 国产极品天堂在线| 亚洲高清免费不卡视频| 成人一区二区视频在线观看| 美女大奶头视频| 色播亚洲综合网| 乱码一卡2卡4卡精品| 国产女主播在线喷水免费视频网站 | 国产在视频线在精品| 免费观看在线日韩| 国产精品一及| 丝袜美腿在线中文| 免费人成在线观看视频色| 欧美+亚洲+日韩+国产| 一级毛片久久久久久久久女| 午夜久久久久精精品| 日本黄大片高清| 九九热线精品视视频播放| 亚洲精品456在线播放app| av天堂中文字幕网| 久久久久久久午夜电影| 日本黄色片子视频| 高清午夜精品一区二区三区 | 成人永久免费在线观看视频| 国产精品,欧美在线| 成人av在线播放网站| 亚洲性久久影院| 少妇丰满av| 免费观看的影片在线观看| 免费一级毛片在线播放高清视频| 亚洲欧美成人综合另类久久久 | 精品国内亚洲2022精品成人| av国产免费在线观看| videossex国产| 成人综合一区亚洲| 91av网一区二区| 国产精品久久视频播放| 少妇的逼水好多| 亚洲国产精品sss在线观看| 18禁在线无遮挡免费观看视频| 精品久久久久久久久久免费视频| 只有这里有精品99| 男人的好看免费观看在线视频| 99久久精品一区二区三区| 欧美日韩在线观看h| 欧美成人a在线观看| 天堂网av新在线| 99国产极品粉嫩在线观看| 免费人成视频x8x8入口观看| 国产在线男女| 菩萨蛮人人尽说江南好唐韦庄 | 国产精品久久久久久久久免| 国产女主播在线喷水免费视频网站 | 日韩人妻高清精品专区| 嫩草影院入口| 国产在线男女| av免费在线看不卡| 干丝袜人妻中文字幕| 国产精品1区2区在线观看.| 床上黄色一级片| 成年免费大片在线观看| 亚洲电影在线观看av| 乱人视频在线观看| 男女那种视频在线观看| 3wmmmm亚洲av在线观看| 午夜激情福利司机影院| 在线天堂最新版资源| 久久久精品94久久精品| 精品99又大又爽又粗少妇毛片| 免费观看在线日韩| 两个人视频免费观看高清| 久久久久久久亚洲中文字幕| 日韩中字成人| 99久久精品国产国产毛片| 日日摸夜夜添夜夜添av毛片| 亚洲美女视频黄频| 桃色一区二区三区在线观看| av在线蜜桃| 国产精品精品国产色婷婷| 日本在线视频免费播放| av在线观看视频网站免费| 国模一区二区三区四区视频| 成人永久免费在线观看视频| 麻豆国产97在线/欧美| 成人亚洲精品av一区二区| 日本-黄色视频高清免费观看| 久久久精品欧美日韩精品| 久久精品国产亚洲av天美| 男女啪啪激烈高潮av片| 春色校园在线视频观看| 不卡视频在线观看欧美| 久久久久久伊人网av| 欧美xxxx性猛交bbbb| 亚洲色图av天堂| 国产精品日韩av在线免费观看| 亚洲av不卡在线观看| av在线蜜桃| 91在线精品国自产拍蜜月| 国内精品久久久久精免费| 亚洲精品久久久久久婷婷小说 | 日韩 亚洲 欧美在线| 九草在线视频观看| 日韩强制内射视频| 黄色欧美视频在线观看| 少妇丰满av| 国产成人aa在线观看| 伦理电影大哥的女人| 国产黄a三级三级三级人| 免费看日本二区| 永久网站在线| 一边摸一边抽搐一进一小说| 天天一区二区日本电影三级| 九色成人免费人妻av| 国产精品久久电影中文字幕| 2021天堂中文幕一二区在线观| 插阴视频在线观看视频| 久久久久久久久久久免费av| 美女xxoo啪啪120秒动态图| 久久久a久久爽久久v久久| 久久99蜜桃精品久久| 国产精华一区二区三区| 欧美最新免费一区二区三区| 亚洲七黄色美女视频| 美女被艹到高潮喷水动态| 成人国产麻豆网| 中文字幕精品亚洲无线码一区| 美女脱内裤让男人舔精品视频 | 99久久精品一区二区三区| 久久精品国产亚洲av香蕉五月| 久久精品国产亚洲网站| 国内少妇人妻偷人精品xxx网站| 嫩草影院入口| 国内精品一区二区在线观看| 一级毛片久久久久久久久女| 国产成人午夜福利电影在线观看| 美女黄网站色视频| 亚洲人成网站在线播| 久久久午夜欧美精品| 国产成人影院久久av| 国产探花极品一区二区| 欧美性猛交╳xxx乱大交人| 国产成人精品一,二区 | 国产午夜精品一二区理论片| 国产亚洲91精品色在线| 久久人妻av系列| 黄片wwwwww| 99热全是精品| 男女啪啪激烈高潮av片| av女优亚洲男人天堂| 三级经典国产精品| 嫩草影院新地址| 99精品在免费线老司机午夜| 国产免费男女视频| 精品久久久久久久久av| 国产精品麻豆人妻色哟哟久久 | 国产极品天堂在线| 国产麻豆成人av免费视频| 99视频精品全部免费 在线| 最近中文字幕高清免费大全6| 精品久久久久久久人妻蜜臀av| 91麻豆精品激情在线观看国产| 婷婷色av中文字幕| av天堂在线播放| 久久久久网色| 亚洲欧美精品自产自拍| 一级毛片我不卡| 内射极品少妇av片p| 亚洲国产欧美人成| 日韩视频在线欧美| 亚洲国产欧洲综合997久久,| 国产 一区 欧美 日韩| 国产 一区精品| 午夜激情欧美在线| 国产乱人视频| 欧美xxxx黑人xx丫x性爽| 中文在线观看免费www的网站| 免费看a级黄色片| 男女视频在线观看网站免费| 成人欧美大片| 国产亚洲91精品色在线| 成人美女网站在线观看视频| 美女黄网站色视频| 国产亚洲av片在线观看秒播厂 | 国产精品一区二区在线观看99 | 在线观看一区二区三区| 成人二区视频| 国产精品福利在线免费观看| 在线免费观看不下载黄p国产| 永久网站在线| 久久午夜亚洲精品久久| 亚洲电影在线观看av| 性插视频无遮挡在线免费观看| 欧美人与善性xxx| 亚洲精品国产成人久久av| 99热这里只有是精品50| 国产成年人精品一区二区| 亚洲第一电影网av| 草草在线视频免费看| 国产精品蜜桃在线观看 | 99在线视频只有这里精品首页| 一级黄色大片毛片| 一边摸一边抽搐一进一小说| 国产老妇女一区| 卡戴珊不雅视频在线播放| 熟妇人妻久久中文字幕3abv| 天美传媒精品一区二区| 成人永久免费在线观看视频| 综合色丁香网| 内射极品少妇av片p| 中国国产av一级| 赤兔流量卡办理| 亚洲精品粉嫩美女一区| a级一级毛片免费在线观看| 久久精品影院6| 亚洲天堂国产精品一区在线| av免费在线看不卡| 十八禁国产超污无遮挡网站| 免费不卡的大黄色大毛片视频在线观看 | 久久6这里有精品| 久久久久久久久久成人| 精品免费久久久久久久清纯| 人妻制服诱惑在线中文字幕| 欧美在线一区亚洲| 亚洲av二区三区四区| 亚洲七黄色美女视频| 国产精品99久久久久久久久| 国产黄片视频在线免费观看| 天天一区二区日本电影三级| 欧美zozozo另类| 精品少妇黑人巨大在线播放 | 国产单亲对白刺激| 五月伊人婷婷丁香| 蜜桃久久精品国产亚洲av| 少妇被粗大猛烈的视频| 亚洲美女搞黄在线观看| 精品一区二区免费观看| 久久人妻av系列| 色吧在线观看| 十八禁国产超污无遮挡网站| 又粗又硬又长又爽又黄的视频 | 久久婷婷人人爽人人干人人爱| 国产精品蜜桃在线观看 | 秋霞在线观看毛片| 欧美日本亚洲视频在线播放| 婷婷色av中文字幕| 久久精品国产自在天天线| 又爽又黄无遮挡网站| 免费看a级黄色片| 欧美色欧美亚洲另类二区| 亚洲激情五月婷婷啪啪| 少妇熟女欧美另类| 最好的美女福利视频网| 国产精品久久久久久久久免| 男女下面进入的视频免费午夜| 校园春色视频在线观看| 男女做爰动态图高潮gif福利片| 成年av动漫网址| 色噜噜av男人的天堂激情| 国产精品久久久久久久久免| 成人国产麻豆网| 欧洲精品卡2卡3卡4卡5卡区| 亚洲av男天堂| 国产日韩欧美在线精品| 变态另类丝袜制服| 美女被艹到高潮喷水动态| 色哟哟·www| 久久这里有精品视频免费| 国产精品日韩av在线免费观看| 麻豆一二三区av精品| 久久久久久国产a免费观看| 国产av麻豆久久久久久久| 有码 亚洲区| 国产在线男女| 免费观看人在逋| 伊人久久精品亚洲午夜| 亚洲aⅴ乱码一区二区在线播放| www.av在线官网国产| 久久韩国三级中文字幕| 成人亚洲欧美一区二区av| 18禁裸乳无遮挡免费网站照片| 成人特级黄色片久久久久久久| 99久久精品国产国产毛片| 寂寞人妻少妇视频99o| 亚洲国产精品成人久久小说 | 国产精品蜜桃在线观看 | 久久久国产成人精品二区| 久久人人精品亚洲av| 久久草成人影院| 26uuu在线亚洲综合色| 亚洲精品乱码久久久v下载方式| 男女做爰动态图高潮gif福利片| av黄色大香蕉| 国产国拍精品亚洲av在线观看| 国产激情偷乱视频一区二区| 欧美性猛交╳xxx乱大交人| 亚洲在久久综合| 一区二区三区四区激情视频 | 国产精品嫩草影院av在线观看| 色哟哟哟哟哟哟| 可以在线观看的亚洲视频| 2022亚洲国产成人精品| 免费看日本二区| 男女下面进入的视频免费午夜| 在线播放国产精品三级| АⅤ资源中文在线天堂| 日韩成人伦理影院| av视频在线观看入口| 欧美最新免费一区二区三区| 精品欧美国产一区二区三| 国内少妇人妻偷人精品xxx网站| 99热全是精品| 少妇被粗大猛烈的视频| 久久久a久久爽久久v久久| 国产国拍精品亚洲av在线观看| 小蜜桃在线观看免费完整版高清| 校园人妻丝袜中文字幕| 婷婷精品国产亚洲av| 成人无遮挡网站| 成人特级黄色片久久久久久久| 国产精品久久久久久精品电影小说 | 日本与韩国留学比较| 亚洲精品456在线播放app| 搡女人真爽免费视频火全软件| 最近最新中文字幕大全电影3| 亚洲av成人av| 91精品一卡2卡3卡4卡| 免费观看人在逋| 国产黄片美女视频| 国产亚洲av嫩草精品影院| 国产精品综合久久久久久久免费| 欧美色欧美亚洲另类二区| 日韩中字成人| 国产精品麻豆人妻色哟哟久久 | 2021天堂中文幕一二区在线观| 久久精品久久久久久噜噜老黄 | 国产一区二区激情短视频| 级片在线观看| avwww免费| 国产精品久久久久久久久免| 日日啪夜夜撸| 国产人妻一区二区三区在| 国产熟女欧美一区二区| 天天躁日日操中文字幕| 国产探花极品一区二区| 激情 狠狠 欧美| 国产欧美日韩精品一区二区| 中文精品一卡2卡3卡4更新| 国产综合懂色| 久久婷婷人人爽人人干人人爱| 亚洲久久久久久中文字幕| eeuss影院久久| 久久久国产成人精品二区| 国产高清激情床上av| 亚洲在线观看片| 亚洲国产欧美人成| 国内精品一区二区在线观看| 亚洲va在线va天堂va国产| 12—13女人毛片做爰片一| 亚洲成a人片在线一区二区| 男女下面进入的视频免费午夜| 日本一本二区三区精品| 青春草亚洲视频在线观看| 成年版毛片免费区| 国产淫片久久久久久久久| 亚洲真实伦在线观看| 国产91av在线免费观看| 天堂√8在线中文| 亚洲欧美成人精品一区二区| 亚洲高清免费不卡视频| 22中文网久久字幕| 晚上一个人看的免费电影| 丰满人妻一区二区三区视频av| 小蜜桃在线观看免费完整版高清| 女人十人毛片免费观看3o分钟| 欧美三级亚洲精品| 韩国av在线不卡| 久久精品夜夜夜夜夜久久蜜豆| 天堂网av新在线| 美女xxoo啪啪120秒动态图| .国产精品久久| 人妻少妇偷人精品九色| 黄片wwwwww| 国产日韩欧美在线精品| 亚洲av一区综合| 波野结衣二区三区在线| 精品人妻熟女av久视频| 少妇丰满av| 日韩大尺度精品在线看网址| 一个人免费在线观看电影| 三级男女做爰猛烈吃奶摸视频| 乱系列少妇在线播放| www日本黄色视频网| 亚洲av成人精品一区久久| 成人漫画全彩无遮挡| 婷婷精品国产亚洲av| 成人性生交大片免费视频hd| 黄片wwwwww| 啦啦啦观看免费观看视频高清| 少妇裸体淫交视频免费看高清| 在线观看免费视频日本深夜| 少妇人妻精品综合一区二区 | 色视频www国产| 久久亚洲精品不卡| 国产激情偷乱视频一区二区| 一个人观看的视频www高清免费观看| 亚洲成人久久性| 欧美成人精品欧美一级黄| 日韩三级伦理在线观看| 亚洲无线观看免费| 午夜老司机福利剧场| 亚洲精品乱码久久久久久按摩| 蜜桃久久精品国产亚洲av| 精品久久国产蜜桃| 色综合色国产| 成人av在线播放网站| 欧美bdsm另类| 欧美最黄视频在线播放免费| 女人十人毛片免费观看3o分钟| av在线蜜桃| 99久国产av精品| 亚洲国产欧洲综合997久久,| 亚洲最大成人中文| 久久久久网色| 在线观看一区二区三区| 成人午夜精彩视频在线观看| 美女内射精品一级片tv| 日韩一本色道免费dvd| 亚洲图色成人| 欧美性猛交黑人性爽| 久久久久久久久久久免费av| 少妇人妻精品综合一区二区 | 极品教师在线视频| 哪里可以看免费的av片| 免费无遮挡裸体视频| 1024手机看黄色片| 哪个播放器可以免费观看大片| 精品久久久噜噜| av在线播放精品| 91久久精品电影网| 高清毛片免费看| 亚洲精品国产成人久久av| 日本成人三级电影网站| 日本黄色片子视频| 亚洲国产欧美在线一区| 欧美激情久久久久久爽电影| 国产成人aa在线观看| 菩萨蛮人人尽说江南好唐韦庄 | .国产精品久久| 99热全是精品| 99热只有精品国产| 国产伦一二天堂av在线观看| 99久久精品国产国产毛片| 哪个播放器可以免费观看大片| 免费观看精品视频网站| 亚洲欧美日韩卡通动漫| 日韩av不卡免费在线播放| 尤物成人国产欧美一区二区三区| 亚洲最大成人av| 亚洲国产高清在线一区二区三| 男人狂女人下面高潮的视频| 波多野结衣高清作品| 国产精品一区二区三区四区免费观看| 亚洲人与动物交配视频| 国产私拍福利视频在线观看| 亚洲在久久综合| 简卡轻食公司| www.色视频.com| 国产伦在线观看视频一区| 免费观看的影片在线观看| 婷婷色综合大香蕉| 麻豆一二三区av精品| 99riav亚洲国产免费| 成人三级黄色视频| 国产探花在线观看一区二区| 丝袜美腿在线中文| 久久久午夜欧美精品| 国产精品麻豆人妻色哟哟久久 | av在线天堂中文字幕| 日本黄大片高清| 蜜桃久久精品国产亚洲av| 我要看日韩黄色一级片| 寂寞人妻少妇视频99o| 免费黄网站久久成人精品| 青春草亚洲视频在线观看| 亚洲中文字幕一区二区三区有码在线看| 蜜桃亚洲精品一区二区三区| 色哟哟·www| 亚洲人成网站在线播放欧美日韩| 成人二区视频| 国国产精品蜜臀av免费| 全区人妻精品视频| 蜜桃亚洲精品一区二区三区| 久久精品综合一区二区三区| 在线观看免费视频日本深夜| 久久草成人影院| 丰满乱子伦码专区| 欧美三级亚洲精品| 内地一区二区视频在线| 青春草亚洲视频在线观看| 悠悠久久av|