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

    DMMIC derivatization-assisted liquid chromatography-mass spectrometry method for metabolite profiling of the glutathione anabolic pathway in esophageal cancer tissues and cells

    2023-12-14 10:54:58LiLiuYuHnLuMinDnWngQunFeiZhoXiuPingChenHngYinChenGuoFengFngZhng
    Journal of Pharmaceutical Analysis 2023年11期

    Li Liu ,Yu-Hn Lu ,Min-Dn Wng ,Qun-Fei Zho ,Xiu-Ping Chen ,*** ,Hng Yin ,Chen-Guo Feng ,**,Fng Zhng ,*

    a The Research Center of Chiral Drugs, Shanghai Frontiers Science Center for TCM Chemical Biology, Innovation Research Institute of Traditional Chinese Medicine, Shanghai University of Traditional Chinese Medicine, Shanghai, 201203, China

    b Department of Thoracic Surgery, Renji Hospital, School of Medicine, Shanghai Jiaotong University, Shanghai, 200127, China

    c School of Public Health, Shanghai University of Traditional Chinese Medicine, Shanghai, 201203, China

    Keywords:Glutathione anabolic pathway Metabolite profiling DMMIC derivatization LC-MS Esophageal squamous cell carcinoma p-Hydroxycinnamaldehyde KYSE-150cell

    ABSTRACT In this work,a new pyrylium derivatization-assisted liquid chromatography-mass spectrometry(LC-MS)method was developed for metabolite profiling of the glutathione anabolic pathway (GAP) in cancer tissues and cells.The pyrylium salt of 6,7-dimethoxy-3-methyl isochromenylium tetrafluoroborate(DMMIC) was used to label the amino group of metabolites,and a reductant of dithiothreitol (DTT) was employed to stabilize the thiol group.By combining DMMIC derivatization with LC-MS,it was feasible to quantify the 13 main metabolites on the GAP in complex biological samples,which had good linearity(R2=0.9981-0.9999),precision(interday precision of 1.6%-19.0%and intraday precision of 1.4%-19.8%)and accuracy (83.4%-115.7%).Moreover,the recovery assessments in tissues (82.5%-107.3%) and in cells(98.1%-118.9%)with GSH-13C2,15N,and Cys-15N demonstrated the reliability of the method in detecting tissues and cells.Following a methodological evaluation,the method was applied successfully to investigate difference in the GAP between the carcinoma and para-carcinoma tissues of esophageal squamous cell carcinoma(ESCC)and the effect of p-hydroxycinnamaldehyde(CMSP)on the GAP in KYSE-150 esophageal cancer cells.The results demonstrate that the developed method provides a promising new tool to elucidate the roles of GAP in physiological and pathological processes,which can contribute to research on drugs and diseases.

    1.Introduction

    Glutathione is an endogenous component of cellular metabolism,and is a tripeptide composed of glycine,cysteine,and glutamic acid [1].Many studies have shown that glutathione is involved in various biological functions due to its anabolic pathway[2-4];in addition,glutathione is associated with the progression of many human diseases,including aging diseases[5,6],cystic fibrosis[7,8],cardiovascular diseases [9],inflammatory diseases [10,11],immune diseases [12],metabolic diseases [13],neurodegenerative diseases[14-16],and cancers[17,18].In particular,the glutathione anabolic pathway (GAP) plays a dual role in the development of cancers[19,20].GAP can remove and detoxify carcinogens,but high expression of GAP may protect cancer cells by conferring resistance to some drugs [21,22].Targeting GAP is a potential strategy to induce cancer cell death [22,23].Therefore,monitoring the metabolite profiling of GAP is of great importance and contributes to research on drugs and diseases.

    There are 13 main metabolites in the GAP (Fig.1),namely,glutamine (Gln),glutamic acid (Glu),glycine (Gly),serine (Ser),methionine(Met),S-adenylyl homocysteine(SAH),O-acetyl-serine(OAS),homocysteine (HCys),cystathionine (Cysta),cysteine (Cys),cysteinylglycine (CG),glutamylcysteine (GC),and glutathione(GSH).In terms of chemical structure,some of these metabolites contain thiol groups while others contain small-molecule oligopeptides;thus,analysis can be difficult because these compounds exhibit chemical instability,are highly polar,are not very abundant and easily undergo enzymatic hydrolysis.Thus,many efforts have been dedicated to achieving efficient,sensitive,and accurate analysis[24].Among them,the enzymatic recycling assay[25]is a classic approach that can quantitatively determine glutathione.However,the amount of enzyme used is large and costly,and its good specificity cannot be used to simultaneously detect other metabolites in the pathway.As a new technology,fluorescent probes have also been reported to be an ideal tool for the real-time monitoring of glutathione [26,27].However,similar to the enzymatic recycling assay,fluorescent probes cannot cover metabolites with different structures in the pathway[28].Compared with these methods,liquid chromatography-mass spectrometry (LC-MS) is a preferred tool for metabolite profiling due to its good selectivity,sensitivity,and high-coverage [29,30].In particular,with the assistance of chemical derivatization [31,32],LC-MS has greatly improved the analysis of compounds that contain unstable functional groups,exhibit large polarity,show no chromatographic retention,and exhibit poor MS ionization.

    Fig.1.Glutathione anabolic pathway (GAP).Gln: glutamine;Ser: L-serine;Gly: L-glycine;GC: glutamylcysteine;SAH: S-adenylyl homocysteine;Glu: L-glutamine;CG: cysteinylglycine;OAS:O-acetyl-serine;Cys: L-cysteine;Cysta:cystathionine;Hcys:homocysteine;Met: L-methionine;GSH:glutathione;GTT:γ-glutamyl transferase;GSS:glutathione synthetase;GCLM: glutamate-cysteine ligase modifier subunit gene;GCLC: glutamate-cysteine ligase catalytic subunit gene.

    In previous studies,it was demonstrated that an effective pyrylium salt [33] named 6,7-dimethoxy-3-methyl isochromenylium tetrafluoroborate(DMMIC)could greatly improve the quantification of amino compounds (ACs) when combined with LC-MS (Fig.2).Considering that the 13 metabolites involved in GAP each contain an amino group,the aim of this study is to exploit the advantages of DMMIC combined with LC-MS to realize metabolite profiling.In addition,to stabilize the metabolites containing the thiol group,a reductant of dithiothreitol (DTT) was added to the sample pretreatment;the power and effect of DTT on the subsequent derivation and detection were investigated.Finally,the developed method was evaluated using a systematic methodology and applied to investigate the difference in the GAP between carcinoma and para-carcinoma tissues of esophageal squamous cell carcinoma(ESCC).Furthermore,the effect ofp-hydroxycinnamaldehyde(CMSP)on the GAP in KYSE-150 esophageal cancer cells was evaluated.

    Fig.2.Reaction of 6,7-dimethoxy-3-methyl isochromenylium tetrafluoroborate(DMMIC) with amines.

    2.Materials and methods

    2.1.Reagents

    The synthesis of [d0]-/[d3]-DMMIC was performed in-house according to a previously reported method [33].L-glycine (Gly),glutamine (Gln),and L-glutamine (Glu) were purchased from Aladdin (Shanghai,China).L-serine (Ser),L-methionine (Met),formic acid,L-cysteine-15N (Cys-15N) and glutathione-13C2,15N(GSH-13C2,15N) sodium salt were purchased from J&K Scientific(Beijing,China).Hydrochloric acid (HCl),L-cysteine (Cys) and DL-homocysteine (HCys) were purchased from Sinopharm Group Co.,Ltd.(Shanghai,China).p-Hydroxycinnamaldehyde(CMSP)was purchased from JiuDing (Shanghai,China).Cysteinylglycine (CG),glutamylcysteine(GC),and glutathione(GSH)were purchased from BiDe (Shanghai,China).O-acetyl-serine (OAS) and cystathionine(Cysta) were purchased from Macklin (Shanghai,China).2-Chlorophenylalanine,dichloromethane,and phosphate buffered solution (PBS) solution were purchased from Adamas Life(Shanghai,China).Pyridine was purchased from Acros Organics(Geel,Belgium).Dithiothreitol (DTT) was procured from Beyotime(Shanghai,China)and the Cell Counting Kit-8(CCK8)was procured from MeiLun Biotechnology (Shanghai,China).Acetonitrile and methanol of HPLC grade were purchased from Merck KGaA(Darmstadt,Germany).Deionized water was produced by a Millipore Direct Q water purification system(El Paso,TX,USA).All other chemicals and solvents used were of analytical grade.Stock solutions of standards were prepared in deionized water with 20 mM DTT at a concentration of 0.1 mol/L.Working solutions were obtained by mixing and diluting the stock solutions.

    2.2.Collection and extraction of tissues and cells

    2.2.1.Tissues

    Tissue samples were collected from the esophagus of 10 patients with esophageal squamous at Shanghai Renji Hospital,conducted with the approval of the corresponding regulatory agencies(Approval number: No.2021-931-06-01).The patients were informed of the objectives of the project.The tissue samples were collected immediately after surgery,frozen in liquid nitrogen,and stored at -80°C until analysis.Before DMMIC derivatization,the frozen tissue samples were washed with ice-cold PBS solution and homogenized at 2 μL/mg in an aqueous solution containing 20 mM DTT and 10 μM 2-chlorophenylalanine.Subsequently,the homogenate was centrifuged for 3 min at 10,000 rpm,and 3 μL/mg extract was added and homogenized again.The supernatant was diluted with 4 vol of methanol to precipitate the protein.After the mixture was frozen overnight at-20°C,the supernatant was transferred for derivatization.

    2.2.2.Cells

    KYSE-150 cells (4 × 105) were seeded in a 6-well plate.Cells were divided into the following groups:the Con and CMSP groups.After incubation overnight,the cells of the CMSP group were treated with 2 mL of 200 μM CMSP solution in culture medium.Blank culture medium was added to the control group.Twenty-four hours later,the culture medium was fully aspirated,and the cells were washed twice with 1 mL of PBS solution.Then,the cells were digested with 500 μL of trypsin and washed with PBS three times.Fifty microliters of 80%methanol aqueous solution containing DTT and 2-chlorophenylalanine was quickly added as the extraction and quenching solvent to the cell precipitates.The samples were repeatedly subjected to freeze-thaw cycles in liquid nitrogen and ultrasonicated 3 times,and then the extracts were stored at-80°C until metabolic profile analysis.

    2.3.DMMIC derivatization

    For DMMIC derivatization,a 20 μL aliquot of the standard solution or cell and tissue extract was added to 10 μL of[d0]-or[d3]-DMMIC reagent(250 mM)and 70 μL of pyridine in methanol(1%).The mixture was vortexed and incubated at 50°C for 1 h at 1000 rpm.Then,the mixture was dried by a stream of nitrogen,and the residue was dissolved in 100 μL of an aqueous solution containing 10 mM DTT and 1% HCl.The [d0]-/[d3]-DMMIC derivatized samples were mixed at a ratio of 1:1 (V/V).The solution was centrifuged for 3 min at 10,000 rpm,and the supernatant was used for LC-MS analysis.

    2.4.LC-MS analysis

    The liquid chromatography quadrupole time-of-flight mass spectrometry(LC-QTOF MS)measurements were performed on an instrument consisting of a 1290 HPLC system and a 6545 UHD quadrupole time of flight mass spectrometer(Agilent Technologies,Santa Clara,CA,USA).LC separation was carried out on a Waters XSelect HSS T3 column (2.1 mm × 150 mm,2.5 μm;Waters,Milford,MA,USA).The mobile phase was composed of 0.4%formic acid in water (A) and MeOH (B).The profile of gradient elution was as follows:0-1 min,90%A;1-2 min,90%-70%A;2-7 min,70%-65%A;7-10 min,65%-5%A;10-13.5 min,5%-5%A;13.5-14 min,5%-90% A;and 13.5-18 min,90% A.The column temperature was maintained at 40°C,the flow rate was 0.3 mL/min,and the sample injection volume was 1 μL.The QTOF MS was operated with a jet stream electrospray ionization (ESI) source in positive mode.The ESI capillary voltage was set to 2500 V,the nozzle voltage was set to 0 V,and the fragment voltage was 150 V.The nitrogen drying gas flow was heated to 325°C at a flow rate of 13 L/min.The sheath gas temperature was 350°C at a flow rate of 12 L/min.The nebulizer was 40 psig.The mass scan range of the primary MS wasm/z100-1000.Before measurements were performed,the QTOF MS instrument was calibrated externally using a series of standards(m/z100 -3200) to ensure that the mass accuracy was less than 2 ppm and the mass resolution was 25,000 (m/z322).

    2.5.Data processing and statistical analysis

    Agilent MassHunter Acquisition B.08.00 (Agilent Technologies)was used for the control of the equipment and acquisition of data.MassHunter Qualitative Analysis Navigator B.08.00 and Mass-Hunter Quantitative Analysis B.08.00 (Agilent Technologies) were used for data processing.Orthogonal partial least squaresdiscriminant analysis (OPLS-DA),boxplot,heatmap and receiver operating characteristic (ROC) curves were generated using MetaboAnalyst 5.0 (http://www.metaboanalyst.ca/).All quantitative data in this study are expressed as the mean ± standard deviation (SD).Differences were evaluated by an unpairedt-test.P<0.05 was recognized as statistically significant.

    3.Results and discussion

    3.1.Optimization of sample pretreatment

    The stability of target metabolites is an important consideration when establishing a method to profile GAP.Cys,HCys,and other thiol metabolites are easily oxidized,while GC,CG,and GSH are vulnerable to hydrolysis by enzymes.Thus,the method used to pretreat thawed tissues was investigated at the beginning of this work;for example,the effect of freeze-drying and the addition of a reductant of DTT were examined [34].As shown in Fig.3,between freeze-drying groups(Figs.3A and B),HCys,CG,GC,and SAH were successfully detected,and the signal of Cys was increased 2-3 orders of magnitude when the samples was homogenized with DTT aqueous solution compared to homogenizationwith water.This indicated that the addition of DTT exerted a positive effect,which reduced the disulfide bonds and ensured that the thiol groups were in a stable state,as previously reported [35,36].However,Cysta and GSH were not detected regardless of whether DTT was added during homogenization.The molecular structure of the two metabolites might undergo irreversible degradation during the freeze-drying process.When thawed tissues were not freeze-dried but subjected to subsequent treatments directly(Figs.3C and D),all 13 target metabolites were successfully detected in the presence of DTT,especially Cysta and GSH.Therefore,directly homogenizing tissue samples with an aqueous solution containing DTT after thawing is the preferred pretreatment method for profiling the GAP.Normally,freeze-drying tissue can effectively eliminate the interference of its moisture on subsequent quantitative results.Considering this factor,an internal standard of 2-chlorophenylalanine was added during homogenization and its signal was used to normalize the target signals among different tissue samples during data processing.

    Fig.3.Chromatograms for the metabolites in thawed tissue pretreated by (A) freeze-drying followed by homogenization with dithiothreitol (DTT) aqueous solution,(B) freezedrying followed by homogenization with water,(C) homogenization with DTT aqueous solution,and (D) homogenization with water.Gln: glutamine;Ser: L-serine;Gly: Lglycine;GC: glutamylcysteine;SAH: S-adenylyl homocysteine;Glu: L-glutamine;CG: cysteinylglycine;OAS: O-acetyl-serine;Cys: L-cysteine;Cysta: cystathionine;Hcys: homocysteine;Met: L-methionine;GSH: glutathione.

    3.2.Optimization of DMMIC derivatization

    Because of the reductant DTT was used in the pretreatment,the effect of its presence on DMMIC derivatization was investigated.Five sulfhydryl-free metabolites were used as the test compounds;these metabolites were dissolved in DTT aqueous solution and water,and then underwent derivatization.Surprisingly,all DMMIC derivatives from the DTT group were higher than those from the water group,and Gln,Ser,Gly,Glu and Met were increased by 33%,31%,23%,25% and 51%,respectively (Fig.4).This result indicated that the presence of DTT exerted a certain positive impact on the derivative reaction.It was speculated that DTT might have a catalytic effect on the reaction,improving the reaction efficiency by forming active intermediates with the DMMIC reagent.

    Fig.4.Comparison of the derivatized metabolites in the 6,7-dimethoxy-3-methyl isochromenylium tetrafluoroborate (DMMIC) reaction using dithiothreitol (DTT) solution and water as solvents.Gln: glutamine;Ser: L-serine;Gly: L-glycine;Glu: Lglutamine;Met: L-methionine.

    Fig.5.The multi-extracted ionchromatograms(mEICs,m/z at molecular masses)of 6,7-dimethoxy-3-methyl isochromenylium tetrafluoroborate(DMMIC)-derivatized metabolites detected by liquid chromatography-mass spectrometry (LC-MS) from (A) esophageal cancer tissue and (B) KYSE-150 cells.Gln: glutamine;Ser: L-serine;Gly: L-glycine;GC: glutamylcysteine;SAH: S-adenylyl homocysteine;Glu: L-glutamine;CG: cysteinylglycine;OAS: O-acetyl-serine;Cys: L-cysteine;Cysta: cystathionine;Hcys: homocysteine;Met: Lmethionine;GSH: glutathione.

    It was previously demonstrated that the following key factors affected the efficiency of the reaction: reaction temperature,reaction time,and the amount of pyridine in the DMMIC derivatization[33].Considering the addition of DTT in this study,the reaction conditions were re-optimized by varying the concentrations of pyridine methanol solution(0,0.5%,1%and 1.5%),reaction times(1,2,3 and 4 h),and reaction temperatures (40,50,60 and 70°C).Overall,the test metabolites tended to increase and then decrease with prolonged reaction time,especially for Gly and Cys.The reaction time in the range evaluated showed no obvious effect on the DMMIC derivative.The addition of pyridine did benefit the reaction,but no effect was observed within 0.5%-1.5% pyridine methanol solution.Finally,a 1%pyridine methanol solution was selected.The reaction temperature and time were maintained at 50°C and 1 h,respectively (Fig.S1).The results obtained indicate that the DMMIC derivatization reaction was achieved under milder conditions with a faster kinetic reaction.After the reaction,the solution was dried with nitrogen gas to remove the pyridine,preventing the ionization of LC-MS detection from incurring interferences.For the residue,an aqueous solution containing 10 mM DTT and 1%hydrochloric acid was the preferred redissolved solvent(Fig.S2),which helped stabilize the DMMIC derivatives in the form of organic salt and replenish the DTT lost during drying with nitrogen.In addition,it was observed that extracting the redissolved solution with dichloromethane helped to remove the liposoluble impurities caused by the biologic matrix and DMMIC derivatization and increased the signals of the derivatized metabolites.

    3.3.Methodology evaluation of DMMIC derivatization-assisted LCMS

    3.3.1.Linearity,limit of detection,and limit of quantification

    Before the methodology evaluation,the composition of the mobile phase was evaluated for the LC-MS detection of DMMICderivatized metabolites.The results showed that 0.4% formic acid water and methanol could elute the derivatized metabolites well(Figs.S3 and S4).As shown in Figs.5 and S5,under the optimal conditions,all test derivatized metabolites in tissues and cells exhibited satisfactory chromatographic behavior with good peak shapes and chromatographic retentions.

    To study the feasibility of the method for linearity,a series of 12 mixed standard solutions of metabolites (SAH was uncommerciable) with different concentrations were prepared by[d0]-DMMIC labeling and detected by LC-MS.The [d3]-DMMIC derivatized standard solutions were used as internal standards(ISs).The calibration curve was constructed by plotting the ratio of the abundance of [d0]-DMMIC-derivatized samples to that of IS at different concentrations.As shown in Fig.6 and Table S1,all metabolites were confirmed to exhibit good linearity (R2=0.9981-0.9999) and accuracy (83.4%-115.7%) in a concentration range of 0.05-1 μM.As the peak signal consists of at least 7 and 15 MS spectra,the concentration corresponding to the signal-to-noise(S/N)ratio ≥3 and 10 peak signals was used as the limit of detection(LOD)and limit of quantification(LOQ),respectively.The LODs and LOQs of the compounds are shown in Table S1.The LOQs were 0.5 nM for Gly,Gln,Ser and Cys,1 nM for OAS;2 nM for CG,5 nM for Glu and HCys,and 10 nM for Met,and the LOQs of GSH,Cysta and GC were 50 nM.Inevitably,SAH was omitted here in the methodology evaluation due to the lack of reference standard,but assessment of the other 12 target metabolites,especially HCys owning a certain structural similarity to SAH,could demonstrate whether the developed method meets the quantitative requirements and is used for metabolite profiling of GAP including the relative quantitative analysis of SAH.

    Fig.6.Calibration curves of the 6,7-dimethoxy-3-methyl isochromenylium tetrafluoroborate (DMMIC)-derivatized metabolites.Gln: glutamine;Ser: L-serine;Gly: L-glycine;GC:glutamylcysteine;Glu: L-glutamine;CG: cysteinylglycine;OAS: O-acetyl-serine;Cys: L-cysteine;Cysta: cystathionine;Hcys: homocysteine;Met: L-methionine;GSH: glutathione.

    3.3.2.Stability,precision,and accuracy

    Additional work,including stability,precision,and accuracy,was performed to validate the method.The stability of the DMMICderivatized metabolites was investigated by using the[d3]-DMMICderivatized metabolites as the ISs of the [d0]-DMMIC-derivatized.As shown in Fig.S6,the tested DMMIC-derivatized metabolites exhibited satisfactory stability over 48 h.The inter-and intra-day precision was evaluated by examining mixed standard solutions at low,medium,and high concentrations.The inter-day precision was 1.6%-19.0%,and the intra-day precision was 1.4%-19.8%(Table S1).

    To verify the accuracy of the metabolic profile analysis method in tissue and cell samples,two recovery assessments were investigated using GSH-13C2,15N and Cys-15N(for details,see supporting information) at low,medium and high concentrations (100,500 nM,and 1 μM).As shown in Table 1 and Fig.S7,the recoveries at different concentrations ranged from 82.5% to 107.3% for tissue and 98.1%-118.9% for cells.These results indicated that the developed method could be used to perform accurate quantitative analysis of metabolites in tissue and cell samples.

    Table 1 Recovery assessment of 6,7-dimethoxy-3-methyl isochromenylium tetrafluoroborate(DMMIC)-derivatized metabolites in tissues and cells.

    3.4.Applications of metabolite profiling in tissues and cells

    3.4.1.Difference in the GAP between the carcinoma and paracarcinoma tissues of ESCC

    Esophageal cancer is among the most common cancers worldwide,causing hundreds of thousands of deaths each year [37].Studies have shown that the concentration of glutathione in esophageal cancer is abnormal[38].To investigate the difference in the GAP between the carcinoma and para-carcinoma tissues of ESCC,thirteen metabolites were detected,and their relative quantification was carried out (Fig.S8).

    In the OPLS-DA model,significant differences were found in the GAP between the carcinoma and para-carcinoma tissues of ESCC with variable importance in projection(VIP)values of 7 metabolites(Glu,Cys,Gly,HCys,Ser,CG,and Met)greater than 1(Fig.7A).The heatmap showed that GAP was highly expressed in carcinoma tissues(Fig.7B).Moreover,the 7 metabolites with significant differences were confirmed using box diagrams(Fig.S9).Except for CG,the expression of the remaining 6 metabolites in cancer tissues was significantly higher thanthatinpara-cancerous tissues.ROCcurve analysis showed that the area under curve (AUC) values of the six metabolites were greater than 0.7,which further indicated that Gly,Glu,Met,Ser,HCys,andCysweresignificantlyup-regulatedincarcinomatissues(Fig.S10).

    Fig.7.Difference in the glutathione anabolic pathway (GAP) between the carcinoma and para-carcinoma tissues of esophageal squamous cell carcinoma (ESCC).(A) Score plot of orthogonal partial least squares-discriminant analysis(OPLS-DA)of the carcinoma and para-carcinoma tissues of esophageal squamous cell carcinoma.(B)Heatmap of the relative abundances of the metabolites in each individual patient and hierarchical clustering.Gln: glutamine;Ser: L-serine;Gly: L-glycine;GC: glutamylcysteine;SAH: S-adenylyl homocysteine;Glu: L-glutamine;CG: cysteinylglycine;OAS: O-acetyl-serine;Cys: L-cysteine;Cysta: cystathionine;Hcys: homocysteine;Met: L-methionine;GSH: glutathione.

    3.4.2.Effect of CMSP on GAP in KYSE-150 esophageal cancer cells

    Momordica cochinchinensis(Lour.)Spreng.,commonly known as Gac fruit,is a perennial dioecious cucurbit plant that originated in South and Southeast Asia and is widely sold for dietary and medicinal purposes.Reportedly,Momordica cochinchinensiscan effectively inhibit the growth of cancer cells,such as gastric cancer,liver cancer,and esophageal cancer cells [39].CMSP,as an important active component ofMomordica cochinchinensis,has been reported to affect NF-κB signaling pathways [40],which are closely related to the expression of glutathione.Thus,we conducted a study on the effect of CMSP intervention on GAP in esophageal cancer cells (Figs.S11 and S12).

    As shown in Fig.8,the OPLS-DA model was used to analyze the administration group and the control group,and an obvious separation trend was found among different groups.The VIP values >1 of GSH and GC were significantly different metabolites in this model.The results of the box diagram showed that the expression levels of GSH and GC in esophageal squamous cells after administration were higher than those in the control group.According to ROC curve analysis,the AUC values of the two compounds were also greater than 0.7,and the AUC value of GC even reached 0.97.Although CMSP can interfere with the production of GSH and GC in cells,the experimental results showed that CMSP exhibits a more significant effect on the expression of GC.

    Fig.8.Effect of CMSP intervention on glutathione anabolic pathway(GAP)in esophageal cancer cells.(A)Orthogonal partial least squares-discriminant analysis(OPLS-DA)scores of the cell administration group and control group.(B) Box plots,and (C) receiver operating characteristic (ROC) curves of the two differential metabolites.CON: control;CMSP: phydroxycinnamaldehyde;GC: glutamylcysteine;GSH: glutathione;AUC: area under curve.

    To verify the accuracy of the assessment,differences in the expression of genes regulating glutathione synthesis and metabolism,including glutamate-cysteine ligase catalytic subunit gene(GCLC),glutamate-cysteine ligase modifier subunit gene (GCLM),glutathione synthetase (GSS),and γ-glutamyl transferase (GTT),were detected by reverse-transcription-polymerase chain reaction(RT-PCR).The GCLC and GCLM genes controlling GC synthesis were significantly up-regulated under CMSP intervention,while GSS and GTT were not significantly changed.This was consistent with the results obtained from the glutathione metabolic profile(Fig.S13).When all the data was integrated,it was observed that CMSP could affect the expression of the GCLC and GCLM genes in esophageal cancer cells and thus significantly affect the concentration of GC.As a key precursor compound of GSH synthesis,GC could also promote the up-regulation of intracellular GSH concentration when the expression of the GSS gene did not change significantly.In addition,the GTT gene,which controls the enzymatic hydrolysis of GSH,was not significantly changed after administration and might also cause GSH to accumulate intracellularly.

    4.Conclusion

    In this study,an effective DMMIC derivatization-assisted LC-MS method was developed for metabolite profiling of GAP in cancer tissues and cells.DMMIC was used as a reagent to label the amino group of metabolites and DTT was used as a reductant to stabilize the thiol group.When DMMIC was combined with LC-MS,the 13 main metabolites on the GAP in complex biological samples could be feasibly quantified.Through systematic methodology evaluation,the developed method was demonstrated to exhibit good linearity(R2=0.9981-0.9999),precision (interday precision of 1.6%-19.0%and intraday precision of 1.4%-19.8%),and accuracy(83.4%-115.7%).Moreover,recovery assessments in tissues (82.5%-107.3%) and in cells(98.1%-118.9%)with GSH-13C2,15N and Cys-15N demonstrated the reliability of the method in the detection of biological samples.

    The method was applied to investigate differences in the GAP between the carcinoma and para-carcinoma tissues of ESCC,as well as the effect of CMSP on the pathway in KYSE-150 esophageal cancer cells.In the former case,the amounts of Glu,Cys,Gly,HCys,Ser,and Met,among the 13 main metabolites of GAP in the carcinoma tissues of ESCC were significantly higher than those in the para-carcinoma tissues,which indicated that upstream glutathione synthesis was activated in the cancer tissue.In the latter case,CMSP affected the expression of GCLC and GCLM genes in KYSE-150 esophageal cancer cells and thus interfered with the synthesis of GC in cells.As a key precursor compound for GSH synthesis,GC can also promote the upregulation of intracellular GSH concentration without significant changes in GSS gene expression.As the glutathione pathway induces cancer cell death,CMSP may trigger the oxidative stress response of cancer cells and lead to cancer cell death.

    In summary,a profiling method to achieve full coverage of the metabolites of GAP was developed.By monitoring the changes in each node in the metabolic pathway,the pathologic mechanisms of diseases including cancer and neurodegenerative diseases can be studied.The developed method is expected to offer a promising new tool to help researchers elucidate the roles of GAP in physiological and pathological processes.

    CRediT author statement

    Li Liu:Conceptualization,Methodology,Investigation,Writing-Original draft preparation;Yu-Han Lu:Conceptualization,Methodology,Investigation,Writing -Original draft preparation;Min-Dan Wang:Methodology,Formal analysis,Investigation;Qun-Fei Zhao:Methodology,Investigation;Xiu-Ping Chen:Conceptualization,Methodology,Investigation,Resources,Writing-Reviewing and Editing;Hang Yin:Conceptualization,Methodology,Investigation,Resources,Writing -Reviewing and Editing;Chen-Guo Feng:Supervision,Conceptualization,Funding acquisition,Writing-Reviewing and Editing;Fang Zhang:Supervision,Conceptualization,Funding acquisition,Writing -Reviewing and Editing.

    Declaration of competing interest

    The authors declare that there are no conflicts of interest.

    Acknowledgments

    We thank the Shanghai Municipal Committee of Science and Technology (Grant Nos.: 20XD1423400,23ZR1460900 and 20DZ2201100),Shanghai Municipal Health Commission/Shanghai Municipal Administration of Traditional Chinese Medicine (Grant No.: ZY(2021-2023)-0501),Shanghai Science and Technology Development Fund from Central Leading Local Government(Grant No.: YDZX20223100001004),National Natural Science Foundation of China (Grant No.: 21672249),and Expenditure Budget Program of Shanghai University of Traditional Chinese Medicine(Grant Nos.:2020LK051,and 2021LK001).

    Appendix A.Supplementary data

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

    国产 一区 欧美 日韩| 插阴视频在线观看视频| 亚洲国产欧美人成| 综合色av麻豆| 搡老乐熟女国产| 综合色丁香网| 欧美激情在线99| 少妇 在线观看| 青春草视频在线免费观看| 三级男女做爰猛烈吃奶摸视频| av国产免费在线观看| 搡老乐熟女国产| 最近2019中文字幕mv第一页| 看黄色毛片网站| 国产精品一区二区性色av| 亚洲国产精品专区欧美| 婷婷色综合www| 别揉我奶头 嗯啊视频| 不卡视频在线观看欧美| 国产午夜精品久久久久久一区二区三区| 一级毛片电影观看| 99热全是精品| 国产精品99久久久久久久久| 黄片wwwwww| 精品少妇久久久久久888优播| 午夜激情福利司机影院| 熟女av电影| 亚洲精品亚洲一区二区| 偷拍熟女少妇极品色| 免费黄色在线免费观看| 亚洲无线观看免费| 汤姆久久久久久久影院中文字幕| 偷拍熟女少妇极品色| 国产精品国产三级国产专区5o| 又黄又爽又刺激的免费视频.| 精品酒店卫生间| 国产 一区精品| 综合色av麻豆| av在线播放精品| 99热6这里只有精品| 久热这里只有精品99| 我的老师免费观看完整版| 黑人高潮一二区| 69av精品久久久久久| 特大巨黑吊av在线直播| 国产高清不卡午夜福利| av又黄又爽大尺度在线免费看| 亚洲精品国产av成人精品| 极品少妇高潮喷水抽搐| 三级国产精品欧美在线观看| 色哟哟·www| 亚洲精品日本国产第一区| 男女那种视频在线观看| 22中文网久久字幕| 18禁在线播放成人免费| 婷婷色av中文字幕| 看免费成人av毛片| 亚洲综合色惰| 国产老妇伦熟女老妇高清| 亚洲av成人精品一区久久| 你懂的网址亚洲精品在线观看| 久久久久久久久大av| 嫩草影院入口| 国产成人精品一,二区| av天堂中文字幕网| 高清av免费在线| 国产精品久久久久久精品电影小说 | 国产亚洲最大av| 丝袜脚勾引网站| 国产欧美日韩精品一区二区| 久久99蜜桃精品久久| 18禁在线播放成人免费| 免费看av在线观看网站| 黄色日韩在线| 日韩av免费高清视频| 丝袜脚勾引网站| 国产在线男女| 久久国产乱子免费精品| 王馨瑶露胸无遮挡在线观看| 一级爰片在线观看| 亚洲精品自拍成人| 女人十人毛片免费观看3o分钟| 欧美极品一区二区三区四区| 在线天堂最新版资源| 免费电影在线观看免费观看| 日韩亚洲欧美综合| 欧美老熟妇乱子伦牲交| 国产精品秋霞免费鲁丝片| videossex国产| 最近最新中文字幕大全电影3| 国产成人精品久久久久久| 777米奇影视久久| 亚洲国产色片| 久久久久久久久久成人| 少妇丰满av| av在线播放精品| 建设人人有责人人尽责人人享有的 | 久久精品人妻少妇| 亚洲成人精品中文字幕电影| 久久精品国产自在天天线| 亚洲综合色惰| 一级片'在线观看视频| 综合色丁香网| 亚洲四区av| 国产中年淑女户外野战色| 内地一区二区视频在线| 亚洲三级黄色毛片| 网址你懂的国产日韩在线| 在线看a的网站| 少妇丰满av| 久久久久久国产a免费观看| 欧美激情在线99| 国产精品一二三区在线看| 成年女人在线观看亚洲视频 | 亚洲av一区综合| 日韩亚洲欧美综合| 亚洲国产精品国产精品| 不卡视频在线观看欧美| 国产免费福利视频在线观看| 91午夜精品亚洲一区二区三区| 精品久久久久久久久亚洲| 最近最新中文字幕免费大全7| 狂野欧美激情性xxxx在线观看| 色播亚洲综合网| 国产成人aa在线观看| 成人无遮挡网站| 日韩亚洲欧美综合| 成人鲁丝片一二三区免费| 欧美极品一区二区三区四区| 国产精品国产三级国产专区5o| 色视频www国产| 成年免费大片在线观看| 成年版毛片免费区| 高清av免费在线| 日韩一区二区三区影片| 看免费成人av毛片| 少妇的逼好多水| 免费在线观看成人毛片| 日韩国内少妇激情av| 日本wwww免费看| 男人和女人高潮做爰伦理| 国产伦精品一区二区三区四那| 亚洲真实伦在线观看| av国产精品久久久久影院| 久久精品国产亚洲av天美| 免费观看在线日韩| 听说在线观看完整版免费高清| 18禁动态无遮挡网站| 欧美区成人在线视频| 欧美精品国产亚洲| 欧美 日韩 精品 国产| 国产男女超爽视频在线观看| 激情五月婷婷亚洲| 久久久久久久国产电影| 国产免费一级a男人的天堂| 亚洲精品国产成人久久av| 免费看av在线观看网站| 三级国产精品欧美在线观看| 久久人人爽av亚洲精品天堂 | 大香蕉久久网| 麻豆国产97在线/欧美| 女人被狂操c到高潮| 网址你懂的国产日韩在线| 亚洲国产欧美在线一区| 我的老师免费观看完整版| 五月玫瑰六月丁香| 成人鲁丝片一二三区免费| 在线a可以看的网站| 国产国拍精品亚洲av在线观看| 久久精品国产亚洲网站| 日韩成人伦理影院| 97超视频在线观看视频| 免费在线观看成人毛片| av.在线天堂| 欧美国产精品一级二级三级 | 午夜激情福利司机影院| 国产熟女欧美一区二区| 欧美国产精品一级二级三级 | 国产乱来视频区| 国产白丝娇喘喷水9色精品| 黑人高潮一二区| 青春草国产在线视频| 久久女婷五月综合色啪小说 | 乱码一卡2卡4卡精品| 在线观看av片永久免费下载| av卡一久久| 成人一区二区视频在线观看| 免费播放大片免费观看视频在线观看| 亚洲精品第二区| 伦理电影大哥的女人| 久久久国产一区二区| 久久99热6这里只有精品| 国产精品.久久久| 久久人人爽人人爽人人片va| 亚洲久久久久久中文字幕| 亚洲欧美精品自产自拍| 欧美zozozo另类| 国产成人精品一,二区| 国产精品无大码| 一本久久精品| 免费看光身美女| 免费不卡的大黄色大毛片视频在线观看| 狂野欧美激情性bbbbbb| 精品国产一区二区三区久久久樱花 | 国产一区二区三区av在线| 18禁在线无遮挡免费观看视频| 国产伦精品一区二区三区四那| 久久精品国产自在天天线| 国产老妇伦熟女老妇高清| 十八禁网站网址无遮挡 | 少妇人妻 视频| 下体分泌物呈黄色| 成人免费观看视频高清| 亚洲av电影在线观看一区二区三区 | 91午夜精品亚洲一区二区三区| 国产免费视频播放在线视频| 日韩成人伦理影院| 777米奇影视久久| 欧美 日韩 精品 国产| 1000部很黄的大片| 午夜老司机福利剧场| 亚洲精品一二三| 亚洲伊人久久精品综合| 成人二区视频| 亚洲精品视频女| 丝袜美腿在线中文| 人妻一区二区av| 一级毛片 在线播放| 日本黄色片子视频| 两个人的视频大全免费| 成人无遮挡网站| 亚洲国产精品999| 高清日韩中文字幕在线| 亚洲人成网站在线播| 交换朋友夫妻互换小说| 在线观看免费高清a一片| 可以在线观看毛片的网站| 亚洲精品成人久久久久久| 在线免费十八禁| 伦理电影大哥的女人| 久久久精品94久久精品| 97在线人人人人妻| av免费观看日本| 国产精品嫩草影院av在线观看| 国产成人freesex在线| 国产男人的电影天堂91| 新久久久久国产一级毛片| 一区二区三区精品91| 最近手机中文字幕大全| 网址你懂的国产日韩在线| 亚洲最大成人av| 免费黄频网站在线观看国产| 久久精品国产自在天天线| 最近2019中文字幕mv第一页| 国产精品一区二区性色av| 亚洲国产精品成人综合色| 午夜视频国产福利| 2018国产大陆天天弄谢| 国产精品麻豆人妻色哟哟久久| 日本午夜av视频| 简卡轻食公司| 亚洲欧美一区二区三区黑人 | 我要看日韩黄色一级片| 美女被艹到高潮喷水动态| 精品久久久久久电影网| 欧美 日韩 精品 国产| 日韩欧美 国产精品| 欧美日韩视频高清一区二区三区二| 国产精品久久久久久精品古装| 国产av码专区亚洲av| 69人妻影院| 久热这里只有精品99| 中国三级夫妇交换| av在线app专区| 精华霜和精华液先用哪个| 久久久久精品性色| 麻豆久久精品国产亚洲av| 91精品伊人久久大香线蕉| 一级毛片aaaaaa免费看小| 99热网站在线观看| 欧美极品一区二区三区四区| 黄片无遮挡物在线观看| a级毛色黄片| 身体一侧抽搐| 嫩草影院入口| 国产黄片美女视频| 午夜福利视频精品| 色综合色国产| 天美传媒精品一区二区| 国产成人一区二区在线| 欧美国产精品一级二级三级 | 亚洲精品乱码久久久v下载方式| 久久99蜜桃精品久久| 乱系列少妇在线播放| 国产精品一区二区三区四区免费观看| 黄色怎么调成土黄色| 波多野结衣巨乳人妻| 一本色道久久久久久精品综合| 国产爽快片一区二区三区| 国产精品国产三级专区第一集| 亚洲av免费高清在线观看| 三级经典国产精品| 久久6这里有精品| 日日摸夜夜添夜夜添av毛片| 精品一区在线观看国产| 午夜免费观看性视频| 亚洲av免费在线观看| 欧美高清性xxxxhd video| 亚洲精品视频女| 国产淫片久久久久久久久| 国产熟女欧美一区二区| 亚洲va在线va天堂va国产| 成年版毛片免费区| av卡一久久| 亚洲欧美日韩东京热| 18禁在线播放成人免费| av在线亚洲专区| 亚洲国产欧美人成| 亚洲av中文字字幕乱码综合| 国产精品久久久久久精品电影| 中文字幕亚洲精品专区| 成人黄色视频免费在线看| 在线免费观看不下载黄p国产| 全区人妻精品视频| 性色avwww在线观看| h日本视频在线播放| 视频中文字幕在线观看| 久久久精品免费免费高清| 一级片'在线观看视频| 日本一本二区三区精品| 97超视频在线观看视频| 国产日韩欧美在线精品| 亚洲av福利一区| 老女人水多毛片| 国产精品久久久久久精品电影| 内地一区二区视频在线| 99久久精品热视频| 中文字幕人妻熟人妻熟丝袜美| 亚洲成人精品中文字幕电影| 免费看日本二区| 国产精品成人在线| 国产精品久久久久久精品古装| 免费看a级黄色片| 国产中年淑女户外野战色| 欧美精品人与动牲交sv欧美| 午夜爱爱视频在线播放| 亚洲激情五月婷婷啪啪| 少妇 在线观看| 久久热精品热| 少妇人妻精品综合一区二区| 免费高清在线观看视频在线观看| 国产69精品久久久久777片| 亚洲精品久久久久久婷婷小说| 国产精品一区www在线观看| 成人毛片60女人毛片免费| 大香蕉97超碰在线| 1000部很黄的大片| 欧美xxxx性猛交bbbb| 在线观看av片永久免费下载| 一级黄片播放器| 中文字幕av成人在线电影| 精品一区二区三区视频在线| 交换朋友夫妻互换小说| 久久久久国产精品人妻一区二区| 亚洲av国产av综合av卡| 韩国高清视频一区二区三区| 亚洲自拍偷在线| 啦啦啦啦在线视频资源| av免费观看日本| 欧美 日韩 精品 国产| av国产免费在线观看| 特大巨黑吊av在线直播| 日本wwww免费看| 亚洲精品一区蜜桃| av免费在线看不卡| 亚洲在线观看片| 少妇人妻久久综合中文| 欧美丝袜亚洲另类| 99久久精品一区二区三区| 中文字幕制服av| 免费人成在线观看视频色| 2022亚洲国产成人精品| 超碰av人人做人人爽久久| 别揉我奶头 嗯啊视频| 国产伦精品一区二区三区视频9| 男女下面进入的视频免费午夜| 亚洲精品日本国产第一区| 内射极品少妇av片p| 精品久久久久久久久av| 别揉我奶头 嗯啊视频| 2021天堂中文幕一二区在线观| 毛片女人毛片| 人妻制服诱惑在线中文字幕| 欧美国产精品一级二级三级 | 国产成人免费无遮挡视频| 少妇高潮的动态图| 日韩欧美精品免费久久| 日本黄大片高清| 人妻少妇偷人精品九色| 亚洲av国产av综合av卡| 韩国av在线不卡| 亚洲真实伦在线观看| 亚洲成人一二三区av| 少妇人妻一区二区三区视频| 亚洲美女搞黄在线观看| 亚洲av.av天堂| 免费黄频网站在线观看国产| 国产男女内射视频| 国产精品人妻久久久久久| 国产精品精品国产色婷婷| 欧美老熟妇乱子伦牲交| 欧美激情在线99| 九草在线视频观看| 18禁裸乳无遮挡免费网站照片| 国产有黄有色有爽视频| 欧美zozozo另类| 人人妻人人看人人澡| 国产成人freesex在线| 看黄色毛片网站| 久久久久久伊人网av| 亚洲高清免费不卡视频| 91精品一卡2卡3卡4卡| 亚洲电影在线观看av| 美女内射精品一级片tv| 国产一区有黄有色的免费视频| 哪个播放器可以免费观看大片| 在线免费十八禁| 高清视频免费观看一区二区| 久久人人爽人人爽人人片va| 99热这里只有是精品50| 久久久久久久久大av| 久久久久久久大尺度免费视频| 欧美亚洲 丝袜 人妻 在线| 亚洲欧美中文字幕日韩二区| 亚洲经典国产精华液单| 国产精品精品国产色婷婷| 国产美女午夜福利| 免费大片黄手机在线观看| 激情五月婷婷亚洲| av黄色大香蕉| 五月伊人婷婷丁香| 亚洲精品影视一区二区三区av| 成人美女网站在线观看视频| 黄片wwwwww| 亚洲国产最新在线播放| 97超碰精品成人国产| 大码成人一级视频| 丝瓜视频免费看黄片| 深夜a级毛片| 少妇熟女欧美另类| 美女脱内裤让男人舔精品视频| 亚洲,一卡二卡三卡| 少妇的逼水好多| 国产av不卡久久| 最近最新中文字幕大全电影3| 国产av不卡久久| 欧美老熟妇乱子伦牲交| 一个人观看的视频www高清免费观看| 在线天堂最新版资源| 国产人妻一区二区三区在| 亚洲真实伦在线观看| 亚洲精品国产av蜜桃| 日韩大片免费观看网站| 小蜜桃在线观看免费完整版高清| 日本爱情动作片www.在线观看| 国产男女超爽视频在线观看| 国产老妇伦熟女老妇高清| 男人舔奶头视频| 免费看av在线观看网站| 男女下面进入的视频免费午夜| av在线老鸭窝| 嫩草影院精品99| 国产精品女同一区二区软件| 国产毛片在线视频| 性插视频无遮挡在线免费观看| 亚洲国产高清在线一区二区三| 国产成人aa在线观看| 亚洲天堂av无毛| 国产久久久一区二区三区| 看十八女毛片水多多多| 亚洲精品日本国产第一区| 免费看a级黄色片| 亚洲精品乱码久久久v下载方式| 国产午夜精品久久久久久一区二区三区| 麻豆乱淫一区二区| 国产精品偷伦视频观看了| 3wmmmm亚洲av在线观看| 国产淫片久久久久久久久| 日本-黄色视频高清免费观看| 亚洲精品影视一区二区三区av| 人妻制服诱惑在线中文字幕| 青青草视频在线视频观看| 国产午夜精品一二区理论片| 夫妻午夜视频| 日韩电影二区| 有码 亚洲区| 久久精品熟女亚洲av麻豆精品| 七月丁香在线播放| 97在线人人人人妻| 国产精品久久久久久精品电影| 极品少妇高潮喷水抽搐| 亚洲国产精品专区欧美| 亚洲精品456在线播放app| 欧美精品一区二区大全| av女优亚洲男人天堂| 成人亚洲精品一区在线观看 | 身体一侧抽搐| 久久久久九九精品影院| 熟女电影av网| 欧美成人一区二区免费高清观看| 99精国产麻豆久久婷婷| 国语对白做爰xxxⅹ性视频网站| 白带黄色成豆腐渣| 国产精品三级大全| 亚洲精品一二三| 精品少妇久久久久久888优播| 中文资源天堂在线| 亚洲内射少妇av| 精品久久久噜噜| 精品久久国产蜜桃| 91狼人影院| av卡一久久| 寂寞人妻少妇视频99o| 狂野欧美激情性bbbbbb| 在线观看av片永久免费下载| 国产在视频线精品| 亚洲电影在线观看av| 国产精品99久久99久久久不卡 | 大片电影免费在线观看免费| 亚洲在线观看片| 亚洲内射少妇av| 欧美一级a爱片免费观看看| 九色成人免费人妻av| 亚洲在线观看片| 青春草亚洲视频在线观看| 国产精品99久久久久久久久| 最近中文字幕高清免费大全6| av天堂中文字幕网| 亚洲四区av| 久久热精品热| 亚洲欧美日韩另类电影网站 | 欧美精品人与动牲交sv欧美| 精品人妻熟女av久视频| 91午夜精品亚洲一区二区三区| 亚洲人与动物交配视频| 黄片无遮挡物在线观看| 免费黄网站久久成人精品| 一区二区av电影网| 精品人妻一区二区三区麻豆| 亚洲av免费在线观看| 午夜免费男女啪啪视频观看| 久久综合国产亚洲精品| 人人妻人人澡人人爽人人夜夜| 一区二区三区免费毛片| 欧美日本视频| 在线观看免费高清a一片| 在线a可以看的网站| 五月天丁香电影| 国产精品偷伦视频观看了| 狠狠精品人妻久久久久久综合| 免费电影在线观看免费观看| 午夜视频国产福利| 天美传媒精品一区二区| 国产成人freesex在线| 69av精品久久久久久| 99久久九九国产精品国产免费| 国产免费一区二区三区四区乱码| 三级男女做爰猛烈吃奶摸视频| 麻豆成人午夜福利视频| 国产午夜精品久久久久久一区二区三区| 亚洲精品第二区| 男女下面进入的视频免费午夜| 国产成人精品一,二区| 搡老乐熟女国产| 内射极品少妇av片p| 男男h啪啪无遮挡| 三级男女做爰猛烈吃奶摸视频| 老司机影院成人| 国产精品无大码| 老女人水多毛片| 国产高清国产精品国产三级 | 成人亚洲精品一区在线观看 | 一本久久精品| 日韩av不卡免费在线播放| 美女视频免费永久观看网站| 一区二区av电影网| 欧美极品一区二区三区四区| 各种免费的搞黄视频| 看黄色毛片网站| 嫩草影院入口| 国产永久视频网站| 啦啦啦啦在线视频资源| 亚洲国产精品专区欧美| 80岁老熟妇乱子伦牲交| 亚洲色图综合在线观看| 国产免费福利视频在线观看| 干丝袜人妻中文字幕| 精品久久久久久电影网| 亚洲在线观看片| 噜噜噜噜噜久久久久久91| 免费大片黄手机在线观看| 久久国内精品自在自线图片| 国产亚洲5aaaaa淫片| 精品人妻一区二区三区麻豆| 精品一区在线观看国产| 日韩,欧美,国产一区二区三区| 久久国内精品自在自线图片| 国产永久视频网站| 中文字幕免费在线视频6| 亚洲内射少妇av| 国产成人精品福利久久| 毛片女人毛片| 22中文网久久字幕| 午夜精品国产一区二区电影 | 日韩 亚洲 欧美在线|