• <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.

    夜夜夜夜夜久久久久| 精品国产一区二区三区四区第35| 又黄又爽又免费观看的视频| 又黄又爽又免费观看的视频| 丰满饥渴人妻一区二区三| 黄色 视频免费看| 热re99久久精品国产66热6| 嫩草影视91久久| 久久人人精品亚洲av| 国产精品久久久av美女十八| 一级a爱片免费观看的视频| 一区二区日韩欧美中文字幕| 欧美在线黄色| 日本 av在线| 国产91精品成人一区二区三区| 天天躁夜夜躁狠狠躁躁| 国产乱人伦免费视频| 国产高清国产精品国产三级| 精品欧美一区二区三区在线| 一级片'在线观看视频| 久久精品国产99精品国产亚洲性色 | 免费在线观看黄色视频的| 久久久久久亚洲精品国产蜜桃av| 亚洲avbb在线观看| 婷婷丁香在线五月| 国产精品一区二区免费欧美| 无限看片的www在线观看| 欧美成人午夜精品| 黄色 视频免费看| 日本欧美视频一区| 精品国产美女av久久久久小说| 国产亚洲精品第一综合不卡| 两人在一起打扑克的视频| 精品久久久精品久久久| 国产99白浆流出| 日本欧美视频一区| 国产一区二区在线av高清观看| 一二三四社区在线视频社区8| 少妇裸体淫交视频免费看高清 | 亚洲中文av在线| 久久久久久久午夜电影 | x7x7x7水蜜桃| 亚洲avbb在线观看| 成人国产一区最新在线观看| 不卡一级毛片| 不卡av一区二区三区| 免费观看精品视频网站| av电影中文网址| 又黄又爽又免费观看的视频| 丝袜在线中文字幕| 夜夜看夜夜爽夜夜摸 | 国产xxxxx性猛交| 女人高潮潮喷娇喘18禁视频| 精品一品国产午夜福利视频| 国产精品二区激情视频| 叶爱在线成人免费视频播放| 国产av又大| 久久草成人影院| 精品一区二区三区av网在线观看| 日韩欧美在线二视频| 怎么达到女性高潮| 亚洲 欧美 日韩 在线 免费| 国产无遮挡羞羞视频在线观看| 国产亚洲精品久久久久5区| 成人三级黄色视频| 男人舔女人的私密视频| 搡老熟女国产l中国老女人| 亚洲av熟女| videosex国产| av在线播放免费不卡| 精品乱码久久久久久99久播| 欧美日韩亚洲综合一区二区三区_| 别揉我奶头~嗯~啊~动态视频| 国产精品av久久久久免费| 久久精品影院6| 亚洲精品粉嫩美女一区| 久久人妻福利社区极品人妻图片| cao死你这个sao货| 波多野结衣一区麻豆| 女人精品久久久久毛片| xxx96com| 国产亚洲精品综合一区在线观看 | 少妇粗大呻吟视频| 一a级毛片在线观看| 日韩三级视频一区二区三区| 最近最新免费中文字幕在线| 亚洲欧美激情综合另类| 精品国产乱子伦一区二区三区| 亚洲av五月六月丁香网| 九色亚洲精品在线播放| 免费一级毛片在线播放高清视频 | 黑人猛操日本美女一级片| 国产精品 欧美亚洲| 国产单亲对白刺激| 国产精品 欧美亚洲| 9191精品国产免费久久| 在线观看www视频免费| 午夜福利一区二区在线看| 日韩成人在线观看一区二区三区| 亚洲精品一二三| 免费女性裸体啪啪无遮挡网站| 欧美人与性动交α欧美精品济南到| 欧洲精品卡2卡3卡4卡5卡区| 亚洲成国产人片在线观看| 色尼玛亚洲综合影院| 国产1区2区3区精品| 精品福利观看| 97碰自拍视频| 美女午夜性视频免费| 欧美日本中文国产一区发布| 亚洲av熟女| 久久这里只有精品19| 国产高清国产精品国产三级| 不卡av一区二区三区| 欧美一区二区精品小视频在线| 欧美在线黄色| 女性生殖器流出的白浆| 桃红色精品国产亚洲av| 婷婷丁香在线五月| 看黄色毛片网站| 天堂俺去俺来也www色官网| 国产精品国产av在线观看| 一级,二级,三级黄色视频| 国产成人啪精品午夜网站| 精品少妇一区二区三区视频日本电影| 丰满的人妻完整版| 操出白浆在线播放| 岛国在线观看网站| 丰满的人妻完整版| 性少妇av在线| 一边摸一边做爽爽视频免费| 露出奶头的视频| 高清毛片免费观看视频网站 | 精品午夜福利视频在线观看一区| 成熟少妇高潮喷水视频| 99久久综合精品五月天人人| 制服人妻中文乱码| 免费在线观看视频国产中文字幕亚洲| 亚洲 欧美 日韩 在线 免费| 精品人妻在线不人妻| 精品久久久久久久毛片微露脸| 国产色视频综合| 亚洲男人的天堂狠狠| 亚洲av成人一区二区三| 国产成人欧美| 亚洲视频免费观看视频| cao死你这个sao货| www国产在线视频色| 午夜日韩欧美国产| 久久婷婷成人综合色麻豆| 亚洲欧美一区二区三区黑人| 最近最新中文字幕大全免费视频| 国产亚洲欧美精品永久| 国产精品98久久久久久宅男小说| 午夜免费观看网址| 99久久精品国产亚洲精品| 超碰成人久久| 19禁男女啪啪无遮挡网站| 欧美丝袜亚洲另类 | 国产又爽黄色视频| 成人手机av| 中文字幕高清在线视频| 亚洲国产看品久久| 99久久国产精品久久久| 欧美日韩av久久| 亚洲视频免费观看视频| 精品无人区乱码1区二区| 欧美乱色亚洲激情| 男女床上黄色一级片免费看| av网站在线播放免费| 亚洲三区欧美一区| 亚洲七黄色美女视频| 久久影院123| 大型黄色视频在线免费观看| 日韩人妻精品一区2区三区| 欧美最黄视频在线播放免费 | 操出白浆在线播放| 久久精品亚洲熟妇少妇任你| 成年人免费黄色播放视频| 久久草成人影院| 久久精品91蜜桃| 国产成人精品在线电影| 午夜久久久在线观看| 欧美日韩黄片免| 黄色视频,在线免费观看| 在线播放国产精品三级| 国产野战对白在线观看| 国产av又大| 不卡av一区二区三区| 99精品欧美一区二区三区四区| 国产成人一区二区三区免费视频网站| 久久精品91蜜桃| 欧美日韩亚洲综合一区二区三区_| 国产人伦9x9x在线观看| av中文乱码字幕在线| 一个人观看的视频www高清免费观看 | av网站在线播放免费| 久久婷婷成人综合色麻豆| 精品免费久久久久久久清纯| 我的亚洲天堂| 亚洲成人免费av在线播放| 久久精品91无色码中文字幕| 欧美人与性动交α欧美软件| 欧美午夜高清在线| 桃色一区二区三区在线观看| 老熟妇仑乱视频hdxx| 亚洲av日韩精品久久久久久密| 免费久久久久久久精品成人欧美视频| 精品人妻1区二区| 电影成人av| 啦啦啦在线免费观看视频4| 国产成人系列免费观看| 成人手机av| 欧美性长视频在线观看| 露出奶头的视频| 级片在线观看| 亚洲国产毛片av蜜桃av| 别揉我奶头~嗯~啊~动态视频| 又黄又爽又免费观看的视频| 久久中文字幕人妻熟女| 久久国产精品人妻蜜桃| 亚洲自拍偷在线| 午夜免费激情av| 男男h啪啪无遮挡| 国产高清激情床上av| 欧美日韩一级在线毛片| 精品国产亚洲在线| 夜夜爽天天搞| 亚洲午夜理论影院| 男女床上黄色一级片免费看| 国产av一区在线观看免费| 亚洲伊人色综图| 久久久久九九精品影院| 美女福利国产在线| 亚洲久久久国产精品| 国产成人欧美在线观看| 国产无遮挡羞羞视频在线观看| 极品教师在线免费播放| 日韩欧美在线二视频| 男人舔女人的私密视频| 男人操女人黄网站| av天堂在线播放| 交换朋友夫妻互换小说| 香蕉国产在线看| 久久久国产一区二区| 麻豆成人av在线观看| 一二三四社区在线视频社区8| 高清av免费在线| 12—13女人毛片做爰片一| 天堂动漫精品| 成人18禁高潮啪啪吃奶动态图| 国产无遮挡羞羞视频在线观看| 久久久久久久精品吃奶| 大码成人一级视频| 久久精品91蜜桃| 亚洲av电影在线进入| 黄片播放在线免费| 亚洲一区二区三区不卡视频| 日韩欧美一区二区三区在线观看| aaaaa片日本免费| 麻豆一二三区av精品| av中文乱码字幕在线| 啪啪无遮挡十八禁网站| 午夜福利影视在线免费观看| 女人高潮潮喷娇喘18禁视频| 男人操女人黄网站| 亚洲国产精品合色在线| 免费观看人在逋| 伦理电影免费视频| 身体一侧抽搐| 中亚洲国语对白在线视频| 正在播放国产对白刺激| 少妇被粗大的猛进出69影院| 免费在线观看亚洲国产| 热re99久久精品国产66热6| 黑丝袜美女国产一区| 久久久久久久久免费视频了| 国产一区在线观看成人免费| 久久久精品国产亚洲av高清涩受| 一二三四在线观看免费中文在| 免费观看人在逋| 久久午夜综合久久蜜桃| 欧美日韩亚洲高清精品| 亚洲国产欧美一区二区综合| 成人精品一区二区免费| 国产欧美日韩一区二区三| 久久精品成人免费网站| 亚洲色图av天堂| 丁香欧美五月| 久久精品国产99精品国产亚洲性色 | 老司机午夜十八禁免费视频| 国产视频一区二区在线看| 久久久精品欧美日韩精品| 操出白浆在线播放| 久久久久久人人人人人| 夜夜爽天天搞| 国产激情久久老熟女| 一区二区三区激情视频| 一边摸一边抽搐一进一小说| 久久精品亚洲熟妇少妇任你| 亚洲中文av在线| 国产在线精品亚洲第一网站| 国产三级在线视频| 99香蕉大伊视频| 久久精品亚洲熟妇少妇任你| 久久人妻熟女aⅴ| 怎么达到女性高潮| 99久久人妻综合| 精品免费久久久久久久清纯| 亚洲第一欧美日韩一区二区三区| 亚洲欧美日韩另类电影网站| 美女午夜性视频免费| 免费在线观看完整版高清| 18禁黄网站禁片午夜丰满| 国产在线精品亚洲第一网站| 免费在线观看视频国产中文字幕亚洲| 国产野战对白在线观看| 亚洲国产欧美一区二区综合| 国产成人啪精品午夜网站| 欧美精品啪啪一区二区三区| 国产蜜桃级精品一区二区三区| 男女下面插进去视频免费观看| 成人亚洲精品av一区二区 | 国产一区二区在线av高清观看| 欧美乱妇无乱码| 咕卡用的链子| 欧美av亚洲av综合av国产av| 69av精品久久久久久| 国产麻豆69| 欧美黄色淫秽网站| 黑人猛操日本美女一级片| 国产99久久九九免费精品| 精品人妻在线不人妻| 极品人妻少妇av视频| 伊人久久大香线蕉亚洲五| av网站免费在线观看视频| www日本在线高清视频| av欧美777| 久久精品国产99精品国产亚洲性色 | www.熟女人妻精品国产| 在线观看一区二区三区激情| 久久久国产成人免费| 亚洲中文字幕日韩| 免费搜索国产男女视频| 97人妻天天添夜夜摸| 午夜福利免费观看在线| 亚洲午夜理论影院| 国产又色又爽无遮挡免费看| 亚洲色图 男人天堂 中文字幕| 黄色毛片三级朝国网站| 日韩免费av在线播放| 久久中文字幕人妻熟女| 免费在线观看亚洲国产| 老司机深夜福利视频在线观看| 欧美午夜高清在线| 欧美日韩一级在线毛片| 国产在线精品亚洲第一网站| 国产欧美日韩精品亚洲av| 精品午夜福利视频在线观看一区| av欧美777| 操美女的视频在线观看| 涩涩av久久男人的天堂| 欧美精品一区二区免费开放| 午夜视频精品福利| 午夜免费激情av| 夜夜看夜夜爽夜夜摸 | xxxhd国产人妻xxx| 亚洲视频免费观看视频| 91老司机精品| 少妇被粗大的猛进出69影院| 9色porny在线观看| 日本 av在线| √禁漫天堂资源中文www| 亚洲全国av大片| 精品熟女少妇八av免费久了| 日本黄色视频三级网站网址| 一级作爱视频免费观看| 亚洲精品中文字幕在线视频| 18禁黄网站禁片午夜丰满| 欧美精品亚洲一区二区| 亚洲国产精品一区二区三区在线| 欧美国产精品va在线观看不卡| 国产精品日韩av在线免费观看 | 国产亚洲精品久久久久久毛片| 亚洲九九香蕉| 久久精品亚洲熟妇少妇任你| 香蕉丝袜av| 欧美性长视频在线观看| 19禁男女啪啪无遮挡网站| 看片在线看免费视频| av欧美777| 色婷婷av一区二区三区视频| 黄色怎么调成土黄色| 极品人妻少妇av视频| 欧美另类亚洲清纯唯美| 亚洲精品国产区一区二| 91麻豆精品激情在线观看国产 | 久久青草综合色| 热99国产精品久久久久久7| 80岁老熟妇乱子伦牲交| 国产麻豆69| 在线观看舔阴道视频| 在线观看一区二区三区激情| 一级作爱视频免费观看| 国产伦一二天堂av在线观看| 国产深夜福利视频在线观看| 麻豆av在线久日| 亚洲色图综合在线观看| 国产不卡一卡二| 欧美精品一区二区免费开放| 国产精品野战在线观看 | 亚洲一区二区三区欧美精品| 怎么达到女性高潮| 亚洲第一av免费看| 久久精品成人免费网站| 亚洲自偷自拍图片 自拍| 婷婷六月久久综合丁香| 性欧美人与动物交配| 啦啦啦免费观看视频1| 很黄的视频免费| 日韩欧美国产一区二区入口| 国产精品免费一区二区三区在线| 黑人巨大精品欧美一区二区蜜桃| 色综合站精品国产| 少妇被粗大的猛进出69影院| 久久久久国内视频| 在线观看免费高清a一片| av视频免费观看在线观看| 极品人妻少妇av视频| 国产成+人综合+亚洲专区| 国产高清激情床上av| 久久久久九九精品影院| 香蕉久久夜色| 亚洲精品美女久久久久99蜜臀| 高清毛片免费观看视频网站 | 在线观看日韩欧美| 欧美激情久久久久久爽电影 | 一级a爱视频在线免费观看| 亚洲成人国产一区在线观看| 精品欧美一区二区三区在线| 亚洲va日本ⅴa欧美va伊人久久| 国产精品永久免费网站| 两个人免费观看高清视频| 国产午夜精品久久久久久| 亚洲成a人片在线一区二区| 欧美一区二区精品小视频在线| 国产免费现黄频在线看| 免费在线观看亚洲国产| 亚洲精品一卡2卡三卡4卡5卡| 亚洲中文av在线| 国产亚洲精品久久久久5区| 国产主播在线观看一区二区| 女警被强在线播放| 麻豆成人av在线观看| 99在线人妻在线中文字幕| 黑人猛操日本美女一级片| 午夜免费成人在线视频| 99香蕉大伊视频| 激情在线观看视频在线高清| 一进一出好大好爽视频| 熟女少妇亚洲综合色aaa.| 精品电影一区二区在线| 亚洲人成电影免费在线| 欧美 亚洲 国产 日韩一| 日本 av在线| 老司机福利观看| av天堂在线播放| 无人区码免费观看不卡| 黄色成人免费大全| 国产三级黄色录像| 久久午夜综合久久蜜桃| 国产xxxxx性猛交| 久久草成人影院| ponron亚洲| av在线天堂中文字幕 | 男人的好看免费观看在线视频 | 亚洲五月婷婷丁香| 国产乱人伦免费视频| 美女高潮到喷水免费观看| www.熟女人妻精品国产| 18禁观看日本| 黑丝袜美女国产一区| 中文字幕最新亚洲高清| 久久人人97超碰香蕉20202| 久久人妻熟女aⅴ| 可以免费在线观看a视频的电影网站| 一边摸一边做爽爽视频免费| 午夜免费观看网址| 欧美成狂野欧美在线观看| 高清av免费在线| 欧美一区二区精品小视频在线| 天堂√8在线中文| 热99re8久久精品国产| 热re99久久国产66热| 啦啦啦免费观看视频1| 久久久精品欧美日韩精品| 日本vs欧美在线观看视频| 亚洲一区二区三区不卡视频| 男女做爰动态图高潮gif福利片 | 欧美精品亚洲一区二区| 一夜夜www| 亚洲黑人精品在线| 亚洲精品国产色婷婷电影| 欧美+亚洲+日韩+国产| 日韩视频一区二区在线观看| 黄色片一级片一级黄色片| 人妻丰满熟妇av一区二区三区| 日韩国内少妇激情av| 亚洲少妇的诱惑av| 国产av精品麻豆| 人人妻人人澡人人看| 丝袜美足系列| 亚洲精品在线观看二区| 久久亚洲真实| 精品日产1卡2卡| 国产精品亚洲一级av第二区| 在线永久观看黄色视频| 美女 人体艺术 gogo| 黑人操中国人逼视频| 国产xxxxx性猛交| 亚洲五月天丁香| 亚洲中文日韩欧美视频| 亚洲欧美日韩无卡精品| 欧美日韩乱码在线| 淫妇啪啪啪对白视频| 一区二区三区激情视频| 99在线人妻在线中文字幕| 无人区码免费观看不卡| av免费在线观看网站| 久久亚洲精品不卡| 久久精品亚洲av国产电影网| xxxhd国产人妻xxx| 巨乳人妻的诱惑在线观看| 久久人妻福利社区极品人妻图片| 欧美一级毛片孕妇| 无限看片的www在线观看| 日韩三级视频一区二区三区| 99热国产这里只有精品6| 国产精品野战在线观看 | a在线观看视频网站| 又黄又粗又硬又大视频| 欧美激情 高清一区二区三区| 国产高清videossex| 99re在线观看精品视频| 国产亚洲欧美精品永久| 亚洲欧美日韩无卡精品| 两性夫妻黄色片| 午夜激情av网站| 国产一区二区三区综合在线观看| 亚洲在线自拍视频| 欧美黑人精品巨大| 免费不卡黄色视频| 91成人精品电影| 在线观看日韩欧美| 久久国产精品影院| 亚洲中文av在线| 老司机午夜十八禁免费视频| 日韩欧美国产一区二区入口| 99精品久久久久人妻精品| 亚洲精品在线观看二区| 日本vs欧美在线观看视频| 成人三级黄色视频| 热99re8久久精品国产| 黄色 视频免费看| 亚洲欧美激情综合另类| 国产成人精品在线电影| 人人妻人人添人人爽欧美一区卜| 免费在线观看完整版高清| 国产97色在线日韩免费| 黑人巨大精品欧美一区二区蜜桃| 欧美日韩视频精品一区| 成人国产一区最新在线观看| 麻豆av在线久日| 亚洲在线自拍视频| 日本精品一区二区三区蜜桃| 国产精品一区二区在线不卡| 国产精品免费一区二区三区在线| svipshipincom国产片| 成人三级黄色视频| 国产麻豆69| 欧美激情久久久久久爽电影 | 国产熟女xx| 香蕉丝袜av| 免费久久久久久久精品成人欧美视频| 黄色成人免费大全| 国产精华一区二区三区| 岛国在线观看网站| 中国美女看黄片| 亚洲av成人一区二区三| 亚洲激情在线av| 日韩大尺度精品在线看网址 | 好男人电影高清在线观看| 成年人黄色毛片网站| ponron亚洲| 日本wwww免费看| 欧美在线黄色| 日韩一卡2卡3卡4卡2021年| 无人区码免费观看不卡| 国产亚洲精品一区二区www| 日本黄色日本黄色录像| 日本撒尿小便嘘嘘汇集6| 99re在线观看精品视频| 亚洲免费av在线视频| 国产又爽黄色视频| 午夜老司机福利片| 久久人人97超碰香蕉20202| 69精品国产乱码久久久| 亚洲少妇的诱惑av| 国产欧美日韩一区二区精品| 一边摸一边做爽爽视频免费| 国产精品成人在线| 麻豆成人av在线观看| 在线观看免费日韩欧美大片| 欧美国产精品va在线观看不卡| 91成年电影在线观看|