亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        Analysis of anticancer compound,indole-3-carbinol,in broccoli using a new ultrasound-assisted dispersive-filter extraction method based on poly(deep eutectic solvent)-graphene oxide nanocomposite

        2022-05-26 07:02:28YnnYunHunhunChenYehongHnFengxiQioHongyunYn
        Journal of Pharmaceutical Analysis 2022年2期

        Ynn Yun,Hunhun Chen,Yehong Hn,Fengxi Qio,Hongyun Yn,*

        aKey Laboratory of Medicinal Chemistry and Molecular Diagnosis of Ministry of Education,College of Pharmaceutical Science,Hebei University,Baoding,Hebei,071002,China

        bKey Laboratory of Public Health Safety of Hebei Province,Institute of Life Science and Green Development,College of Public Health,Hebei University,Baoding,Hebei,071002,China

        cCollege of Biochemistry,Baoding University,Baoding,Hebei,071000,China

        Keywords:

        Indole-3-carbinol

        Broccoli

        Poly(deep eutectic solvent)-graphene oxide

        Ultrasound-assisted dispersive-filter

        extraction

        Rapid analysis

        A B S T R A C T

        Indole-3-carbinol(I3C),an important anticancer compound found in broccoli,has attracted considerable attention.The rapid extraction and accurate analysis of I3C in the pharmaceutical industry in broccoli is challenging as I3C is unstable at low pH and high temperature.In this study,a rapid,accurate,and lowcost ultrasound-assisted dispersive-filter extraction(UADFE)technique based on poly(deep eutectic solvent)-graphene oxide(PDES-GO)adsorbent was developed for the isolation and analysis of I3C in broccoli for the first time.PDES-GO with multiple adsorption interactions and a fast mass transfer rate was synthesized to accelerate adsorption and desorption.UADFE was developed by combining dispersive solid-phase extraction(DSPE)and filter solid-phase extraction(FSPE)to realize rapid extraction and separation.Based on the above two strategies,the proposed PDES-GO-UADFE method coupled with high-performance liquid chromatography(HPLC)allowed the rapid(15-16 min),accurate(84.3%-96.4%),and low-cost(adsorbent:3.00 mg)analysis of I3C in broccoli and was superior to solid-phase extraction,DSPE,and FSPE methods.The proposed method showed remarkable linearity(r=0.9998;range:0.0840-48.0 μg/g),low limit of quantification(0.0840μg/g),and high precision(relative standard deviation≤5.6%).Therefore,the PDES-GO-UADFE-HPLC method shows significant potential in the field of pharmaceutical analysis for the separation and analysis of anti-cancer compounds in complex plant samples.

        1.Introduction

        Indole-3-carbinol(I3C)has received significant attention in the field of natural product chemistry owing to its anticancer properties[1,2].I3C is naturally found in cruciferous plants,particularly broccoli[3].When the broccoli tissues are damaged(such as during mastication or homogenization),glucobrassicin,an important glucosinolate,comes in contact with myrosinase,an endogenous enzyme present in cruciferous plants,to release I3C[4].I3C can interact with the WW domain containing E3 ubiquitin protein ligase 1,which affects the function of the tumor suppressor gene phosphatase and tensin homologue to afford anticancer effects[5].I3C performs anticancer activities against lung,nasopharyngeal,liver,ovarian,and other types of cancer[6,7].Thus,it is important to develop an analytical method to extract and analyze I3C in broccoli.

        It has been reported that I3C is unstable under many conditions such as low pH and high temperature.Under acidic conditions,I3C is rapidly converted into a series of acid condensation compounds[8-10].The instability of I3C limits its accurate analysis in broccoli.Notably,the sample preparation process accounts for 50%-70% of the entire analysis time[11].In addition to considering the experimental conditions that are likely to cause I3C instability during its analysis,decreasing the time for sample preparation can effectively reduce the effect of I3C instability on the analytical results;therefore,developing a rapid sample preparation method is necessary.

        Most sample preparation methods such as solid-phase extraction(SPE)and derivatization methods rely primarily on the extraction method and the adsorbent[12-15].Therefore,the rapid extraction and accurate determination of I3C in broccoli can be achieved considering these two aspects[16,17].Several extraction methods such as SPE,dispersive solid-phase extraction(DSPE),and matrix solid-phase dispersion(MSPD)have been developed[18,19].SPE affords high recovery and remarkable reproducibility,but the adsorbent amount required is large(tens to hundreds of milligrams),and it often requires an activation step(pretreatment with methanol and water)before loading the sample solution[20,21].In DSPE,the adsorbent is uniformly dispersed to ensure sufficient contact with the sample solution,washing solvent,and elution solvent,such that adsorption and desorption occur rapidly[22,23].However,each step in DSPE requires additional centrifugation,which limits the sample preparation process[24].In MSPD,extraction,filtration,and purification steps are completed in one step,but large amounts of the sample and adsorbent are lost during the transfer of the ground semidry mixture from the mortar to the SPE cartridge[25].Therefore,a new,rapid,and lowcost extraction method with a low sample loss is required to increase the speed and decrease the cost of the sample preparation method,which meets the requirements for I3C analysis in broccoli.

        In addition to the extraction method,adsorbents with high mass transfer rates and high adsorption capacities can accelerate the adsorption and desorption,thereby decreasing the sample preparation time[26,27].Conventional adsorbents such as C8,C18,silica,and Al2O3cannot meet the requirements of sample preparation because of a single adsorption interaction,low mass transfer rate,and low adsorption capacity.Therefore,several studies have focused on developing new adsorbents.Graphene oxide(GO),which has a two-dimensional sheet structure and a large specific surface area,contains many oxygen-containing functional groups that can be easily modified with other functional groups to improve its adsorption capacity[28,29].A deep eutectic solvent(DES)is a new type of ionic liquid that can be easily designed and synthesized and has multiple functional groups[30-32].DES has been used to modify various materials to increase the adsorption interactions and improve the adsorption capacity of the material[33-35].Therefore,it is proposed that grafting a thin layer of poly(deep eutectic solvent)(PDES)on the GO surface can not only enrich the adsorption interactions and increase the capacity for adsorbing the target analyte,but also afford a high mass transfer rate because of the thin PDES layer,which can increase the speed and extraction efficiency of the sample preparation process for the extraction and analysis of I3C in broccoli.

        In this study,a new ultrasound-assisted dispersive-filter extraction(UADFE)method based on the PDES-GO adsorbent was proposed for the rapid,accurate,and low-cost analysis of I3C in broccoli for the first time.The PDES-GO-UADFE method combines the advantages of PDES-GO(high mass transfer rate and high adsorption capacity)and UADFE(rapid extraction and separation).The extraction parameters were optimized,and the proposed method was validated in terms of linearity,accuracy,and precision.Finally,the developed method was applied to analyze I3C in different varieties and parts of the broccoli samples,and it was also used to determine the effect of high-temperature processing on the I3C content in broccoli.

        2.Experimental

        2.1.Chemicals and reagents

        Choline chloride was purchased from Guangfu Chemical Co.,Ltd.(Tianjin,China).Methacrylic acid was purchased from Sigma Aldrich(Shanghai,China).Acetonitrile,methanol,ammonia water(25%),and dichloromethane were purchased from Kermel Chemical Co.,Ltd.(Tianjin,China).Silica,NH2,strong cation exchange adsorbent,and C18were purchased from Bonna-Agela Technologies(Tianjin,China).Reversible addition-fragmentation chain transfer agent functionalized graphene oxide(GO-RAFT)was synthesized in this laboratory[36].Standard stock solution(1.00 mg/mL)of I3C was prepared in methanol.

        2.2.Synthesis and characterization of adsorbent

        Choline chloride and methacrylic acid were mixed in a 1:2 M ratio,and the mixture was stirred in a water bath maintained at 80°C until a transparent and uniform solution(DES)was formed,which was used for the subsequent synthesis of PDES-GO.GO-RAFT(300 mg)was dispersed in acetonitrile(60 mL)and sonicated for 1 h.DES(3.00 g)and 2,2'-azobis(2-methylpropionitrile)(75.0 mg)were then added to the mixture and reacted at 60°C for 24 h under N2atmosphere to prepare PDES-GO,which was then washed with methanol and water.Thereafter,the material was freeze-dried under vacuum.

        The I3C adsorption performance of PDES-GO was evaluated via a static adsorption experiment.After the addition of PDES-GO(2.00 mg)and different concentrations of standard solution(10.0,30.0,50.0,80.0,100,and 200μg/mL,2.0 mL),the centrifuge tubes were placed in a constant temperature oscillator(300 r/min,25°C)for 12 h.The supernatant was then analyzed using highperformance liquid chromatography(HPLC).

        2.3.Sample preparation

        Broccoli samples were purchased from the local markets(Baoding,China).For the sample preparation,100 g of broccoli and 100 g of phosphate buffer solution(PBS,pH 7.4,0.1 M)were mixed and homogenized for 5 min.Then,25.0 g of the homogenized mixture was weighed and placed in a constant temperature oscillator at 25°C for 3 h to release I3C.Thereafter,dichloromethane(10 mL)was added to the mixture,followed by magnetic stirring for 10 min.The mixture was centrifuged(10,000 r/min,5 min),and the upper solvent layer was transferred into an eggplant-shaped bottle.The extraction process was repeated twice.Subsequently,the organic phase(dichloromethane)was evaporated to dryness(30°C)and then reconstituted in 5.0 mL of PBS with ultrasonication(thrice).Finally,the reconstituted solution was filtered and transferred into a volumetric flask(25 mL)and diluted up to the mark with PBS.

        The procedure of the PDES-GO-UADFE method is shown in Fig.1.PDES-GO(3.00 mg)and the sample solution(0.5 mL)were mixed and sonicated for 10 min to extract I3C.The mixture was then drawn into a syringe,and a syringe filter(0.45μm)was connected to rapidly separate the adsorbent and the sample solution.The syringe filter was unplugged,and 0.5 mL of the washing solvent(PBS)was drawn into the syringe.The same syringe filter was connected to the syringe to wash out the impurities on the adsorbent(~1 min).Thereafter,1.0 mL of the elution solvent(ammonia water-acetonitrile,1:9,V/V)was drawn into the syringe,and the syringe filter and the needle were connected.Subsequently,the elution solvent was repeatedly drawn up and down into the syringe three times to elute sufficient I3C from the adsorbent(~2 min).This elution process was repeated once with the elution solvent(0.5 mL;~2 min).Finally,the eluent and PBS were mixed in a 1:1 ratio by vortexing for 10 s,followed by HPLC analysis.

        Fig.1.Schematic of poly(deep eutectic solvent)-graphene oxide (PDES-GO)-ultrasound-assisted dispersive-filter extraction (UADFE)-high-performance liquid chromatography(HPLC)method.

        2.4.Instruments and analytical conditions

        The morphology of PDES-GO was examined by scanning electron microscopy(SEM,Phenom-World BV,Eindhoven,The Netherlands).Fourier-transform infrared(FTIR)spectroscopy was performed using a Thermo Scientific Nicolet iS10 FTIR spectrometer(Thermo Fisher Scientific,Waltham,MA,USA).A Thermo Ulti-Mate3000 DGLC HPLC system(Thermo Fisher Scientific,Waltham,MA,USA)with a Chromeleon 7.2 workstation,diode array detector,and chromatographic column(Accucore C18,100 mm×4.6 mm,2.6μm)was employed for the analysis of I3C.The mobile phase consisted of acetonitrile-water(40:60,V/V)at a flow rate of 1.0 mL/min.The analysis wavelength of the detector was 218 nm,and the injection volume was 20μL.

        2.5.Comparison of the proposed method with other sample preparation methods

        The established extraction method(UADFE)was compared to other sample preparation methods(DSPE,filter solid-phase extraction(FSPE),and commercial SPE)[37-39].In the extraction process,these sample preparation methods mainly included the following steps:activation,sample loading,washing,and elution.The detailed parameters of the procedure are provided in Table 1[37-39].

        Table 1Procedure of different sample preparation methods.

        3.Results and discussion

        3.1.Synthesis and characterization of adsorbent

        The composition and structure of PDES-GO were determined using FTIR and SEM.As shown in Fig.2A,the peaks at 1500 and 1600 cm-1correspond to the C=C bonds,which comprise the GO skeleton,and the peaks at 3400,1715,and 1250 cm-1correspond to the O-H,C=O,and C-O-C bonds of GO,respectively.The peaks at 1100,1060,and 1000 cm-1were attributed to the Si-O-C,C=S,and-C-S-C-bonds,respectively,indicating that the RAFT agent was grafted onto the GO sheet.The peak intensity of C=O increased in the spectrum of PDES-GO,indicating that GO-RAFT was successfully modified with DES,as DES contains a carboxyl group.In terms of the morphology,the prepared PDES-GO has a threedimensional skeleton with a porous morphology,based on the SEM data(Fig.2B).The modification of the RAFT agent and PDES on the GO sheets provides support for GO,such that the PDES-GO sheets are not stacked closely.Around the hole,the edge of the PDES-GO sheets with a thin lamellar structure can be clearly observed,which is conducive to achieving the rapid adsorption and desorption of I3C.

        The I3C adsorption performance of PDES-GO was investigated via the static adsorption experiments using a series of I3C standard solutions(10.0-200μg/mL).As shown in Fig.2C,a high I3C adsorption with PDES-GO(~80.0μg/mg)is achieved when the adsorption dynamic equilibrium was attained,which was mainly due to multiple adsorption interactions between the functional groups of PDES-GO and I3C.Thereafter,a series of commercial adsorbents(silica,NH2,C18,and SCX)with different adsorption interactions were selected to further examine the multiple adsorption interactions of PDES-GO.Silica and NH2can form hydrogen bonding;C18is capable of hydrophobic interactions,and SCX can undergo ion exchange.As shown in Fig.2D,commercial adsorbents with a single adsorption interaction exhibited poor I3C adsorption performance,which limited their use in the extraction and isolation of I3C.Graphene has a higher adsorption capacity than the commercial adsorbents because graphene is capable of hydrophobic interactions and π-π conjugation for I3C adsorption.PDES-GO exhibits multiple adsorption interactions(π-π conjugation,hydrogen bonding,hydrophobic interactions,and electrostatic adsorption)after modification with PDES and affords the highest I3C adsorption(83.8μg/mg)among these adsorbents,which demonstrates the excellent adsorption performance of PDES-GO in comparison to the performances of other adsorbents for I3C extraction.

        Fig.2.(A)Fourier-transform infrared(FTIR)spectra and(B)scanning electron microscopy(SEM)data of poly(deep eutectic solvent)-graphene oxide(PDES-GO);(C)static I3C adsorption of PDES-GO;(D)I3C adsorption amounts of different adsorbents.GO-RAFT:reversible addition-fragmentation chain transfer agent functionalized GO.

        3.2.Optimization of sample preparation procedure

        As I3C in broccoli does not exist in a prototype form,it is necessary to release I3C through an autolysis step before analysis.In the autolysis step,broccoli was homogenized and placed in a constant temperature oscillator at 25°C for 3 h[40].During autolysis,glucobrassicin in broccoli was contacted with myrosinase in cruciferous plants to release I3C.After autolysis,the homogenized mixture was extracted with dichloromethane,followed by rotary evaporation and reconstitution to afford a sample solution containing I3C,which required further extraction and purification.Owing to the instability of I3C,the optimization of a series of extraction parameters(adsorbent amount,ultrasonication time,washing solvent,and elution solvent)in PDES-GO-UADFE was performed to afford rapid extraction and accurate analysis of I3C in broccoli.

        3.2.1.Optimizationofadsorbentdosage

        The volume of the sample solution was set to 0.5 mL,and the amount of adsorbent(1.00-5.00 mg)was determined because the adsorption performance of the adsorbent directly affected the extraction efficiency of the analyte(Fig.3A).When the adsorbent amount was 3.00 mg,the I3C loss ratio was<5%.A small amount of adsorbent was needed mainly because of the multiple interactions of PDES-GO,which could effectively adsorb I3C through π-π conjugation,hydrophobic interactions,electrostatic adsorption,and hydrogen bonding with the functional groups of I3C(benzene ring,-NH-,and-OH).

        3.2.2.Optimizationofultrasonicationtime

        As UADFE was a time-dependent extraction procedure,the ultrasonication time for extraction was investigated.As shown in Fig.3B,suitable extraction results(loss ratio of I3C<5%)were obtained within 10 min of ultrasonication,owing to the thin layer of PDES-GO sheet that facilitated the rapid adsorption and desorption of I3C.Ultrasound-assisted extraction mode facilitated sufficient dispersion and contact of the adsorbent and sample solution for I3C extraction,which improved the extraction efficiency and mass transfer rate of the adsorbent,thereby decreasing the sample preparation time.Therefore,10 min of ultrasonication was sufficient to afford the rapid extraction of I3C in the sample solution.

        Fig.3.Optimization of(A)adsorbent amount;(B)ultrasonication time;(C)types of washing solvents(1:water;2:phosphate buffer solution(pH 7.4,0.1 M);3:methanol-water(1:9,V/V);4:methanol-water(3:7,V/V);5:acetonitrile-water(1:9,V/V);6:acetonitrile-water(3:7,V/V));(D)types of elution solvents(1:methanol;2:ammonia water-methanol(1:9,V/V);3:acetonitrile;4:ammonia water-acetonitrile(1:9,V/V);5:acetonitrile-water(8:2,V/V);(E)volumes of ammonia water-acetonitrile(1:9,V/V).(F)Chromatograms of broccoli sample before(a)and after(b)pretreatment with PDES-GO-UADFE method,and I3C standard solution(c).

        3.2.3.Optimizationofwashingsolvent

        In the washing step,six types of washing solvents(0.5 mL)were investigated(Fig.3C).When the washing solvent contained 10%-30% of organic solvent,a high loss of I3C(6.4%-28.2%)occurred.When water or PBS was used as the washing solvent,only a small I3C loss(≤4.1%)was observed.However,water typically dissolves a certain amount of CO2,affording a pH of<7,which is not conducive to the stability of I3C;therefore,PBS(0.5 mL)was selected as the washing solvent.

        3.2.4.Optimizationofelutionsolvent

        A series of elution solvents was then examined to achieve complete I3C elution.As shown in Fig.3D,the addition of a small amount(10%)of ammonia water to methanol or acetonitrile improved I3C elution because ammonia water can destroy the interactions between I3C and the adsorbent,which include hydrogen bonding and electrostatic adsorption.After the addition of ammonia water,the basic elution solvent was also beneficial for I3C stability.Upon further investigation of the elution solvent volume(Fig.3E),it was determined that 1.5 mL of ammonia wateracetonitrile(1:9,V/V)solution was sufficient to elute I3C from the adsorbent with high recovery(>90%).

        3.3.Method validation

        The PDES-GO-UADFE-HPLC method for I3C analysis was validated by determining the detection limit(LOD),quantification limit(LOQ),linearity,accuracy,and precision.The LOD and LOQ values were 0.0250 and 0.0840μg/g,respectively,which were determined based on the signal-to-noise(S/N)ratios of 3 and 10,respectively.The calibration curves were constructed by plotting the chromatographic peak area versus a series of analyte concentrations.Sufficient linearity in the range of 0.0840-48.0μg/g was obtained with a coefficient of determination(r)of 0.9998.Accuracy was assessed at three spiking levels of 1.00,10.0,and 30.0μg/g(n=3),and the recoveries were 84.3%-96.4% with relative standard deviations(RSDs)of≤5.0%.The method precision was indicated in terms of repeatability and reproducibility,which were calculated by the extraction and analysis of I3C from the broccoli matrix at the spiking level of 10.0μg/g using the methods in one(intraday,n=6)and three consecutive days(interday,n=3),respectively.The intraday and interday precisions were expressed as RSDs,which had values≤5.6%.A wide range(0.0840-48.0μg/g),good linearity(r=0.9998),high sensitivity(LOD:0.0250μg/g),high accuracy(84.3%-96.4%),and precision(RSD≤5.6%)of the method showed the application potential of the proposed method for I3C analysis in the broccoli samples.

        3.4.Analysis of real samples

        The developed method was applied to the extraction and analysis of I3C in broccoli.The broccoli solution obtained after autolysis and preliminary solvent extraction contained many interfering substances,which affected the accurate analysis of I3C(Fig.3F(a)).After pretreatment using the PDES-GO-UADFE method,the interference peaks near I3C disappeared and the I3C peak could be clearly identified(Fig.3F(b)),indicating the purification effect of the proposed method.

        The developed method was applied for I3C analysis in different varieties of broccoli samples,and the results are shown in Table 2.The I3C contents detected in different varieties of broccoli were 6.58 μg/g(sample#1)and 9.11 μg/g(sample#4).The results indicated that broccoli sample#4 had a higher I3C content than broccoli sample#1.This indicates that the proposed method can be used to identify desirable broccoli varieties with high I3C contents.

        Table 2Application of the proposed PDES-GO-UADFE-HPLC method for I3C analysis in different varieties of broccoli samples.

        Furthermore,the developed method was applied for the analysis of I3C in different edible parts of the broccoli samples.The edible parts of broccoli mainly include the flower bud and the flower stem;the I3C content in these two parts was determined and compared.As shown in Table 2(samples#2 and#3),the result indicates that the I3C content in the flower bud(5.43-7.21μg/g)was significantly higher than that in the flower stem(2.29-2.60μg/g),providing detailed distribution data of I3C in different parts of broccoli.

        The developed method was also applied to investigate the effect of high temperature(boiling)on the I3C content in broccoli.The I3C content in broccoli sample#4 was determined before and after processing with boiling water for 10 min,which decreased significantly from 9.11 to 2.68μg/g.Therefore,it is recommended to reduce the time required for high-temperature processing,which is beneficial for retaining more I3C in broccoli.

        3.5.Comparison of UADFE with other sample preparation methods

        The UADFE method was compared to three common sample preparation methods,namely,FSPE,DSPE,and SPE,to show the advantages of the proposed method.The general operation procedures of FSPE,DSPE,and SPE were based on the literature references[37-39].In comparison,the types and volumes of solvents used in each step were identical,and only the extraction method was changed for suitable comparison.The results are shown in Table 3.The FSPE procedure was rapid,but the recovery(21.4%-25.2%)was the lowest among the four methods.Owing to the small adsorbent amount(3.00 mg),the adsorbent could not be easily and uniformly spread in the filter;therefore,the adsorbent and the sample solution did not sufficiently contact the FSPE.In DSPE,PDES-GO is hydrophilic and light,which was easily lost upon discarding the sample solution or washing solvent after centrifugation,resulting in a low I3C recovery(51.8%-60.2%).SPE,the most common sample preparation method,had a high recovery of I3C.However,the adsorbent consumption was high(500 mg of C18),and the method was more time-consuming(42-52 min)than the other methods.In the UADFE method,the DSPE mode was adopted to afford rapid and sufficient adsorption and desorption,and the FSPE mode was employed to achieve rapid separation of the adsorbent and solution with a high recovery(84.3%-96.4%)and low adsorbent consumption(3.00 mg).Therefore,the UADFE method was better than the other three methods in terms of cost,recovery,and time consumption,indicating that it is more suitable for I3C analysis in the broccoli samples and can meet the requirements of speed,accuracy,and low cost for pharmaceutical analysis.

        Table 3Comparison of UADFE method with other sample preparation methods.

        4.Conclusions

        In conclusion,a rapid,accurate,and low-cost PDES-GO-UADFE method was developed for the first time,which was coupled to HPLC for the analysis of I3C in broccoli.In the PDES-GO-UADFE method,PDES-GO with multiple adsorption interactions,high adsorption capacity,and fast mass transfer was synthesized to afford the rapid extraction of I3C,which was superior to that obtained using the commercial adsorbents.The rapid UADFE method with low adsorbent consumption was designed by combining the advantages of DSPE and FSPE to realize rapid adsorption,desorption,and separation.The proposed method was successfully applied to the extraction and analysis of I3C in broccoli,thereby affording a new analytical method for the analysis of anti-cancer compounds in the field of pharmaceutical analysis.

        CRediT author statement

        Yanan Yuan:Validation,Writing-Original draft preparation;HuanhuanChen:Investigation;YehongHan:Visualization,Methodology;Fengxia Qiao:Data curation;Hongyuan Yan:Supervision,Writing-Reviewing and Editing.

        Declaration of competing interest

        The authors declare that there are no conflicts of interest.

        Acknowledgments

        This work is supported by the National Natural Science Foundation of China(Grant Nos.:82073605 and 81803287),the Natural Science Foundation of Hebei Province(Grant Nos.:B2018201270 and H2021201002),the Talent Engineering Training Foundation of Hebei Province(Grant No.:A201802002),the Research Project of Hebei Provincial Administration of Traditional Chinese Medicine(Grant No.:2021174),Post-graduate's Innovation Fund Project of Hebei University(Grant No.:HBU2021ss006),and the Outstanding Doctoral Cultivation Project of Hebei University(Grant No.:YB201703).

        AV在线毛片| 午夜三级a三级三点| 极品熟妇大蝴蝶20p| 久久久久久人妻一区精品| 一区二区三区观看视频在线| 日韩久久无码免费毛片软件| 又白又嫩毛又多15p| 久久尤物av天堂日日综合| 日韩人妻美乳中文字幕在线| 欧美性猛交99久久久久99按摩| 波多野结衣乳巨码无在线| 性一交一乱一伦一视频一二三区| 亚洲av手机在线一区| 欧美丰满少妇xxxx性| a级国产乱理论片在线观看| 国产成人cao在线| 蜜桃夜夜爽天天爽三区麻豆av| 亚洲av成人片色在线观看高潮| 欧美老妇与禽交| 亚洲免费不卡av网站| 在线播放草猛免费视频| 国产精品狼人久久久久影院| 久久久久久久98亚洲精品| 国产一区二区三区在线观看免费版| 久久久久高潮综合影院| 人妻忍着娇喘被中进中出视频| 99精品国产闺蜜国产在线闺蜜| 日本免费看一区二区三区| 人妻洗澡被强公日日澡电影| 亚洲 欧美 国产 日韩 精品| 中文熟女av一区二区| 不卡一区二区视频日本| 国产特级毛片aaaaaaa高清 | 无码一区二区三区亚洲人妻 | 麻豆国产原创视频在线播放| 久久精品国产亚洲AV无码不| 熟女少妇av一区二区三区 | 1精品啪国产在线观看免费牛牛| 精品黑人一区二区三区久久hd| 亚洲av综合av一区| 亚州综合激情另类久久久|