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

    Trace determination and characterization of ginsenosides in rat plasma through magnetic dispersive solid-phase extraction based on core-shell polydopamine-coated magnetic nanoparticles

    2020-02-28 11:11:18NingningZhoShuLiuJunpengXingZifengPiFengruiSongZhiqingLiu
    Journal of Pharmaceutical Analysis 2020年1期

    Ningning Zho ,Shu Liu ,Junpeng Xing ,Zifeng Pi ,Fengrui Song ,Zhiqing Liu ,c

    a National Center of Mass Spectrometry in Changchun and Jilin Province Key Laboratory of Chinese Medicine Chemistry,Changchun Institute of Applied Chemistry,Chinese Academy of Sciences,Changchun 130022,China

    b University of Science and Technology of China,Hefei 230029,China

    c State Key Laboratory of Electroanalytiacl Chemistry,Changchun Institute of Applied Chemistry,Chinese Academy of Sciences,Changchun 130022,China

    Keywords:

    ABSTRACT

    1.Introduction

    Many traditional Chinese medicines(TCMs),which are homologous to food and medicine,have been widely used in daily food and disease treatment[1].However,their pharmacodynamic material basis or food nutrients remain unclear.The matrix effect of complex biological samples seriously limits the research on active ingredients and metabolites of TCMs[2,3].Thus,trace bioactive constituents and metabolites from biological samples with highly abundant endogenous interference should be enriched to better investigate the pharmacodynamic material basis.At present,the ultra-high-performance liquid chromatography coupled with highresolution mass spectrometry(UPLC-MS)has been widely used to enrich and detect the complex biological samples[4,5].In addition,new acquisition modes that can obtain good data of all precursors and fragment ions for the determination of targeted molecules are developed by the development of data-dependent acquisition(DDA)and data-independent acquisition(MSE)[6,7].Moreover,the analytical sensitivity and capacity are remarkably enhanced[8].Thus,these superiorities have provided an advanced analytical platform to gain comprehensive information for the separation and characterization of metabolites.

    Despite the presence of advanced analytical instrumental platforms,there is an urgent need to develop pretreatment methods in conjunction with UPLC-MS for better analysis.At present,traditional solvent extraction methods[9-13],liquid-liquid extraction(LLE)and solid-phase extraction(SPE)[14]are widely used to extract ginsenosides.But these methods require high temperature,long hours,and lots of efforts and solvent.And SPE method needs expensive commercial cartridges and fussy procedures.Although our previously proposed eco-friendly deep eutectic solvent(DES)[2]for extraction of ginsenosides showed high efficiency and solubility,sample concentration is still limited by negligible volatility of DES.The development of magnetic solid-phase extraction(MSPE)and magnetic dispersion solid-phase extraction(MDSPE)[15]has provided a new potential value for the efficient extraction of ginsenosides,but it is a challenge to synthesize functional adsorbents.Table S1 summarizes the advantages and drawbacks of extraction techniques for ginsenosides.

    Recently,magnetic nanoparticles(MNPs)have attracted considerable attention due to their outstanding properties,such as good dispersion,controllable core-shell structure,fast and efficient separation,sustainable usage and no need to pack columns.MNPs have been widely used as adsorbents for enriching small molecules[16,17].However,MNPs cannot form other interactions except for hydrogen bonds between most TCMs or food bioactive constituents due to the lack of specific functional groups.Moreover,the syntheses of materials require complex procedures and harmful solvents in the preparation process,and few are applied to complex biological samples.

    Dopamine(DA),a small molecule that mimics the adhesive proteins of mussels,has been confirmed to play a unique role in the interfacial adhesion of numerous material surfaces by selfpolymerization under sufficient mild conditions(weak alkaline medium)[18-20].For example,self-polymerization of DA has been used for the extraction of estrogenic compounds from cow,goat,sheep and human milk[21].Unfortunately,the reported core-shellbased polydopamine NPs(PDA NPs)still have some limitations.The first limitation is that PDA possesses strong interfacial adhesion,which results in blocking cavities of buried substance in molecular imprinted materials(MIPs).Moreover,few studies report the detection of targeted molecules from complex matrix samples.

    In this study,a novel strategy was developed to enrich and detect ginsenosides from complex biological samples by introducing UNIFITWlibraries assisted with Fe3O4@SiO2@PDA NPs as adsorbents of MDSPE.The strategy incorporated the advantages of the three approaches:(1)MDSPE with Fe3O4@SiO2@PDA NPs for fast and efficiently capturing trace bioactive constituents from biological samples without packing column;(2)UPLC-MS technology for fast and accurate detection of trace constituents;(3)Supplemental UNIFI libraries for fast,high-throughput and automatic match of data and identification of compounds.This strategy achieved a promote platform including pre-treatment,sample analysis based on MS and data analysis.For experiments,the conditions were statistically optimized by response surface methodology(RSM)with total ginsenosides to obtain the best extraction process.Subsequently,the strategy was applied for the fast detection and determination of metabolites from rat plasma,which eliminated the interference of several endogenous substances.Moreover,the synergetic strategy was accomplished in compliance to 9 of the 12 green chemistry principles with green analytical chemistry's(GAC)[22].The proposed method was proven to be a promising strategy to improve the enrichment capability and sensitivity of MS for tracing the targeted molecules in biological samples.

    2.Experimental

    2.1.Chemicals

    White ginseng(WG),total ginsenosides(extract of WG)and all the standards of ginsenosides Rb1,Rb2,Rb3,Rc,Rd,Re,Rg1,PPD,Rf,Rg2,Rg3,Rh1,Rh2 and compound K were obtained from the College of Pharmacy,Jilin University.Dioscin as an internal standard(IS)was purchased from National Institutes for Food and Drug Control(Beijing,China).HPLC-grade acetonitrile and methanol were acquired from Fisher Scientific(Loughborough,UK).Ultrapure water was prepared by a Milli-Q water(18.2 MΩ)purification system(Milford,MA,USA).Tris(hydroxyl methyl amino methane)(Tris),dopamine hydrochloride(DA·HCl),heparin,ammonium sulfate(NH4)2SO4,aminopropyl triethoxy silane(APTES),tetraethoxysilane(TEOS)and Fe3O4nanoparticles(20-30 nm)were obtained from Aladdin Industrial Corporation(Shanghai,China).HCl,acetic acid and methanol were provided by Sinopharm Chemical Reagent Co(Shenyang,China).A 5 mM Tris-HCl buffer solution(pH 7-9)was prepared by us.

    2.2.Instrumentations

    The morphology and Energy Dispersive X-Ray(EDX)spectroscopies of MNPs were obtained by the scanning electron microscopy(SEM,JSM-6460LV,Japan).Fourier transform infrared spectroscopy(FT-IR,IFS 66 V/S,Germany)was performed to identify the materials with KBr.Other characteristics were measured by vibrating sample magnetometry(VSM,MPMS XL,USA)at room temperature and a D8 advance X-ray diffractometer(XRD,D8,Germany)with Cu Kαradiation(λ=1.5406?).The UPLC-MS/MS analyses were performed by a Waters ACQUITY UPLC system(Waters Corp.,Milford,USA)coupled to a triple quadrupole(TQ)mass spectrometer equipped with the electrospray ionization(ESI)source in multiple reactions monitoring(MRM)mode.A Q-TOF SYNAPT G2-S High Definition Mass Spectrometer coupled to ESI(Waters Crop.,Manchester,UK)was adopted for detection of ginsenosides metabolites before and after treatment with Fe3O4@SiO2@PDA NPs.The C18 column(50 mm×2.1 mm,1.9μm;Thermo scientific)was employed in the entire experiments.MassLynxV4.1 and Marker-Lynx Application Manager(Waters Corporation,Milford,USA)and UNIFI 1.8(Waters,USA)were carried out to process the collected data.And all detailed steps of mass spectrometry analyses were described in Supplementary data(Methods section),including UPLC-MS/MS analysis,Q-TOF analysis and mass spectral data analysis with UNIFITW.

    2.3.Synthesis of the Fe3O4@SiO2,Fe3O4@SiO2@APTES and Fe3O4@SiO2@PDA NPs

    The Fe3O4@SiO2NPs were synthesized by the previous method[23].The Fe3O4@SiO2@APTES NPs were synthesized through the modified method[24].The synthesis strategy of Fe3O4@SiO2@PDA NPs based on the modified method[25]through one-pot synthesis is as follows:briefly,the Fe3O4NPs(100 mg)were dispersed into an ethanol solution(95% ethanol:water,v/v).Then,ammonium hydroxide was added into the mixture solvents(pH=8.5).Ultrasound 10 min later,0.2 mL of TEOS was added dropwise into the system for 6 h at 40℃.After that,90 mg of DA·HCl was added into the bottle,followed by stirring for 1 h.Subsequently,the obtained product was isolated by the external magnetic field and washed by water and ethanol.Eventually,the product was dried under vacuum at 40℃ for 24 h.

    2.4.Animals

    Six male rats,weighting 200±20 g,were obtained by Jilin Medical University(Jilin,China).After a week of adjustable feeding,they were administered WG(7 g/kg)once a day for 5 days.The rats were sacrificed 2 h after the last gavage,and the blood was collected and centrifuged.The supernatants were stored at-80℃.All animals’experiments and the procedures followed the Guide for the Care and the Use of Laboratory Animals of Jilin University.

    Scheme 1.Schematic presentation of the proposed strategy for capture and characterization of ginsenosides by Fe3O4@SiO2@PDA NPs based on UPLC-MS method.

    2.5.Enrichment of ginsenosides from total ginsenosides

    During MNPs screening and tailoring,5 mg of the MNPs were mixed with 1 mL of the total ginsenosides in water(500μg/mL).The mixture was vortexed and shaken for 2 h.Later,the MNPs were isolated by external magnetic field,followed by ultrasound-assisted cleaning with 1 mL of methanol-acetic acid(9:1,v/v)and then direct mass spectrometry analysis.

    The adsorption capacity of Fe3O4@SiO2@PDA NPs for ginsenosides was investigated at different concentrations of total ginsenosides while other conditions remained unchanged.Meanwhile,the dynamic adsorption capacity was evaluated by varied time from 10 min to 360 min.Desorption time was estimated at varied time from 10 min to 120 min.The adsorption amounts(Q,mg/g)of MNPs were calculated according to Eq(1)[26].

    where,Co(μg/mL),Ce(μg/mL),V(mL)and m(g)represent the initial and equilibrium concentration of ginsenosides,the volume of the total ginsenosides solutions and the mass of the MNPs,respectively.

    To investigate the maximum adsorption properties of Fe3O4@-SiO2@PDA NPs,the adsorption isotherm was calculated by using a series of total ginsenosides solutions with different concentrations.The kinetic tests were calculated by varied different adsorption time.The Langmuir and Freundlich isotherm models,and the pseudo-first-order and pseudo-second-order models were adopted to evaluate the equilibrium and kinetic data,respectively[27].

    2.6.Experimental design and statistical analysis for total ginsenosides

    The screening experiments,including synthetic species of magnetic materials,content of dopamine,pH of the solvent system,concentrations of the total ginsenosides,enrichment time and desorption time were carried out.According to the above obtained result of single factor experiments,the appropriate variable range in extraction process(shown in Table S2)was employed for experimental design using RSM.The Box-behnken design(BBD)combined with RSM was performed as follows:pH(7-9),concentrations of the total ginsenosides(100-700μg/mL),enrichment time(30-80 min),and elution time(10-80 min)at levels(-1,0,1).Here,the reason for choosing such pH range was that the extraction efficiency was better at weak alkalinity for ginsenosides[15].The extraction recoveries of ginsenosides in 29 experiments(three replicates)were estimated by BBD in Table S3.

    Design-Expert 8.0.6 software[28]was operated to evaluate the statistical significance in design experiment.Based on analysis of variance(ANOVA),the quadratic multiple regression models were built to investigate model quality according to the square of correlation coefficient(R2)and the lack of fit which evaluated the failure of the model.

    2.7.Real samples extraction procedure

    Fig.1.SEM images of Fe3O4,Fe3O4@SiO2,Fe3O4@SiO2@APTES and Fe3O4@SiO2@PDA NPs.

    Fig.2.The EDX spectrum and element compositions of Fe3O4@SiO2@PDA NPs.

    The real plasma samples after lavage WG were analyzed by our synergetic strategy in Scheme 1.The 200μL plasma samples were divided into conventional methanol group and Fe3O4@SiO2@PDA NPs enrichment group(six replicates).First,all plasma samples were dealt with 3 folds methanol.After centrifugation,the supernatants were dried under nitrogen in methanol group.However,in the PDA group,5 mg Fe3O4@SiO2@PDA NPs were added into 1 mL Tris-HCl buffer system containing the supernatants of the same volume(final pH=8.78).After vortexing for 2 min,the system was shaken for 80 min.Later,the obtained supernatant was removed by external magnetic field.And the MNPs were eluted with 1 mL of methanol-acetic acid(9:1,v/v)by ultrasound for 68.25 min.Finally,the eluted solvents were also dried under nitrogen.All both dried samples were re-dissolved with 50μL methanol for TOF analysis.

    2.8.Method validation and measurement uncertainty

    A series of working solutions of Rb1,Rb2,Rb3,Rc,Rd,Re,Rg1,PPD,Rf,Rg2,Rg3,Rh1,Rh2 and compound K were obtained by stepwise dilution of the stock solutions in methanol.The calibration was created in the form of y=ax+b.The limits of detection(LODs)and limits of quantification(LOQs)were recognized at a signal-to-noise ratio(S/N)of 3 and 10 times.For quality control(QC)samples,fourteen ginsenosides were operated at low,medium,and high concentrations in WG powder and plasma,respectively.The IS was prepared at a concentration of 0.1μg/mL in each sample.Intra-and inter-day accuracy and precision were assessed by QC samples at concentration levels of low,middle and high for 14 ginsenosides in 3 days.The relative standard deviation(RSD)and relative error(RE)represented the precision and accuracy,respectively.To investigate the matrix effect,the analyses of total ginsenosides and drug-free plasma samples spiked with the 14 analytes were carried out.To further estimate the extraction recovery,the samples of WG powder and drug-free plasma samples spiked with the 14 analytes were directly analyzed using our proposed method.

    Proper quality assurance/quality control(QA/QC)procedures are vital for the measurement results,which provide reliable analytical information.According to Konieczka and Namie′snik[29],amount of sample,analyte concentration,calibration step,repeatability of results and recovery determination(trueness)of the uncertainty components were estimated.The uncertainty of results was calculated based on Karimi et al.[30]and Borecka et al.[31].Expanded uncertainty was determined by the Eq(2).

    The meaning of each element in the Eq(2)is as follows:U-expanded uncertainty(%),k-coverage factor(k=2 for P=95%),ur(sample)-relative standard uncertainty of sample,ur(cal)-relative standard uncertainty of calibration step,ur(true)-relative standard uncertainty of trueness,ur(rep)-relative standard uncertainty of repeatability,and ur(MDL)-relative standard uncertainty of MDL determination.Meanwhile,ur(cal),ur(true),ur(rep)and ur(MDL)were calculated by the following Eqs(3)-(6),respectively.

    The sensitivity of the method is expressed by method detection limit(MDL)and method quantitation limit(MQL)for the entire procedure.And the MDL and MQL were calculated by the Eq(7)and Eq(8).

    Fig.3.(A)X-Ray diffraction patterns of Fe3O4,Fe3O4@SiO2,Fe3O4@SiO2@APTES and Fe3O4@SiO2@PDA NPs.(B)FT-IR spectra of all the NPs.(C)VSM magnetization curves of all the NPs.(D)VSM magnetization curves of synthetic Fe3O4@SiO2@PDA NPs at different reaction times.

    3.Results and discussions

    where SDxyis the standard deviation of the calibration curve,p is the repetition number of measurements for a given sample,m is the total number of standard samples used for plotting the calibration curve,Xsampleis the concentration of sample,Xmis the mean value of all concentrations of a standard solution for which the measurement was made to plot a calibration curve,Xiis the calculated concentration,b is the slope of the calibration curve,RSDRis the RSD derived from the recovery,RSDresultsis the RSD in real samples and Cmis the mean concentration of analyte.The analytical result is expressed as mean±expanded uncertainty.

    3.1.Characterization of the magnetic nanoparticles

    3.1.1.SEM and EDX characterizations

    All MNPs were characterized through SEM.As shown in Fig.1,the diameters of Fe3O4,Fe3O4@SiO2,Fe3O4@SiO2@PDA and Fe3O4@SiO2@APTES NPs were about 25 nm,130 nm,180 nm,and 220 nm,respectively.The results revealed that the size of particles increased with the layer-by-layer modification on the Fe3O4core,and the order is expressed as follows:Fe3O4NPs<Fe3O4@SiO2 NPs<Fe3O4@SiO2@PDA NPs<Fe3O4@SiO2@APTES NPs.In addition,the enrichment capability of different MNPs for ginsenosides was also estimated through MS analysis(MRM conditions,calibration curves and analytical results are shown in Tables S4 and S5 and Fig.S1).The order of enrichment capability was in consistent with the order of MNPs size.Although the specific surface areas of Fe3O4and Fe3O4@SiO2NPs were larger than the other two,the enrichment capabilities of bare MNPs and hydroxyl-containing Fe3O4@-SiO2NPs were weaker than the other two amino-containing MNPs.It is obvious that the effect of amino groups was greater than hydroxyl groups,and the effect of functional groups was relatively more important than specific surface area in this experiment.Moreover,the higher enrichment capability for Fe3O4@SiO2@PDA NPs compared with Fe3O4@SiO2@APTES NPs was obtained due to two reasons:(1)PDA possessed multiple recognition sites,such as amino and hydroxyl groups;(2)Fe3O4@SiO2@PDA NPs with small size had relatively high surface areas.

    Fig.4.RSM response surface plots of the model for extraction of ginsenosides.

    Furthermore,the EDX spectrum and main elemental composition analysis for Fe3O4@SiO2@PDA NPs are presented in Fig.2.The contents of N elements were further validated the existence of PDA.All the results suggested that PDA was successfully grafted to the Fe3O4@SiO2NPs.

    3.1.2.XRD characterization

    The XRD patterns are shown in Fig.3A.It was obviously observed that six major characteristic diffraction peaks of Fe3O4NPs appeared at 2θ=30.2°(220),35.5°(311),43.2°(400),53.4°(422),57.2°(511),and 62.9°(440).Although the Fe3O4NPs were gradually modified,the spectra of modified MNPs were similar to the XRD spectrum of Fe3O4NPs at the six major peaks,revealing that the modification had not been destroyed on the crystal structure of Fe3O4particles.The peaks attributed to SiO2and PDA were not observed,because SiO2or PDA coating on the surface of Fe3O4mainly existed as amorphous forms[32].

    3.1.3.FT-IR characterization

    The Fe3O4,Fe3O4@SiO2,Fe3O4@SiO2@APTES and Fe3O4@-SiO2@PDA NPs were further validated by FT-IR(Fig.3B).The typical peaks located at approximately 576 cm-1and 1090 cm-1were attributed to the stretching vibration of Fe-O and Si-O bonds,respectively[33].The peak at approximately 1626 cm-1was assigned to the N-H vibration in APTES[34].The peaks for PDA at approximately 1621 cm-1and 1508 cm-1 were assigned to C=C stretching vibrations in the aromatic ring and N-H bending vibration,respectively[35].Moreover,in the entire spectrum,a broad band at approximately 3420 cm-1corresponded to the-NH or-OH stretching vibration[32].These results reaffirmed the formation of PDA shell on Fe3O4@SiO2@PDA NPs.

    3.1.4.VSM analysis

    As observed in Fig.3C,the saturation magnetizations of Fe3O4and Fe3O4@SiO2NPs were 53.08 emu/g and 42.40 emu/g,which reduced to 33.07 and 37.71 emu/g after coating with APTES and PDA,respectively.Although the saturation magnetization declined clearly,the MNPs still met the experimental requirements.The saturation magnetization of Fe3O4@SiO2@PDA NPs was stronger than that of Fe3O4@SiO2@APTES,which proved its efficient separation in the adsorption of ginsenosides.As shown in Fig.3D,saturation magnetization reduced with the increase of polymerization time.And the intensity was the highest when the polymerization time was 60 min.This optical condition demonstrated the existence of critical values for polymerization time.

    3.2.Effect of content of DA

    The content of DA was estimated for the synthesis of Fe3O4@-SiO2@PDA NPs to improve their enrichment capability(Q).As shown in Fig.S2A,Q was the strongest when DA·HCl was 90 mg.Excess amounts of DA·HCl caused the decrease of extraction capacity for ginsenosides.The results showed that superfluous DA·HCl promoted self-aggregation of DA,or the shell of MNPs thickened and resulted in the decrease of specific surface area[36].

    3.3.Effect of reaction time

    The polymerization time was also estimated for the synthesis of Fe3O4@SiO2@PDA NPs.As displayed in Fig.S2B,Fe3O4@SiO2@PDA NPs exhibited the highest extraction efficiencies of ginsenosides when the polymerization time was 60 min,which resulted from the increasement of recognition sites,a stepwise suitable thickness andthe above-mentioned strong saturation magnetization.However,the enrichment effect became relatively low with long polymerization time.The low enrichment effect was caused by the reduction of magnetism and specific surface area due to the thickened shell of the MNPs as time went by,and this finding was consistent with the characterization data in Fig.3D.

    Table 1Calculated values of relative standard uncertainties,combined standard uncertainties and expanded uncertainties for the determination of 14 analytes in real plasma.

    3.4.Optimization for ginsenosides extraction using RSM

    Single-factor experiments were implemented to obtain the suitable value ranges for statistical analysis(Fig.S3).In Fig.S3A,when pH ranged from 7 to 9,the values of adsorption capacities increased.It is confirmed that the adsorption capacities for ginsenosides were better at weak alkalinity for the following two reasons:(1)ginsenosides formed more hydrogen bonds with the multiple recognition sites of PDA at weak alkalinity;(2)hydrophobic interaction was better at weak alkalinity.For the isotherm test,Langmuir model(R2=0.94)was more successful and suitable than Freundlich model(R2=0.89)at 10-1000μg/mL(Total ginsenosides solutions),and the fitting Q was greater by Langmuir model(Fig.S3B).For the adsorption kinetic test(Fig.S3C),the Q was coincided with the pseudo-first-order model,and the pseudo-firstorder model(R2=0.99)was more suitable than pseudo-secondorder model(R2=0.86).As shown in Fig.S3D,the desorption capacity remained unchanged after about 60 min.All those results suggested that the single-factor experiment could provide certain ranges for RSM analysis.Furthermore,Design-Expert 8.0.6 was used to evaluate the statistical relevance in the suitable range(Fig.4).Second-order polynomial quadratic equations between extraction yields(Y)and coded variables(A,B,C and D)were calculated and expressed as follows:

    The data of ANOVA analysis and 3D response plots indicated the relevant variables for Y of ginsenosides and the interaction effect between the variables(Table S6).In the model,extraction concentrations(B),enrichment time(C),and elution time(D)exhibited significant difference at p<0.05.pH displayed an insignificant efficiency(p=0.0856).The results indicated that B,C and D presented remarkable effects,and no obvious interactions with each other were observed during the enrichment process.A predicted optimal Y of 29.63 mg/g was obtained when the extraction system pH was 8.78,concentration was 576.32μg/mL,enrichment time was 80 min,and elution time was 68.25 min.Meanwhile,the highest Yof 28.32±0.82 mg/g was achieved and the result reflected the repeatability of the model obtained.Generally,a good methodological result was obtained.In addition,the recovery exceeded 80%after recycling Fe3O4@SiO2@PDA NPs 6 times(data are shown in Fig.S4),and this finding revealed that the MNPs were efficient and sustainable for the enrichment of ginsenosides.

    3.5.Application to plasma samples and metabolite identification

    From BBD experiments,the best extraction process was optimized and the recovery was extremely close to predicted values.The LODs and LOQs of analytes were obtained through the verification of methodology,and good linearities were higher than 0.99 and wide linear ranges were obtained.In Table S7,the accuracy and precision of all the analytes were evaluated.All values conformed to methodological standard.In Table S8,mean accuracies of extraction recoveries and matrix effects of all the analytes ranged from 85%to 115% at three QC levels(n=6).The high extraction recovery and less matrix effect of 14 analytes indicated that the Fe3O4@SiO2@PDA NPs extraction method was viable.The method was feasible when ginsenosides were enriched from blank plasma spiked with standard mixtures of ginsenosides(in methanol)by the Fe3O4@-SiO2@PDA NPs.In addition,the QA/QC procedures were successfully applied to the measurement and analytical results.Finally,uncertainties of all measurement and analytical results were calculated and the final expanded uncertainty results are shown in Table 1.To further evaluate whether the proposed method has high sensitivity for Q-TOF MS analysis in real complex samples,the proposed approach was implemented to the rat plasmaadministered WG,and UNIFI?libraries were supplemented to automatically assist the separation and characterization of compounds.The base peak ion(BPI)chromatograms of conventional methanol group and Fe3O4@SiO2@PDA NPs enrichment group are summarized in Fig.5.As shown in BPI chromatograms,ion suppressions from plasma had a severe effect on the mass spectrum of ginsenosides,although the signal intensity was slightly high for the conventional methanol method.

    Fig.5.The BPI chromatograms acquired by conventional methanol extraction method(A)and Fe3O4@SiO2@PDA NPs enrichment method(B).

    Compound 3(recognized as ginsenosides Re)was taken as examples to illustrate the structural identification of ginsenosides based on mass spectral fragmentation pathways[15]and UNIFI?libraries.As shown in Fig.6,characteristic ions of compound 3 were m/z 945.5379[M-H]-and 991.5471[M+HCOOH-H]-in the negative ion mode.For the high-energy function,the main fragmentation ions were at m/z 799.4739,783.4918,637.4239,475.3763,391.2855 and 161.0451,which corresponded to[M-Rha]-,[M-Glc]-,[M-Glc-Rha]-,[M-Glc-Rha-Glc]-,[M-Glc-Rha-Glc-C6H12]-and[C6H9O5]-,respectively.These fragmentation ions were identical to those of Re in the UNIFI?libraries.Moreover,the liquid retention time of compound 3 was the same as the Re in UNIFI?libraries.So,compound 3 was identified to ginsenoside Re through constant comparison and confirmation.All the prototypes and metabolites in vivo were matched through the above characterization and confirmation steps,the results are listed in Table 2.The combination of MSEand UNIFI showed a huge potential in the separation and determination of compounds.

    3.6.Comparison of the proposed method with other methods

    Table S9 shows the distribution information of each ginsenoside in the biological samples detected by different enrichment methods.Surprisingly,23 compounds were detected and determined using our proposed method,whereas only 8 ginsenosides were determined by the methanol method.For the proposed enrichment method,the ion signal enhancement effect of targeted compounds was obvious.Moreover,the proposed strategy can rapidly separate and identify many targets in the spectrum,which avoided several endogenous substances.This finding proved that Fe3O4@SiO2@PDA NPs had a huge potential for the enrichment of ginsenosides from a wide range of samples.

    Fig.6.The major fragmentation pathways of ginsenoside Re.

    Table 2The major prototypes and metabolites of enrichment and identification using Fe3O4@SiO2@PDA NPs extraction method in the biological samples in the negative ion mode.

    Table 3Comparison of adsorption capability of pre-concentration methods for determination of ginsenosides.

    As shown in Table 3,the proposed strategy was compared with previous extraction methods for the determination of ginsenosides[2,9-15].The proposed strategy had several advantages over normal methods.Compared with conventional solvent extraction and LLE,MDSPE did not require centrifugation and can perform rapid separation.Moreover,the adsorption capability of Fe3O4@-SiO2@PDA NPs for ginsenosides was the highest(28.32 mg/g in total ginseng and 24.56 mg/g in plasma).All these results demonstrated that the proposed method was a fast,sensitive and reusable technique that could be used for the enrichment and determination of ginsenosides in biological samples.

    As for green aspects,the above-mentioned techniques for the enrichment of ginsenosides are comprehensively assessed in Table S10.The proposed strategy was accomplished in compliance to 9 of the 12 GAC principles,which stood for reduction of the use of chemical substances and the recycle of them.All other methods met no more requirments of the principles than the one used in this work.

    4.Conclusions

    In summary,Fe3O4@SiO2@PDA NPs with multiple recognition sites in PDA were synthesized through a green one-pot procedure.It is easy to synthesize and there is no need to pack columns for sample analysis.Importantly,the materials could recycle at least 6 times,conforming to the GAC principles.Moreover,the fast separation and characterization of 23 ginsenoside prototypes and metabolites from plasma were achieved by combining UPLC-MS with UNIFI libraries.In comparison,conventional methanol method can only detect 8 ginsenosides from the same plasma samples.This new synergetic strategy for the accurate characterization of trace ginsenosides using an easy,fast,economical,automated and ecofriendly way shows potential applications.The integration of facile and one-pot materials,advanced analytical platforms and automatically high-throughput post data-processing softwares will promote the research on the pharmacodynamic material basis of related food and TCMs.

    久久精品国产亚洲av涩爱| 成人亚洲精品一区在线观看| 成年动漫av网址| kizo精华| 18禁国产床啪视频网站| 久热这里只有精品99| 黄片小视频在线播放| a级片在线免费高清观看视频| 亚洲人成网站在线观看播放| 国产又色又爽无遮挡免| 久久婷婷青草| 国产精品香港三级国产av潘金莲 | 日韩不卡一区二区三区视频在线| 91aial.com中文字幕在线观看| 哪个播放器可以免费观看大片| 免费人妻精品一区二区三区视频| 91国产中文字幕| 波多野结衣一区麻豆| 免费播放大片免费观看视频在线观看| 亚洲色图 男人天堂 中文字幕| 日韩av在线免费看完整版不卡| 国产极品天堂在线| 最新的欧美精品一区二区| 亚洲欧美一区二区三区黑人 | 午夜免费鲁丝| 国产精品久久久久成人av| 卡戴珊不雅视频在线播放| 卡戴珊不雅视频在线播放| 精品国产国语对白av| 在线观看美女被高潮喷水网站| 国产亚洲午夜精品一区二区久久| 两个人看的免费小视频| 人人澡人人妻人| 欧美中文综合在线视频| 高清不卡的av网站| 午夜激情av网站| 中文字幕精品免费在线观看视频| 国产精品久久久久久av不卡| 色婷婷久久久亚洲欧美| 美女福利国产在线| av国产精品久久久久影院| 9191精品国产免费久久| 色播在线永久视频| 精品卡一卡二卡四卡免费| 99九九在线精品视频| 韩国av在线不卡| 国产亚洲一区二区精品| 999精品在线视频| 国产淫语在线视频| 日韩精品免费视频一区二区三区| 宅男免费午夜| 精品视频人人做人人爽| 午夜免费男女啪啪视频观看| 人人妻人人爽人人添夜夜欢视频| 黄片小视频在线播放| 成人亚洲精品一区在线观看| 国产一区二区三区综合在线观看| 精品少妇黑人巨大在线播放| 熟女电影av网| 日产精品乱码卡一卡2卡三| 一本色道久久久久久精品综合| 青春草亚洲视频在线观看| 精品少妇内射三级| 国产黄色视频一区二区在线观看| 日韩制服丝袜自拍偷拍| 精品视频人人做人人爽| 黑丝袜美女国产一区| 国产黄色免费在线视频| 亚洲伊人久久精品综合| 成人影院久久| 91久久精品国产一区二区三区| 99久久综合免费| 97精品久久久久久久久久精品| tube8黄色片| 18禁国产床啪视频网站| 亚洲精品一二三| 久久精品国产亚洲av涩爱| 熟女电影av网| 汤姆久久久久久久影院中文字幕| 日本av免费视频播放| 看非洲黑人一级黄片| www.精华液| 校园人妻丝袜中文字幕| 国产人伦9x9x在线观看 | 十八禁高潮呻吟视频| 亚洲人成网站在线观看播放| 亚洲av电影在线进入| 纵有疾风起免费观看全集完整版| 1024视频免费在线观看| 日本91视频免费播放| 嫩草影院入口| 伊人久久大香线蕉亚洲五| 性高湖久久久久久久久免费观看| 免费观看a级毛片全部| 国产日韩一区二区三区精品不卡| 18禁动态无遮挡网站| 国产在线免费精品| 免费不卡的大黄色大毛片视频在线观看| 欧美人与性动交α欧美软件| 18禁国产床啪视频网站| 亚洲天堂av无毛| 午夜福利视频精品| 精品国产一区二区三区四区第35| 成人毛片60女人毛片免费| www.精华液| 一边摸一边做爽爽视频免费| 人人妻人人澡人人爽人人夜夜| 色婷婷av一区二区三区视频| 欧美日韩精品网址| 国产一区二区三区综合在线观看| 亚洲欧美成人综合另类久久久| 精品卡一卡二卡四卡免费| 免费黄频网站在线观看国产| 男的添女的下面高潮视频| 久久久精品94久久精品| 女的被弄到高潮叫床怎么办| 天堂中文最新版在线下载| 自拍欧美九色日韩亚洲蝌蚪91| a级毛片黄视频| 亚洲视频免费观看视频| 国产欧美日韩综合在线一区二区| 亚洲婷婷狠狠爱综合网| 午夜福利乱码中文字幕| 99香蕉大伊视频| 国产精品久久久久久av不卡| 精品人妻在线不人妻| 久久精品亚洲av国产电影网| av国产久精品久网站免费入址| 国产免费又黄又爽又色| 丝袜人妻中文字幕| 国产成人av激情在线播放| 老熟女久久久| av网站免费在线观看视频| 在线观看国产h片| 国产精品久久久久久精品古装| 国产欧美亚洲国产| 大片电影免费在线观看免费| 国产亚洲一区二区精品| av.在线天堂| 国产精品 欧美亚洲| 韩国av在线不卡| 久久鲁丝午夜福利片| 桃花免费在线播放| 热re99久久精品国产66热6| 久久久久国产网址| 各种免费的搞黄视频| 久久久久精品性色| 高清视频免费观看一区二区| 最近最新中文字幕免费大全7| 国产片特级美女逼逼视频| 熟女av电影| 人人妻人人澡人人看| 我要看黄色一级片免费的| 精品人妻在线不人妻| 哪个播放器可以免费观看大片| 亚洲欧美一区二区三区久久| 国产淫语在线视频| 亚洲欧美成人综合另类久久久| 99久久精品国产国产毛片| 午夜免费观看性视频| 制服诱惑二区| 亚洲av免费高清在线观看| 国产精品国产av在线观看| 一区在线观看完整版| 一级毛片我不卡| 国产淫语在线视频| 亚洲国产av新网站| 久久人人爽人人片av| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 91精品伊人久久大香线蕉| 午夜影院在线不卡| 婷婷色综合大香蕉| 国产精品熟女久久久久浪| 熟女av电影| 亚洲国产av影院在线观看| a级毛片黄视频| 亚洲在久久综合| 日韩大片免费观看网站| 晚上一个人看的免费电影| 精品国产一区二区三区久久久樱花| 夫妻午夜视频| 91成人精品电影| 91午夜精品亚洲一区二区三区| 高清不卡的av网站| 日韩在线高清观看一区二区三区| 最近最新中文字幕免费大全7| 麻豆乱淫一区二区| av.在线天堂| 两个人看的免费小视频| 午夜福利视频在线观看免费| 免费人妻精品一区二区三区视频| 中文乱码字字幕精品一区二区三区| 尾随美女入室| 日本vs欧美在线观看视频| 国产精品免费视频内射| 国产有黄有色有爽视频| 精品少妇久久久久久888优播| 国产亚洲av片在线观看秒播厂| 日韩av不卡免费在线播放| 久久99热这里只频精品6学生| 午夜福利视频在线观看免费| 久久精品久久久久久噜噜老黄| av女优亚洲男人天堂| 99re6热这里在线精品视频| 成人国语在线视频| 在线亚洲精品国产二区图片欧美| 亚洲av中文av极速乱| 国产精品久久久久久av不卡| 亚洲人成网站在线观看播放| av又黄又爽大尺度在线免费看| 久久久精品区二区三区| 精品人妻偷拍中文字幕| 大片免费播放器 马上看| 少妇被粗大的猛进出69影院| 国产一区亚洲一区在线观看| 午夜精品国产一区二区电影| 男女啪啪激烈高潮av片| 99九九在线精品视频| 日韩在线高清观看一区二区三区| av网站在线播放免费| 欧美 日韩 精品 国产| 午夜福利在线观看免费完整高清在| 亚洲情色 制服丝袜| 国产一区有黄有色的免费视频| 欧美日韩成人在线一区二区| 在线 av 中文字幕| 亚洲人成77777在线视频| 精品一区在线观看国产| 欧美精品高潮呻吟av久久| 成人手机av| 免费日韩欧美在线观看| 亚洲经典国产精华液单| 青青草视频在线视频观看| 免费少妇av软件| 国产 一区精品| 成年人午夜在线观看视频| 日韩,欧美,国产一区二区三区| 国产精品人妻久久久影院| 亚洲成人一二三区av| 边亲边吃奶的免费视频| videossex国产| 亚洲欧美成人精品一区二区| 欧美精品亚洲一区二区| 亚洲成av片中文字幕在线观看 | 亚洲精品中文字幕在线视频| 寂寞人妻少妇视频99o| xxxhd国产人妻xxx| 熟女av电影| 亚洲国产欧美在线一区| 亚洲av在线观看美女高潮| 国产黄色免费在线视频| 久久国内精品自在自线图片| 久久久国产精品麻豆| 伦理电影免费视频| 国产乱人偷精品视频| 欧美精品一区二区免费开放| 亚洲欧美日韩另类电影网站| 国产福利在线免费观看视频| 精品久久久精品久久久| 国产精品一区二区在线观看99| 亚洲精品国产色婷婷电影| av在线app专区| 一区在线观看完整版| 国产精品99久久99久久久不卡 | 亚洲成国产人片在线观看| 亚洲精品中文字幕在线视频| 亚洲精品,欧美精品| 亚洲四区av| videosex国产| xxxhd国产人妻xxx| 国产精品二区激情视频| 一边亲一边摸免费视频| av在线老鸭窝| 日本vs欧美在线观看视频| 国产一级毛片在线| 欧美日韩成人在线一区二区| 久久精品夜色国产| 成人毛片a级毛片在线播放| 亚洲成人手机| 国产成人一区二区在线| 亚洲精品国产色婷婷电影| 国产日韩欧美在线精品| 精品一区二区免费观看| 三上悠亚av全集在线观看| 亚洲一区二区三区欧美精品| 亚洲欧美精品综合一区二区三区 | 91在线精品国自产拍蜜月| 久久精品aⅴ一区二区三区四区 | 精品亚洲乱码少妇综合久久| 免费在线观看黄色视频的| 免费黄色在线免费观看| 亚洲av欧美aⅴ国产| 亚洲欧美一区二区三区久久| 久久久久久久久久久免费av| 日韩成人av中文字幕在线观看| 两个人免费观看高清视频| 成人国产av品久久久| 秋霞伦理黄片| 自拍欧美九色日韩亚洲蝌蚪91| 精品亚洲乱码少妇综合久久| 精品一区二区免费观看| 老熟女久久久| 欧美bdsm另类| 欧美 日韩 精品 国产| 免费人妻精品一区二区三区视频| 三级国产精品片| 电影成人av| 人妻人人澡人人爽人人| 欧美老熟妇乱子伦牲交| 成人国语在线视频| tube8黄色片| 国产黄频视频在线观看| 午夜91福利影院| 日本黄色日本黄色录像| 久久99热这里只频精品6学生| 亚洲欧美成人精品一区二区| 欧美国产精品一级二级三级| 久热这里只有精品99| xxx大片免费视频| 十八禁高潮呻吟视频| 成人国产麻豆网| 欧美日韩亚洲高清精品| 亚洲天堂av无毛| 一区二区三区激情视频| 91aial.com中文字幕在线观看| 亚洲精品成人av观看孕妇| 精品人妻偷拍中文字幕| 久久精品夜色国产| 亚洲精品久久久久久婷婷小说| 国产片特级美女逼逼视频| 日本-黄色视频高清免费观看| 亚洲国产欧美网| 女性生殖器流出的白浆| xxxhd国产人妻xxx| 日本黄色日本黄色录像| 一二三四中文在线观看免费高清| 高清视频免费观看一区二区| 黄色怎么调成土黄色| 两个人免费观看高清视频| 波多野结衣一区麻豆| 交换朋友夫妻互换小说| 捣出白浆h1v1| 中文天堂在线官网| 国产日韩一区二区三区精品不卡| 国产无遮挡羞羞视频在线观看| 伦精品一区二区三区| 五月伊人婷婷丁香| 午夜老司机福利剧场| 国产淫语在线视频| 成人手机av| 老司机亚洲免费影院| www.av在线官网国产| av天堂久久9| 国产野战对白在线观看| 国产精品偷伦视频观看了| 不卡av一区二区三区| 精品少妇黑人巨大在线播放| 爱豆传媒免费全集在线观看| 久久精品国产亚洲av天美| 边亲边吃奶的免费视频| 国产亚洲午夜精品一区二区久久| 久久久久久久大尺度免费视频| 男女下面插进去视频免费观看| 精品第一国产精品| 久热这里只有精品99| 国产麻豆69| 丝袜脚勾引网站| 国产男女超爽视频在线观看| 男女高潮啪啪啪动态图| 中文字幕制服av| 欧美精品人与动牲交sv欧美| 在线观看美女被高潮喷水网站| 欧美日韩亚洲高清精品| 国产av精品麻豆| 国产视频首页在线观看| 久久久久久人人人人人| 国产av码专区亚洲av| 伦理电影免费视频| 国产爽快片一区二区三区| 交换朋友夫妻互换小说| 欧美国产精品一级二级三级| 日日啪夜夜爽| 欧美精品一区二区大全| 美女xxoo啪啪120秒动态图| 精品99又大又爽又粗少妇毛片| 亚洲av电影在线进入| xxx大片免费视频| 久久久精品区二区三区| 亚洲经典国产精华液单| 亚洲综合色惰| 水蜜桃什么品种好| 亚洲av电影在线进入| 欧美中文综合在线视频| 成年动漫av网址| 免费av中文字幕在线| 成人18禁高潮啪啪吃奶动态图| 亚洲第一青青草原| av国产精品久久久久影院| 亚洲中文av在线| 久久韩国三级中文字幕| 你懂的网址亚洲精品在线观看| 国产精品久久久久久av不卡| av.在线天堂| 久久久久网色| 亚洲精品视频女| 性色av一级| 99香蕉大伊视频| 黑人巨大精品欧美一区二区蜜桃| av片东京热男人的天堂| 国产日韩欧美在线精品| 男女国产视频网站| 国产在线视频一区二区| 亚洲欧美一区二区三区久久| 大话2 男鬼变身卡| 久久久久网色| 一区二区三区精品91| 亚洲精品在线美女| 亚洲美女视频黄频| 丝袜喷水一区| 丰满乱子伦码专区| 亚洲久久久国产精品| 日日爽夜夜爽网站| 国产有黄有色有爽视频| 欧美精品亚洲一区二区| 国产精品一国产av| 久久人人97超碰香蕉20202| 亚洲国产av影院在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 丰满乱子伦码专区| 亚洲国产日韩一区二区| 最近手机中文字幕大全| 亚洲激情五月婷婷啪啪| 亚洲精品中文字幕在线视频| 少妇人妻久久综合中文| 免费少妇av软件| 欧美少妇被猛烈插入视频| 国产伦理片在线播放av一区| 国产精品久久久久久av不卡| 国产免费福利视频在线观看| 97精品久久久久久久久久精品| 91aial.com中文字幕在线观看| 国产爽快片一区二区三区| 日韩大片免费观看网站| 少妇的逼水好多| 黄色怎么调成土黄色| 亚洲人成77777在线视频| 国产精品嫩草影院av在线观看| 亚洲精品自拍成人| 久久久久久人人人人人| 少妇的丰满在线观看| 亚洲精品久久成人aⅴ小说| 亚洲国产精品一区二区三区在线| 丝袜脚勾引网站| 国产精品偷伦视频观看了| 狠狠精品人妻久久久久久综合| 性少妇av在线| 80岁老熟妇乱子伦牲交| 高清不卡的av网站| 大话2 男鬼变身卡| 久久精品久久久久久久性| 欧美日韩精品成人综合77777| 亚洲图色成人| 2018国产大陆天天弄谢| 99国产精品免费福利视频| 少妇猛男粗大的猛烈进出视频| 黄频高清免费视频| 97精品久久久久久久久久精品| 哪个播放器可以免费观看大片| 极品少妇高潮喷水抽搐| 夫妻午夜视频| 国产一区二区激情短视频 | 日本黄色日本黄色录像| 高清欧美精品videossex| 一级,二级,三级黄色视频| 亚洲国产欧美日韩在线播放| 一本—道久久a久久精品蜜桃钙片| 成人午夜精彩视频在线观看| 一区二区三区四区激情视频| 你懂的网址亚洲精品在线观看| 一级a爱视频在线免费观看| 国产熟女午夜一区二区三区| 国产日韩欧美亚洲二区| 伦理电影免费视频| 亚洲经典国产精华液单| 97人妻天天添夜夜摸| 亚洲美女黄色视频免费看| 97在线人人人人妻| 国产亚洲最大av| 亚洲,欧美,日韩| 91精品国产国语对白视频| 老汉色av国产亚洲站长工具| 欧美日韩视频高清一区二区三区二| 免费大片黄手机在线观看| 好男人视频免费观看在线| 免费高清在线观看日韩| 一本久久精品| 母亲3免费完整高清在线观看 | 丝袜人妻中文字幕| 亚洲伊人色综图| 国产精品欧美亚洲77777| 国产探花极品一区二区| 日韩免费高清中文字幕av| 这个男人来自地球电影免费观看 | 波野结衣二区三区在线| 国产精品免费视频内射| 人人妻人人澡人人看| 国产淫语在线视频| 亚洲精品av麻豆狂野| 肉色欧美久久久久久久蜜桃| 欧美成人午夜精品| 亚洲av.av天堂| 国产不卡av网站在线观看| 国产成人免费无遮挡视频| 永久免费av网站大全| 免费看av在线观看网站| av免费在线看不卡| 亚洲国产精品999| 一区二区三区乱码不卡18| 亚洲激情五月婷婷啪啪| 久久精品国产a三级三级三级| √禁漫天堂资源中文www| 不卡视频在线观看欧美| 热99国产精品久久久久久7| 91aial.com中文字幕在线观看| 国产亚洲av片在线观看秒播厂| 中国国产av一级| 国产精品av久久久久免费| 人妻系列 视频| 日日撸夜夜添| 亚洲精华国产精华液的使用体验| av卡一久久| 美女大奶头黄色视频| 纯流量卡能插随身wifi吗| 又粗又硬又长又爽又黄的视频| 国产在线视频一区二区| 亚洲国产色片| av福利片在线| 日韩一区二区视频免费看| 久久97久久精品| 亚洲av欧美aⅴ国产| 男人舔女人的私密视频| xxx大片免费视频| 国产在线视频一区二区| 亚洲欧美中文字幕日韩二区| 国产极品粉嫩免费观看在线| 成年女人在线观看亚洲视频| 日日爽夜夜爽网站| 久久免费观看电影| 日日啪夜夜爽| 久久综合国产亚洲精品| 男女免费视频国产| 午夜精品国产一区二区电影| 男女高潮啪啪啪动态图| 人体艺术视频欧美日本| 久久久精品区二区三区| 大码成人一级视频| 大香蕉久久网| 亚洲精品第二区| 精品一区二区免费观看| 亚洲国产精品999| 亚洲四区av| 最黄视频免费看| 日韩一区二区视频免费看| 丝袜美足系列| 婷婷色麻豆天堂久久| 高清在线视频一区二区三区| 久久综合国产亚洲精品| 满18在线观看网站| 九九爱精品视频在线观看| 欧美国产精品一级二级三级| 国产av精品麻豆| 婷婷色av中文字幕| 在线观看三级黄色| 女性生殖器流出的白浆| 一边摸一边做爽爽视频免费| 久久久久视频综合| 日本黄色日本黄色录像| 亚洲欧美一区二区三区国产| 国产成人免费观看mmmm| 久久久久精品人妻al黑| 国产成人aa在线观看| 亚洲av.av天堂| 亚洲精品美女久久久久99蜜臀 | 丰满少妇做爰视频| 午夜av观看不卡| 日韩视频在线欧美| 黄色毛片三级朝国网站| 2022亚洲国产成人精品| 亚洲欧美色中文字幕在线| 看免费av毛片| 伊人亚洲综合成人网| 新久久久久国产一级毛片| 校园人妻丝袜中文字幕| a级毛片在线看网站| 在线观看人妻少妇| 精品少妇久久久久久888优播| 一级片免费观看大全| 亚洲av在线观看美女高潮| 国产综合精华液| 国产福利在线免费观看视频| 国产深夜福利视频在线观看| 久久久久久久久免费视频了| 中文天堂在线官网| 丁香六月天网| 精品少妇久久久久久888优播| 国产精品 国内视频| 老女人水多毛片| 免费观看性生交大片5| 亚洲精品自拍成人| 国产日韩欧美在线精品| 嫩草影院入口| 国产精品一国产av| 欧美变态另类bdsm刘玥| 亚洲四区av| 菩萨蛮人人尽说江南好唐韦庄|