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

    Tartary buckwheat (Fagopyrum tataricum (L.) Gaertn) protein-derived antioxidant peptides: mechanisms of action and structure-activity relationship in Caco-2 cell models

    2022-07-11 05:49:08YimingZhouXunmingSheZhidongChenYunWeiYingXioXioliZhou

    Yiming Zhou, Xunming She, Zhidong Chen, Yun Wei, Ying Xio, Xioli Zhou,b,*

    a School of Perfume and Aroma Technology, Shanghai Institute of Technology, Shanghai 201418, China

    b University Think Tank of Shanghai Municipality, Institute of Beautiful China and Ecological Civilization, Shanghai 201418, China

    ABSTRACT

    Excessive reactive oxygen species (ROS) can cause oxidative damage and lead to various metabolic disease.Tartary buckwheat (Fagopyrum tataricum (L.) Gaertn) is a new kind of protein-rich functional food, the protein in which has been proved to have good antioxidant capacity. In this study, in order to further explore the antioxidant mechanism of Tartary buckwheat protein, 4 peptides (CR-8, LR-8, GK-10 and SR-12) were isolated and identified from it. H2O2 was used to induce oxidative damage to Caco-2 cells to evaluate antioxidant capacity of these peptides. The results of superoxide dismutase (SOD), total antioxidant capacity (T-AOC)and mitochondrial membrane potential etc. showed that these peptides have superior antioxidant capacity.CR-8 has the best antioxidant capacity. In order to further clarify the antioxidant mechanism of CR-8,metabolomics was used to analyze related metabolites and metabolic pathways. The results showed that after CR-8 intervention, the content of metabolites such as L-acetyl carnitine has increased. This indicated that CR-8 can improve the antioxidant capacity of damaged cells by intervening in multiple metabolic pathways. This also revealed the anti-oxidant mechanism of tartary buckwheat protein. In conclusion, it provided a theoretical basis for further studying the activity of tartary buckwheat portein and utilizing buckwheat resources.

    Keywords:

    Antioxidant protective

    Tartary buckwheat peptides

    Metabolic mechanim

    Caco-2 cells

    1. Introduction

    Reactive oxygen species (ROS) are oxygen-derived small molecules [1] including hydrogen peroxide and the hydroxyl radical; they are generated from normal metabolic processes. ROS are chemically reactive molecules that have essential functions in living organisms [2]. A moderate increase in ROS can promote cell proliferation and differentiation. ROS at low levels can act as cell signaling molecules by reversibly oxidizing thiol groups and function [3]. However, high levels of ROS can harm cells by nonspecifically damaging macromolecules including DNA, proteins,and lipids [4,5]. ROS can lead to multiple metabolic disease like cardiovascular disease [6], type-2 diabetes [7], Alzheimer’s disease [8], and certain types of cancer [9].

    Maintaining ROS homeostasis is crucial for normal cell growth and survival. Cells control ROS levels by balancing ROS generation with ROS elimination via ROS-scavenging systems such as superoxide dismutase, glutathione peroxidase, peroxiredoxins,glutaredoxin, thioredoxin and catalase [5]. Food also contains many biologically active substances that scavenge free radicals [10].Currently, the enzymatic hydrolysis of peptides from proteins has become a research hotspot, including plant proteins and animal proteins. The molecular weight of peptides are between amino acids and proteins and have rich biological activities including antioxidant,cholesterol-lowering, and antibacterial effects [11].

    Tartary buckwheat is a functional food and its protein has excellent biological activity. Compared with other proteins, it has a unique amino acid composition [12]. In the previous study, we have isolated 13S globulin from Tartary buckwheat and analyzed its structural characteristics. And we also investigated the protective effect of acidic subunit enzymatic hydrolysates of Tartary buckwheat 13S globulin against oxidative stress induced by H2O2[13]. In this study, to further evaluate the antioxidant effect of these peptides, we used Caco-2 cells to construct an oxidative stress model. At the same time, we used metabolomics to deeply analyze the in fluence of these peptides on related metabolic pathways.

    PyMOL is one of the most widely used molecular visualization systems: It is a simple and powerful tool for sequence and structure analysis and homology modeling [14] and is suitable for drug development [15] including the evaluation of binding sites and interaction forces, protein preparation, and molecular mechanics [16].

    Metabolomics can analyze the changes in the body’s metabolites under the effects of endogenous or exogenous stimuli and interference such as diseases and drug treatments [17]. It can analyze changes in cell endogenous metabolites, find biomarkers that can reflect the state of the body, and help people understand related metabolic pathways [18]. Currently, liquid chromatography mass spectrometry (LC-MS), gas chromatography mass spectrometry (GC-MS)and capillary electrophoresis mass spectrometry (CE-LC) are among the most commonly used tools in chromatographic analysis of the metabonomics [19]. This experiment used high performance liquid chromatography mass spectrometry (HPLC-MS) to analyze the changes in the metabolism of Caco-2 cells after the intervention of Tartary buckwheat peptides to evaluate the effect of Tartary buckwheat peptides on the cells.

    2. Materials and methods

    2.1 Materials

    The human colon adenocarcinoma Caco-2 cell line was obtained from Chinese Academy of Sciences. Dulbecco’s Modified Eagle Medium (DMEM), nonessential amino acid (NEAA),Hanks’ balanced salt solution (HBSS) were products of Gibco Life Technologies. Malondialdehyde (MDA) determination kit,catalase (CAT) determination kit, superoxide dismutase (SOD)determination kit and total antioxidant capacity (T-AOC)determination kit were purchased from the Beyotime Institute of Biotechnology (Shanghai, China). Synthesized peptides,namely, CTGFVAVR, LRGENDQR, GFIVQAQDLK, and SFFLAGQSQQGR, with purities higher than 98% were provided by GL Biochem Ltd. (Shanghai, China). All chemicals were of chromatographic or analytical quality and purchased from Sinopharm Chemical Reagent (Shanghai, China) unless mentioned otherwise.

    2.2 Structure and property analysis of Tartary buckwheat peptides

    The identified amino acid sequences were input into PyMOL to simulate the three-dimensional structure of Tartary buckwheat peptides. Then the relationship between structure and antioxidant function was analyzed.

    2.3 Cell culture and treatments

    Caco-2 cells were obtained from the Chinese Academy of Sciences. All cells were routinely cultured in DMEM medium supplemented with 10% fetal bovine serum (FBS). The cytoprotective effect of peptides against H2O2-induced oxidative damage was evaluated in parallel experiments. Caco-2 cells were incubated with synthesized peptides as indicated above followed by stimulation with 600 μmol/L H2O2for 4 h. To determine intracellular enzymatic viability, Caco-2 cells (5 × 104viable cells/well) were seeded onto 96-well black, clear-bottom plates and allowed to attach for 24 h.The cells were then stimulated with 600 μmol/L H2O2for an extra 4 h followed by incubation with synthesized peptides (500 μmol/L).According to the above processing method, the cells were divided into three main groups. They were control group which referred to the Caco-2 cells without treatment. The model group referred to Caco-2 cells under oxidative stress which were induced by H2O2.The intervention group referred to the Caco-2 cells under oxidative stressed which were intervened by Tartary buckwheat peptides. The intervention group were further divided into the CR-8 group, LR-8 group, GK-10 group and SR-12 group.

    2.4 3-(4,5-Dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) assay

    After incubation, the medium was removed and then 100 μL of serum-free medium containing MTT solution at afinal concentration of 5 mg/mL was added to each well. Cells were then incubated for another 4 h in a humidified atmosphere of 95% air with 5% CO2at 37 °C. Then, the medium was removed and 100 mL of DMSO was added into each well to dissolve the purple formazan crystals.Finally, the absorbance of each well was measured at 490 nm using a microplate reader. Cell viability was expressed as the percentage of the non-treated control, which was set to 100% [20].

    2.5 The measurement of apoptosis and necrosis

    Referred to Ormerod et al. [21] for measurement of apoptosis and necrosis with some adjustment. The cells to be tested were washed and centrifuged, and mixed with fluorescent dyes Hoechst 33342 and Propidium Iodide (PI) in at a volume of 1:1; they were then added to the cells for staining at 4 °C. The cells were collected by enzymatic digestion after the intervention, and a double staining method was used to determine the degradation and necrosis of the cells. The dark wavelength was used for 20–30 min. Flow cytometry analysis included two channels red and blue; there were 5 000 events.

    2.6 Effect of Tartary buckwheat peptides on cell growth cycle

    To judge the effect of Tartary buckwheat antioxidant peptides on cell oxidative stress, cells (2 × 105) were collected after the intervention. The Caco-2 cells were digested and blown evenly from a 12-well plate, washed with PBS, and centrifuged for 5 min. Then 70% precooled ethanol was added and fixed at 4 °C for overnight. And then, the cells were centrifuged to remove the supernatant and added 500 μL RNase A (50 μg/mL) at 37 °C for 30 min to digest the RNA in the cells. After PBS washing and adding PI dye (50 μg/mL), the cells were stained at 4 °C for 30 min.

    2.7 The measurement of cells membrane potential

    After digesting and collecting the cells, we added 1 mL of JC-1 staining working solution and mixed well. The sample was incubated for 20 min at 37 °C in a cell incubator. After incubation,the supernatant was aspirated and washed twice with JC-1 staining buffer; 2 mL of cell culture solution was added mixed well, and the Caco-2 cells were observed with through an inverted fluorescence microscope.

    2.8 Determination of T-AOC, CAT, SOD, and MDA in cells

    Caco-2 cells were cultured in 24-well plates. After the cells adhered to the wall and spread over the bottom, the samples were interfered for 24 h. Then cells were then collected and crushed with an ultrasonic pulverizer. The T-AOC, MDA, SOD, and CAT activity were evaluated using an intracellular determination kit. The relevant tests were carried out in accordance with the instructions provided by the kit.

    2.9 Data processing and metabolites identification of metabolomics

    2.9.1 Sample preparation

    The intracellular metabolites were extracted according refer to Li et al. [22]. Samples were reconstituted in 200 μL MeOH/H2O (85:15,V/V); 100 μL was diluted to 50% (V/V) MeOH in water for non-targeted metabolomics study followed by 30 μL of sample pooled together from each participant as a quality control (QC) sample. The other 100 μL was used for targeted metabolomics. All samples were detected within 24 h of reconstitution.

    2.9.2 HPLC-MS conditional selection

    Referred to Ma et al. [23] for chromatographic separation and massed spectrometry detection conditions with some adjustment.HPLC-MS was used in this experiment and the related conditions were optimized.

    The optimal conditions of HPLC were set as followed: the metabolite was fractionated using a C18column with a flow rate of 0.3 mL/min. Mobile phase A was 0.1% formic acid in HPLC water,mobile phase B was 0.1% formic acid in ACN and used 40 °C as column temperature. The solvent gradient changed according to the following conditions: 0–2 min, 95% A; 2?12 min, 95%?5% A;12?15 min, 5% A; 15?17 min, 5% A.

    The optimal conditions of MS were set as followed: positive ion mode, heater temperature: 300 °C; capillary temperature: 350 °C;electrospray voltage: 3.0 kV; exhaust gas flow rate: 1 arb; sheath gas flow rate: 45 arb; auxiliary gas flow rate: 15 arb; S-Lens RF Level,30%. Negative ion mode: heater temperature: 300 °C; capillary temperature: 350 °C; electrospray voltage: 3.2 kV; exhaust gas flow rate: 1 arb; sheath gas flow rate: 45 arb; auxiliary gas flow rate: 15 arb;S-Lens RF Level, 60%.

    2.10 Statistical analysis

    Data were expressed as means ± standard deviation (SD). All of the statistical analyses were carried out using SPSS 20.0 software(SPSS Inc., Chicago, IL, USA). Multiple group comparisons were performed using a one-way analysis of variance (ANOVA). The level of statistical significance was set at P < 0.05.

    The mass spectrometry data were converted into GDF format by Chroma TOF software and imported into SIMCA-P software,version 11.0 for principal component analysis (PCA), partial least-squares discrimination analysis (PLS-DA), orthogonal partial leastsquares-discrimination analysis (OPLS-DA) and other multivariate statistical analysis.

    3. Results

    3.1 Structure and properties of Tartary buckwheat peptides

    The 4 peptides consist of 8–12 amino acids and have a molecular weight of about 800–1 400 Da. The properties of amino acids can determine the activity of peptides [24], and 4 peptides have been identified as containing amino acids with antioxidant capacity.Therefore, these 4 peptides may have antioxidant properties, but the antioxidant activity needs to be further verified. These Tartary buckwheat peptides all have random coils and β-sheets, and some peptides have α-helices. These structures greatly affect the properties of the peptides. The sequence of antioxidant peptides included a high hydrophobic moiety. This moiety can scavenge lipid-derived radicals due to interaction with lipid molecules [25]. Small peptides were extraordinarily potent in preventing cell death due to pro-oxidants [26].The peptides contained 27 amino acids (3 037.49 Da) and showed antioxidant activity that positively correlated with concentration; there was no antimicrobial activity. The small size and hydrophobic amino residues contributed to the high activity [24,27] (Fig. 1).

    Fig. 1 Three-dimensional simulated structure of Tartary buckwheat peptides.

    Theoretical isoelectric point, hydrophobicity, and net charge of each peptide was calculated using ExPASyf (http://ExPASyf.com).

    3.2 Determination of cell viability and oxidative damage model

    3.2.1 Cell viability

    The toxicity of Tartary buckwheat peptides to cell were evaluated including the protective effect of Tartary buckwheat peptides on oxidative damage to Caco-2 cells. The results were shown in Fig. 2.The cell viability measured by the MTT method showed that at a concentration of 500 μg/mL, the cell viability of Tartary buckwheat peptides intervention was not significantly different from that of the control group. Cells can neither show inhibitory viability nor promote cell growth, i.e., Tartary buckwheat peptides do not impact oxidation.The peptides hydrolyzed by alkaline protease were relatively safe for the cell event at 500 μg/mL.

    Fig. 2 Construction of cell viability and cell damage model. (A) The effect of Tartary buckwheat peptides content on cell viability. (B) The effect of different concentrations of H2O2 and intervention time on Caco-2 cells.Different lowercase letters indicate significant differences (P < 0.05).

    3.2.2 Construction of oxidative damage model

    H2O2can induce oxidative stress and damage to cells. In order to determine the amount of H2O2used and the intervention time.Under the conditions of different concentration gradients of H2O2and intervention time, a cell viability of about 50% was determined as the best intervention condition. The experimental results were shown in the Fig. 2. The survival rate of Caco-2 cells was related to the concentration of H2O2and the intervention time. With the extension of the intervention time, the cell viability showed a downward trend.After the intervention for 4 h, the cell viability under different H2O2concentrations tended to be flat. When the intervention concentration of H2O2was 800 μmol/L, the cell viability could reach 55% relative to the control group. Therefore, the best result was when the intervention time was 4 h, the H2O2concentration was 800 μmol/L.

    3.3 Effects of Tartary buckwheat peptides on Caco-2 cells

    3.3.1 Analysis of cell apoptosis and necrosis

    The apoptosis and necrosis of cells can be distinguished by flow cytometry (Fig. 3A). The apoptosis and necrosis rate of the cells treated with H2O2increased significantly, which indicated that the cells were damaged and gradually necrotized under the effect of H2O2.Long term oxidative stress state would lead to irreversible damage of most organelles: ROS damaged double-stranded DNA and caused apoptosis [28]. The intervention of Tartary buckwheat peptides can effectively reduce the proportion of apoptotic cells, which showed that Tartary buckwheat peptides could act on the early apoptotic state of cells and improve the viability of cells; however, this could not reach the level of normal cells, i.e., the control group, in which the effect of CR-8 peptide was the best, and that of GK-10 was second.Compared with the cells in the early stage of apoptosis, Tartary buckwheat peptides had no obvious effect on the cells in the late stage of necrosis.

    3.3.2 Analysis of cell cycle

    The effect of Tartary buckwheat peptides on the Caco-2 cell cycle were shown in Fig. 3B. Compared with the control group, the Caco-2 cells in the G1 phase of the model group were significantly reduced.This was due to the concentration of the nucleus and fragmentation of the DNA, resulting in a decrease in the fluorescence intensity of the cell DNA. H2O2would induce oxidative fragmentation of DNA,which would eventually lead to cell death [29]. At the same time, the proportion of cells in the S phase increases, indicating that the cells were blocked at this stage. After the intervention of Tartary buckwheat peptides, the proportion of cells in the G1 phase increased, and the S phase decreased. These changes suggested that the cells resumed their normal growth cycle. The smaller molecular weight meant that the peptides could more easily penetrate the cell membrane.

    3.3.3 Mitochondrial membrane potential

    Under oxidative stress, the mitochondrial membrane potential decreases [30]. After the intervention of Tartary buckwheat peptides,the mitochondrial membrane potential of most cells was significantly improved (Fig. 3C). The cell mitochondrial viability after CR-8 intervention was the highest. The cells were covered with clear red fluorescence. SR-12 had the worst effect, and the red fluorescence distribution was scattered and weak. Changes in mitochondrial viability affected mitochondrial dynamics leading to a loss in membrane potential, increased production of reactive oxygen species,and depletion of ATP. This led to severely impaired mitochondrial function, mitochondrial dysfunction-mediated apoptosis [31].Buckwheat peptides could significantly improve the vitality of damaged mitochondria and gradually restore the cells to normal.

    Fig 3 Effects of Tartary buckwheat polypeptide on Caco-2 cells. (A) The effects of Tartary buckwheat peptides on apoptosis and necrosis rate of Caco-2 cells under oxidative stress. It shows the distribution of apoptosis and necrosis of different groups of Caco-2 cells. Numbers adjacent to the outline indicate the percentage of cells in each area. (B) The effects of Tartary buckwheat peptides on Caco-2 cell cycle. Numbers indicate the percentage of cells in each stage.(C) The effects of Tartary buckwheat peptides on mitochondrial membrane potential of Caco-2 cells. The distribution of red fluorescence represents the viability of cell mitochondria. (1) Control group; (2) H2O2 treated group; (3–6) Tartary buckwheat peptides (CR-8, LR-8, GK-10, SR-12) treated before the same H2O2 stimulation. In thefigures, red circles were used to circle the area where the membrane potential changes more obviously.

    Fig 3 (Continued)

    3.4 Cytoprotective effect of Tartary buckwheat peptides against H2O2-induced oxidative stress in Caco-2 cells

    The antioxidant capacity of cells is composed of enzymatic and non-enzymatic systems. The total antioxidant capacity can re flect the health of the cells [32]. The total antioxidant capacity of Caco-2 cells decreased significantly (P < 0.05) under the induction of hydrogen peroxide (Fig. 4A). The total antioxidant capacity of Caco-2 cells could be improved by Tartary buckwheat peptides intervention, and CR-8 had the best effect: it was 51.38% higher than the H2O2group,but 33.53% lower than the control group. The molecular weight and length of Tartary buckwheat peptides could affect its biological activity, but there seemed to be no definite relationship between activity and weight/length. A relatively small molecular weight can more easily act intracellularly with better antioxidant activity.However, the relatively long peptide (SR-12) also had a better effect,which showed that the amino acid composition of peptides had a large in fluence on antioxidant activity.

    The CAT activity (Fig. 4B) in the model group decreased by 45.21%, compared with the control group, and H2O2decreased the activity of antioxidant enzymes. All four peptides could significantly increase the activity of CAT enzyme in injured cells. Especially, after the CR-8 intervention, the CAT activity increased by 64.66%, and the activity of SOD (Fig. 4C) was similar to CAT. The activity of the antioxidant enzyme in the CR-8 group increased by 64.66% compared with the model group, but the activity of the antioxidant enzyme in the Tartary buckwheat peptides group was still far from that in the control group.

    Fig. 4 Cytoprotective effect of Tartary buckwheat peptides against H2O2-induced oxidative stress in Caco-2 cells. The bar graph shows the impact of Tartary buckwheat peptides on different groups from four different aspects. (A) The effect of Tartary buckwheat peptides on the total antioxidant capacity of Caco-2 cells under oxidative stress. (B) Effect of buckwheat peptides on CAT enzyme activity of Caco-2 cells under oxidative stress. (C) Effect of Tartary buckwheat peptides on SOD activity in Caco-2 cells under oxidative stress. (D) The effect of Tartary buckwheat peptides on MDA content in Caco-2 cells under oxidative stress.Different lowercase letters indicate significant differences (P < 0.05).

    The MDA content (Fig. 4 D) in the model group was significantly higher than the control group (47.33%, P < 0.05) perhaps due to the peroxidation of unsaturated fatty acids in the cell biofilm by H2O2, and the production and dissociation of a high concentration of MDA in the cells. This enhanced lipid peroxidation damages the cell membrane structure. Tartary buckwheat peptides could significantly reduce the content of MDA in cells and gradually returned them to the level of MDA in the control group. This showed that the lipid peroxidation in Caco-2 cells had been reduced, and the function of the cell membrane had been restored. This result may be related to the hydrophobicity of peptides, i.e., the peptides with a higher hydrophobicity were more likely to act on the cell membrane composed of lipids to remove accumulated peroxides, prevented the peroxidation of unsaturated fatty acids, and maintained the normal structure and function of the cell membrane.

    Based on the above experimental results, Tartary buckwheat peptides could improve cell activity by improving cell mitochondrial activity and avoid early cell apoptosis. In addition, Tartary buckwheat peptides could significantly improve the antioxidant capacity of Caco-2 cells under oxidative stress. However, different Tartary buckwheat peptides had different effects. Among them, CR-8 had the best effect.

    3.5 Multivariate statistical analysis of cell metabolomics

    In order to further explore the protective mechanism of Tartary buckwheat peptides on Caco-2 cells under oxidative stress. We selected the CR-8 group with the best experimental results as the samples for metabolomics experiments, and conducted in-depth identification and analysis of its metabolites. The differences among the samples from different cell groups were examined using principal components analysis (PCA). After mass spectrometry analysis, the PCA scores of the positive ion mode and the negative ion mode were shown in Fig. 5. In order to further differentiate the data and highlight the differences between groups, the data was analyzed by PLS-DA and OPLS-DA. At the same time, the permutation test was used to verify the reliability of the model with the intercept of the fitted regression line on the Y axis being less than 0.

    Fig. 5 The mean PCA, PLS-DA, Replacement verification plot and OPLS-DA of the control, the H2O2 and CR-8 groups. (A) Score graph in negative ion mode.(A1) PCA score plot (R2X = 0.607, Q2 = 0.439). (A2) PLS-DA score plot (R2X = 0.53, R2Y = 0.985, Q2 = 0.931). (A3) Replacement verification plot, intercept of the regression Q2 = -0.260 773. (A4) OPLS-DA score graph. (B) Score graph in positive ion mode. (B1) PCA score plot (R2X = 0.729, Q2 = 0.594). (B2) PLS-DA scoreplot (R2X = 0.729, R2Y = 0.99, Q2 = 0.948). (B3) Replacement verification plot, intercept of the regression Q2 = -0.260 773. (B4) OPLS-DA score graph.

    In the negative ion mode, a total of 3 035 metabolites were obtained. Three principal components are obtained by PCA. In the positive ion mode, a total of 1 383 metabolites were obtained, and 5 main components were obtained by PCA. Through PLS-DA, OPLS-DA and Permutation test, we found that the three groups were more distinct in the negative ion mode, while the samples in the group were relatively concentrated. However, the quality control group clustered better in the positive ion mode. However, the experimental group was not obvious, and there was a cross-stacking phenomenon in the control group, model group and the intervention group, indicating that there was a certain similarity in the metabolites in the positive ion mode. It was further confirmed by PLS-DA and OPLS-DA.

    Based on the above analysis, although the PCA score chart in the negative ion mode showed more obvious clustering, the cumulative principal component score was lower than that in the positive ion mode,and the modelfit was worse than that in the positive ion mode. Therefore,we chose the metabolites in positive ion mode for further analysis.

    The control group and the intervention group were well distinguished, and the model group was more scattered. Compared with the control group, the model group was closer to the intervention group, indicating that the two groups of metabolites detected in the positive ion mode were similar to a certain degree, which may be due to the same metabolic pathways, resulting in similar scatter scores.After H2O2intervention, the cell metabolites were significantly different and the points were relatively scattered, indicating that the metabolites were complex, and the relatively concentrated scatter score chart indicated that the internal clustering was better. This may be different from the pathway of cellular metabolism that peptides participate in. It showed that CR-8 could affect multiple cellular metabolic pathways to a certain extent.

    3.6 Identification of differential metabolites

    In this study an S-plot map of OPLS-DA was used to screen the differential metabolites under the conditions of |P(corr)| > 0.5 and|P| > 0.1. We used the FC value (fold change) to observe the change in differential metabolites between groups. The qualitative analysis of metabolites was based on the comparison ofm/zand HMDB. As it showed in Fig. 6, the control group and the model group involved relatively few differential metabolites, while the difference between the model group and the CR-8 group increased significantly. This may indicate that Tartary buckwheat peptide CR-8 can interfere with multiple metabolic pathways to achieve antioxidant effects.

    Fig. 6 Metabolic difference analysis between control group, model group and CR-8 group. (a) OPLS-DA score plot of the control group vs. model group.(b) OPLS-DA S-plot of the control group vs. model group. (c) OPLS-DA score plot of the model group vs. CR-8 group. (d) OPLS-DA S-plot of the model group vs. CR-8 group.

    The major potential metabolites in the positive ion models were screened in samples of cells (Table 1). The analyzed metabolic differences mainly included: L-acetylcarnitine; L-tryptophan;acetylspermidine; lysophosphatidic acid (20:4); lysophosphatidic acid (20:3); lysophosphatidic acid (16:0); pantothenic acid;glycerophosphoethanolamine; L-histidine; L-alanine. These metabolites involve multiple metabolic pathways. Table 1 showed the logarithm of the ratio of the model group to the control group and the intervention group to the model group (base 2). A positive number means an increase, a negative number means a decrease. Among them, L-histidinol and L-aspartic acid levels increased. It showed that the change of free amino acid content could directly affect cell growth and proliferation, and also confirmed the experimental results of apoptosis and necrosis and cell cycle.

    Table 1The effect of CR-8 on the marker differential metabolites and metabolic pathway of Caco-2 cells under oxidative stress.

    3.7 Metabolic pathway analysis

    Some metabolic pathways have significant effects on the oxidative stress of Caco-2 cells: fatty acid metabolism; glycerol phospholipid metabolism; tryptophan metabolism; pantothenic acid and coenzyme A biosynthesis; histidine metabolism; and alanine, aspartic acid and glutamate metabolism (Fig. 7).

    Fig. 7 Tartary buckwheat peptides involved in metabolic pathways.Summary of pathway analysis using MetaboAnalyst 3.5. Each point represents one metabolic pathway; dot size and shade of color positively correlate with the impact of the metabolic pathway.

    Among them, glycerol phospholipid metabolism had the highest contribution value, and glycerol phospholipid has an important in fluence on cell membrane function. This also proved that CR-8 can improve Caco-2 cells which under oxidative stress by interfering with glycerol phospholipids and other related metabolic pathways.

    4. Discussion

    A comparison of the structure of Tartary buckwheat peptides suggest that Val can mprove the CAT activity in the cell when the N-terminal and ASP are at the C-terminal of the peptides these conclusions are consistent with the results whereby hydrophobic amino acids such as Leu, Pro, Phe, and Val play an important role in the antioxidant activity of the peptides [26,33]. ROS involved in cell growth and metabolism [34]. ROS accumulation can damage organelles including mitochondria, endoplasmic reticulum, and cell membrane [35]. Caco-2 cell necrosis during oxidative stress, may be altered by Tartary buckwheat peptides. These peptides affect cells with mild damage but not cells with severe damage.

    The peptides sequence has a stronger antioxidant activity at about 8 amino acid [24]. Hydrophilic amino acids also affect antioxidant activity [36]. For example, the total antioxidant capacity of SR-12 is next to that of CR-8, which may be related to the high content of glutamine [37]. However, CR-8 had strong free radical scavenging ability, which was related to the N-terminal cysteine (Cys). Since tryptophan had lost H, it could show stronger oxygen free radical binding ability.

    L-acetylcarnitine is mainly involved in fatty acid metabolism [38].Its main function is to transfer fatty acids from cytoplasm to mitochondria for β-oxidation and to store energy in the form of acyl carnitine Disorders of lipid metabolism could cause abnormal TCA cycle intermediate products [39], elevated free fatty acids [40], and changed in amino acid content [41] may directly affect cell growth and proliferation. Compared with the control group, the carnitine content of the model group decreased. This may be caused by H2O2inducing oxidative damage to the cells, which made the cell fatty acid metabolism disorder. After CR-8 intervention, the content of L-acetylcarnitine increased, lipid metabolism gradually returned to normal, and cell oxidative damage was improved. Cell cycle,apoptosis and necrosis experiments can also prove it. In fatty acid metabolism, fatty acid synthase is a key enzyme in the process of fatty acid synthesis, which catalyzes the synthesis of endogenous long fatty acids [42]. The increase in the expression of fatty acid synthase will significantly increase the deposition of triglycerides in the body.

    Lysophosphatidic acid (LPC) is involved in glycerol phospholipid metabolism [43]. It is mainly regulated by lysophosphatidylcholine transferase and phospholipase in the blood [44]. The PLA2/LPC pathway is one of the most important metabolic pathways in the metabolism of glycerophospholipids. And the metabolism of glycerophosphatidyl plays an important role in lipid metabolism [45].Compared with the control group, the LPC content of the model group decreased, which was caused by severe oxidative damage.The increase in LPC content in the CR-8 group indicated that the antioxidant capacity of cells and the lipid oxidation situation has been improved. Because PLA2 not only maintains cell membrane phospholipid homeostasis, it also plays an important role in regulating oxidative damage and inflammation under physiological and pathological conditions [37]. This proved that Tartary buckwheat peptides can regulate lipid metabolism to a certain extent, and protect cells under oxidative damage.

    The content of alanine, histidine and tryptophan changed, and the amino acid content of the intervention group increased compared with the model group. This may indicate that CR-8 can interfere with the amino acid metabolism pathway of the cell and improve cell viability.The results of cell cycle experiments can also prove this.

    Therefore, Tartary buckwheat polypeptide represented by CR-8 can increase the activity of Caco-2 cells under oxidative stress. The main principle is to improve cell viability by improving the state of organelles. Because different kinds of Tartary buckwheat polypeptides have different effects and can only partially restore cells. This may be related to the intervention concentration of Tartary buckwheat polypeptide, so further research is needed.

    In conclusion, this study provided new evidence that Tartary buckwheat peptides has an effect on cells with mild damage from ROS. The CR-8 peptide has the best effect and contains a peptides sequence with 8 amino acids. We analyzed the related metabolites of lipid oxidation, and the intervention of Tartary buckwheat peptides CR-8 partially improved lipid metabolism and restored cells to a normal state. At the same time, we analyzed the metabolic pathways related to amino acids, because the Tartary buckwheat peptides participates in the related metabolic pathways. These pathways improve cell viability and promote cell growth. These results may influence cell reparation and affect multiple metabolic diseases.Therefore, this study will provide a research direction for the development of Tartary buckwheat in food science.

    Con flicts of interest

    The authors confirm that there are no con flicts of interest in this manuscript.

    Acknowledgements

    The authors thank Shanghai Natural Science Foundation(20ZR1455800), the National Science Foundation of China(31871805), China Agriculture Research System (CARS-08-D2)and Shanghai Municipal Education Commission (Plateau Discipline Construction Program).

    七月丁香在线播放| 中文字幕最新亚洲高清| 丝袜美腿诱惑在线| 国产爽快片一区二区三区| 不卡av一区二区三区| 丰满少妇做爰视频| 国产免费视频播放在线视频| 大码成人一级视频| 一级片'在线观看视频| 一级毛片电影观看| 亚洲伊人久久精品综合| 操美女的视频在线观看| 欧美黄色淫秽网站| 国产亚洲欧美精品永久| 亚洲欧洲日产国产| 欧美日韩视频高清一区二区三区二| 啦啦啦在线免费观看视频4| 久久ye,这里只有精品| 精品久久蜜臀av无| av一本久久久久| 18禁裸乳无遮挡动漫免费视频| 国产激情久久老熟女| 多毛熟女@视频| 激情五月婷婷亚洲| 夫妻午夜视频| 高潮久久久久久久久久久不卡| 国产成人a∨麻豆精品| 美国免费a级毛片| a 毛片基地| 少妇 在线观看| 好男人电影高清在线观看| 国产精品香港三级国产av潘金莲 | 免费观看av网站的网址| 1024视频免费在线观看| 亚洲色图 男人天堂 中文字幕| 又大又爽又粗| 黄色片一级片一级黄色片| 亚洲国产av影院在线观看| 精品一区二区三区四区五区乱码 | 国产精品久久久av美女十八| 国产无遮挡羞羞视频在线观看| 亚洲av综合色区一区| 三上悠亚av全集在线观看| 人人妻,人人澡人人爽秒播 | 伊人亚洲综合成人网| 99香蕉大伊视频| 日韩伦理黄色片| 18禁裸乳无遮挡动漫免费视频| 久久ye,这里只有精品| 久久久久国产一级毛片高清牌| 少妇人妻久久综合中文| 成人国语在线视频| 亚洲精品久久久久久婷婷小说| 18禁黄网站禁片午夜丰满| 午夜老司机福利片| 国产一区亚洲一区在线观看| 2018国产大陆天天弄谢| 中国国产av一级| 一级毛片女人18水好多 | 久久精品人人爽人人爽视色| 亚洲一码二码三码区别大吗| 欧美日韩一级在线毛片| 国产欧美日韩综合在线一区二区| 视频区图区小说| 老司机午夜十八禁免费视频| 成人黄色视频免费在线看| 久久久久视频综合| 国产亚洲欧美精品永久| 国产成人av激情在线播放| 国产成人a∨麻豆精品| 日本欧美国产在线视频| 丝袜在线中文字幕| 婷婷色综合大香蕉| 午夜视频精品福利| 亚洲欧美日韩另类电影网站| 不卡av一区二区三区| 久久女婷五月综合色啪小说| 亚洲国产毛片av蜜桃av| 91成人精品电影| 蜜桃在线观看..| 性高湖久久久久久久久免费观看| 人体艺术视频欧美日本| 老鸭窝网址在线观看| 国产色视频综合| 亚洲伊人久久精品综合| 精品久久蜜臀av无| 亚洲精品久久午夜乱码| 国产精品熟女久久久久浪| 建设人人有责人人尽责人人享有的| 欧美+亚洲+日韩+国产| 日本91视频免费播放| 高清视频免费观看一区二区| 免费高清在线观看日韩| 精品免费久久久久久久清纯 | 日本a在线网址| 久久精品成人免费网站| 国产一区二区在线观看av| svipshipincom国产片| 多毛熟女@视频| 午夜福利影视在线免费观看| a级毛片黄视频| 亚洲国产欧美日韩在线播放| 精品国产一区二区三区四区第35| 亚洲中文av在线| 欧美精品一区二区大全| e午夜精品久久久久久久| 人人妻人人添人人爽欧美一区卜| 亚洲自偷自拍图片 自拍| 亚洲欧美精品综合一区二区三区| 国产深夜福利视频在线观看| 国产午夜精品一二区理论片| www.999成人在线观看| 亚洲av日韩在线播放| 免费在线观看黄色视频的| 亚洲精品国产一区二区精华液| 国产免费一区二区三区四区乱码| 国产日韩欧美视频二区| 久久国产精品男人的天堂亚洲| 久久久久精品人妻al黑| 日韩av不卡免费在线播放| 国产成人免费观看mmmm| 9色porny在线观看| 99国产精品一区二区三区| 一级黄片播放器| 亚洲少妇的诱惑av| 亚洲人成网站在线观看播放| 精品免费久久久久久久清纯 | 丁香六月欧美| 欧美精品啪啪一区二区三区 | 国产在视频线精品| 免费在线观看视频国产中文字幕亚洲 | 日本av免费视频播放| 国产97色在线日韩免费| 午夜老司机福利片| 国产精品欧美亚洲77777| 99久久精品国产亚洲精品| 汤姆久久久久久久影院中文字幕| 老司机午夜十八禁免费视频| 亚洲精品久久午夜乱码| 日韩 亚洲 欧美在线| av天堂在线播放| 啦啦啦 在线观看视频| 亚洲精品久久久久久婷婷小说| 国产精品 欧美亚洲| 美国免费a级毛片| 纯流量卡能插随身wifi吗| 人妻 亚洲 视频| 国产熟女欧美一区二区| 十八禁人妻一区二区| 大码成人一级视频| 日日摸夜夜添夜夜爱| 欧美 亚洲 国产 日韩一| 国产一级毛片在线| 日韩免费高清中文字幕av| 伦理电影免费视频| 久久久精品区二区三区| 国产av一区二区精品久久| 免费观看av网站的网址| 国产精品秋霞免费鲁丝片| 丁香六月欧美| 天堂俺去俺来也www色官网| 成人午夜精彩视频在线观看| 亚洲国产精品一区二区三区在线| 成年人黄色毛片网站| 国产精品 欧美亚洲| 夜夜骑夜夜射夜夜干| 亚洲精品国产av蜜桃| 亚洲久久久国产精品| 国产黄色免费在线视频| 大片电影免费在线观看免费| 一本久久精品| 嫩草影视91久久| 国产成人精品在线电影| 大型av网站在线播放| 日韩一卡2卡3卡4卡2021年| 日日摸夜夜添夜夜爱| 日韩中文字幕欧美一区二区 | 免费在线观看黄色视频的| 国产精品免费大片| 人人妻人人澡人人爽人人夜夜| 999精品在线视频| 最近最新中文字幕大全免费视频 | 1024香蕉在线观看| 午夜久久久在线观看| 永久免费av网站大全| 欧美av亚洲av综合av国产av| 男女边摸边吃奶| 久久毛片免费看一区二区三区| 亚洲欧美日韩高清在线视频 | 欧美97在线视频| 亚洲欧洲精品一区二区精品久久久| 亚洲中文日韩欧美视频| 久久精品人人爽人人爽视色| 亚洲,一卡二卡三卡| av福利片在线| 一区二区日韩欧美中文字幕| 50天的宝宝边吃奶边哭怎么回事| 欧美精品一区二区免费开放| 国产野战对白在线观看| 91精品国产国语对白视频| 日本午夜av视频| 国产精品九九99| 婷婷成人精品国产| 啦啦啦视频在线资源免费观看| 18禁黄网站禁片午夜丰满| 国产午夜精品一二区理论片| 亚洲欧美清纯卡通| 99久久99久久久精品蜜桃| 国产精品国产三级专区第一集| 久久精品人人爽人人爽视色| 一个人免费看片子| 亚洲精品一区蜜桃| avwww免费| 亚洲国产精品999| 久久人妻熟女aⅴ| 在线观看免费日韩欧美大片| 免费在线观看黄色视频的| 国产成人一区二区三区免费视频网站 | 国产精品一区二区在线观看99| 日本a在线网址| 亚洲黑人精品在线| 亚洲国产欧美一区二区综合| 黄频高清免费视频| 欧美日韩综合久久久久久| 丝袜人妻中文字幕| 久久久亚洲精品成人影院| 亚洲一区二区三区欧美精品| 欧美精品一区二区大全| 两人在一起打扑克的视频| av又黄又爽大尺度在线免费看| 777米奇影视久久| 美女扒开内裤让男人捅视频| 久久久久精品国产欧美久久久 | 夫妻性生交免费视频一级片| 人人妻人人爽人人添夜夜欢视频| 国产伦理片在线播放av一区| 宅男免费午夜| videosex国产| 国产av一区二区精品久久| 亚洲午夜精品一区,二区,三区| 建设人人有责人人尽责人人享有的| 欧美日韩黄片免| 激情视频va一区二区三区| av电影中文网址| 国产亚洲欧美精品永久| 成人手机av| 狠狠婷婷综合久久久久久88av| 久久久久久久久久久久大奶| 香蕉丝袜av| 男的添女的下面高潮视频| 欧美成狂野欧美在线观看| 亚洲天堂av无毛| 亚洲人成电影观看| 国产成人91sexporn| 18禁黄网站禁片午夜丰满| 亚洲五月婷婷丁香| 黄色片一级片一级黄色片| 色精品久久人妻99蜜桃| 老司机午夜十八禁免费视频| 91麻豆av在线| 亚洲欧洲日产国产| 99国产精品99久久久久| 九色亚洲精品在线播放| 天天躁夜夜躁狠狠躁躁| av欧美777| 大香蕉久久网| 久久天堂一区二区三区四区| 亚洲中文字幕日韩| 国产xxxxx性猛交| 欧美老熟妇乱子伦牲交| 一边摸一边抽搐一进一出视频| 欧美日韩av久久| a级毛片黄视频| 又紧又爽又黄一区二区| 免费在线观看视频国产中文字幕亚洲 | 97精品久久久久久久久久精品| 国产精品久久久久成人av| 欧美中文综合在线视频| 黄色毛片三级朝国网站| 国产成人a∨麻豆精品| 午夜久久久在线观看| av网站免费在线观看视频| 亚洲av在线观看美女高潮| 两个人免费观看高清视频| 国产野战对白在线观看| 日本欧美视频一区| 免费在线观看影片大全网站 | 97在线人人人人妻| 80岁老熟妇乱子伦牲交| 99久久人妻综合| 精品国产超薄肉色丝袜足j| 国产真人三级小视频在线观看| 亚洲成人免费av在线播放| 免费在线观看完整版高清| h视频一区二区三区| 成年av动漫网址| 午夜免费成人在线视频| 亚洲伊人色综图| 久久久久久久国产电影| 人妻人人澡人人爽人人| 激情五月婷婷亚洲| 欧美亚洲 丝袜 人妻 在线| 国产福利在线免费观看视频| 久久性视频一级片| 国产极品粉嫩免费观看在线| 久久精品久久久久久噜噜老黄| 各种免费的搞黄视频| 国产人伦9x9x在线观看| 黑人欧美特级aaaaaa片| 啦啦啦在线观看免费高清www| 国产97色在线日韩免费| 国产99久久九九免费精品| 欧美精品啪啪一区二区三区 | 菩萨蛮人人尽说江南好唐韦庄| av有码第一页| 后天国语完整版免费观看| 亚洲五月婷婷丁香| 一级片免费观看大全| 国产精品成人在线| 大型av网站在线播放| av视频免费观看在线观看| 老司机影院毛片| 日韩人妻精品一区2区三区| 亚洲久久久国产精品| 少妇裸体淫交视频免费看高清 | 色94色欧美一区二区| 国产色视频综合| 国产精品 欧美亚洲| 欧美日韩视频高清一区二区三区二| 老司机深夜福利视频在线观看 | 在线观看国产h片| 91九色精品人成在线观看| 天天躁日日躁夜夜躁夜夜| 女性被躁到高潮视频| 国产成人欧美| 亚洲av美国av| 大话2 男鬼变身卡| 亚洲国产av新网站| 啦啦啦中文免费视频观看日本| 免费观看a级毛片全部| bbb黄色大片| 亚洲精品久久午夜乱码| 狠狠精品人妻久久久久久综合| 国产午夜精品一二区理论片| 久久久久精品人妻al黑| 亚洲国产精品一区二区三区在线| 久久久亚洲精品成人影院| 亚洲精品国产色婷婷电影| 亚洲五月色婷婷综合| 丝袜喷水一区| 国产av精品麻豆| 麻豆av在线久日| 国产成人啪精品午夜网站| 日韩免费高清中文字幕av| 精品亚洲成a人片在线观看| 欧美老熟妇乱子伦牲交| 欧美日韩成人在线一区二区| 色播在线永久视频| 嫩草影视91久久| 欧美日韩综合久久久久久| 一区在线观看完整版| 色播在线永久视频| 天天躁夜夜躁狠狠躁躁| 日本91视频免费播放| 叶爱在线成人免费视频播放| 黄色片一级片一级黄色片| 精品高清国产在线一区| 一本大道久久a久久精品| 日韩 亚洲 欧美在线| www.自偷自拍.com| 欧美精品一区二区大全| 18禁裸乳无遮挡动漫免费视频| 国产成人影院久久av| 十分钟在线观看高清视频www| 国产成人av教育| 久久中文字幕一级| 嫩草影视91久久| 国产欧美亚洲国产| 制服诱惑二区| 人人妻人人添人人爽欧美一区卜| 国产男女内射视频| 亚洲欧美清纯卡通| 99国产精品免费福利视频| 在线观看免费午夜福利视频| 欧美日韩av久久| 天天影视国产精品| 大话2 男鬼变身卡| av一本久久久久| 国产精品亚洲av一区麻豆| 久久综合国产亚洲精品| 久久久久久久大尺度免费视频| 国产在线免费精品| 如日韩欧美国产精品一区二区三区| 国产极品粉嫩免费观看在线| 一本综合久久免费| 无遮挡黄片免费观看| 国产一卡二卡三卡精品| 亚洲色图 男人天堂 中文字幕| 亚洲国产精品一区三区| 丝袜喷水一区| 亚洲国产精品一区二区三区在线| www.精华液| 国产精品久久久久久精品电影小说| 亚洲成av片中文字幕在线观看| 日韩伦理黄色片| 国产在视频线精品| 国产亚洲精品第一综合不卡| 成年美女黄网站色视频大全免费| 一本色道久久久久久精品综合| 热99国产精品久久久久久7| 欧美+亚洲+日韩+国产| 啦啦啦 在线观看视频| 国产精品秋霞免费鲁丝片| 欧美 亚洲 国产 日韩一| 在线观看人妻少妇| 成人黄色视频免费在线看| 色播在线永久视频| 乱人伦中国视频| 国产熟女欧美一区二区| 亚洲精品久久午夜乱码| 精品亚洲乱码少妇综合久久| 永久免费av网站大全| 亚洲国产中文字幕在线视频| 亚洲中文字幕日韩| 90打野战视频偷拍视频| 国产成人免费观看mmmm| 啦啦啦中文免费视频观看日本| 一区福利在线观看| 日本欧美视频一区| 深夜精品福利| 在线精品无人区一区二区三| 欧美国产精品一级二级三级| 国产在线视频一区二区| 乱人伦中国视频| 亚洲成人手机| 日日爽夜夜爽网站| 日韩,欧美,国产一区二区三区| 三上悠亚av全集在线观看| 美女脱内裤让男人舔精品视频| 亚洲综合色网址| 美女福利国产在线| 免费黄频网站在线观看国产| 在线观看免费日韩欧美大片| 国产免费视频播放在线视频| 久久国产精品男人的天堂亚洲| 大话2 男鬼变身卡| 久久久国产一区二区| 男人爽女人下面视频在线观看| 欧美黑人精品巨大| 天堂中文最新版在线下载| 国产亚洲午夜精品一区二区久久| 一本大道久久a久久精品| 国产伦理片在线播放av一区| av片东京热男人的天堂| 老司机影院成人| 熟女av电影| 波多野结衣一区麻豆| 国产在视频线精品| 日韩欧美一区视频在线观看| 日韩 欧美 亚洲 中文字幕| 国产精品偷伦视频观看了| 亚洲人成电影免费在线| 国产精品久久久久成人av| 人妻一区二区av| 黄色片一级片一级黄色片| 下体分泌物呈黄色| 深夜精品福利| 免费人妻精品一区二区三区视频| 一本久久精品| 亚洲av美国av| 国产97色在线日韩免费| 在现免费观看毛片| 国产在线视频一区二区| 精品一区二区三区四区五区乱码 | 色精品久久人妻99蜜桃| 日本av手机在线免费观看| 又紧又爽又黄一区二区| 亚洲国产日韩一区二区| 自拍欧美九色日韩亚洲蝌蚪91| 不卡av一区二区三区| 欧美日韩黄片免| 最新的欧美精品一区二区| 高清不卡的av网站| 在线观看www视频免费| 大陆偷拍与自拍| 在线观看免费高清a一片| 丝袜喷水一区| 香蕉国产在线看| 国产成人a∨麻豆精品| 男人添女人高潮全过程视频| 中文欧美无线码| 国产免费一区二区三区四区乱码| 少妇精品久久久久久久| 成人手机av| 又黄又粗又硬又大视频| 国产欧美日韩精品亚洲av| 国产成人av教育| 国产爽快片一区二区三区| 中国美女看黄片| 亚洲av电影在线观看一区二区三区| 纯流量卡能插随身wifi吗| 狠狠婷婷综合久久久久久88av| 男女高潮啪啪啪动态图| 亚洲 国产 在线| 精品一区二区三区四区五区乱码 | 超碰成人久久| 新久久久久国产一级毛片| 国产欧美亚洲国产| 大码成人一级视频| 久久久久久久国产电影| 成人免费观看视频高清| 男女国产视频网站| 中文字幕制服av| 男女边摸边吃奶| 2018国产大陆天天弄谢| 久久亚洲精品不卡| 女人精品久久久久毛片| 亚洲欧美中文字幕日韩二区| 国产精品.久久久| 高潮久久久久久久久久久不卡| 热re99久久精品国产66热6| 午夜福利免费观看在线| 国产激情久久老熟女| 少妇粗大呻吟视频| 国产欧美日韩一区二区三 | 99久久综合免费| 午夜影院在线不卡| 视频区图区小说| 亚洲av片天天在线观看| 熟女av电影| 自拍欧美九色日韩亚洲蝌蚪91| 男人舔女人的私密视频| 这个男人来自地球电影免费观看| 黄色毛片三级朝国网站| 国产精品国产av在线观看| bbb黄色大片| 亚洲欧美日韩另类电影网站| 丝袜美腿诱惑在线| 日日夜夜操网爽| 免费看不卡的av| 在线观看国产h片| 欧美人与善性xxx| 校园人妻丝袜中文字幕| 亚洲男人天堂网一区| 日日爽夜夜爽网站| 久久精品人人爽人人爽视色| netflix在线观看网站| 亚洲人成77777在线视频| 91九色精品人成在线观看| 色网站视频免费| 又黄又粗又硬又大视频| 满18在线观看网站| 2021少妇久久久久久久久久久| 99国产精品一区二区蜜桃av | 人人妻,人人澡人人爽秒播 | 男女免费视频国产| 手机成人av网站| 黑丝袜美女国产一区| 欧美少妇被猛烈插入视频| 十八禁高潮呻吟视频| 狠狠婷婷综合久久久久久88av| 欧美久久黑人一区二区| 精品亚洲成a人片在线观看| 国产色视频综合| 色综合欧美亚洲国产小说| 国产xxxxx性猛交| 精品久久久久久电影网| 涩涩av久久男人的天堂| 在线观看免费日韩欧美大片| 中文乱码字字幕精品一区二区三区| 你懂的网址亚洲精品在线观看| 久久热在线av| 亚洲成人国产一区在线观看 | a级毛片黄视频| 欧美激情 高清一区二区三区| 一本—道久久a久久精品蜜桃钙片| 国产亚洲一区二区精品| 亚洲国产欧美一区二区综合| 久久午夜综合久久蜜桃| 韩国高清视频一区二区三区| 日韩人妻精品一区2区三区| 老司机亚洲免费影院| 成年人免费黄色播放视频| 操美女的视频在线观看| a 毛片基地| 免费观看av网站的网址| 午夜日韩欧美国产| 久久久欧美国产精品| 国产99久久九九免费精品| 99re6热这里在线精品视频| 欧美另类一区| 麻豆乱淫一区二区| 久久精品成人免费网站| e午夜精品久久久久久久| 亚洲伊人色综图| 在线观看免费日韩欧美大片| 天堂俺去俺来也www色官网| 看十八女毛片水多多多| 91老司机精品| 亚洲精品久久成人aⅴ小说| 久久久久久亚洲精品国产蜜桃av| 久热这里只有精品99| 亚洲国产精品一区三区| 亚洲中文av在线| 啦啦啦在线免费观看视频4| 亚洲欧美一区二区三区黑人| 久久国产精品影院| 欧美成狂野欧美在线观看| 午夜免费成人在线视频| 亚洲欧洲日产国产| 久久精品国产亚洲av涩爱| 一级毛片我不卡| 男女边摸边吃奶| 久久人人97超碰香蕉20202| 久久久国产一区二区|