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

    Characterization of multiple chemical components of GuiLingJi by UHPLC-MS and 1H NMR analysis

    2022-07-22 03:52:30JinghoShiXioxiGoAirongZhngXuemeiQinGunhuDu
    Journal of Pharmaceutical Analysis 2022年3期

    Jingho Shi,Xioxi Go,Airong Zhng,Xuemei Qin,*,Gunhu Du,d,**

    aModern Research Center for Traditional Chinese Medicine,Shanxi University,Taiyuan,030006,China

    bSchool of Traditional Chinese Materia Medica and Food Engineering,Shanxi University of Chinese Medicine,Jinzhong,Shanxi,030619,China

    cShanxi Guangyuyuan Chinese Medicine Co.,Ltd.,Jinzhong,Shanxi,030800,China

    dInstitute of Materia Medica,Chinese Academy of Medical Sciences&Peking Union Medical College,Beijing,100050,China

    ABSTRACT

    GuiLingJi(GLJ),a classic traditional Chinese medicine(TCM)formula,is composed of over 20 herbs,according to the Pharmacopeia of the People's Republic of China.Owing to its various activities,GLJ has been used in clinical settings for more than 400 years in China.However,the ambiguous chemical material basis limits the development of studies on the quality control and pharmacological mechanisms of GLJ.Therefore,comprehensive characterization of the multiple chemical components of GLJ is of great significance for the modernization of this formula.Given the great variety of herbs in GLJ,both UHPLCMS and1H NMR techniques were employed in this study.In addition,solvent extraction with different polarities was used to eliminate signal interference and the concentration of trace components.A variety of MS analytic methods were also used,including implementation of a self-built compound database,diagnostic ion filtering,mass defect filtering,and Compound Discoverer 3.0 analysis software.Based on the above strategies,a total of 150 compounds were identified,including 5 amino acids,13 phenolic acids and glycosides,11 coumarins,72 flavones,20 triterpenoid and triterpenoid saponins,23 fatty acids,and 6 other compounds.Moreover,13 compounds were identified by1H NMR spectroscopy.The UHPLC-MS and1H NMR results supported and complemented each other.This strategy provides a rapid approach to analyzing and identifying the chemical composition of Chinese herbal prescriptions.The current study provides basis for further research on the quality control and pharmacological mechanism of GLJ.

    Keywords:

    GuiLingJi

    Chemical component

    UHPLC-MS

    1H NMR

    Traditional Chinese medicine

    1.Introduction

    Traditional Chinese medicine(TCM)formulas have been used in clinical settings for thousands of years in China.However,although the clinical efficacy of TCM formulas is definite,their ambiguous effective substances pose a main obstacle to the modernization of these formulas[1].TCM formulas apply a multicomponent and multitarget approach in the treatment of diseases,which is different from Western medicines,and the pharmacologic mechanisms have been increasingly clarified.Therefore,comprehensive characterization and identification of the multiple components of TCM formulas is of particular importance.

    Recently,UHPLC-MS has been increasingly used for rapid chemical profiling of complicated medicinal herbs.This technique is especially suitable for the identification of unknown organic compounds because of its high resolution,high sensitivity,high selectivity,and accurate mass values[2-4].With the MS/MS technique,accurate mass values can help elucidate the molecular weight and candidate chemical formula.Moreover,both elemental composition and fragment ions can be obtained,which significantly reduces the number of possible structures of putative compounds and enables the efficiency analysis of complex samples.The diagnostic ion and mass defect filtering strategy minimizes ubiquitous matrix interference,extracts the MS features of some potentially active compounds,and rapidly identifies a series of compounds with the same skeleton[5].However,MS-based approaches can only identify the components by comparison with those of reference compounds,so they are only suitable for molecules that are easily ionized.NMR spectroscopy is one of the most definitive identification methods and shows non-selectivity.However,the lower sensitivity and the overlaps of signals in the complex matrix limitits applications in identifying the chemical components of TCM formulas.Combining NMR and UHPLC-MS methods to identify the matrix chemical components of TCM formulas can facilitate the final characterization[6,7]and provide a new strategy for analyzing the complex systems of TCM formulas.

    GuiLingJi(GLJ)is a well-known and classical TCM formula that can be traced back hundreds of years ago to the Ming and Qing dynasties,where it was used as an elixir.Modern pharmacological studies on GLJ indicate that it can significantly improve kidney-yang deficiency syndrome and testicular dysfunction[8,9],promote spermatogenesis in oligospermia rats by increasing the density and amount of sperms[10],and improve the mitochondrial function of sperms in idiopathic asthenospermia and the level of sperm viability[11].Additionally,GLJ can improve the impairment of learning and memory functions caused by aging[12,13].GLJ has various medical efficacies based on its multiple components.According to the Chinese Pharmacopeia(2020 edition),GLJ is one of the confidential products of national brands of TCM with 20 opened herbs,including Ginsen Radix et Rhizoma Rubar,Cervi Cornu Pantotrichum,Hippocampus,Lycii Fructus,Caryophylli Flos,Cyathulae Radix,Cynomorii Herba,Rehmanniae Radix Praeparata,Psoraleae Fructus,Cuscutae Semen,Eucommiae Cortex,Cistanches Herba,Glycyrrhizae Radix Et Rhizoma,Asparagi Radix,Epimrdii Folium,and Amomi Fructus,among others[14].Each herb has its own unique processing method.For instance,Cervi Cornu Pantotrichum is fried with mature vinegar,Lycii Fructus is fried with honey,and Cistanches Herba is fried with yellow wine.More than 20 kinds of single herbs are processed and mixed with powder.Then,different from traditional decocting or extraction methods,GLJ is made of the above mixed powder by means of alchemy.For such an exclusive and complicated TCM formula possessing precise curative effects,few studies have investigated the systematic characterization of chemical components in GLJ.Therefore,clarifying the material basis of GLJ will provide useful chemical information for quality control and mechanism-of-action research.

    In this study,an enhanced identification strategy was developed for comprehensive characterization of the multiple chemical components of GLJ.First,the chemical components of the 16 herbs in the prescription were summarized to establish a selfbuilt database of internal compounds to predict the chemical components of GLJ.Second,the methanol extract of GLJ was partitioned into different polar parts and characterized with UHPLC-MS techniques.A series of methods,such as summarizing the fragmentation regularity,diagnostic ion filtering,and mass defect filtering and inference by Compound Discoverer 3.0(CD)software,was used to analyze the mass spectrometric data.Moreover, chemical compounds of GLJ were identified by comparing the retention time and MS/MSndata with those of individual herbs.Furthermore,two polar fractions of GLJ obtained by two-phase extraction were characterized by1H NMR spectroscopy.Finally,a series of methods was combined to explore the complicated chemical composition of GLJ.The framework of this study is displayed in Fig.1.

    2.Experimental

    2.1.Reagents and materials

    GLJ(batch No.:20170726)was generously provided by Shanxi Guangyuyuan Traditional Chinese Medicine Co.,Ltd.(Jinzhong,China).Standard references included kaempferol(≥98%,batch No.:MUST-17032911),psoralen(≥98%,batch No.:MUST-16092810),baohuoside I(≥98%,batch No.:MUST-19042203),isobavachalcone(≥98%,batch No.:MUST-16102615),neobavaisoflavone(≥98%,batch No.:MUST-16092810),bavachinin A(≥98%,batch No.:MUST-18101705),oleanolic acid(≥98%,batch No.:MUST-17042706),ginsenoside Re(≥98%,batch No.:180829),and ginsenoside Rg1(≥98%,batch No.:180917),which were obtained from Chengdu Must Biotechnology Co.,Ltd.(Chengdu,China)and Shanghai Ronghe Pharmaceutical Technology Development Co.,Ltd.(Shanghai,China).Methanol,petroleum ether,ethyl acetate,and nbutanol of analytical grade were purchased from Tianjin Four Fine Chemicals Co.,Ltd.(Tianjin,China).Ultrapure water was prepared with a Milli-Q water purification system(Millipore,Bedford,MA,USA).Methanol of LC-MS grade and formic acid of HPLC grade were obtained from Thermo Fisher Scientific Inc.(Thermo Fisher Scientific,Waltham,MA,USA).Deuterated chloroform(CDCl3,99.8%D)containing tetramethylsilane(TMS,0.03%,m/V)and methanol-d4(99.8%D)were obtained from J&K Scientific Ltd.(Beijing,China).

    2.2.Sample preparation

    The reference standards(0.1-1.0 mg/mL)were accurately weighed and then mixed and dissolved in methanol of LC-MS grade.An accurately weighed 1.00 g sample of GLJ powder was ultrasonically extracted with 50 mL of methanol at room temperature for 30 min.The extracts were filtered and concentrated to be dried.One part of the residue was transferred to a 5 mL volumetric flask and dissolved in LC-MS grade methanol to volume as the methanol extract of GLJ(GLJ-M).The 16 herbs were extracted under the same process as GLJ-M to identify the source of ingredients.Meanwhile,the other part of the residue was suspended in 30 mL of water and successively treated with petroleum ether,ethyl acetate,and n-butanol.Then,the petroleum ether fraction(GLJ-P),ethyl acetate fraction (GLJ-E)and n-butanol fraction (GLJ-B)were concentrated to dryness.The residue was dissolved in LC-MS grade methanol,and the solution was completely transferred into a 2 mL volumetric flask.All solutions were filtered through a nylon syringe filter(0.22μm)for UHPLC-MS analysis.

    Fig.1.Summary diagram of the combined strategy of UHPLC-MS and1H NMR analysis.GLJ:GuiLingJi;GLJ-M:the methanol extract of GLJ;GLJ-P:the petroleum ether fraction;GLJ-E:ethyl acetate fraction;GLJ-B:the n-butanol fraction;GLJ-C:the chloroform fraction;GLJ-W/M:the aqueous methanol fraction.

    For1H NMR analysis,1.00 g of GLJ powder was exactly weighed and ultrasonically extracted with 24 mL of extraction solvent composed of water,methanol,and chloroform at a volume ratio of 1:1:2 at room temperature for 30 min.The extracts were centrifuged at 3,500 r/min for 20 min.Then,the two layers were separated and concentrated to dryness,and the residues were obtained.The chloroform fraction(GLJ-C)was dissolved in 800μL of CDCl3,and the aqueous methanol fraction(GLJ-W/M)was dissolved in 800μL of CD3OD.Then,the resulting mixture was centrifuged for 15 min at 13,000 r/min.Finally,600μL of supernatant was transferred to an oven-dried 5 mm NMR tube for1H NMR analysis.

    2.3.UHPLC-MS conditions

    UHPLC-MS analysis was conducted on a Thermo-Fisher UHPLC system coupled to a Q Exactive Orbitrap high-resolution mass spectrometer(UHPLC-QE HRMS;Thermo Fisher Scientific,Waltham,MA,USA).Chromatographic separation was carried out on a Waters Acquity UPLC HSS T3 column(2.1mm×100mm,1.8μm;part No.:186003539)maintained at 30°C.The mobile phase of solvent A(water with 0.1% formic acid)and solvent B(acetonitrile)was operated under the following program:0-10 min,5%-30%B;10-15 min,30%-52%B;15-18 min,52%B;18-26 min,52%-87%B,26-30 min,87%B;30-35min,87%-100%B;and 35-38 min,100%-5%B.The sample injection volume was 3μL,and the flow rate was 0.2mL/min.

    The mass spectrometer was fitted with an electrospray ionization source,and mass detection was carried out in both positive and negative ion modes with the following settings:capillary temperature,320°C;heater temperature,300°C;sheath gas velocity,35 arb;auxiliary gas flow rate,10 arb;full scan range,m/z 100-1500;and MS2scan range,m/z 50-1500.Fast data-dependent acquisition was applied to trigger MS/MS acquisition of precursor ions with an intensity threshold of 5000.Full scan spectra were acquired at a scan time of 0.3s,followed by MS/MS spectra acquisition at a scan time of 0.1 s.The top 15 most intense ions were selected to perform MS/MS acquisition.The collision energy was set to a ramp of 25-60 V for MS/MS acquisition.

    2.4.1H NMR conditions

    1H NMR spectra were measured at 600.13 MHz using an Avance 600Bruker spectrometer (Bruker Biospin, Rheinstetten, Germany) equipped with a Bruker 5 mm double resonance BBI probe.CD3OD and CDCl3were used for internal lock purposes.Spectra of GLJ-C and GLJ-W/M were acquired using the zg30 sequence.Sixty-four scans were recorded,resulting in a total acquisition time per sample of 5min.

    2.5.Data processing

    The UHPLC-QE HRMS raw data were automatically processed using an Xcalibur workstation(Thermo Fisher Scientific,Waltham,MA,USA).Then,the processed data were imported into CD software(Thermo Fisher Scientific,Waltham,MA,USA)for peak extraction and normalization.The peak width at 5% height and the peak-to-peak baseline noise were automatically determined,and the processing parameters were set as follows:retention time range,0-38 min;mass range,100-1500Da;mass tolerance,5 ppm;retention time tolerance,0.05min;S/N threshold,1.5;negative adducts,[M-H]-,[M+HCOO]-;and positive adducts,[M+H]+and[M+Na]+.The compounds characterized by CD software were identified by searching the mz Cloud,Masslist,Chemspider,and PubChem databases.

    The NMR spectra were processed using MestReNova(version 6.1.1,Mestrelab Research,Santiago de Compostella,Spain)to manually correct the baseline and phase and calibrated to TMS at δ0.00 for CDCl3.All spectra were referenced to the residual signal of CD3OD at δ3.31.

    3.Results and discussion

    3.1.Establishment of compound databases and analysis of fragmentation rules

    According to Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform(https://tcmsp-e.com/tcmsp.php),Web of Science(https://www.isiknowledge.com/),and CNKI(http://www.cnki.net/),a database of compounds containing compound names,molecular formulas,accurate molecular masses,and chemical structures was established.The established compound database included 16 of the medicinal herbs in GLJ without some mineral and animal drugs,as the chemical compositions of these mineral and animal drugs had not been studied.In addition,mass spectrum fragment ions of some compounds were obtained by Massbank(http://massbank.eu/)and the literature,and the1H NMR chemical shifts of some primary metabolites were obtained by HMBD(https://hmdb.ca/).Finally,a total of 1246 compounds including 358 flavonoids,179 triterpenes and their glycosides,143 phenylpropanes and their glycosides,123 steroids and their glycosides,98 organic acid,54 saccharides,20 amino acids,and 271 other compounds(alkaloids,tannins,polyphenols,unsaturated hydrocarbons,etc.)were summarized.

    The fragmentation rules of the compounds were also studied,taking flavonoids and ginsenoside as examples.Representative authentic standards of flavonoids,including kaempferol,baohuoside I,isobavachalcone,neobavaisoflavone and bavachinin A,were subjected to analysis by a UHPLC-QE HRMS system in both positive and negative ion modes.Fragmentation patterns of each standard were summarized based on assigned fragment ions.The prominent product ions were selected as diagnostic ions.The prime diagnostic ions or neutral loss fragments of these compounds were confirmed by MS/MS spectra analysis of 7 standards as well as by the self-built compound database.Because the characteristic fragmentation pattern of flavonoids indicated a retro-Diels-Alder(RDA)reaction,significant product ions of the RDA reaction at m/z 119.04922,135.04453,133.02904,and 151.00322 were set as diagnostic ions.The possible fragment pathways of flavonoids in negative ion mode are shown in Fig.2A.

    From the fragment ions of the reference standard including ginsenosides Re and Rg1and the reported literature[15-17],fragmentation rules could be found.The protopanaxadiol,protopanaxatriol,and oleanolic acid types of ginsenosides were found at m/z 459,475,and 455,corresponding to successive losses of sugar moieties,respectively.Characteristic neutral losses of these sugar residues were 162 Da(hexose),146 Da(deoxyhexose),and 132 Da(fructose).The lost sugar units,e.g.,glucose,rhamnose,and xylose,produced common highly abundant ions at m/z 101,113,119,161,and 179.Based on the above fragmentation regularity,prominent product ions such as m/z 455,459,and 475 were selected as diagnostic ions of the analogous skeletons,and m/z 101,113,119,161,and 179 were selected as diagnostic ions of the sugar residues.The possible fragment pathways of ginsenosides in negative ion mode are shown in Fig.2B,with ginsenoside Rg1as an example.

    Fig.2.Possible fragment pathways and diagnostic ions of(A)flavonoids and(B)ginsenoside.

    3.2.UHPLC-QE HRMS analysis and identification of the main constituents in GLJ

    Prior to real sample analysis,extraction methods and chromatography conditions were fully evaluated and optimized.To achieve better chromatographic separation and mass spectrometric detection,different extraction solvents including 100% methanol,75% methanol,and 100% water were explored.The extracts with 100% methanol showed higher peak responses but poor peak separations;thus,solvents with different polarities were used.Various compositions of mobile phases such as methanol-water,acetonitrile-water,and acetonitrile-water with formic acid were explored.Acetonitrile-water with formic acid was chosen as the mobile phase due to its more symmetrical peak shape and ionization efficacy in both positive and negative ion modes.

    Methodological verification was carried out to monitor the stability of the LC-MS system.Ten ions with different retention time and m/z values from the base peak intensity chromatograph of 6 GLJ-M samples were extracted.The RSDs for retention time were 4.88×10-4%-3.77×10-3%,the m/z values were 3.64×10-6%-4.85×10-5% and the relative peak areas were 5.02%-12.94%(Table S1),which demonstrated the favorable stability of the system and the reliability of the data analysis.

    A total of 150 compounds were tentatively characterized,including 5 amino acids,13 phenolic acids and glycosides,11 coumarin,72 flavone,20 triterpenoid and triterpenoid saponins,23 fatty acids and 6 others.Among them,9 components were definitively identified by comparing their reference standards,while the others were tentatively identified by comparing the accurate masses,ESI-MSnfragment ions,and molecular weights with those in the self-built compound database.CD software provided some references including matched data and fitting compounds.The total ion chromatograms of GLJ-M,GLJ-P,GLJ-E,and GLJ-B in both negative and positive ion modes are shown in Fig.3.All the compounds that we identified in GLJ were already known;no new compounds were identified.The corresponding chromatographic and mass spectral data of 150 compounds are shown in Tables S2-S13.The source of ingredients was identified by comparing the accurate masses and retention time with the 16 herbs.Some amino acids and fatty acids may be mainly derived from multiflavor animal medicines without exclusive sources.The exact peaks with mass ranges and ion fragments of standard references and GLJ are shown in(Figs.S1-S9).

    3.2.1.Identification of amino acids in GLJ

    Five amino acids were identified in GLJ(Table S2).Amino acids might lose COOH in the reaction of α-cracking.For instance,compound 3 was detected at 2.42min,which indicated that it is a polar compound.The molecular formula of compound 3 was determined to be C9H11NO3according to the result that the typical MS/MS spectra of the protonated ions at m/z 182.08119[M+H]+produced the ions at m/z 136.07578 and 119.04925 due to successive losses of COOH and NH2(Fig.4).Compound 3 was identified as tyrosine by comparing the fragment ions with those in the literature[18].Other amino acids were inferred by CD and further identified by comparing the MS/MSndata with those in the literature[18].The chemical structures of amino acids in GLJ are shown in Fig.S10.

    3.2.2.Identification of phenolic acids and glycosides in GLJ

    Thirteen phenolic acids and glycosides were identified in GLJ in negative ion mode(Table S3).The characteristic fragmentation pattern of phenolic acids was neutral losses of CO2,OH or COOH.For instance,compound 2 was deduced as gallic acid by CD software.Its retention time was 2.10 min,and the molecular formula was determined to be C7H6O5.The typical MS/MS spectra of the deprotonated ions at m/z 169.01378[M-H]-produced ions at m/z 125.02368[M-H-CO2]-and 107.01310[M-H-CO2-OH]-.The MS/MS2spectra and fragmentation patterns are shown in Fig.5.Compound 2 was identified as gallic acid by comparing the fragment ions with those in the literature[6].Other phenolic acids were identified by comparing the exact mass,molecular formulas and fragment ions with those in a self-built compound database.

    Five phenylpropanoid glycosides were identified in GLJ.For example,compound 13 displayed an[M-H]-ion at m/z 487.14587 and generated fragment ions at m/z 179.03442 and 161.02377.The fragment ion at m/z 179.03442 was a caffeoyl species and was generated by neutral losses of 162 Da(hexose)and 146 Da(deoxyhexose).The fragment ions at m/z 161.02377 confirmed the existence of hexose(Fig.6).Accordingly,compound 13 was identified as a phenylpropanoid glycoside named cistanoside F.Then,compounds 21,23,29,and 45 were filtered by taking the fragment ions m/z 179 as diagnostic ions or the mass defect of 162,146 Da.These compounds were identified as echinacoside,caffeic acid 3-glucoside,acteoside,and ferulic acid-4-β-glucoside by comparing the fragment ions and molecular formulas with the literature[19].The chemical structures of phenolic acids and glycosides in GLJ are shown in Fig.S11.

    3.2.3.Identification of coumarins in GLJ

    Eleven coumarins were identified in GLJ(Table S4).The main MS/MS behavior of coumarins was neutral loss of 28Da(CO)and 44Da(CO2).For example,compound 72 was identified by comparing the retention time and ion fragments with the reference compound.Protonated ions at m/z 187.03894[M+H]+were observed,which further generated the ions at m/z 159.04366[M+H-CO]-,m/z 143.04915[M-H-CO2]-,m/z 131.04919[M-H-CO-CO]-and 115.05441[M+H-CO2-CO]-.The mass spectra and proposed major fragmentation with structures are shown in Fig.7.Similarly,compound 90 displayed a fragment ion at m/z367.11911[M-H]-.Then,it generated ions at m/z 352.09702[M-H-CH3]-and m/z 297.04062[M-H-C5H10]-,which provided methyl and isopentane groups.The characteristic fragment ions at m/z 324.10007[M-H-CH3-CO]-and m/z 269.04532[M-H-C5H10-CO]-were explained by the dissociation of the lactone ring.Compound 90 was identified as isoglycycoumarin based on the above MS/MSndata.Thus,compounds 43,56,74,81,87,89,106,110,and 119 were tentatively identified by comparing the MS/MSndata with those in the literature[20-22].The chemical structures of coumarins in GLJ are shown in Fig.S12.

    Fig.3.The total ion chromatograms of GLJ-M,GLJ-P,GLJ-E,GLJ-B in(A)the positive and(B)the negative ion modes.

    Fig.4.The MS/MS spectra and proposed fragmentation for tyrosine in positive ion mode.

    3.2.4.Identification of flavones in GLJ

    A total of 72 compounds were characterized as flavonoids,including flavonols,flavones,chalcones,isoflavones,flavanones,and flavone glycosides.These flavonoids with different structural skeletons had different characteristic fragment ions.

    3.2.4.1.Identification of flavonols.Eleven flavonols were identified in GLJ(Table S5).Compound 71 was determined to be kaempferol by comparing the retention time and ion fragments with the reference compound.Compound 71 displayed deprotonated ions at 285.04077[M-H]-and produced ions at m/z 257.04501 by the neutral loss of CO.A series of ions at m/z 151.00322 and m/z 133.02904 were generated by RDA cleavage.Similarly,the mass spectra of compound 5 showed an abundant[M-H]-ion at m/z 301.03488,which generated a product ion at m/z 273.04028 by neutral loss of CO,a product ion at m/z 178.99866 by loss of the B-ring and CH,and a product ion at m/z 151.00308 by the RDA reaction.Compound 5 was identified as quercetin by comparing the fragment ions with those in the literature[23].Other flavonols were filtered by taking the fragment ions at m/z 285,255,151,301,and 178 as diagnostic ions.Then,compounds 10,16,52,60,73,94,95,96,and 104 were further identified by comparing the MS/MSndata with the literature[23,24].The chemical structures of flavonols identified in GLJ are shown in Fig.S13.

    Fig.5.The MS/MS spectra and proposed fragmentation for gallic acid in negative ion mode.

    Fig.6.The MS/MS spectra and proposed fragmentation for cistanoside F in negative ion mode.

    3.2.4.2.Identification of flavones.Ten flavones were identified in GLJ(Table S6).For example,the retention time of compound 55 was 14.99min,and its molecular formula was C15H10O6.Compound 55 displayed a[M-H]-ion at m/z 285.04077.In addition,the fragment ions at m/z 257.04565 and m/z 241.35442 were generated by successive loss of CO and O,respectively,and the fragment ions at m/z 151.00298 and m/z 133.02856 were generated by RDA cleavage.The relative abundance of m/z 133.02856 was higher than that of m/z 151.00298,which is different from the fragment ion abundance of kaempferol.Compound 55 was identified as luteolin by comparing the MS/MSndata with those in the literature[25].Compounds 69,80,and 100 were filtered by taking the fragment ion m/z 151 as the diagnostic ion and further identified by comparing the MS/MSndata with those in the literature[23,26,27].Other flavones were tentatively characterized by comparing the retention time and ion fragments with those in the database of internal compounds.The chemical structures of flavones in GLJ are shown in Fig.S14.

    3.2.4.3.Identification of chalcone.Eight chalcones were identified in GLJ(Table S7).Compound 103 was determined to be isobavachalcone by comparing the retention time and ion fragments with those of the reference compound.Compound103 displayed an[M-H]-ion at m/z 323.12918 and a fragment ion at m/z 268.07498[M-H-C4H7]-,indicating the existence of an isopentane group.The fragment ions at m/z 203.07111 and m/z 119.04950 were generated from the α-broken keto groups.Compounds 93,112 and 122 were filtered by taking the fragment ion m/z 119 as the diagnostic ion and further identified by comparing the MS/MSndata with the literature[21,28].Other chalcones were tentatively characterized by comparing the exact mass and ion fragments with those in the self-built compound database.The chemical structures of chalcones in GLJ are shown in Fig.S15.

    Fig.7.The MS/MS spectra and proposed fragmentation for psoralen in positive ion mode.

    3.2.4.4.Characterization of isoflavones.Twenty isoflavones were identified in GLJ(Table S8).Compound 97 was determined to be neobavaisoflavone by comparing the retention time and ion fragments with those of the reference compound.Compound 97 displayed an[M+H]+ion at m/z 323.12775 and a fragment ion at m/z 267.06509[M+H-C4H7]-,indicating the existence of an isopentene group.The fragment ion at m/z 239.06970[M+H-C4H7-CO]-was observed by neutral loss of CO successively,and the fragment ion at m/z 137.02341 was generated from RDA cleavage.Other isoflavones were filtered by taking the fragment ion m/z 151 as the diagnostic ion and further identified by comparing the MS/MSndata recorded in the literature[24,29-31].The chemical structures of isoflavones contained in GLJ are shown in Fig.S16.

    3.2.4.5.Characterization of flavanone.Seven flavanones were identified in GLJ(Table S9).Compound 121 was determined to be bavachinin A by comparing the retention time and ion fragments with those of the reference compound.Compound 121 has a precursor ion at m/z 339.15903[M+H]+.The product fragment ions at m/z 271.09631[M-C5H9]+were generated by the loss of the isopentane group. The fragment ions at m/z 219.10146[M+H-C8H8O]+,m/z 151.03885[M+H-C5O9-C8H8O]+,and m/z 119.04922[M+H-C13H14O3]-were the characteristic ions of RDA cleavage.The other 6 compounds were filtered by taking the fragment ions m/z 135 and 119 as diagnostic ions and identified by comparing the MS/MSndata in the literature[20,21,28].The chemical structures of flavanone in GLJ are shown in Fig.S17.

    3.2.4.6.Characterization of flavone glycosides.Sixteen flavone glycosides were identified in GLJ(Table S10).Flavone glycosides showed characteristic neutral losses of 176 Da(glucuronic acid,GlcA),162 Da(glucosyl,Glu),146 Da(rhamnosyl,Rha),and 132 Da(arabinosyl,Ara).The aglycon structure might lose CO,CO2,CHO and so on during RDA cleavage.For example,compound 91 was determined to be baohuoside I by comparing the retention time and ion fragments with those of the reference compound.Compound 91 has a precursor ion at m/z 513.17688[M-H]-.The fragment ion at m/z 366.11108 was generated by characteristic neutral losses of 146 Da(rhamnose);then,the fragment ions at m/z 351.08771 and m/z 323.09189 were generated by successive losses of CH3and CO.Other flavone glycosides were filtered by using 176,162,146,and 132 Da as mass defects;these compounds were tentatively identified by comparing the MS/MSndata with those in the literature[32-37].The chemical structures of flavone glycosides in GLJ are shown in Fig.S18.

    3.2.5.Identification of triterpenes and triterpenoid saponins in GLJ

    Six triterpenes and 14 triterpenoid saponins were identified in GLJ(Table S11).Compound 142 was determined to be oleanolic acid by comparing the retention time and ion fragments with those of the reference compound.The deprotonated ion at m/z 455.35361 produced ions at m/z 407.32980[M-H-CO2-4H]-and 391.25125[M-H-CO2-2H-H2O]-.The fragment ion at m/z247.28014[M-H-C14H24O]-was the characteristic ion of RDA cleavage.Fragment ions such as[M-H-CO2]-and[M-H-HCOOH]-were regarded as the characteristic ions of triterpenes.Compounds 65,116,129,130,and 132 were tentatively characterized as triterpene acids based on the above fragment rules[27,38].

    The triterpenoid saponins in GLJ were mainly ginsenosides.Ginsenosides,such as ginsenoside Rg1,ginsenoside Rb,ginsenoside Re,and ginsenoside Rd,which have relatively higher contents in ginseng,were detected in the GLJ-M and GLJ-E fractions.However,many ginsenosides such as ginsenoside Ro,ginsenoside Rf,and ginsenoside Rc with relatively low contents could be detected only in the n-butanol fraction.As an example,compound 38 was determined to be ginsenoside Re by comparing the retention time and ion fragments with those of the reference compounds.Typical MS/MS spectra of the deprotonated ions of the standard reference samples at m/z 991.54944[M-H+HCOOH]-were assigned to ginsenoside Re.The deprotonated ion at m/z 945.54315 produced ions at m/z 783.48767,637.43170,and 475.37979 due to successive losses of glucose,rhamnose,and glucose.The structures of the other 12 ginsenosides were screened out by the diagnostic ion strategy and further identified by comparing the MS/MSndata with those in the literature[39,40].The chemical structures of triterpenes and triterpenoid saponins in GLJ are shown in Fig.S19.

    3.2.6.Identification of fatty acids in GLJ

    Twenty-three fatty acids were identified in GLJ(Table S12).Fatty acids were mainly detected in the GLJ-P fraction.These compounds could be rapidly identified by characteristic ions including[M-H-CO]-,[M-H-CO2]-,[M-H-H2O]-,[M-H-COOH]-,[M-H-COOH]-,[M-H-CH]-,[M-H-CH2]-,and[M-H-CnH2n]-.Take compound 138 as an example.It had a retention time of 29.67min,and the molecular formula was C18H30O2.Compound 138 had a precursor ion at m/z 277.21759[M-H]-and generated product ions m/z 259.20721 and m/z 233.22884 by neutral loss of H2O and CO2,respectively.The fragment ions at m/z 84.41108,m/z 71.01276,and m/z 59.01290 were the characteristic ions of aliphatic and unsaturated aliphatic compounds.Compound 138 could be identified as α-linolenic acid or γ-linolenic acid by comparing the MS/MSndata with MassBank.Other fatty acids were inferred by CD software and then identified by comparing the MS/MSndata with MassBank data.

    3.3.1H NMR analysis of different polar parts of GLJ

    To obtain a better separation degree and more compound information,various extraction conditions were investigated,including solvent extraction with different polarities and twophase extraction.Eventually,the two-phase extraction with water/methanol and chloroform was selected because the1H NMR spectra obtained from this extraction showed a higher peak response value and better separation degree and covered all distinct regions such as organic acids,amino acids,and phenolic regions.

    A total of 13 compounds,namely,6 organic acids or their ester,2 amino acids,2 organic bases,and 3 sugars,were characterized by1H NMR spectroscopy.Moreover,1H NMR analysis supported the existence of fatty acids and triglycerides.

    3.3.1.Characterization of chloroform fraction

    The1H NMR spectrum of GLJ-C showed strong signals in the region from δH 0.5 to δH 3.2 ppm,where most signals were derived from amino acids and organic acids(Fig.8A).The terminal methyl(δ0.880,t,J= 7.2 Hz),methylene(δ2.346,t,J= 7.2 Hz),allylic CH2(δ2.015,m),diallylic CH2(δ2.771 t,J= 6.6 Hz),other hydrocarbon chain(δ1.213-1.328),and olefinic(δ5.348,m)protons confirmed the existence of saturated and unsaturated fatty acids or their ester or amides in the chloroform fraction.The1H NMR signals of GLJ-C also showed signals at δ4.146(dd,J=6,12 Hz)and δ4.293(dd,J=4.2,12 Hz),which represented the methylene group of the glycerol backbone from triglycerides[6,41].

    3.3.2.Characterization of the aqueous methanol fraction of GLJ

    The1H NMR spectrum of GLJ-W/M can be divided into three distinct regions(Fig.8B).Organic acids and amino acids including threonine(δ 1.340,d,J= 7.2Hz),acetic acid(δ 1.929s),succinic acid(δ2.454,s),succinate(δ2.432,s),creatine(δ 3.059s;δ 3.960,s),and glycine(δ3.689,s)were identified in the corresponding characterized regions.Other signals(δ3.216,s and δ3.278,s;δ3.920,s)were identified as belonging to choline and betaine,respectively.Lactate(δ 4.134,q,J=7.2Hz),raffinose(δ4.964,d,J=4.2Hz;δ5.443,d,J=3.6 Hz),and sucrose(δ5.422,d,J= 4.2Hz)were identified in the carbohydrate region,and formic acid(δ8.462,s)was identified in the phenolic region.The chemical shifts and coupling constants of all the identified compounds are summarized in Table 1.

    Table 1Chemical shift assignments in different fractions of GuiLingJi(GLJ).

    Fig.8.1H NMR spectra of(A)GLJ-C and(B)GLJ-W/M.Numbers 1-15 represent the characteristic group peaks of compounds 1-15 in Table 1.

    3.4.Chemical composition of GLJ

    An overview of UHPLC-QEHRMS and1H NMR analyses combined with a polar division strategy for comprehensive characterization of the chemical composition of GLJ is shown in Fig.9.To further clarify the distribution of the 150 compounds in different polar parts,a relative content heat map and Venn diagram were employed.Thus,only 82 compounds could be identified in GLJ-M;many components in the methanol extraction failed to be detected under the current test conditions due to serious signal interference.Moreover,77 compounds were identified in GLJ-P,114 compounds in GLJ-E,and 46 compounds in GLJ-B.More compounds were identified by using the method of polar division,and a total of 150 compounds were rapidly identified in all polar parts of GLJ by UHPLC-QE HRMS.In addition,13 compounds including organic acids,amino acids,organic bases,and sugars were rapidly identified by1H NMR spectroscopy.The1H NMR spectra also definitively confirmed the existence of fatty acids,which were identified by UHPLC-QE HRMS.

    3.5.Pharmacological activities of the main chemical components in GLJ

    This study demonstrated that GLJ consists of multiple chemical components,which mainly include flavonoids and saponins,ginsenosides,coumarins,phenolic acids and so on.These chemical compositions have various pharmacological activities.For example,multiple pharmacological actions of ginsenosides have been observed to have anti-inflammatory activity and anti-central nervous system disorder effects,enhancing energy and sexuality[42].In particular,ginsenosides Re and Rp1possess anti-inflammatory properties by suppressing the NF-κB signaling pathway[43].Ginsenoside Rb1has an antifatigue effect by suppressing skeletal muscle oxidative stress and improving energy metabolism[44].Ginsenosides Rb1and Rg1may significantly ameliorate learning and memory impairment by increasing the acetycholine level and inhibiting the decrease in 5-hydroxytryptamine in the brain[45].In addition,ginsenosides Rb1and Rg1have been shown to improve testicular function in spermatogenesis by upregulating follicle stimulating hormone levels and luteinizing hormone levels[46].Ginsenoside Rb1also has a potent antifatigue effect caused by its suppression of skeletal muscle oxidative stress and improves energy metabolism[44].A large number of flavonoids were identified in GLJ.Flavonoids can reverse age-related declines in neurocognitive performance by interacting with the cellular and molecular architecture of the brain responsible for memory and reducing neuronal loss due to neurodegenerative processes[47].For example,epimedium flavonoids have been reported to improve learning-memory impairment and decrease the production of β-amyloid in the brains of amyloid-β precursor protein(APP)transgenic mice.These compounds might also improve cognitive impairment and white matter lesions induced by chronic cerebral hypoperfusion by inhibiting the Lingo-1/Fyn/ROCK signaling pathway and activating the BDNF/TrkB,NRG-1/ErbB4 and PI3K/Akt/CREB pathways in white matter[48].Icariin,as the major ingredient of Epimrdii Folium,could improve learning and memory functions in APP/presenilin-1 transgenic mice,possibly by stimulating the NO/cGMP signaling pathway and coordinating the induction of nitric oxide synthase(NOS)isoforms[49].Furthermore,it was reported to improve spatial learning and memory abilities as it causes a decreased expressions of tumor necrosis factor-α,interleukin-1,and cyclooxygenase-2 in the hippocampus[50].In addition,phenolic acids such as caffeic acid have beneficial and/or protective effects in improving reproductive function,preventing sperm damage,ameliorating sperm function,and improving testicular function[51].Coumarins have shown inhibitory effects against infection from various viruses such as human immunodeficiency virus,flu virus,enterovirus 71,and coxsackievirus A16[52].In conclusion,the synergistic combination of various compounds could achieve the multiple functions of GLJ.However,the current quality control standards of GLJ only include the content determination of two ginsenosides,which cannot fully reflect the quality of this formulation.The results of this study suggest that further attention should be paid to some flavonoids,which will help to find quality markers and improve quality standards.

    Fig.9.The overview of the UHPLC-QE HRMS and1H NMR for comprehensive characterization of chemical composition of GLJ.

    4.Conclusion

    In this study,a combined approach of UHPLC-QE HRMS and1H NMR techniques coupled with a polarity partitioning strategy was used to comprehensively characterize multiple chemical components of GLJ.The compound database,CD software,diagnostic ion filtering,and mass defect filtering strategy were employed to analyze the mass spectrum.Finally,a total of 150 compounds were rapidly identified in GLJ by UHPLC-QE HRMS,of which 9 compounds were unambiguously identified by comparing the retention time and MS/MS data with those of reference compounds.In addition,13 compounds were rapidly identified by1H NMR spectroscopy.These two techniques could potentially complement each other.This analysis strategy can be applied to other TCM formulas for comprehensive characterization and identification of chemical components.

    CRediT author statement

    Jingchao Shi:Data curation,Writing-Original draft preparation;Xiaoxia Gao:Investigation,Writing-Reviewing and Editing;Airong Zhang:Data curation,Funding acquisition;Xuemei Qin:Supervision;Guanhua Du:Project administration.

    Declaration of competing interest

    The authors declare that there are no conflicts of interest.

    Acknowledgments

    We acknowledge financial support from the Shanxi Provincial Science and Technology Key Research and Development Plan(Grant Nos.:201603D3113006 and 201903D311012).We also thank the Scientific Instrument Center of Shanxi University for the UHPLC-QE HRMS and1H NMR analyses.

    Appendix A.Supplementary data

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

    国产伦精品一区二区三区视频9 | 国产精品久久久久久久久免 | 亚洲av第一区精品v没综合| 国产视频一区二区在线看| 午夜免费男女啪啪视频观看 | 国产免费av片在线观看野外av| 亚洲欧美日韩无卡精品| 国产精品久久久久久亚洲av鲁大| 三级国产精品欧美在线观看| 一区二区三区激情视频| 特级一级黄色大片| 免费一级毛片在线播放高清视频| 一级毛片女人18水好多| 亚洲在线观看片| 18+在线观看网站| 午夜两性在线视频| 国产午夜精品久久久久久一区二区三区 | 美女cb高潮喷水在线观看| 少妇的逼好多水| 午夜福利18| 99视频精品全部免费 在线| 免费看光身美女| 国产午夜精品论理片| 国产精品亚洲美女久久久| 啦啦啦韩国在线观看视频| 法律面前人人平等表现在哪些方面| 最近在线观看免费完整版| 天天躁日日操中文字幕| or卡值多少钱| 国产成人福利小说| 小说图片视频综合网站| 精品不卡国产一区二区三区| 欧美不卡视频在线免费观看| 国产午夜精品论理片| 搞女人的毛片| 免费人成在线观看视频色| 校园春色视频在线观看| 在线观看午夜福利视频| x7x7x7水蜜桃| 一区二区三区国产精品乱码| 久久久久久九九精品二区国产| 草草在线视频免费看| 亚洲av日韩精品久久久久久密| 一区二区三区免费毛片| 桃色一区二区三区在线观看| 女人被狂操c到高潮| 欧美最黄视频在线播放免费| 激情在线观看视频在线高清| 亚洲精品一区av在线观看| 国产成年人精品一区二区| 色综合站精品国产| 午夜免费激情av| 中出人妻视频一区二区| 国产精品野战在线观看| 亚洲专区国产一区二区| 国内揄拍国产精品人妻在线| 最近最新中文字幕大全免费视频| 少妇的丰满在线观看| 国内精品久久久久精免费| av福利片在线观看| 国产不卡一卡二| 嫩草影院精品99| 最近在线观看免费完整版| av黄色大香蕉| av黄色大香蕉| 小蜜桃在线观看免费完整版高清| 不卡一级毛片| 国产成人系列免费观看| 国产黄片美女视频| 久久久久久久精品吃奶| 美女黄网站色视频| 男人舔奶头视频| 成人国产综合亚洲| 久久久精品大字幕| 成人18禁在线播放| aaaaa片日本免费| 欧美zozozo另类| 免费高清视频大片| 身体一侧抽搐| 免费在线观看成人毛片| 亚洲五月婷婷丁香| 亚洲欧美一区二区三区黑人| 99热只有精品国产| 国产精品久久久久久精品电影| 国产精品久久久久久精品电影| 国产精品免费一区二区三区在线| 国产精品av视频在线免费观看| 午夜免费观看网址| 亚洲国产欧洲综合997久久,| 国产一区二区激情短视频| 国产精品三级大全| 露出奶头的视频| 欧美成人性av电影在线观看| 国产精品久久视频播放| 久久久久国内视频| 久久久色成人| 欧美日韩亚洲国产一区二区在线观看| 淫妇啪啪啪对白视频| 亚洲av五月六月丁香网| 国产高清激情床上av| 欧美性猛交╳xxx乱大交人| 久9热在线精品视频| 国产av一区在线观看免费| 久久精品国产亚洲av涩爱 | 99久久综合精品五月天人人| 在线播放无遮挡| 成人高潮视频无遮挡免费网站| 久久欧美精品欧美久久欧美| 精品久久久久久久毛片微露脸| 国产精品美女特级片免费视频播放器| 国产日本99.免费观看| av天堂在线播放| 亚洲第一欧美日韩一区二区三区| 无人区码免费观看不卡| 国产成人影院久久av| 欧美成人a在线观看| 欧美成人a在线观看| 午夜福利视频1000在线观看| 一级作爱视频免费观看| 国产成人aa在线观看| 亚洲第一欧美日韩一区二区三区| 日韩欧美在线二视频| 亚洲人成网站高清观看| 日日夜夜操网爽| 免费av观看视频| 午夜免费男女啪啪视频观看 | 高清毛片免费观看视频网站| 欧美成狂野欧美在线观看| 精品无人区乱码1区二区| 欧美另类亚洲清纯唯美| 中文资源天堂在线| 又黄又爽又免费观看的视频| 全区人妻精品视频| av在线蜜桃| 精品电影一区二区在线| 俄罗斯特黄特色一大片| 老司机福利观看| 啦啦啦观看免费观看视频高清| 亚洲电影在线观看av| 亚洲激情在线av| 日本 欧美在线| 亚洲精品乱码久久久v下载方式 | 精品久久久久久久毛片微露脸| 在线观看日韩欧美| 久久久久久久久中文| 日韩精品青青久久久久久| 哪里可以看免费的av片| 毛片女人毛片| av专区在线播放| 日韩欧美 国产精品| 国产国拍精品亚洲av在线观看 | 亚洲国产日韩欧美精品在线观看 | 欧美日韩一级在线毛片| a级一级毛片免费在线观看| 久久久久久人人人人人| 国产成人福利小说| 国产中年淑女户外野战色| 久久久久久久久中文| 国产精品久久久久久亚洲av鲁大| 天美传媒精品一区二区| 欧美日韩亚洲国产一区二区在线观看| 久久香蕉精品热| 搡女人真爽免费视频火全软件 | 国产高潮美女av| 美女高潮喷水抽搐中文字幕| 成人性生交大片免费视频hd| 亚洲无线观看免费| 在线观看一区二区三区| 午夜免费激情av| 国产色婷婷99| av中文乱码字幕在线| 一个人免费在线观看的高清视频| 狂野欧美白嫩少妇大欣赏| 真人做人爱边吃奶动态| 亚洲精品日韩av片在线观看 | 国产成人啪精品午夜网站| 伊人久久精品亚洲午夜| 日韩 欧美 亚洲 中文字幕| 亚洲aⅴ乱码一区二区在线播放| 可以在线观看的亚洲视频| 可以在线观看毛片的网站| 不卡一级毛片| 欧美日韩中文字幕国产精品一区二区三区| 蜜桃亚洲精品一区二区三区| 麻豆成人av在线观看| 国产伦精品一区二区三区视频9 | 天堂影院成人在线观看| 99久久久亚洲精品蜜臀av| 国产欧美日韩精品亚洲av| 99在线视频只有这里精品首页| 日日夜夜操网爽| 久久精品人妻少妇| 日本免费a在线| 欧美极品一区二区三区四区| 午夜福利高清视频| 免费搜索国产男女视频| 一本精品99久久精品77| 国产午夜精品久久久久久一区二区三区 | 91久久精品国产一区二区成人 | 精品不卡国产一区二区三区| 国产亚洲欧美在线一区二区| 亚洲欧美精品综合久久99| 免费av观看视频| 美女免费视频网站| 高清毛片免费观看视频网站| 欧美日本视频| 欧美午夜高清在线| e午夜精品久久久久久久| 少妇高潮的动态图| 亚洲av免费在线观看| 黄片大片在线免费观看| 男人和女人高潮做爰伦理| 亚洲av免费在线观看| 国产午夜精品久久久久久一区二区三区 | 欧美黄色片欧美黄色片| 午夜亚洲福利在线播放| 免费在线观看成人毛片| 99热精品在线国产| 久久国产精品影院| 国产一区在线观看成人免费| 哪里可以看免费的av片| 一本久久中文字幕| 亚洲人成网站在线播| 18+在线观看网站| 日韩免费av在线播放| 国产不卡一卡二| 男女之事视频高清在线观看| 美女 人体艺术 gogo| 国产精品久久久久久久久免 | 九色成人免费人妻av| 最近最新免费中文字幕在线| 黄片大片在线免费观看| 99热6这里只有精品| 亚洲第一电影网av| 9191精品国产免费久久| 久久香蕉精品热| h日本视频在线播放| 国产高清视频在线播放一区| 少妇丰满av| 欧美性猛交╳xxx乱大交人| 嫁个100分男人电影在线观看| 女人被狂操c到高潮| 精品久久久久久久毛片微露脸| 老鸭窝网址在线观看| 亚洲人成伊人成综合网2020| 热99re8久久精品国产| 国产在视频线在精品| 亚洲一区二区三区不卡视频| 手机成人av网站| 男人的好看免费观看在线视频| 欧美成人a在线观看| av女优亚洲男人天堂| 欧美日韩亚洲国产一区二区在线观看| 婷婷丁香在线五月| 无人区码免费观看不卡| www国产在线视频色| 淫妇啪啪啪对白视频| 色综合婷婷激情| 少妇高潮的动态图| 国产三级中文精品| 老汉色av国产亚洲站长工具| 国产精品久久久久久人妻精品电影| 国产中年淑女户外野战色| 久久久国产成人精品二区| 18禁裸乳无遮挡免费网站照片| 久久亚洲真实| 99精品久久久久人妻精品| 午夜精品久久久久久毛片777| 久久久国产成人精品二区| 日日夜夜操网爽| 精品福利观看| 午夜免费成人在线视频| 哪里可以看免费的av片| 国内精品美女久久久久久| 成年女人看的毛片在线观看| 在线观看免费视频日本深夜| 欧美高清成人免费视频www| 亚洲无线观看免费| 此物有八面人人有两片| 色综合亚洲欧美另类图片| 最近最新中文字幕大全免费视频| 国产精品久久久久久久久免 | 精品一区二区三区人妻视频| 亚洲精品日韩av片在线观看 | 免费大片18禁| 午夜福利在线观看吧| 国产日本99.免费观看| 中文在线观看免费www的网站| 国产亚洲精品综合一区在线观看| 又紧又爽又黄一区二区| 少妇人妻一区二区三区视频| 午夜福利免费观看在线| 久99久视频精品免费| www国产在线视频色| 日日摸夜夜添夜夜添小说| 日本黄色视频三级网站网址| 一区二区三区高清视频在线| 男人舔女人下体高潮全视频| 久久精品综合一区二区三区| 久久性视频一级片| 一卡2卡三卡四卡精品乱码亚洲| 在线十欧美十亚洲十日本专区| 99久久久亚洲精品蜜臀av| 欧美黄色片欧美黄色片| 国模一区二区三区四区视频| 亚洲成a人片在线一区二区| 中文在线观看免费www的网站| 国语自产精品视频在线第100页| 午夜福利高清视频| 久久人人精品亚洲av| 亚洲av成人av| 毛片女人毛片| 熟女少妇亚洲综合色aaa.| 91九色精品人成在线观看| 一区二区三区高清视频在线| 黄色女人牲交| 两个人视频免费观看高清| 国产成人aa在线观看| 国产又黄又爽又无遮挡在线| 欧美不卡视频在线免费观看| 午夜福利18| 国产黄a三级三级三级人| 国模一区二区三区四区视频| 中出人妻视频一区二区| 亚洲国产欧美网| 色哟哟哟哟哟哟| 99热只有精品国产| 色精品久久人妻99蜜桃| 成年女人毛片免费观看观看9| 9191精品国产免费久久| 在线观看av片永久免费下载| 日韩欧美三级三区| 久久精品91无色码中文字幕| 校园春色视频在线观看| 国产熟女xx| 国内少妇人妻偷人精品xxx网站| 岛国视频午夜一区免费看| 桃红色精品国产亚洲av| svipshipincom国产片| 欧美绝顶高潮抽搐喷水| 成人av一区二区三区在线看| 91麻豆精品激情在线观看国产| 久久6这里有精品| 在线免费观看不下载黄p国产 | 久久精品国产亚洲av涩爱 | 亚洲乱码一区二区免费版| 精品免费久久久久久久清纯| 欧美中文综合在线视频| 91在线观看av| 欧美又色又爽又黄视频| 一本精品99久久精品77| av天堂中文字幕网| 日韩 欧美 亚洲 中文字幕| 国产成人影院久久av| 亚洲欧美日韩高清在线视频| 国产乱人视频| 免费大片18禁| 99热这里只有精品一区| 成人高潮视频无遮挡免费网站| 欧美av亚洲av综合av国产av| 听说在线观看完整版免费高清| 中文字幕精品亚洲无线码一区| 免费av观看视频| 亚洲欧美日韩无卡精品| 日本一二三区视频观看| 中文字幕av在线有码专区| 热99在线观看视频| 亚洲成人久久爱视频| 一个人看的www免费观看视频| 亚洲国产欧洲综合997久久,| 日本黄色片子视频| 成人特级黄色片久久久久久久| 亚洲性夜色夜夜综合| 久久久久精品国产欧美久久久| 久久久久久九九精品二区国产| 亚洲aⅴ乱码一区二区在线播放| 啦啦啦观看免费观看视频高清| 色哟哟哟哟哟哟| 国产成人福利小说| 操出白浆在线播放| 国产精品99久久久久久久久| 男女下面进入的视频免费午夜| 欧美成人一区二区免费高清观看| 亚洲人成网站高清观看| 12—13女人毛片做爰片一| 99国产综合亚洲精品| 色视频www国产| 国产精品国产高清国产av| 国产蜜桃级精品一区二区三区| 在线观看午夜福利视频| 久久久精品大字幕| 成人国产一区最新在线观看| 丰满人妻一区二区三区视频av | 日韩欧美免费精品| 亚洲av免费高清在线观看| 免费在线观看日本一区| 波多野结衣高清无吗| 久久中文看片网| 中亚洲国语对白在线视频| 久久人妻av系列| 国内少妇人妻偷人精品xxx网站| 免费电影在线观看免费观看| 18禁美女被吸乳视频| 黄片大片在线免费观看| www国产在线视频色| av天堂中文字幕网| 麻豆国产97在线/欧美| 久久草成人影院| 黄色成人免费大全| 99在线人妻在线中文字幕| 欧美丝袜亚洲另类 | 午夜福利欧美成人| 国产亚洲av嫩草精品影院| 欧美高清成人免费视频www| 国产野战对白在线观看| 日本免费a在线| 99热这里只有精品一区| 国语自产精品视频在线第100页| 一区福利在线观看| 久久久色成人| 国产成+人综合+亚洲专区| 最新中文字幕久久久久| 男人的好看免费观看在线视频| 看片在线看免费视频| 人妻夜夜爽99麻豆av| 亚洲精品色激情综合| 1024手机看黄色片| 成人精品一区二区免费| av中文乱码字幕在线| 岛国视频午夜一区免费看| 午夜激情福利司机影院| 搞女人的毛片| 久久婷婷人人爽人人干人人爱| 91久久精品国产一区二区成人 | 亚洲国产日韩欧美精品在线观看 | 亚洲最大成人手机在线| 国产色婷婷99| 国产野战对白在线观看| 欧美一区二区国产精品久久精品| 两人在一起打扑克的视频| 黄色成人免费大全| 国产精品,欧美在线| 国产成人av教育| 成人永久免费在线观看视频| 天堂√8在线中文| 午夜福利在线观看吧| 美女被艹到高潮喷水动态| 色在线成人网| 久久亚洲真实| 精品国内亚洲2022精品成人| 精品人妻1区二区| 男女视频在线观看网站免费| av专区在线播放| 色噜噜av男人的天堂激情| 18禁裸乳无遮挡免费网站照片| 日本一本二区三区精品| 久久国产乱子伦精品免费另类| 日本一二三区视频观看| 伊人久久精品亚洲午夜| 国内揄拍国产精品人妻在线| 久9热在线精品视频| 亚洲欧美日韩东京热| 亚洲国产欧洲综合997久久,| 99久久精品国产亚洲精品| 动漫黄色视频在线观看| 69人妻影院| 啦啦啦韩国在线观看视频| av在线蜜桃| 欧美性猛交╳xxx乱大交人| 亚洲男人的天堂狠狠| 99国产精品一区二区三区| 美女黄网站色视频| 人妻久久中文字幕网| 1000部很黄的大片| 嫩草影院精品99| 免费无遮挡裸体视频| 99精品欧美一区二区三区四区| 内地一区二区视频在线| 亚洲专区中文字幕在线| 国产精品 国内视频| 午夜免费男女啪啪视频观看 | 国产av麻豆久久久久久久| 久久婷婷人人爽人人干人人爱| 男人的好看免费观看在线视频| 欧美日本视频| 国产精品,欧美在线| 国产成人欧美在线观看| 精品人妻偷拍中文字幕| 一区二区三区免费毛片| 久久精品国产综合久久久| 国产精品,欧美在线| 首页视频小说图片口味搜索| 国产极品精品免费视频能看的| 超碰av人人做人人爽久久 | 免费看美女性在线毛片视频| 在线播放国产精品三级| 精品电影一区二区在线| 免费人成视频x8x8入口观看| 中文字幕久久专区| 久久亚洲精品不卡| 国产麻豆成人av免费视频| 老熟妇乱子伦视频在线观看| 床上黄色一级片| 老司机福利观看| 日韩欧美 国产精品| 岛国在线免费视频观看| 热99re8久久精品国产| 精品久久久久久久久久久久久| 悠悠久久av| 极品教师在线免费播放| 欧美3d第一页| 久久香蕉精品热| 最近最新免费中文字幕在线| 特级一级黄色大片| 99久久久亚洲精品蜜臀av| 国内毛片毛片毛片毛片毛片| 级片在线观看| 一区福利在线观看| 亚洲人成网站在线播| avwww免费| 日韩欧美在线二视频| 黄色女人牲交| 在线播放国产精品三级| 亚洲激情在线av| 国产aⅴ精品一区二区三区波| 日韩欧美三级三区| 黑人欧美特级aaaaaa片| 国产精品一区二区三区四区久久| 国产成人av激情在线播放| 天堂√8在线中文| 日日摸夜夜添夜夜添小说| 久久久久久大精品| 91在线观看av| 午夜亚洲福利在线播放| 成人国产综合亚洲| 最后的刺客免费高清国语| 国内精品美女久久久久久| 欧美中文日本在线观看视频| 天堂动漫精品| 日本熟妇午夜| 成人av在线播放网站| 国产精品香港三级国产av潘金莲| a级一级毛片免费在线观看| 母亲3免费完整高清在线观看| 免费在线观看影片大全网站| 欧美另类亚洲清纯唯美| 色综合站精品国产| 国产精品99久久99久久久不卡| 黄色片一级片一级黄色片| 99久久精品热视频| 亚洲国产日韩欧美精品在线观看 | 国产午夜福利久久久久久| 夜夜夜夜夜久久久久| 亚洲av美国av| 内射极品少妇av片p| 日韩亚洲欧美综合| 男女下面进入的视频免费午夜| 真人做人爱边吃奶动态| 免费看十八禁软件| 村上凉子中文字幕在线| 欧美一级a爱片免费观看看| 制服丝袜大香蕉在线| 麻豆久久精品国产亚洲av| 一卡2卡三卡四卡精品乱码亚洲| 成年女人永久免费观看视频| 亚洲专区国产一区二区| 桃色一区二区三区在线观看| 日本精品一区二区三区蜜桃| 久久久久国产精品人妻aⅴ院| 操出白浆在线播放| 色综合亚洲欧美另类图片| 香蕉久久夜色| 欧美色视频一区免费| 久久精品夜夜夜夜夜久久蜜豆| 在线国产一区二区在线| 男插女下体视频免费在线播放| 亚洲精品日韩av片在线观看 | 18禁国产床啪视频网站| 国产免费男女视频| 亚洲国产色片| 久久精品夜夜夜夜夜久久蜜豆| 亚洲国产欧美网| 免费观看的影片在线观看| 国产真人三级小视频在线观看| 亚洲av中文字字幕乱码综合| 午夜老司机福利剧场| 国产三级在线视频| 看免费av毛片| 国产精品久久视频播放| 村上凉子中文字幕在线| 啦啦啦韩国在线观看视频| 久久天躁狠狠躁夜夜2o2o| 亚洲乱码一区二区免费版| 免费观看精品视频网站| 亚洲精华国产精华精| 午夜福利18| 99精品在免费线老司机午夜| 中国美女看黄片| 少妇人妻一区二区三区视频| av专区在线播放| 国产成人影院久久av| av国产免费在线观看| 久久久久久久久中文| 毛片女人毛片| 精品国产三级普通话版| 亚洲无线观看免费| 成人av在线播放网站| 久久久久久久久大av| 在线十欧美十亚洲十日本专区| 热99re8久久精品国产| 欧美黑人欧美精品刺激| 欧美激情久久久久久爽电影| 久久久久性生活片| 国产av不卡久久| 在线观看免费视频日本深夜|