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

    成年美女黄网站色视频大全免费 | 在线观看免费视频网站a站| 97精品久久久久久久久久精品| 高清不卡的av网站| 在线观看三级黄色| 亚洲av欧美aⅴ国产| 精品人妻熟女毛片av久久网站| 精品亚洲成国产av| 成人二区视频| 国产精品嫩草影院av在线观看| 国产日韩一区二区三区精品不卡 | 有码 亚洲区| 一级毛片黄色毛片免费观看视频| 国产成人午夜福利电影在线观看| 亚洲国产毛片av蜜桃av| 亚洲成人一二三区av| 九色亚洲精品在线播放| 国产片特级美女逼逼视频| 国产成人免费无遮挡视频| √禁漫天堂资源中文www| 日本午夜av视频| 久久免费观看电影| 99精国产麻豆久久婷婷| 免费大片黄手机在线观看| 成人午夜精彩视频在线观看| 午夜精品国产一区二区电影| 亚洲欧美一区二区三区国产| 国语对白做爰xxxⅹ性视频网站| 久久人人爽av亚洲精品天堂| 色哟哟·www| 久久久久视频综合| 欧美精品人与动牲交sv欧美| 看十八女毛片水多多多| 人人妻人人爽人人添夜夜欢视频| 丰满饥渴人妻一区二区三| a 毛片基地| 一本大道久久a久久精品| 91久久精品国产一区二区成人| 成人18禁高潮啪啪吃奶动态图 | 美女国产视频在线观看| 日本色播在线视频| 亚洲精品,欧美精品| 日韩一本色道免费dvd| 欧美三级亚洲精品| 欧美激情极品国产一区二区三区 | 亚洲色图综合在线观看| 国产有黄有色有爽视频| 伊人久久精品亚洲午夜| 久久久欧美国产精品| 夜夜爽夜夜爽视频| av在线老鸭窝| 亚洲国产精品国产精品| 建设人人有责人人尽责人人享有的| 国产精品久久久久久精品电影小说| 国产免费视频播放在线视频| 人人妻人人添人人爽欧美一区卜| 国产国拍精品亚洲av在线观看| 亚洲精品色激情综合| 中国三级夫妇交换| 亚洲伊人久久精品综合| 亚洲av成人精品一二三区| 欧美bdsm另类| 国产永久视频网站| 91久久精品电影网| 日韩精品有码人妻一区| 97超视频在线观看视频| av天堂久久9| 美女内射精品一级片tv| 午夜激情av网站| 国产精品蜜桃在线观看| 国产日韩欧美在线精品| 国产成人一区二区在线| 亚洲欧美一区二区三区黑人 | 日韩一区二区三区影片| 99久久人妻综合| 久久久精品免费免费高清| 国产精品一区二区三区四区免费观看| 亚洲精品一区蜜桃| 国产欧美日韩一区二区三区在线 | 午夜免费观看性视频| 亚洲丝袜综合中文字幕| 高清毛片免费看| 欧美成人午夜免费资源| 999精品在线视频| 女的被弄到高潮叫床怎么办| 999精品在线视频| 蜜桃国产av成人99| 飞空精品影院首页| 国产精品久久久久久精品电影小说| 中文乱码字字幕精品一区二区三区| 人人澡人人妻人| 狂野欧美激情性xxxx在线观看| 国产国拍精品亚洲av在线观看| 2018国产大陆天天弄谢| 国产亚洲欧美精品永久| 成人国产麻豆网| 国产国语露脸激情在线看| 国产熟女午夜一区二区三区 | 亚洲,一卡二卡三卡| 久久久久人妻精品一区果冻| 精品久久国产蜜桃| 精品久久蜜臀av无| 在线观看三级黄色| 97超视频在线观看视频| 久久久国产一区二区| 国产成人精品婷婷| 一级a做视频免费观看| 欧美日韩精品成人综合77777| 一区二区三区精品91| 只有这里有精品99| 亚洲精品自拍成人| 国产黄色免费在线视频| 这个男人来自地球电影免费观看 | 国产日韩欧美亚洲二区| 99热全是精品| 欧美3d第一页| 性色av一级| 国产乱来视频区| 一边摸一边做爽爽视频免费| 少妇熟女欧美另类| 一本久久精品| 一本久久精品| 国产精品久久久久久精品电影小说| 欧美精品亚洲一区二区| 国产精品99久久99久久久不卡 | 免费高清在线观看日韩| 日本色播在线视频| 18禁在线播放成人免费| 午夜激情av网站| 亚洲不卡免费看| 免费黄频网站在线观看国产| 精品亚洲乱码少妇综合久久| 大又大粗又爽又黄少妇毛片口| 秋霞伦理黄片| 精品一品国产午夜福利视频| 久久久精品免费免费高清| 大话2 男鬼变身卡| 中国三级夫妇交换| 国产在视频线精品| av黄色大香蕉| .国产精品久久| 国国产精品蜜臀av免费| 亚洲五月色婷婷综合| 亚洲三级黄色毛片| 永久网站在线| 大片电影免费在线观看免费| 99久久精品国产国产毛片| 欧美 亚洲 国产 日韩一| 欧美成人午夜免费资源| 亚洲国产日韩一区二区| 91国产中文字幕| 大片电影免费在线观看免费| 狂野欧美激情性bbbbbb| 成年女人在线观看亚洲视频| 边亲边吃奶的免费视频| 国产亚洲午夜精品一区二区久久| 三级国产精品片| 国产精品蜜桃在线观看| 女性生殖器流出的白浆| 中文欧美无线码| 国产片内射在线| 午夜老司机福利剧场| 王馨瑶露胸无遮挡在线观看| 女性被躁到高潮视频| av女优亚洲男人天堂| 秋霞在线观看毛片| 精品亚洲乱码少妇综合久久| 极品少妇高潮喷水抽搐| 国产日韩欧美亚洲二区| 秋霞伦理黄片| 九草在线视频观看| 国产熟女欧美一区二区| 成年美女黄网站色视频大全免费 | 国产探花极品一区二区| 亚洲经典国产精华液单| 国产成人精品无人区| 91国产中文字幕| 狂野欧美激情性xxxx在线观看| 亚洲欧美中文字幕日韩二区| 亚洲精品成人av观看孕妇| 国产乱来视频区| 免费看不卡的av| 高清黄色对白视频在线免费看| 美女国产视频在线观看| 嫩草影院入口| 三上悠亚av全集在线观看| 久久鲁丝午夜福利片| 91aial.com中文字幕在线观看| 国产男女内射视频| 美女cb高潮喷水在线观看| 欧美成人午夜免费资源| 岛国毛片在线播放| 美女主播在线视频| 青青草视频在线视频观看| 韩国av在线不卡| 亚洲人成网站在线观看播放| 国产精品 国内视频| 一级片'在线观看视频| 夜夜骑夜夜射夜夜干| 亚洲精品一区蜜桃| 日韩熟女老妇一区二区性免费视频| 精品国产乱码久久久久久小说| 久久综合国产亚洲精品| 最近2019中文字幕mv第一页| 久久这里有精品视频免费| 国产成人精品久久久久久| 亚州av有码| 免费观看的影片在线观看| 国产乱人偷精品视频| 国产精品.久久久| 3wmmmm亚洲av在线观看| 在线观看美女被高潮喷水网站| 99热全是精品| 国产精品人妻久久久久久| 美女大奶头黄色视频| 久久人人爽人人爽人人片va| 中文字幕精品免费在线观看视频 | 在线亚洲精品国产二区图片欧美 | 久久99热这里只频精品6学生| 国产视频内射| 国产精品不卡视频一区二区| 久久久久久久久久久久大奶| 国产精品免费大片| a 毛片基地| 成人午夜精彩视频在线观看| 91午夜精品亚洲一区二区三区| 多毛熟女@视频| 欧美日韩av久久| 久久精品熟女亚洲av麻豆精品| 纵有疾风起免费观看全集完整版| 大码成人一级视频| 99久久精品国产国产毛片| av在线老鸭窝| 男女免费视频国产| 一级片'在线观看视频| 下体分泌物呈黄色| 久久久精品免费免费高清| 亚洲精品视频女| 国产亚洲最大av| 99热全是精品| 国产伦理片在线播放av一区| 欧美日韩成人在线一区二区| av天堂久久9| 啦啦啦啦在线视频资源| 人人妻人人爽人人添夜夜欢视频| 欧美日韩在线观看h| 人人澡人人妻人| 亚洲国产欧美日韩在线播放| 国产精品一区二区三区四区免费观看| 国产国语露脸激情在线看| 王馨瑶露胸无遮挡在线观看| 日本91视频免费播放| 大话2 男鬼变身卡| 国产亚洲午夜精品一区二区久久| 麻豆乱淫一区二区| 亚洲国产精品999| 美女大奶头黄色视频| 视频区图区小说| 国语对白做爰xxxⅹ性视频网站| 国产成人精品一,二区| 伊人亚洲综合成人网| 精品久久国产蜜桃| 少妇人妻精品综合一区二区| 久久久欧美国产精品| 亚洲美女黄色视频免费看| 一区二区三区四区激情视频| 自线自在国产av| 人体艺术视频欧美日本| 一本—道久久a久久精品蜜桃钙片| 一本大道久久a久久精品| 国产成人精品一,二区| 天天操日日干夜夜撸| 国产视频首页在线观看| a级毛色黄片| 国产成人freesex在线| 99热国产这里只有精品6| 国产日韩一区二区三区精品不卡 | 伦理电影大哥的女人| 色婷婷av一区二区三区视频| 一本久久精品| 成年女人在线观看亚洲视频| 精品国产一区二区三区久久久樱花| 成年美女黄网站色视频大全免费 | 中文字幕人妻熟人妻熟丝袜美| 肉色欧美久久久久久久蜜桃| 大片电影免费在线观看免费| 亚洲五月色婷婷综合| 国产成人freesex在线| 九九久久精品国产亚洲av麻豆| 午夜老司机福利剧场| 成人综合一区亚洲| 2022亚洲国产成人精品| 人人妻人人添人人爽欧美一区卜| 高清毛片免费看| 天堂中文最新版在线下载| 欧美国产精品一级二级三级| 成人二区视频| .国产精品久久| 国产成人精品婷婷| a级毛色黄片| videosex国产| 中文精品一卡2卡3卡4更新| 97精品久久久久久久久久精品| 人妻一区二区av| 亚洲美女视频黄频| 欧美 亚洲 国产 日韩一| 午夜免费男女啪啪视频观看| 亚洲精品,欧美精品| 国产高清有码在线观看视频| 精品亚洲成a人片在线观看| 中文天堂在线官网| 黄色怎么调成土黄色| 久久精品国产a三级三级三级| 91久久精品电影网| 日本av手机在线免费观看| 99久久人妻综合| 边亲边吃奶的免费视频| 青青草视频在线视频观看| 亚洲精品一区蜜桃| 交换朋友夫妻互换小说| 亚洲情色 制服丝袜| 丝瓜视频免费看黄片| 三级国产精品欧美在线观看| 亚洲内射少妇av| 久久精品久久久久久久性| 亚洲中文av在线| 欧美精品国产亚洲| 欧美3d第一页| 99视频精品全部免费 在线| 黑人猛操日本美女一级片| 欧美日韩亚洲高清精品| 亚洲第一区二区三区不卡| 丰满乱子伦码专区| 欧美日韩亚洲高清精品| 精品熟女少妇av免费看| 91久久精品国产一区二区三区| 亚洲国产精品国产精品| 男女啪啪激烈高潮av片| av在线播放精品| 久久国产精品男人的天堂亚洲 | 人人妻人人澡人人看| 亚洲av免费高清在线观看| 99久久精品一区二区三区| 91国产中文字幕| 午夜老司机福利剧场| 国产精品人妻久久久久久| 亚洲欧洲国产日韩| 成年人午夜在线观看视频| 91精品一卡2卡3卡4卡| 国产极品粉嫩免费观看在线 | 99久久精品一区二区三区| av免费观看日本| 免费观看的影片在线观看| 尾随美女入室| 人人妻人人澡人人看| 国产精品一区www在线观看| 99re6热这里在线精品视频| 97超视频在线观看视频| 少妇 在线观看| 黄色一级大片看看| 午夜激情久久久久久久| 韩国高清视频一区二区三区| 人妻制服诱惑在线中文字幕| 国产精品人妻久久久影院| 国产无遮挡羞羞视频在线观看| 大陆偷拍与自拍| 制服人妻中文乱码| av线在线观看网站| 高清黄色对白视频在线免费看| 综合色丁香网| 欧美97在线视频| 最新中文字幕久久久久| 九色亚洲精品在线播放| 亚洲av二区三区四区| 久久久久人妻精品一区果冻| 国产乱人偷精品视频| 亚洲精品,欧美精品| 国产精品久久久久成人av| 色5月婷婷丁香| 一本—道久久a久久精品蜜桃钙片| 妹子高潮喷水视频| 午夜福利在线观看免费完整高清在| 欧美性感艳星| 天天影视国产精品| av免费观看日本| 熟女人妻精品中文字幕| 国产在线免费精品| 激情五月婷婷亚洲| 伊人久久国产一区二区| 这个男人来自地球电影免费观看 | 97在线人人人人妻| 欧美日韩成人在线一区二区| 亚洲成人手机| 国产一区有黄有色的免费视频| 欧美日韩视频精品一区| 久久97久久精品| 国产淫语在线视频| 亚洲美女视频黄频| 国产精品一二三区在线看| 九九在线视频观看精品| 色视频在线一区二区三区| 在线观看一区二区三区激情| 国产精品偷伦视频观看了| 国产亚洲精品久久久com| 欧美+日韩+精品| 亚洲综合色网址| 丝瓜视频免费看黄片| 国产精品偷伦视频观看了| 在线 av 中文字幕| 亚洲av二区三区四区| 丰满迷人的少妇在线观看| 天天影视国产精品| 亚洲国产av新网站| 亚洲成人一二三区av| 成人毛片a级毛片在线播放| 美女福利国产在线| a级片在线免费高清观看视频| 亚洲一级一片aⅴ在线观看| 18禁动态无遮挡网站| 国产精品国产三级专区第一集| 亚洲四区av| 大香蕉久久网| 中国国产av一级| 亚洲第一av免费看| 精品一区二区免费观看| 国产精品一二三区在线看| 久久99一区二区三区| 久久韩国三级中文字幕| 久久av网站| 大片电影免费在线观看免费| 伊人久久国产一区二区| 哪个播放器可以免费观看大片| 麻豆乱淫一区二区| 亚洲精品aⅴ在线观看| 亚洲成人手机| 国产成人freesex在线| 成人影院久久| 欧美97在线视频| 99久久中文字幕三级久久日本| 午夜精品国产一区二区电影| 日本与韩国留学比较| 最近的中文字幕免费完整| 九九在线视频观看精品| 亚洲精品av麻豆狂野| 18+在线观看网站| 国产精品人妻久久久影院| 亚洲图色成人| 日韩伦理黄色片| 三级国产精品片| 男女边摸边吃奶| 亚洲欧美中文字幕日韩二区| av.在线天堂| 亚洲高清免费不卡视频| 日韩,欧美,国产一区二区三区| 在线观看国产h片| 一区二区日韩欧美中文字幕 | 国产午夜精品久久久久久一区二区三区| 国产 精品1| 国产免费福利视频在线观看| 亚洲成人手机| 成年人免费黄色播放视频| 制服人妻中文乱码| 少妇的逼水好多| 一级a做视频免费观看| 亚洲第一av免费看| 我的女老师完整版在线观看| 亚洲美女搞黄在线观看| 国产精品无大码| 久久精品夜色国产| 亚洲av欧美aⅴ国产| 国产亚洲精品久久久com| 大陆偷拍与自拍| 大片免费播放器 马上看| 2018国产大陆天天弄谢| 晚上一个人看的免费电影| 三级国产精品欧美在线观看| 日韩免费高清中文字幕av| 成人免费观看视频高清| 久久青草综合色| 欧美激情国产日韩精品一区| 搡女人真爽免费视频火全软件| 国产精品久久久久久久久免| 国产极品天堂在线| 亚洲人与动物交配视频| 亚洲国产精品一区二区三区在线| 卡戴珊不雅视频在线播放| 亚洲精品美女久久av网站| 综合色丁香网| 91精品三级在线观看| 在线观看国产h片| kizo精华| 成人18禁高潮啪啪吃奶动态图 | 国产男女超爽视频在线观看| 午夜福利视频精品| 国产成人freesex在线| 校园人妻丝袜中文字幕| 国产成人一区二区在线| 久久人人爽av亚洲精品天堂| kizo精华| 搡老乐熟女国产| 蜜桃在线观看..| 亚洲欧洲精品一区二区精品久久久 | 日韩欧美精品免费久久| 伦理电影大哥的女人| 男男h啪啪无遮挡| 国产成人精品福利久久| 久久精品久久精品一区二区三区| 人妻一区二区av| 青春草亚洲视频在线观看| 午夜视频国产福利| 黑人猛操日本美女一级片| 午夜免费观看性视频| 免费久久久久久久精品成人欧美视频 | 搡老乐熟女国产| 熟女av电影| 9色porny在线观看| 亚洲无线观看免费| 色视频在线一区二区三区| 欧美激情国产日韩精品一区| 黄片播放在线免费| 考比视频在线观看| 成人影院久久| 久久av网站| 校园人妻丝袜中文字幕| 亚洲精品美女久久av网站| 久久国产亚洲av麻豆专区| 黄色一级大片看看| 美女福利国产在线| 老司机影院毛片| 亚洲国产最新在线播放| 国产av码专区亚洲av| 高清欧美精品videossex| 一级毛片黄色毛片免费观看视频| 亚洲四区av| 天天操日日干夜夜撸| 99热6这里只有精品| 大香蕉久久成人网| 国产一级毛片在线| 观看av在线不卡| 女的被弄到高潮叫床怎么办| 日韩精品有码人妻一区| 亚洲精品一区蜜桃| 一本色道久久久久久精品综合| 亚洲av成人精品一区久久| 午夜激情久久久久久久| 考比视频在线观看| 制服人妻中文乱码| 久久久久国产网址| 一区二区三区四区激情视频| 黄色一级大片看看| 高清午夜精品一区二区三区| 99九九在线精品视频| 欧美日韩一区二区视频在线观看视频在线| 一本大道久久a久久精品| 满18在线观看网站| 简卡轻食公司| 欧美日韩成人在线一区二区| 婷婷色av中文字幕| 免费观看a级毛片全部| 亚洲美女搞黄在线观看| 七月丁香在线播放| 日韩视频在线欧美| 午夜久久久在线观看| 亚洲精华国产精华液的使用体验| 18禁在线播放成人免费| 欧美日本中文国产一区发布| 国产精品久久久久久av不卡| 国产精品久久久久久久电影| 免费日韩欧美在线观看| 男女啪啪激烈高潮av片| 啦啦啦啦在线视频资源| 久久久久精品久久久久真实原创| 97在线人人人人妻| 在线观看三级黄色| 蜜桃在线观看..| 国产高清有码在线观看视频| 亚洲av二区三区四区| 免费观看的影片在线观看| 22中文网久久字幕| 日本黄色片子视频| 国产极品粉嫩免费观看在线 | 国产成人精品无人区| 狂野欧美白嫩少妇大欣赏| 男人操女人黄网站| 国产深夜福利视频在线观看| 亚洲伊人久久精品综合| 国产日韩欧美在线精品| 97在线视频观看| 亚洲性久久影院| 精品少妇内射三级| 肉色欧美久久久久久久蜜桃| 精品人妻熟女av久视频| 亚洲美女黄色视频免费看| 久久精品国产a三级三级三级| 啦啦啦在线观看免费高清www| 在线播放无遮挡| 91精品国产国语对白视频| 国产有黄有色有爽视频| 最近2019中文字幕mv第一页| 哪个播放器可以免费观看大片| 国产乱人偷精品视频| 日韩av免费高清视频| 欧美人与善性xxx| 午夜91福利影院| 国产女主播在线喷水免费视频网站| 国产精品一二三区在线看| 国产国语露脸激情在线看| 国产成人精品一,二区| 国产精品一区二区在线不卡| 国产成人freesex在线| 亚洲综合精品二区| 日韩av在线免费看完整版不卡| 欧美丝袜亚洲另类| 亚洲av成人精品一二三区| 免费人成在线观看视频色|