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

    Taxifolin stability:In silico prediction and in vitro degradation with HPLC-UV/UPLC-ESI-MS monitoring

    2021-05-12 06:30:50FernnCristinStengerMourCrmemosSntosMhoFveroReisorferPulAngeliGriCoutoMurizioRiiVlirCehinelFilhoTigoBonominiLouisSnjoTniMriBellBresolin
    Journal of Pharmaceutical Analysis 2021年2期

    Fernn Cristin Stenger Mour,Crmem Lúi os Sntos Mho,Fvero Reisorfer Pul,Ang′eli Gri Couto,Murizio Rii,Vlir Cehinel-Filho,Tigo J.Bonomini,Louis P.Snjo,Tni Mri Bell′e Bresolin,*

    aPharmaceutical Sciences Graduate Program,School of Health Sciences,Course of Pharmacy,Universidade do Vale do Itajaí,Itajaí,SC,88302-202,Brazil

    bLaboratory of Research and Development of Drugs,Pharmaceutical Sciences Graduate Program,School of Health Sciences,Course of Pharmacy,Universidade Federal do Pampa,Uruguaiana,RS,97500-970,Brazil

    cDepartment of Pharmaceutical Science,Universit`a degli Studi di Perugia,Perugia,PG,06123,Italy

    dChemistry Department,Universidade Federal de Santa Catarina,Trindade,Florian′opolis,SC,88040-900,Brazil

    Keywords:

    Dihydroquercetin

    In silico stability prediction

    Forced degradation

    A B S T R A C T

    Taxifolin has a plethora of therapeutic activities and is currently isolated from the stem bark of the tree Larix gmelinni(Dahurian larch).It is a flavonoid of high commercial interest for its use in supplements or in antioxidant-rich functional foods.However,its poor stability and low bioavailability hinder the use of flavonoid in nutritional or pharmaceutical formulations.In this work,taxifolin isolated from the seeds of Mimusops balata,was evaluated by in silico stability prediction studies and in vitro forced degradation studies(acid and alkaline hydrolysis,oxidation,visible/UV radiation,dry/humid heating)monitored by high performance liquid chromatography with ultraviolet detection(HPLC-UV)and ultrahigh performance liquid chromatography-electrospray ionization-mass spectrometry(UPLC-ESI-MS).The in silico stability prediction studies indicated the most susceptible regions in the molecule to nucleophilic and electrophilic attacks,as well as the sites susceptible to oxidation.The in vitro forced degradation tests were in agreement with the in silico stability prediction,indicating that taxifolin is extremely unstable(class 1)under alkaline hydrolysis.In addition,taxifolin thermal degradation was increased by humidity.On the other hand,with respect to photosensitivity,taxifolin can be classified as class 4(stable).Moreover,the alkaline degradation products were characterized by UPLC-ESI-MS/MS as dimers of taxifolin.These results enabled an understanding of the intrinsic lability of taxifolin,contributing to the development of stability-indicating methods,and of appropriate drug release systems,with the aims of preserving its stability and improving its bioavailability.

    1.Introduction

    Taxifolin,also called dihydroquercetin,was first isolated from Pseudotigusa taxifolia[1].It belongs to the class of dihydroflavonoids,which are present as a yellow pigment in many edible plants.Lavitol?,an enriched extract of the stem bark of Larix gmelinni,contains about 90% taxifolin.This natural product has been marketed and used as a food ingredient in the USA since 2009[2].It has been recently authorized as a novel food ingredient in Europe[3].Taxifolin has been reported in more than 700 articles published in the literature.

    Taxifolin(3,5,7,3′,4′-pentahydroxyflavanone or dihydroquercetin)is a potent antioxidant,comparable to alpha-tocopherol,whose mechanism of action consists of lipid radical scavenging[4].Taxifolin presents several pharmacological properties,such as attenuating diabetic nephropathy,reducing sugar,uric acid and creatinine in human blood[5],decreasing the accumulation ofβ-amyloid and preventing memory deficits[6],protecting against oxidative stress[7],promoting osteogenic differentiation in human bone marrow mesenchymal stem cells[8],reducing blood viscosity and dilating the bloodvessels,andreducing arterialhypertension[9].Italsohelps prevent diabetic cardiomyopathy[10]and protects against alcoholic liver steatosis[11].Its in vivo gastroprotective effects were also demonstrated by our research group,with a similar effect to that of omeprazol,inhibiting 41% of the pump effect[12].Due to its relevant pharmacological activities,this potential phytodrug was previously incorporated into gastro adhesive microparticles for the treatmentof gastric disorders,by our research group[13].

    However,despite its high clinical application potential,there is no sufficient data on its stability,a key quality statement for a drug.Moreover,there are few stability studies carried out with natural products[14].Previous studies revealed that taxifolin polymerizes when submitted to electrolysis in neutral solutions(pH 7.0)[15].This substance is also considered to be highly slightly soluble in water,and its absolute bioavailability after oral administration of lipid solution was only 36%[16].Some authors argue that the results for absorption profile and the parameters of taxifolin vary quite considerably,suggesting that the bioavailability of the compound depends onthe source[3].This compound also appears tobe degraded by the intestinal microflora[17].Based on the above,the investigation of taxifolin’s stability behavior is an important step in developing a new active pharmaceutical ingredient(API)and in stimulating the development of a potential new drug release system.Moreover,the detection of potential degradation products requires stability-indicating methods[18]and the application of forced degradation tests of the API candidate,in both solid and solution forms.The stress condition should be carefully selected,aiming to generate potential degradation products,which are likely to be formed under realistic storage conditions[19].

    Therefore,taxifolin stability was investigated using in silico stability prediction studies and in vitro forced degradation studies,monitoring the formation of degradation products by high performance liquid chromatography with ultraviolet detection(HPLCUV)and ultra performance liquid chromatography coupled to electrospray ionization tandem high resolution mass spectrometry(UPLC-ESI-HRMS/MS).

    2.Experimental

    2.1.Chemicals and reagents

    Methanol and acetonitrile of liquid chromatography with ultraviolet detection(HPLC-UV)and liquid chromatography-mass spectrometry(LC-MS)grade were obtained from PanReac,Castellar del Vall`es,Barcelona,Spain.Phosphoric acid of analytical grade was purchased from Din^amica,Sao Paulo,Brazil.Class 1 water was obtained by ultrapurification(Direct-Q?,Merck KGaA,Darmstadt,Germany).All the solutions were filtered through a regenerated cellulose 0.45 and 0.22 μm membranefilter(Macherey-Nagel,Düren,Germany)prior to injection into HPLC and UPLC,respectively.Taxifolin used in this experiment was isolated from the seed of Mimusops balata(from Itajaí,Santa Catarina,Brazil)with 99.4% purity(estimated by HPLC),and characterized by1H and13C NMR,infrared spectroscopy,HPLC and MS(Supplementary material).

    2.2.In silico stability prediction studies

    Computational analyses were performed using Spartan 08 version 116.2?for Windows(Wave function,Inc.,Irvine,CA,USA)and all initial structures were constructed using fragments of atoms and structural fragments by the molecular editor.Geometric optimization was carried out using the Merck Molecular Strength Field(MMFF94)followed by the Austin Model.The structure of taxifolin(Fig.1A and B)was subjected to conformational analysis.The increment of the torsion angle was 30°in a range of 0-360°,using systematic analysis,by the functional density theory(DFT)method to B3LYP/6,311G*(d,f).The lowest energy conformer calculated was re-optimized using the same method.This structure was used to determine the number of electrons at natural atomic population analysis(NPA)using single point energy at the same level of theory of geometry optimization.For these data,Fukui function(FF)derivatives values,positive(f-j,Eq.(1))for electrophilic attack,negative(f+j,Eq.(2))for nucleophilic attack and f0j(Eq.(3))for radical attack we calculated as follows:

    Fig.1.Chemical structure of taxifolin(A);3D visualization in tube(B);Molecular electronic potential(MEP)(C),beyond the van der Waals isosurface 0.002 eV using Spartan for Windows 08.Color scheme:blue positive to red negative electrostatic potentials values(-40.000-75.000 kcal/mol).

    The possible regions of the molecule susceptible to degradation,under acid and alkaline conditions,were evaluated by the Fukui functions[20].These functions indicate the susceptibility of the electronic density to deform at a given position upon accepting or donating electrons[21].In this case,qjis the number of electrons(evaluated from NPA)at the j th atomic site in the neutral(N),anionic(N+1)or cationic(N-1)chemical species on the reference molecule.

    The dual descriptorΔf(r)of local reactivity,which allows us to obtain a preferred locus for nucleophilic attacks(Δf(r)>0)and a preferred electrophilic attack site(Δf(r)<0)in the system at point r,was calculated using the Eq.(4)[22,23].

    The main mechanism of auto-oxidation observed in the photolytic degradation reaction of the hydrogen atom abstraction energy(EHabstraction)was calculated for the bond of the hydrogen atom with the carbon atom,and bond dissociation energy(BDE)in the molecule,using the Eq.(5)[24,25].

    2.3.HPLC-UV methodology

    The chromatographic profile,purityand extractive content were analyzed by HPLC(Prominence LC-20AT LC,Shimadzu,Tokyo,Japan)equipped with a binary pump,a photo diode array detector(SPD-M20A,Shimadzu,Tokyo,Japan),an auto-sampler(SIL-20AHT,Shimadzu,Tokyo,Japan)in-line degasser(DGU-20A5,Shimadzu,Tokyo,Japan),a column oven(CTO-10AS V,Shimadzu,Tokyo,Japan),a 20μL stainless steel loop and the software Class VP(version 6.14).

    The analyses were carried out with a column Kinetex?F5(Phenomenex,Torrance,California,USA;150 mm×4.6 mm,2.6 μm)conditioned at 30°C,with detection at 288 nm.The best gradient elution(0.6 mL/min)was acidified water pH 3.0 with 0.1%(V/V)phosphoric acid(A),acetonitrile(B)and methanol(C)according to the following gradient A:B:C(V/V/V):85:10:5 to 80:15:5(0-5 min);80:15:5 to 70:25:5(5-15 min);70:25:5 to 50:45:5(15-25 min);and 50:45:5 to 85:10:5(25-30 min),maintaining these conditions until 35 min.

    The HPLC-UV methodology was validated for linearity,limits of detection(LOD)and quantification(LOQ),robustness and selectivity.To access the linearity,taxifolin was dissolved in methanol in the range of 5.0-100.0μg/mL,in triplicate,followed by linear regression studies,using Excel software.The LOQ and LOD were determined based on the standard deviation of the response and the slope[26].Selectivity of the method for taxifolin was provenfirst by accessing the resolution(R)between the taxifolin peak and the nearest peak(impurity)in the original sample.Secondly,the methodology was applied to quantify taxifolin among the degradation products after the stress tests,measuring the purity peak of taxifolin through the photodiode array(PDA)detector.To check the robustness of the method,the flow(±0.06 mL/min),temperature(± 1°C)and storage time of the sample solution(3,8,12 and 24 h)were varied measuring the relative standard deviation(RSD)of the analyte content.The taxifolin content was calculated using the regression linear equation.

    2.4.Forced degradation studies

    Firstly,an aliquot of taxifolin(1 mg)was dissolved in 1 mL of methanol.Then,the volume was completed,to 5 mL,with the respective solvent of each forced degradation test.Hydrolytic stress studies were performed in acid medium(1 M HCl)and the mixture was maintained for 15,30,60,240 min and 24 h.Taxifolin was also treated with an alkaline condition,at 1mM NaOH,during 15,30,60,240 min and 24 h.Oxidative tests were carried out in 30% H2O2for 15,30,60,240 min and 24 h.All the stress experiments were carried out at 25°C.The final concentration of taxifolin in the respective degradation media was 200μg/mL.Acid and alkaline samples were neutralized with NaOH and HCl,respectively,and prior to HPLC analysis,the sample solutions were diluted with methanol to 100μg/mL.This concentration was chosen in order to improve the appearance of the degradation products in the HPLCUV chromatograms.For LC-MS analysis,taxifolin was submitted to 0.01 M NaOH degradation at room temperature,followed by immediate neutralizationwith equimolar HCl.The taxifolinpowder was stored in a photostability chamber(model 91423,ZEM,Minas Gerais,Brazil)under combined visible light(1.2 and 2.4 million lx·h)and ultraviolet-A(UVA)(200 and 400 Wh/m2irradiation).The negative control was protected with aluminium foil coating.The light intensity was monitored by a Luximeter(model MLM-1011,Minipa,S~ao Paulo,Brazil).The samples were also stored in an oven at 40°C and at 40°C/75% relative humidity(RH)(into a desiccator with saturated solution of NaCl)for 30 days[27].After exposure to light,heat and humidity,sample solutions of 100μg/mL in methanol were analyzed by HPLC.All stressed samples were assayed by comparison with non-degraded reference standards in the same concentration,in methanol,to calculate the percentage of degradation.

    2.5.UPLC-ESI-HRMS/MS analysis

    LC-MS analyses were performed on an Acquity UPLC system class H(Waters,St.Milford,MA,USA)composed of a PDA detector,sample manager and a quaternary solvent manager as well as an Acquity UPLC BEH C18column(St.Milford,MA,USA;130 ?,1.0 mm × 50 mm,particle size 1.7μm).Temperatures of 40°C and 20°C were set for the column and sample tray,respectively.A volume of 3μL was injected for each sample,and the separation was obtained in a gradient condition(0-5 min:95%A(water:formic acid,99.9:0.1,V/V)and 5%B(acetonitrile);5-8 min:20%A;8-10 min:50%A;10-13 min:55%A;13-16 min:90%A;and 16-20 min:95%A),at a flow rate of 0.3 mL/min.

    Mass detection was conducted on a Xevo G2-S QTof(Waters,St.Milford,MA,USA)with an electrospray probe operating in negative ionization mode;nebulizer gas:nitrogen,cone gas flow 60 L/h;desolvation gas flow 900 L/h,sampling cone 40 V,source offset 80 V;collision gas,argon.Lockspray reference sample was leucine encephalin with reference mass at m/z 554.2615(ESI-).Temperatures of 300°C and 120°C were used for the desolvation and for the cone,respectively,while the capillary voltage was 3 kV.The collision energy was 30 eV.Data were acquired in a range of 100-1500 Da,at a scan time of 1.0 s during 20 min,and were processed with Mass Lynx V4.1(Waters,St.Milford,MA,USA).

    3.Results and discussion

    The isolated and purified taxifolin was characterized by Fourier transform infrared spectroscopy spectra(Fig.S1),mass spectrometry data(Fig.S2)and1H NMR spectral data(Fig.S3)and then,subjected to the stress tests.

    3.1.LC-UV method validation

    The method was linear in the range of 5.0-100.0μg/mL(y=113027x-67233,r2=0.9989).The LOD and LOQ were 4.6 and 15.9μg/mL,respectively.The chromatographic conditions showed selectivity with high resolution(R=4.2)between taxifolin and the nearest peak(impurity identified as a taxifolin isomer),showing a high purity peak index of 0.999973,through the PDA detector.The method was robust for changes in temperature(RSD of 0.25%),mobile phase flow(RSD of 0.47%),and storage time(RSD of 0.47%),with P>0.05,showing no statistical difference in the taxifolin assay in the small and deliberate modifications of the methodology,compared to the standard condition.

    3.2.In silico stability prediction study

    Taxifolin stability was first investigated by in silico prediction studies and after by in vitro forced degradation tests;the samples were monitored by HPLC-UV and UPLC-ESI-HRMS/MS.

    In silico stability prediction studies were performed to further understand the chemical reactivity by predicting the most likely position of structural hydrolysis.The chemical structure(Fig.1A)in the tube model(Fig.1B)and the map of electrostatic potential charges(MEP)(Fig.1C)of taxifolin are shown in Fig.1.The MEP plot calculating on the van der Waals surface and focusing on the negative isopotential surfaces might be used to describe the occurrence of the electronic conjugation in the molecule.The most negative potential regions were demonstrated in red to oxygen atoms.The benzene ring showed no coplanarity with the chromen-4-ene ring,but the oxygen atoms from chromen-4-ene and hydroxyl(O10)moieties and benzene attached to position C2 showed a mesomeric effect.The center of the structure presented extended distribution of negative charges,which became more positive as the potential increased to hydrogen attached to oxygen atoms from hydroxyl moieties.

    As shown in Table 1,the major values of the Fukui function(f(r))are those with higher reactivity[28],whose sites may be related to hydrolysis reactions during stress degradation studies.Table 1 also shows that positions C2,C4,and C7 were more susceptible to undergoing nucleophilic attack.According to the dual descriptorΔf(r)results,the following order was observed:C4>C7>O10>C8,indicating preferably the chromen-4-one ring and C2′(Fig.1)in benzene ring sites as the most reactive(Δf(r)>0)into the system at point r.All molecular regions,especially C4 and C7,were shown to be the most probable reactive sites to alkaline or acid hydrolysis(Table 1,Fig.1).

    Table 1 Values of the neutral,positive and negative populational analysis(NPA),electrophilic f-and nucleophilic f+condensed Fukui functions andΔf(r)of the atoms of the taxifolin molecule calculated with the DFT/B3LYP and the 6.311G*(d,f)basis set considering Eqs.(1)-(3).

    The Fukui function f0and BDE,used to estimate the hydrogen abstraction energies,indicate the auto-oxidation of taxifolin,as shown in Table 2.Fukui function showed a high value for the hydrogen atom of C2′in the chemical structure of taxifolin(Fig.1),indicating the main region of susceptibility in the oxidation process.However,this hydrogen atom is stabilized by the resonance effect,and may not easily receive an electron in auto-oxidation,which suggests that the mechanism of electron transfer is not related to the chemical structure oxidation of taxifolin.

    BDE is a measure of the bond strength in a chemical bond,and may be defined as the standard energy(or even the enthalpy)change when a bond is broken by a reaction[29].BDE is an indicator of primary site(s)of the auto-oxidation of organic compounds and drugs[25].Low BDE values of hydrogen atoms from the taxifolin structure(Table 2)are more liable to be abstracted from the auto-oxidation reaction.The energies of hydrogen atoms attached to C2 and C3 from the chromen-4-ene ring and O7′and O8′from hydroxylmoieties attached tothe para-position of the benzene ring(Fig.1)reflect the high oxidative susceptibility that characterizes the instability of this region in the oxidative environment.On the other hand,the hydrogen atom from C8 showed a high BDE calculated value,indicating that this molecular position moiety is more stable in the studied conditions.The mesomeric effect observed in this structural region mayhinder the abstraction of this hydrogen.This procedure can be used for auto-oxidation and also for predicting photolytic degradation[30].

    Table 2 Fukui function values for radical attack and bond dissociation energies of taxifolin hydrogen.

    The Fukui function was previously used to estimate the antioxidant mechanism of taxifolin,compared to quercetin[31].However,the authors studied only the hydrogen atoms attached to oxygen atoms,and showed that in oxidant medium,taxifolin exhibited more susceptibility to hydrogen loss at O3′and O4′,which corresponds to the O7′and O8′position denominated in the present work(Fig.1A).Thus,these findings are in agreement with the results demonstrated here,considering the hydroxyl moieties.

    On the other hand,forced degradation studies are recommended as part of stability studies for new drugs because the molecule is exposed to extreme conditions during its manufacturing,storage,and administration.Thus,to simulate the appearance of potential degradation products,stress tests are useful and allow us to establish the degradation pathways and the intrinsic stability of the molecule,and to validate a stabilityindicating method.The Guideline ICH Q1A(R2)[27]states that the molecule should be exposed to high temperatures,humidity,oxidation,photolysis,and susceptibility to hydrolysis across a range of pH values,in solution or suspension.

    3.3.Forced degradation tests

    Taxifolin was submitted to forced degradation studies,selecting the condition which provides about 10%-20% of degradation[27].

    The taxifolin degradations,for each stress condition,were as follows:20.2%(1 M HCl,30 min),16.3%(1 mM NaOH,15 min),11.7%[30% H2O2,24 h],9.8%(dry heat,at 40°C,30 days),23.1%(humid heat,at 40°C and 75% RH,30 days)and 9.0%(photolysis at 2.4 million lx·h and at 400 Wh/m2of visible and UVA radiation,respectively).Taxifolin was shown to be extremely unstable(class 1)[32]under alkaline hydrolysis.Taxifolin degradation was increased by humidity,compared with dry heat,probably due to the hydrolysis process in the former condition.The stability of a semi-solid formulation containing taxifolinwas also evaluated[33].After 12 weeks,at 40°C,only 3% of the initial amount was found,which is in agreement with the results of the present study,reinforcing the thermolability of this compound,especially in the presence of humidity.Under oxidative conditions,taxifolin could be classified as class 4(stable),with a relative low photosensitivity[32].Thus,according to the above results,taxifolin needs to be protected against these stress conditions with a proper formulation,to avoid intestinal delivery in particular.

    The chromatographic profile of taxifolin showed changes after being subjected to certain degradation conditions,with the appearance of supplementary peaks(Fig.2).Taxifolin used for the stress studies showed low levels of impurity(0.6%)(Fig.2A).In this original sample,four impurity peaks eluted after the major taxifolin peak(Fig.2A1).The acidic hydrolysis tests with 1 M HCl showed supplementary peaks with earlier elution,decreasing from the major peak,as well as the disappearance of impurity peaks,which probably suffered acid hydrolysis(Fig.2B).Alkaline tests were performed starting at 1 M NaOH,followed by 0.1 M NaOH and in both conditions;taxifolin was completely decomposed,as evidenced by the UV profile of the peak with a retention time of 18 min.After decreasing the concentration to 0.01 M NaOH,75% taxifolin degradation was observed,and as the contact time increased,this degradation was completed(Fig.S4).

    Thus,the chosen conditions for LC-UV analysis,for alkaline degradation,were 1 mM NaOH for 15 min(Fig.2C).Oxidation with 30% H2O2resulted in a decrease in taxifolin and a change of chromatographic pro file,with many earlier eluted peaks(Fig.2D).Samples were also submitted to dry heat(40°C)(Fig.2E)and humid heat(40°C/75% RH)for 30 days(Fig.2F),without important chromatographic pro file changes.Finally,taxifolin was submitted to combined daylight fluorescent lamp and UVA(Fig.2G),showing photostability.

    UPLC-ESI-HRMS-MS analysis was used to characterize the degradation products of taxifolin present in samples obtained from various conditions of degradation including alkaline medium(Fig.3).

    The samples were more sensitive to negative ionization mode using the base peak ion as the mode of acquisition of the chromatograms[34].The first chromatogram obtained from the taxifolin sample(Fig.3A)displayed seven peaks,with the most intense one at 6.47 min,with a shoulder at 6.62 min,both with the base peak ion m/z 303.0515 corresponding to[C15H12O7-H]-calculated for m/z 303.0505(Δ3.37 ppm).The second peak at 6.62 min exhibited the same molecular ion,suggesting a taxifolin diastereomer.This suggestion was supported by fragment ions observed intheMS/MSspectraofbothcompoundsatm/z285.04[C15H12O7-H2O]-and 125.02(corresponding to phloroglucinol:ring A).The peak at 7.54 min with the pseudo-molecular ion m/z 301.0367[C15H10O7-H]-(Δ6.22 ppm)was assigned as quercetin based on the fragment ion observed at m/z 151.01 formed by the Retro-Diels-Alder rearrangement of ring C.This molecule is an oxidized form of taxifolin.The other peaks at 8.79,10.73,11.17,and 14.77 min did not give fragment ions that would enable their characterization.However,the last metabolite at 14.77 min(m/z 281.2484[C18H34O2-H]-Δ1.23 ppm)was assigned as oleic acid derivative.

    In addition to the metabolites already present in the taxifolin sample,acid hydrolysis produced two new metabolites.The first was observed at 0.46 min(molecular ion m/z 197.79)and the second at 7.10 min(molecular ion m/z 723.51).The latter peak also appeared in the alkaline,oxidative and dry thermal degradation of taxifolin(Table 3).Diagnostic analysis of the MS/MS data of this peak was not conclusive.

    Table 3 Degradation products of taxifolin after exposition to different stress conditions,monitored by UPLC-ESI-MS.

    Sample oxidative test showed in its LC-MS data features similar to those of acid degradation,with the exception of an additional peak at 0.49 min with m/z 162.88.The dry thermal test afforded a sample for which LC-MS analysis gave two peaks,at 12.46 min and 13.34 min,with m/z 325.18 and m/z 339.20,respectively.The structure of m/z 325.18 could not be assigned;whereas based on a literature search,m/z 339.20[C22H28O3-H]-(Δ4.66 ppm)is related to a steroid.These metabolites are probably contaminants of the starting material,and appear during the degradation process,which decreases the concentration of the main compounds.

    The sample from humid thermal degradation afforded in its LCMS data peaks at 12.57 and 14.47 minwith the same mass value m/z 439.25.Because of the lack of fragment ions on the MS/MS spectrum,the structure of this metabolite could not be assigned.However,a literature search revealed that m/z 439.25 is related to a steroid.

    In the UV photolysis sample,a peak was observed at 8.82 min with m/z 221.1180[C13H18O3-H]-(Δ1.04 ppm),like in visible photolysis.A search using this elemental composition led to a structure related to norsesquiterpene.An additional peak was also observed at 12.60 min(m/z 421.2301)with no relation to taxifolin.Some peaks that appear in the water conditions are probably degradation products of quercetin,because in acid,alkaline,oxidative,and dry thermic degradation,the peak at 7.54 min,for which the molecular ion is m/z 301.04,disappeared.In addition,these samples,when analyzed by HPLC,did not show a significant decrease in the taxifolin peak(Fig.2).

    The purest taxifolin sample(99%)was submitted to alkaline hydrolysis,Specifically at 0.01 M NaOH for LC-MS analysis,and neutralized to afford,in its LC-MS analysis,a peak at 7.25 minflanked with m/z 603.0787[C30H20O14-H]-(calc.m/z 603.0775,Δ 2.02 ppm)(Fig.3B),suggesting the dimerization of taxifolin.The same chromatogram also showed the presence of a substance at 4.86 min with m/z 319.0456[C15H12O8-H]-(calc.m/z 319.0454,Δ 0.65 ppm),different from the taxifolin chemical composition(m/z 303[C15H12O7-H]-)by 16 Da.This substance may be a product formed from the oxidation of taxifolin by oxygen(O2)(Fig.4).

    As ring B is more susceptible to oxidation than ring A[35],it is suggested that alkaline degradation involves the auto-oxidation promoted by the radical peroxide formation.After losing a molecule of water and forming a ketone group,the aromatic system of B ring was rebuilt,resulting in the ion m/z 319.0456(Fig.4A).UPLCESI-MS2data of this product showed fragment ions at m/z 193,195,153 and 163(Fig.S5)supported this hypothesis and indicated that oxidation occurred on B ring of the flavanone.

    The second peak with m/z 337.0568[C15H14O9-H]-(calc.m/z 337.0560,Δ2.50 ppm),observed at 6.18 min(Fig.3B),differed from m/z 319.0456 by 18 Da,corresponding to an H2O molecule(Fig.4A).It was suggestedthat this compound was a resultof the oxidation of ring B of taxifolin and the opening of ring C.The MS/MS spectra of this second product showed a base peak ion of m/z 125 and a radical anion m/z 152.The fragmentation proposal in Fig.S5 revealed that m/z 125 could be related to A and B rings,whereas m/z 163 was produced by the sigma bond cleavage at theα-carbon atom.The latter fragmentation also indicated that ring C was opened.

    The third and fourth peaks,with m/z 605.0964[C30H22O14-H]-(calc.m/z 605.0931,Δ5.40 ppm),at 6.99 and 7.43,min,respectively,were indicative of dimeric products of taxifolin.It was suggested that these dimers were formed by 2 mol of taxifolin in the alkaline medium,which reacted with one mol of O2to produce 2 mol of taxifolin anion radical and H2O2.Furthermore,radical delocalization led to the formationof two dimeric flavanones(Fig.4B).Similar oxidative dimerization was electrochemically produced by Chernikov et al.[15].The MS/MS spectrum of both compounds(m/z 605.0964)showed similar fragments at m/z 393,m/z 259 and m/z 217.The ion m/z 393 was obtained by the elimination of one chromenone moiety and H2O.While the loss of two chromenone units produced m/z 217,the product ion m/z 259 was obtained from the opening of C ring and the loss of a chromenone unit(Fig.S6).

    Fig.2.Chromatograms of pristine taxifolin(A and A1)and of taxifolin degraded with 1 M HCl,30 min(B),1 mM NaOH,15 min(C),30% H2O2,24 h(D),dry heat,30 days(E),humid heat,30 days(F)and photolysis,2.4 million lx?h(G).

    Fig.3.LC-MS of pristine taxifolin(A)and the alkaline forced degradation(0.01 M NaOH immediately neutralized with equimolar HCl)of taxifolin(B).

    The product at m/z 603.0787 differs from at m/z 605.0964 by 2 amu,suggesting oxidation by O2(Fig.4C).In addition,the oxidation could produce an ortho-quinone dimeric flavanone,thus like a quercetin molecule,which has a similar structure to that of taxifolin,differing only in the double bound between C2 and C3.Quercitin also shows susceptibility to pH,temperature and storage conditions,generating the same ortho-quinone,which can interact with glutatione and may explain its high antioxidant effect[36].The latter,when analyzed by MS/MS,generated ions at m/z 409,391,381,193,and 177,demonstrating that oxidation occurred on the B ring but not on the C ring(Fig.S6).

    Fig.4.Hypothetical alkaline degradation mechanism of taxifolin with 0.01 M NaOH leading to m/z 319.0456 and m/z 337.0568(A);m/z 605.0964(B);m/z 603.0787(C).

    4.Conclusions

    Thus,in silico stability prediction analysis suggested that the taxifolin molecule is susceptible to nucleophilic attack in C2′,C4 and C7,as confirmed by the in vitro alkaline degradation,suggesting that the molecular ionwith m/z 603.0787 and m/z 605.0964 could be taxifolin dimers via nucleophilic attack in C2′for both molecules.The degradation products formed are also susceptible to nucleophilic attack.Other carbons are also previewed to be susceptible to nucleophilic attack,with less probability,but in vitro forced degradation studies proved that secondary degradation could occur.In the m/z 319 peak fragmentation,the nucleophilic attack occurred in C4 and O12(m/z 167),in C2,C3,and O12(m/z 153),in C2 and O12(m/z 193)and in C2 and O12(m/z 195),which corroborates the in silico stability prediction for susceptibility in C2,C4 and O12.In the m/z 337 fragmentations,the nucleophilic attack occurred in C2′,C4(m/z 125)and in C2,O12 and C4(m/z 125),in O1,O12(m/z 163),in C2 and C3(m/z 152),in agreement with in silico stability prediction for C2′,C2,C4,O12.Thus,the in silico stability prediction study should be predictive of the reactivitysusceptibility of taxifolin in the studied conditions,contributing to the understanding of taxifolin intrinsic lability and guiding the development of appropriate release systems for this phytodrug.

    Declaration of competing interest

    The authors declare that there are no conflicts of interest.

    Acknowledgments

    This work was partially supported by CAPES(PVE,Grant No.88887.116106/2016-00)(Coordena?~ao de Aperfei?oamento de Pessoal de Nível Superior),Brazil,which provided financial support in the form of a doctoral’s degree scholarship to Stenger,F.C.and financial support(Science Program Without Borders-Researcher Special Visitor-PVE),and CNPq(Conselho Nacional de Desenvolvimento Científico e Tecnolo′gico),Edital Universal(Grant No.88887.122964/2016-00).

    Appendix A.Supplementary data

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

    欧美人与善性xxx| 一级毛片 在线播放| 男人舔奶头视频| 久久久久久久久大av| 中国国产av一级| 春色校园在线视频观看| 久久久久久久大尺度免费视频| 成人欧美大片| av福利片在线观看| 少妇被粗大猛烈的视频| 欧美一区二区亚洲| 日韩 亚洲 欧美在线| 直男gayav资源| 久久久久久久久久成人| 国产精品成人在线| 久久ye,这里只有精品| 在线 av 中文字幕| 久久精品熟女亚洲av麻豆精品| 国产成人午夜福利电影在线观看| 成人亚洲精品av一区二区| 午夜福利高清视频| 人体艺术视频欧美日本| 亚洲激情五月婷婷啪啪| 国产爱豆传媒在线观看| 国产极品天堂在线| 1000部很黄的大片| 日韩av不卡免费在线播放| 亚洲成人精品中文字幕电影| 搡女人真爽免费视频火全软件| 免费不卡的大黄色大毛片视频在线观看| 激情五月婷婷亚洲| 国产免费一区二区三区四区乱码| 国产精品女同一区二区软件| 特大巨黑吊av在线直播| 一级毛片久久久久久久久女| 99视频精品全部免费 在线| 久久久久久久久久人人人人人人| 99视频精品全部免费 在线| 国产爱豆传媒在线观看| 亚洲av国产av综合av卡| 亚洲精品成人久久久久久| 免费观看av网站的网址| 国产成人精品久久久久久| 成人无遮挡网站| 亚洲最大成人av| 日韩三级伦理在线观看| 亚洲伊人久久精品综合| 日本猛色少妇xxxxx猛交久久| 日日撸夜夜添| 亚州av有码| 亚洲熟女精品中文字幕| 日本熟妇午夜| 一个人看视频在线观看www免费| 不卡视频在线观看欧美| 日韩欧美 国产精品| 国产黄色免费在线视频| 免费观看a级毛片全部| 久久人人爽人人爽人人片va| 综合色av麻豆| 国产成人精品婷婷| 黑人高潮一二区| 国产黄色免费在线视频| 菩萨蛮人人尽说江南好唐韦庄| 全区人妻精品视频| 插逼视频在线观看| 日韩av不卡免费在线播放| 亚洲精品乱码久久久久久按摩| 日产精品乱码卡一卡2卡三| 亚洲自拍偷在线| 亚洲av日韩在线播放| 大陆偷拍与自拍| 人体艺术视频欧美日本| 肉色欧美久久久久久久蜜桃 | 日本与韩国留学比较| 国产毛片在线视频| 日韩av免费高清视频| 国产伦在线观看视频一区| 少妇人妻 视频| 1000部很黄的大片| 日本-黄色视频高清免费观看| 国产精品久久久久久久电影| 国产午夜精品久久久久久一区二区三区| 亚洲激情五月婷婷啪啪| 亚洲av国产av综合av卡| 精品人妻偷拍中文字幕| 久久热精品热| 成人黄色视频免费在线看| 欧美bdsm另类| 久久99蜜桃精品久久| 哪个播放器可以免费观看大片| 国产高潮美女av| 亚洲精品乱码久久久久久按摩| 午夜福利在线在线| 少妇丰满av| 亚洲欧美成人综合另类久久久| 七月丁香在线播放| 国产视频首页在线观看| 国产亚洲最大av| 亚洲精品国产av成人精品| 熟妇人妻不卡中文字幕| 国产在线男女| 精品少妇黑人巨大在线播放| 久久久久久久久久久免费av| 在线天堂最新版资源| 中文字幕av成人在线电影| 搡女人真爽免费视频火全软件| 视频中文字幕在线观看| 99久久中文字幕三级久久日本| 中文精品一卡2卡3卡4更新| 日日撸夜夜添| 欧美日韩视频高清一区二区三区二| 在线免费十八禁| 韩国av在线不卡| 九九爱精品视频在线观看| 97超碰精品成人国产| 国产精品.久久久| 久久久久久久久久成人| 日韩欧美 国产精品| 亚洲内射少妇av| 水蜜桃什么品种好| 精品久久久久久久人妻蜜臀av| h日本视频在线播放| 国产精品一区二区在线观看99| 亚洲av日韩在线播放| 久久99蜜桃精品久久| 日本午夜av视频| 中国三级夫妇交换| 亚洲av一区综合| 一区二区三区乱码不卡18| 乱系列少妇在线播放| 成人鲁丝片一二三区免费| 国产高清三级在线| 最近中文字幕2019免费版| 国产精品福利在线免费观看| 国产 一区精品| 免费观看在线日韩| 在线观看av片永久免费下载| 欧美日韩在线观看h| 一级片'在线观看视频| av在线播放精品| 在线免费十八禁| 亚洲图色成人| 欧美一级a爱片免费观看看| 亚洲色图av天堂| 国产精品成人在线| 乱系列少妇在线播放| 高清日韩中文字幕在线| 久久久久精品久久久久真实原创| 欧美国产精品一级二级三级 | 国产一区二区亚洲精品在线观看| 国产在线一区二区三区精| 日韩一区二区三区影片| 国产精品一二三区在线看| 直男gayav资源| 看黄色毛片网站| 看十八女毛片水多多多| www.av在线官网国产| 国产男女超爽视频在线观看| 一个人看视频在线观看www免费| 色综合色国产| 欧美精品人与动牲交sv欧美| 街头女战士在线观看网站| 国产av不卡久久| 亚洲av免费高清在线观看| 久久精品国产亚洲av天美| 亚洲,欧美,日韩| 欧美日韩视频精品一区| 国产男女内射视频| 亚洲精品久久午夜乱码| 久久精品夜色国产| 欧美97在线视频| 亚洲av不卡在线观看| 纵有疾风起免费观看全集完整版| 久久午夜福利片| 好男人在线观看高清免费视频| 免费黄频网站在线观看国产| 天天一区二区日本电影三级| 天天躁夜夜躁狠狠久久av| 2022亚洲国产成人精品| 国产精品av视频在线免费观看| 黄片wwwwww| 99久久人妻综合| 国产高清有码在线观看视频| 女人十人毛片免费观看3o分钟| 色婷婷久久久亚洲欧美| 国产老妇伦熟女老妇高清| 看免费成人av毛片| 成人无遮挡网站| 最近手机中文字幕大全| 久久久久九九精品影院| 伊人久久国产一区二区| h日本视频在线播放| 精品国产三级普通话版| 国产成人a∨麻豆精品| 老师上课跳d突然被开到最大视频| 精品国产露脸久久av麻豆| 日韩人妻高清精品专区| 王馨瑶露胸无遮挡在线观看| 午夜视频国产福利| 最近最新中文字幕免费大全7| 最近2019中文字幕mv第一页| 三级男女做爰猛烈吃奶摸视频| 男女边吃奶边做爰视频| 97在线人人人人妻| 国产午夜精品久久久久久一区二区三区| 久久午夜福利片| 另类亚洲欧美激情| av免费在线看不卡| 蜜桃亚洲精品一区二区三区| 国产av码专区亚洲av| 嫩草影院精品99| 亚洲av在线观看美女高潮| 亚洲av欧美aⅴ国产| 久久精品久久精品一区二区三区| 亚洲人成网站在线观看播放| 国产精品国产三级国产av玫瑰| 99热这里只有是精品在线观看| 大陆偷拍与自拍| 中国国产av一级| videos熟女内射| 蜜臀久久99精品久久宅男| 毛片一级片免费看久久久久| 国产欧美日韩精品一区二区| 99久久中文字幕三级久久日本| 天堂俺去俺来也www色官网| 久久久久久久久久人人人人人人| 嫩草影院新地址| 99热这里只有是精品50| 国产真实伦视频高清在线观看| 亚洲婷婷狠狠爱综合网| 天天一区二区日本电影三级| 搡老乐熟女国产| 精品亚洲乱码少妇综合久久| 一区二区三区精品91| 久久韩国三级中文字幕| 亚洲在久久综合| 啦啦啦在线观看免费高清www| 久久精品久久久久久久性| 亚洲电影在线观看av| 99九九线精品视频在线观看视频| 日韩电影二区| 国产精品国产三级国产av玫瑰| 欧美3d第一页| 亚洲一区二区三区欧美精品 | 精品酒店卫生间| 国国产精品蜜臀av免费| 欧美一级a爱片免费观看看| 3wmmmm亚洲av在线观看| 亚洲精品第二区| 99精国产麻豆久久婷婷| 亚洲欧美日韩卡通动漫| 在线播放无遮挡| 最近最新中文字幕免费大全7| 亚洲成人一二三区av| 1000部很黄的大片| 男插女下体视频免费在线播放| 精品人妻一区二区三区麻豆| 久久精品国产自在天天线| 少妇人妻精品综合一区二区| 天堂中文最新版在线下载 | 别揉我奶头 嗯啊视频| 2022亚洲国产成人精品| 国产片特级美女逼逼视频| 成人高潮视频无遮挡免费网站| 亚洲国产欧美人成| 成人美女网站在线观看视频| 免费电影在线观看免费观看| 亚洲精品乱码久久久v下载方式| 欧美精品一区二区大全| 高清欧美精品videossex| 国产一区二区三区综合在线观看 | 国产淫片久久久久久久久| 国产欧美亚洲国产| 观看美女的网站| 秋霞伦理黄片| xxx大片免费视频| 99热这里只有精品一区| 五月伊人婷婷丁香| 亚洲精品成人av观看孕妇| 国产精品久久久久久久电影| 国产免费一区二区三区四区乱码| 十八禁网站网址无遮挡 | 99久久中文字幕三级久久日本| 国产精品伦人一区二区| 久久99蜜桃精品久久| 亚洲国产色片| 麻豆久久精品国产亚洲av| 国产成人免费无遮挡视频| 中国美白少妇内射xxxbb| av又黄又爽大尺度在线免费看| 蜜桃久久精品国产亚洲av| 汤姆久久久久久久影院中文字幕| av女优亚洲男人天堂| 欧美区成人在线视频| 啦啦啦中文免费视频观看日本| 成人午夜精彩视频在线观看| 男人舔奶头视频| 国产一区二区亚洲精品在线观看| 精品少妇久久久久久888优播| 热99国产精品久久久久久7| 成人午夜精彩视频在线观看| 精品人妻偷拍中文字幕| 97超视频在线观看视频| 久久6这里有精品| 晚上一个人看的免费电影| 国产探花极品一区二区| 最新中文字幕久久久久| 亚洲欧美清纯卡通| 可以在线观看毛片的网站| 一个人看的www免费观看视频| 亚州av有码| 51国产日韩欧美| www.av在线官网国产| 亚洲av电影在线观看一区二区三区 | 亚洲性久久影院| 中国三级夫妇交换| 男女那种视频在线观看| 99热网站在线观看| 亚洲国产精品国产精品| 视频区图区小说| 国产69精品久久久久777片| 国产乱人视频| 精品人妻熟女av久视频| 中文字幕av成人在线电影| 亚洲精品日本国产第一区| 一级毛片黄色毛片免费观看视频| 一区二区三区免费毛片| 波野结衣二区三区在线| 91精品伊人久久大香线蕉| 久久久久国产网址| 亚洲aⅴ乱码一区二区在线播放| 特级一级黄色大片| 观看美女的网站| 日日摸夜夜添夜夜添av毛片| 久久久精品欧美日韩精品| 人妻 亚洲 视频| 久久人人爽人人爽人人片va| 在线观看一区二区三区激情| 精品午夜福利在线看| 国产 一区 欧美 日韩| 国产伦精品一区二区三区视频9| 欧美 日韩 精品 国产| 欧美日韩在线观看h| 五月天丁香电影| 亚洲最大成人中文| 国产真实伦视频高清在线观看| 九草在线视频观看| 亚洲欧洲日产国产| 亚洲精品一区蜜桃| 日本-黄色视频高清免费观看| 日韩不卡一区二区三区视频在线| 色视频www国产| 永久网站在线| 噜噜噜噜噜久久久久久91| 国产综合懂色| 国产爱豆传媒在线观看| 一区二区三区免费毛片| 尾随美女入室| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产精品一区www在线观看| 亚洲精华国产精华液的使用体验| 国内揄拍国产精品人妻在线| 国产中年淑女户外野战色| 国产精品一区二区在线观看99| 色哟哟·www| 国产91av在线免费观看| 两个人的视频大全免费| 亚洲精品久久午夜乱码| 国产精品麻豆人妻色哟哟久久| 日本色播在线视频| 午夜爱爱视频在线播放| 国产白丝娇喘喷水9色精品| 久久精品综合一区二区三区| 热re99久久精品国产66热6| freevideosex欧美| 女人被狂操c到高潮| 久久99精品国语久久久| 91精品国产九色| 亚洲欧美精品专区久久| 久久久久久久大尺度免费视频| 3wmmmm亚洲av在线观看| 国产在线一区二区三区精| 国产精品久久久久久精品电影| 欧美成人午夜免费资源| 99久久中文字幕三级久久日本| 老女人水多毛片| 一本色道久久久久久精品综合| 涩涩av久久男人的天堂| 欧美精品国产亚洲| 亚洲人成网站在线观看播放| 在线观看美女被高潮喷水网站| 国产毛片在线视频| 国产av国产精品国产| 色综合色国产| 一本久久精品| 边亲边吃奶的免费视频| 久久久久性生活片| 精品久久久久久久久亚洲| 欧美极品一区二区三区四区| 亚洲av.av天堂| 国产日韩欧美亚洲二区| 大片电影免费在线观看免费| 国产又色又爽无遮挡免| 美女xxoo啪啪120秒动态图| 欧美老熟妇乱子伦牲交| 69人妻影院| 看十八女毛片水多多多| 99久久精品国产国产毛片| av在线蜜桃| av女优亚洲男人天堂| 91午夜精品亚洲一区二区三区| 韩国av在线不卡| 深夜a级毛片| 欧美成人精品欧美一级黄| 男女下面进入的视频免费午夜| 在线亚洲精品国产二区图片欧美 | 免费在线观看成人毛片| 国产精品无大码| 国产亚洲一区二区精品| 视频区图区小说| 九色成人免费人妻av| 可以在线观看毛片的网站| 综合色av麻豆| 久久久久国产网址| 亚洲综合精品二区| 亚洲经典国产精华液单| 亚洲,一卡二卡三卡| av天堂中文字幕网| 日韩亚洲欧美综合| 国产精品国产三级专区第一集| 婷婷色av中文字幕| 久久精品国产自在天天线| 精品亚洲乱码少妇综合久久| 国产精品久久久久久精品电影| 国产精品麻豆人妻色哟哟久久| 九色成人免费人妻av| 91精品一卡2卡3卡4卡| 国产黄片视频在线免费观看| 人人妻人人澡人人爽人人夜夜| 男插女下体视频免费在线播放| 边亲边吃奶的免费视频| 美女高潮的动态| 哪个播放器可以免费观看大片| 国产久久久一区二区三区| 一本色道久久久久久精品综合| 成人欧美大片| 亚洲精品久久午夜乱码| 久久精品国产亚洲av天美| 看非洲黑人一级黄片| 中文在线观看免费www的网站| 黑人高潮一二区| 亚洲丝袜综合中文字幕| 午夜亚洲福利在线播放| 国产欧美亚洲国产| 日本一二三区视频观看| 亚洲精品乱码久久久久久按摩| 欧美丝袜亚洲另类| eeuss影院久久| 亚洲精品一区蜜桃| 菩萨蛮人人尽说江南好唐韦庄| 亚洲国产成人一精品久久久| 一区二区三区精品91| 亚洲aⅴ乱码一区二区在线播放| 成人无遮挡网站| 国产探花极品一区二区| 一级毛片电影观看| 欧美97在线视频| 国产精品偷伦视频观看了| 国内揄拍国产精品人妻在线| 丝袜美腿在线中文| 国产伦理片在线播放av一区| 97在线视频观看| 极品少妇高潮喷水抽搐| 免费看a级黄色片| 久久久久九九精品影院| 黄色一级大片看看| 青青草视频在线视频观看| 欧美变态另类bdsm刘玥| 男女边吃奶边做爰视频| 中文字幕久久专区| 久久久久久久久大av| 内地一区二区视频在线| 精品人妻视频免费看| 少妇熟女欧美另类| 亚洲激情五月婷婷啪啪| 熟女人妻精品中文字幕| 国产高潮美女av| 亚洲国产色片| 人妻夜夜爽99麻豆av| 一本一本综合久久| 夜夜看夜夜爽夜夜摸| 日韩免费高清中文字幕av| av天堂中文字幕网| 好男人视频免费观看在线| 波野结衣二区三区在线| 欧美成人a在线观看| 女人久久www免费人成看片| 亚州av有码| 久久久精品欧美日韩精品| 中文在线观看免费www的网站| 国内揄拍国产精品人妻在线| 高清在线视频一区二区三区| 国产成人91sexporn| 晚上一个人看的免费电影| 黄片wwwwww| 熟妇人妻不卡中文字幕| 亚洲怡红院男人天堂| 97在线人人人人妻| 国产精品一及| 国产日韩欧美亚洲二区| 亚洲一级一片aⅴ在线观看| av在线老鸭窝| 777米奇影视久久| 久久精品熟女亚洲av麻豆精品| 国产成人福利小说| 午夜视频国产福利| 青春草亚洲视频在线观看| 亚洲无线观看免费| 国产男女内射视频| 国产成人aa在线观看| 女人被狂操c到高潮| 国产一区二区三区综合在线观看 | 国产又色又爽无遮挡免| 亚洲综合色惰| 国产一区有黄有色的免费视频| 亚洲自偷自拍三级| 特大巨黑吊av在线直播| 成人黄色视频免费在线看| 男人狂女人下面高潮的视频| 黄片无遮挡物在线观看| 欧美极品一区二区三区四区| 国产伦精品一区二区三区视频9| 青春草国产在线视频| 亚洲一区二区三区欧美精品 | 大香蕉久久网| 国产黄片视频在线免费观看| 三级经典国产精品| 大香蕉97超碰在线| 人体艺术视频欧美日本| 大香蕉97超碰在线| 成年女人看的毛片在线观看| 大香蕉97超碰在线| 22中文网久久字幕| 午夜福利视频1000在线观看| av网站免费在线观看视频| 亚洲国产色片| 亚洲精品aⅴ在线观看| 国产人妻一区二区三区在| 制服丝袜香蕉在线| 在线 av 中文字幕| 中文资源天堂在线| 91精品一卡2卡3卡4卡| 久久久久久国产a免费观看| 一级爰片在线观看| 精品一区在线观看国产| 国产精品99久久久久久久久| 免费观看性生交大片5| 三级经典国产精品| 色播亚洲综合网| 在线a可以看的网站| 色婷婷久久久亚洲欧美| 天天躁夜夜躁狠狠久久av| 国产视频内射| 国产一区二区在线观看日韩| 国产乱人偷精品视频| 中文精品一卡2卡3卡4更新| 少妇的逼水好多| 观看美女的网站| 成人高潮视频无遮挡免费网站| 久久精品久久精品一区二区三区| 久久久久久久久久人人人人人人| 99久国产av精品国产电影| 免费av不卡在线播放| av播播在线观看一区| 性色av一级| 91午夜精品亚洲一区二区三区| 亚洲国产精品成人久久小说| 寂寞人妻少妇视频99o| 丝袜脚勾引网站| 2022亚洲国产成人精品| 国产精品一二三区在线看| 午夜福利在线在线| 美女视频免费永久观看网站| 久久久成人免费电影| 精品久久久噜噜| 男女国产视频网站| 十八禁网站网址无遮挡 | 一级毛片久久久久久久久女| 国产 一区精品| 久久久久久久大尺度免费视频| 精品一区二区免费观看| 一级爰片在线观看| 日韩强制内射视频| 国产精品一区二区三区四区免费观看| 欧美成人一区二区免费高清观看| 两个人的视频大全免费| 性色av一级| 少妇人妻久久综合中文| 男人添女人高潮全过程视频| 99九九线精品视频在线观看视频| 成年版毛片免费区| 一本一本综合久久| 一区二区三区精品91| 亚洲图色成人| 国产成人福利小说| 一个人看视频在线观看www免费| 丰满少妇做爰视频| 日韩人妻高清精品专区| 晚上一个人看的免费电影| 丰满人妻一区二区三区视频av| 国产在线男女| 亚洲色图综合在线观看| av.在线天堂| 国产精品国产三级国产专区5o| 亚洲欧洲国产日韩|