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

    Synthesis and antifungal activities of phenylhydrazono derivatives as the novel fructose-1,6-bisphosphate aldolase inhibitors

    2020-06-18 00:30:14GUOLiGUOYanrongRENYanliangHANXinya

    GUO Li,GUO Yanrong,REN Yanliang,HAN Xinya

    (1.Hubei Ecological and Environmental Monitoring Center Station,Wuhan 430072,China;2.School of Chemistry and Chemical Engineering,Anhui University of Technology,Maanshan,Anhui 243002,China;3.Laboratory of Pesticide &Chemical Biology,Key Laboratary of Ministry of Education,College of Chemistry,Central China Normal University,Wuhan 430079,China)

    Abstract:Fructose-1,6-bisphosphate aldolase from C.albicans (CaFBA-II) is an attractive new target for the discovery of drugs to combat invasive fungal infection.Our previous study demonstrated that β-arylhydrazono-α,γ-dicarbonylderivatives exhibit moderate inhibitory against CaFBA-II.Herein,several new phenylhydrazono derivatives were found to potently inhibit CaFBA-II.The experimental results show that the NO2 group in the R2 position is favorable,but COOH,SO3H and halogen atom are unfavorable for their CaFBA-II inhibitory activities.Especially,when the R3 was substituted by sulfonic acid group,compound 2g possessed the highest inhibitory activity (IC50=200 nmol·L-1).Furthermore,compound 2g were selected as representative molecule,the binding mode of 2g and the surrounding residues in the active site of CaFBA-II were elucidated by jointly using DOX methods and molecular dynamics (MD) simulations.Notably,antifungal experiments demonstrate that most of our designed compounds possess moderate inhibitory activities (MIC80=16-64 μg · mL-1) against C.albicans.The present results suggest that the phenylhydrazono derivatives can be used as the lead compounds of novel drugs against fungal pathogens of humans in the future.

    Key words:novel inhibitors;rational drug design;phenylhydrazono derivatives;CaFBA-II;antifungal activities

    Bacterial and fungal infections are an important public health problem of major concern[1].C.albicansis the primary systemic fungal affecting humans.Azoles are the most well-known antifungal drugs all over the word to combatC.albicansinfections.However,due to extensive use of these drugs in humans,the resistance of pathogenic microorganisms to current antibiotics is an important public health problem[2-3].Fructose-1,6-bisphosphate aldolase (FBA) is an essential enzyme in GNG pathway,and occurs in two distinct classes in the biosphere,and catalyzes the reversible aldol condensation of dihydroxyacetonephosphate (DHAP) and glyceraldehyde 3-phosphate (GAP) to fructose-1,6-bisphosphate (FBP)[4-6].Class I FBA (FBA-I),which is present in higher organisms (plants and animals) and some prokaryotes,forms a Schiffbase intermediate between the keto substrate (FBP or DHAP) and a lysine residue of the active site.Class II FBA (FBA-II) requires a divalentmetal ion (usually zinc or cobalt ion) to polarize the keto carbonyl group of the substrate (FBP or DHAP) and to stabilize the enediolate intermediate formed during catalysis.Because of their occurrence in many pathogenic microbes (bacteria,yeasts,parasites) and their absence in animals,F(xiàn)BA-II has been regarded as a particularly attractive new target and a lot of the attention was focused on the drug design of the enzyme[5,7-10].In previous work,we discovered several potent FBA-II inhibitors againstC.albicans(CaFBA-II)[11-12].Herein,a positive progress was achieved in this field,several new phenylhydrazono derivatives were found to exhibit potent CaFBA-II inhibitory activities and antifungal effects againstC.albicans.Furthermore,the probable binding-modes between the hit compounds and CaFBA-II were also analyzed by jointly using the DOX method[13](i.e.,molecular docking,ONIOM3 and XO),and molecular dynamics (MD) simulations and enzymatic assays.

    1 Material and Methods

    1.1 Chemistry

    All synthetic work was performed in air and at room temperature.All solvents and chemicals were obtained from commercial sources and used as received.The proton (1H) and carbon (13C) NMR spectra were recorded on a Varian Mercury-Plus 400 or 600 MHz spectrometer in DMSO-d6or acetone-d6.Chemical shifts are given in parts per million (ppm) with TMS as the internal reference.Analytical thin-layer chromatography (TLC) was carried out on Merck silica gel F-254 plates.Flash chromatography purifications were performed on Merck silica gel 60 (230-400 mesh) as the stationary phase,and MeOH and CH2Cl2were used as eluents.High-resolution mass spectrometry (HRMS) analysis was performed using a hydrid IT-TOF mass spectrometer with ESI interface (Shimadzu,Kyoto,Japan) and Agilent 6224 Accurate-Mass time-of-flight mass spectrometer with ESI interface (Agilent Technologies,Waldbronn,Germany).If necessary,the purity was determined by high performance liquid chromatography (HPLC).The purity of novel synthesized compounds was 96% or higher,as illustrated in supporting information.

    General procedure for the preparation of phenylhydrazono (1,2a-g).A suspension of the substituted aniline (5 mmol) in NaOH solution (0.5 mol · L-1,10 mL) were added NaNO2(5 mmol) and 36% HCl (5 mL) at a temperature below 5 ℃.The resulting diazonium solution was used directly without purification.

    A solution of NaOH (0.5 mol · L-1,5 mL) in ethanol (5 mL) was added a portion of theβ-diketone (5 mmol).The solution was cooled in an ice bath to 0 ℃,and then added with the diazonium solution mentioned above under vigorous stirring for 1 h.The mixture was filtered,and the solid product was washed with water,then recrystallized in ethanol to produce1,2a-g.Notably,this approach did engender important safety concerns due to the highly explosive diazonium salt which may decomposed intensively under a relatively high temperature while releasing heat and outgassing.In addition,the raw material sodium nitrite decomposes at 175 ℃,which may easily cause the combustion or explosion of organics in the reaction system,and the combustion or explosion may also be caused by strong oxidants.

    1,1,1,5,5,5-Hexafluoro-3-(2-(2-hydroxy-5-nitrophenyl) hydrazono)pentane-2,4-dione (1).Dark red powder,yield 75%.1H NMR (400 MHz,DMSO-d6)δ:8.98 (d,J=9.7 Hz,1H),8.05 (dd,J=11.3,1.9 Hz,1H),6.59 (d,J=10.1 Hz,1H).13C NMR (101 MHz,DMSO-d6)δ:174.75,132.48,132.13,132.10,129.88,129.84,124.08,90.84.HRMS (ESI):m/zcalculated for C11H5F6N3O5+K+[M+K]+411.9765.Found:412.0574.

    Ethyl 4,4,4-trifluoro-2-(2-(2-hydroxy-5-nitrophenyl) hydrazono)-3-oxobutanoate (2a).Dark red powder,yield 75%.1H NMR (400 MHz,DMSO-d6)δ:15.19 (s,1H),8.49 (s,1H),8.17 (d,J=8.5 Hz,1H),7.22 (d,J=8.9 Hz,1H),4.41 (d,J=6.7 Hz,2H),1.38 (t,J=7.0 Hz,3H).13C NMR (151 MHz,DMSO-d6)δ162.55,153.09,140.64,129.33,123.26,122.99,116.44,110.55,62.13,14.25.HRMS (ESI):m/zcalculated for C12H10F3N3O6-H+[M-H]-348.0443.Found:348.0453.

    Ethyl 4,4,4-trifluoro-2-(2-(2-hydroxy-5-carbonylphenyl) hydrazono)-3-oxobutanoate (2b).Brown powder,yield 45%.1H NMR (400 MHz,DMSO-d6)δ:1H NMR (600 MHz,DMSO-d6)δ12.83 (s,1H),8.10 (s,1H),7.72 (d,J=8.4 Hz,1H),7.05 (d,J=8.5 Hz,1H),4.33 (q,J=7.1 Hz,2H),1.32 (t,J=7.1 Hz,3H).13C NMR (151 MHz,DMSO-d6)δ:176.88,167.23,165.70,163.62,152.66,138.38,128.82,123.85,122.58,116.45,116.20,89.06,61.49,14.36.HRMS (ESI):m/zcalculated for C13H11F3N2O6-H+[M-H]-347.0491.Found:347.0489.

    Ethyl 4,4,4-trifluoro-2-(2-(2-hydroxy-5-sulfophenyl) hydrazono)-3-oxobutanoate (2c).Dark powder,yield 35%.1H NMR (400 MHz,DMSO-d6)δ:15.47 (s,1H),11.79 (s,1H),8.02 (d,J=2.0 Hz,1H),7.46 (dd,J=8.4,2.0 Hz,1H),6.95 (d,J=8.5 Hz,1H),4.39 (q,J=7.1 Hz,2H),1.35 (t,J=7.1 Hz,3H).13C NMR (151 MHz,DMSO-d6)δ:162.55,150.44,141.22,133.09,132.30,128.00,119.44,116.59,61.58,14.73.HRMS (ESI):m/zcalculated for C12H11F3N2O7S[M] 384.0234.Found:384.0000.

    Ethyl 4,4,4-trifluoro-2-(2-(5-chloro-2-hydroxyphenyl) hydrazono)-3-oxobutanoate (2d).Dark red powder,yield 40%.1H NMR (600 MHz,DMSO-d6)δ:13.34 (s,1H),11.47 (s,1H),11.14 (s,1H),7.37 (s,1H),7.29 (d,J=8.4 Hz,1H),7.17 (d,J=8.6 Hz,1H),7.04 (t,J=9.0 Hz,1H),4.34 (dd,J=13.8,6.7 Hz,2H),1.32 (t,J=7.1 Hz,3H).13C NMR (151 MHz,DMSO-d6)δ:162.34,149.62,134.45,133.92,129.78,124.13,119.46,119.19,61.58,14.60.HRMS (ESI):m/zcalculated for C12H10ClF3N2O4-H+[M-H]-337.0203.Found:337.0208.

    Ethyl 4,4,4-trifluoro-2-(2-(2-hydroxyphenyl)hydrazono)-3-oxobutanoate (2e).Red powder,yield 58%.1H NMR (600 MHz,DMSO-d6)δ:13.55 (s,1H),10.91 (s,1H),7.47 (d,J=8.0 Hz,1H),7.14 (t,J=7.7 Hz,1H),7.05 (d,J=8.0 Hz,1H),6.97 (t,J=7.6 Hz,1H),4.34 (q,J=7.0 Hz,2H),1.32 (t,J=7.1 Hz,3H).13C NMR (151 MHz,CDCl3)δ:170.61,163.82,146.33,127.76,127.49,121.46,117.43,117.18,109.99,61.97,14.06.HRMS (ESI):m/zcalculated for C12H11F3N2O4-H+[M-H]-303.0592.Found:303.0604.

    Ethyl 4,4,4-trifluoro-2-(2-(3-chloro-2-hydroxy-5-nitrophenyl)hydrazono)-3-oxobutanoate (2f).Black powder,yield 43%.1H NMR (600 MHz,DMSO-d6)δ:8.73 (s,1H),6.63 (d,J=2.2 Hz,1H),6.53 (d,J=8.1 Hz,1H),6.39 (dd,J=8.1,2.2 Hz,1H),4.37 (s,2H),1.23 (d,J=4.4 Hz,3H),1.18 (s,9H).13C NMR (151 MHz,DMSO-d6)δ:166.25,162.84,136.52,129.75,125.21,124.51,124.27,122.85,120.84,109.45,61.44,14.41.HRMS (ESI):m/zcalculated for C12H9ClF3N3O6-H+[M-H]-382.0054.Found:382.0062.

    Ethyl 4,4,4-trifluoro-2-(2-(2-hydroxy-5-nitro-3-sulfophenyl) hydrazono)-3-oxobutanoate (2g).Dark red powder,yield 50%.1H NMR (400 MHz,DMSO-d6)δ:15.13 (s,1H),13.23 (s,1H),12.31 - 12.09 (m,1H),8.50 (d,J=2.7 Hz,1H),8.29 (d,J=2.7 Hz,1H),4.47 - 4.34 (m,2H),1.36 (dt,J=25.0,7.1 Hz,3H).13C NMR (101 MHz,DMSO-d6)δ:167.29,155.22,145.86,137.89,137.08,132.73,124.10,121.34,66.31,19.46.HRMS (ESI):m/zcalculated for C12H10F3N3O9S-H+[M-H]-428.0012.Found:428.0281

    1.2 Enzyme inhibition activity of CaFBA-II

    To evaluate the inhibitory activity of hit compounds screened and synthesized in the present study,the half maximal inhibitory concentration (IC50) values of hit compounds were determined at the CaFBA-II recombinant protein level by linear regression analysis using the program package origin.The gene sequence of FBA-II fromC.albicanswas obtained from NCBI (XM_717597.1)[14].CaFBA-II was expressed inE.coliBL21(DE3) cells as described previously[15].Commercial preparations from Sigma of glycerol 3-phosphate dehydrogenase (GPDH) from rabbit muscle,triosephosphate isomerase (TIM) from rabbit muscle,and Ra-FBA-I were used.Ra-FBA-I and CaFBA-II activity and its inhibition were determined in reactions using an NADH linked enzymatic assay.The activity of the cleavage reaction was assayed in Tris-K-acetate buffer (0.1 mol·L-1,pH 7.4,with 0.2 mol·L-1potassium acetate for FBA-II),using fructose bisphosphate and inhibitor at the appropriate concentrations,NADH (0.3 mmol·L-1),GPDH (0.26 U),TIM (1 U) and aldolase (4 mU) in a cuvette to give a final volume of 0.42 mL.The decrease in absorbance of NADH at 340 nm was monitored on a spectrophotometer over 5 min.

    To determine thecorresponding inhibitor constants (IC50),the initial rate data taken at a saturating substrate concentration,a fixed effector concentration,and systematically varied inhibitor concentrations were fit to a Hill equation,V=V0-(V0-V∞)/((I0.5/I)n+1)[16],whereV,V0,andV∞are the velocity,maximum velocity (atI=0),and the limiting velocity (atIsaturating),respectively;n is the Hill coefficient associated with the inhibitor;andIC50is the inhibitor concentration at 50% inhibition.All kinetic data were fit to the growth/sigmoidal model in origin 7.0 software.The inhibition assay was run under a concentration of compounds2fand2gillustrated in Fig.1.

    Fig.1 The inhibition assay was run in the concentration range (0-50 μmol·L-1) of hit compounds 2f and 2g (The inhibition of DMSO was taken as blank control,for all hit compounds dissolved in DMSO)

    1.3 Inhibitory Assays on Candida albicans

    Invitroantifungal activity was measured as the minimum inhibitory concentration (MIC) that achieved 80% inhibition of the C.albicans using the broth microdilution method,as given in document M27-A3[17]and M38-A2[18]from the National Committee for Clinical Laboratory Standards (CLSI).TheC.albicansstrain (SC5314) were obtained from the Second Military Medical University.The test compounds were dissolved in DMSO and serially diluted in growth media.The microtiter plates were incubated at 35 ℃ for 48 h for the Candida species.DMSO was used as a growth control and fluconazole (Pfizer) and ketoconazole (Sigma) were used as controls against the tested fungi.

    Inhibitory effect of hit compounds on the growth of weed roots and shoots

    The inhibitory effect of hit compounds on the Growth of Weed roots and shoots was executed as described previously[19].A mixture of 4 g agar powder and 0.5 L distilled water was heated to melt and then cooled down about 45 ℃.The 20 μL solution containing 1.5 mg · mL-1,4.5 mg · mL-1,15 mg · mL-1compounds and 280 μL melting agar were mixed in 96-well plate.The seeds ofCynodondactylonL.,and Arabidopsis thaliana were put on the surface of the agar.The 96-well plates were kept in (25 ± 1) ℃,12 h in the light and 12 h in the dark alternatively,and 50%-55% relative humidity for seven days.Then the lengths of roots and shoots,and the inhibitory rate was determined.The commercial herbicide acetocblor,glyphosate,and meostrione were used as the positive control.

    2 Results and Discussions

    2.1 Synthesis and optimization of phenylhydrazono derivatives

    Our previous work[11]suggested that the hydroxyl group of phenylhydrazonoderivatives is essential to the binding of hit compounds and its target,because it could coordinate with the Zn(II) ion of CaFBA-II.To obtain more potent hit compounds,several novel phenylhydrazono derivatives with hydroxyl group were further synthesized,as illustrated in Tab.1,together with their CaFBA-II inhibitory activities and antifungal activities.

    Tab.1 The inhibitory activities of the hit compounds against CaFBA-II,together with the inhibitions against C.albicans

    Compounds1,2a-gwere synthesized according to procedure described previously[20],as illustrated in Scheme 1.

    Scheme 1 Synthetic Route of the hit compounds1and2a-g

    Reagents and conditions:a) NaNO2/HCl,0-5 ℃;b) CF3COCH2COR1,NaOH,0-5 ℃

    It is shown in our previous study that the methyl group onα,γ-dicarbonyl substructure are unfavorable to the inhibitory activities of hit compounds,thus we attempt to synthesize the new compounds with trifluoromethyl or ethyl ester groups inR1position.As listed in Tab.1,the nitro group in theR2position of the benzene ring is favorable for the inhibitory activities of hit compounds (compounds1and2a).When the nitro group was changed to a carboxyl group or a sulfonic acid group,such as2band2c,the corresponding inhibition activities against CaFBA-II decreased approximately 10-fold compared with that of2a.Notably,when the Cl and SO3H were introduced to theR3,compounds2fand2gexhibit potent CaFBA-II inhibitory activities,with theIC50of 1.5 μmol·L-1and 0.2 μmol·L-1.These results indicate that a hydrophilic group in theR3position is favorable.

    The ribbon represents the helix and sheet of protein,the light blue stick represents the hit compound 2g,and the surrounding important residues were representedby the creamy white stick.The green dashed line represents the important interactions of the ligands and target enzymeFig.2 The proposed binding-mode of compound 2g into the active site of CaFBA-II optimized by MD simulations

    The binding mode of2gand CaFBA-II were further predicted by using DOX strategy[13,21],and optimized by molecular dynamics (MD) simulations (Fig.2 and Fig.3) with the PMEMD module in the AMBER 12 package.Additionally,using FBA-II fromE.coli(PDB ID:1B57) as template,the plausible homology model of CaFBA-II were performed by SWISSMODEL server,which has been documented previously[11].As shown in Fig.1,the sulfonic acid and nitro group on the benzene ring of2gcan form hydrogen bonds with Asn224 and Ser268,respectively.The hydroxyl group on the benzene ring of2gcould coordinate with Zn(II) ion.The hydrogen-bonds between carbonyl group and Asn36 could also be observed in Fig.1.However,when the sulfonic acid group inR3position was substituted by H and Cl,the hydrogen bond between hit compound and Asn224 was lost,thus the inhibitory activity of compound2a(3.2 μmol·L-1and2f(1.5 μmol·L-1) will decrease approximately 10-fold compared with that of2g(0.2 μmol·L-1).This to some extend explain why2gpossess the highest inhibitory activity against CaFBA-II.

    Fig.3 Plots of the molecular dynamic (MD) simulation time vs root-mean-square deviation (RMSD) of the all atoms of CaFBA-II in complex with 2g

    According to the binding mode in Fig.2,the carboxyl group or the sulfonic acid group inR2position can also form hydrogen-bond with Ser268,however,our results show that theIC50values of2b(34 μmol·L-1) and2c(25 μmol·L-1) are 10-fold larger than that of2a(3.2 μmol·L-1).In order to explain this case,the solvation energies of2a,2band2cwere further calculated at wb97xd/6-31G*-smd theory level by using Gaussian 09 software package,which are of 61.45 kJ· mol-1,75.81 kJ · mol-1and 71.66 kJ· mol-1,respectively.Our current calculation results reveal that the larger salvation effect counteracts the binding of hit compounds into CaFBA-II,thus compounds2band2cexhibit weak inhibitory activities than2a,In addition,when theR2position was substituted by chlorine,the inhibitory activity of2ddecreased approximately 3-fold compared with that of2a.This is likely due to the weak interaction (halogen bond) of2dcompared to that (hydrogen bond) of2a.Taken together,the nitro-group in theR2position is favorable.In the future,the compounds with more different substituents inR3position will be designed and synthesized to further insight the structure-activity relationships (SAR) of these hit compounds.

    Additionally,the compound (2g) with highest inhibitory activity against CaFBA-II were also tested against a representative of mammalian FBA-I,isozyme A from rabbit muscle.Our experimental results show that the inhibitory activity of2gagainst rabbit muscle FBA-I(Ra-FBA-I) is larger than 200 μmol·L-1,thus the inhibitory selectivity of this compound against FBA-I and FBA-II is up to 1,000.The selectivity for FBA-II validates that the compound2gexhibit pure competitive inhibition patterns.

    2.2 Antifungal evaluation

    All of our present compounds were assayed using a standard procedure for the inhibition of growth of cultivatedC.albicans.[17-18].The results show that compounds2eand2fexhibit moderate inhibitory activities againstC.albicanswith theMIC80values of 32 μg · mL-1.In comparison,compounds1,2a,2b,2dand2gshow more potent inhibitory activities (16 μg · mL-1) againstC.albicans.The dose responses of2dand2gin the inhibition of growth of cultivatedC.albicanswere illustrated in Fig.4.To our knownledge,several FBA-II inhibitors have been reported previously,however,few of them demonstrate inhibitory activities against microorganisms so far.This is the first report of phenylhydrazonoderivatives with antifungal effects againstC.albicans.In addition,the inhibitions of compound2gagainstC.neoforman(MIC80=16 μg · mL-1),Synechocystissp.PCC 6803 (EC50=12 μg · mL-1),Cynodondactylon L (MIC80>1000 μg · mL-1) and Arabidopsis thaliana (MIC80>1000 μg · mL-1) were tested.The results show that compound2gcould inhibit the microorganisms with FBA-II (e.g.,C.neoformanandSynechocystissp.PCC 6803),but couldn’t inhibit the organisms without FBA-II (e.g.,CynodondactylonL andArabidopsisthaliana).Thus,compound2gshows selectivity in itsinvitroactivities.

    Fig.4 Dose response for compounds 2d and 2g in the inhibition of growth of cultivated C.albicans

    3 Conclusion

    Several new phenylhydrazono derivatives were discovered to potently inhibit CaFBA-II.The current experimental results show that compound1,2a,2fand2gexhibit potent inhibitory activities against CaFBA-II (IC50<5.0 μmol·L-1).Especially,compound2gpossessed the highest inhibitory activity (IC50=200 nmol·L-1).Furthermore,compound2gwere selected as the representative molecule,the binding mode of2gand the surrounding residues in the active site of CaFBA-II were elucidated by jointly using DOX methods and molecular dynamics (MD) simulations.The structure-activity relationships (SAR) analysis of the designed compounds can be consistently explained according to our proposed binding-modes of2gin CaFBA-II.These results further confirm and aid our understanding of the binding mode of inhibitors targeting CaFBA-II.The binding mode of2gand CaFBA-II show that the hydrophilic group in theR2position is favorable own to the hydrogen-bond between hit compound and Ser268 of CaFBA-II,but COOH,SO3H inR2position is unfavorable for their CaFBA-II inhibitory because of their larger salvation effect.It should be noticed that compound2gshows high selectivity (up to 1,000) for CaFBA-II compared with Ra-FBA-I.Notably,invitrotest results demonstrate that most of our designed compounds possess moderate inhibitory activities (MIC80=16-64 μg · mL-1) againstC.albicans.The present enzymatic and in vitro experiments suggest that the phenylhydrazono derivatives can be used as the lead compounds of novel drugs against fungal pathogens of humans in the future.

    国产精品 国内视频| 国产精品综合久久久久久久免费| 亚洲五月婷婷丁香| 日本a在线网址| 国产一区二区激情短视频| 在线观看免费视频日本深夜| 亚洲欧美一区二区三区黑人| 国产成人精品无人区| 久久亚洲精品不卡| 久久久久国产精品人妻aⅴ院| 中文资源天堂在线| 真人做人爱边吃奶动态| 精品一区二区三区av网在线观看| 国产伦一二天堂av在线观看| av在线蜜桃| 丰满的人妻完整版| 这个男人来自地球电影免费观看| 日韩欧美国产在线观看| 九九在线视频观看精品| 久久精品aⅴ一区二区三区四区| 国产精品女同一区二区软件 | 99久久成人亚洲精品观看| 成人永久免费在线观看视频| 99国产极品粉嫩在线观看| 亚洲精品粉嫩美女一区| 国产真实乱freesex| 国产精品久久电影中文字幕| 中文字幕久久专区| 免费在线观看影片大全网站| 九九在线视频观看精品| 亚洲色图av天堂| 精品不卡国产一区二区三区| 黄色日韩在线| 亚洲国产日韩欧美精品在线观看 | 无限看片的www在线观看| 99国产极品粉嫩在线观看| 久久久久久久午夜电影| 亚洲欧美激情综合另类| 国产午夜精品论理片| www.自偷自拍.com| 亚洲精品粉嫩美女一区| 一个人免费在线观看电影 | 无限看片的www在线观看| 9191精品国产免费久久| 成人午夜高清在线视频| 一级黄色大片毛片| 老汉色∧v一级毛片| 久久香蕉精品热| 欧美日本亚洲视频在线播放| 国产又色又爽无遮挡免费看| 欧美日韩中文字幕国产精品一区二区三区| 国产av不卡久久| 国产精品久久久人人做人人爽| 欧美黑人欧美精品刺激| 啦啦啦免费观看视频1| 亚洲自偷自拍图片 自拍| 99久久精品热视频| 久久久久久久久久黄片| 99在线人妻在线中文字幕| 国产亚洲欧美在线一区二区| 亚洲激情在线av| 一二三四在线观看免费中文在| 99久久无色码亚洲精品果冻| 欧美黑人巨大hd| tocl精华| 久久久久久久久免费视频了| 欧美日韩黄片免| 在线免费观看的www视频| 又大又爽又粗| 窝窝影院91人妻| 欧美一区二区精品小视频在线| 波多野结衣高清作品| 国产精华一区二区三区| 欧美日韩一级在线毛片| 男人的好看免费观看在线视频| 黄片大片在线免费观看| 丰满的人妻完整版| 久久人人精品亚洲av| 97人妻精品一区二区三区麻豆| 久久精品aⅴ一区二区三区四区| 国产又色又爽无遮挡免费看| 欧美午夜高清在线| 日韩欧美一区二区三区在线观看| 淫秽高清视频在线观看| 首页视频小说图片口味搜索| 激情在线观看视频在线高清| 草草在线视频免费看| 久久亚洲真实| 真实男女啪啪啪动态图| 欧美三级亚洲精品| 男女下面进入的视频免费午夜| 欧美黑人欧美精品刺激| 国产精品自产拍在线观看55亚洲| 久久久久久久精品吃奶| 日韩有码中文字幕| 婷婷丁香在线五月| 人人妻人人看人人澡| 亚洲精品在线美女| 999久久久国产精品视频| 99国产极品粉嫩在线观看| 亚洲性夜色夜夜综合| 99久久国产精品久久久| 久久久久精品国产欧美久久久| 亚洲国产欧美网| 欧美日韩综合久久久久久 | 高潮久久久久久久久久久不卡| 亚洲电影在线观看av| avwww免费| 国产一区二区三区在线臀色熟女| 啪啪无遮挡十八禁网站| 18禁黄网站禁片午夜丰满| 亚洲一区高清亚洲精品| 99久久精品一区二区三区| 一区二区三区高清视频在线| 免费观看人在逋| 99久国产av精品| 日韩精品中文字幕看吧| 日韩欧美一区二区三区在线观看| 桃红色精品国产亚洲av| 国产av麻豆久久久久久久| 国产亚洲精品久久久久久毛片| 午夜福利视频1000在线观看| 好男人电影高清在线观看| 国产精品亚洲一级av第二区| 搡老熟女国产l中国老女人| 一边摸一边抽搐一进一小说| a级毛片a级免费在线| 18美女黄网站色大片免费观看| 欧美成狂野欧美在线观看| 久久久国产成人免费| 中出人妻视频一区二区| 亚洲av成人精品一区久久| 国产伦在线观看视频一区| 男人舔女人的私密视频| 床上黄色一级片| 午夜精品久久久久久毛片777| 怎么达到女性高潮| 一二三四在线观看免费中文在| www.自偷自拍.com| 亚洲av免费在线观看| av视频在线观看入口| 中出人妻视频一区二区| 欧美高清成人免费视频www| 小说图片视频综合网站| 97人妻精品一区二区三区麻豆| 舔av片在线| 国产成人aa在线观看| 一卡2卡三卡四卡精品乱码亚洲| 窝窝影院91人妻| 免费在线观看视频国产中文字幕亚洲| 999精品在线视频| 亚洲精品在线美女| 欧美日本视频| 真人做人爱边吃奶动态| 精品一区二区三区av网在线观看| 天堂影院成人在线观看| 亚洲专区字幕在线| 欧美乱色亚洲激情| 欧美日韩国产亚洲二区| 看免费av毛片| 国产伦精品一区二区三区四那| 国产精品一区二区三区四区免费观看 | 综合色av麻豆| 白带黄色成豆腐渣| 亚洲精品在线观看二区| АⅤ资源中文在线天堂| 国产精品久久视频播放| 少妇的丰满在线观看| 欧美性猛交黑人性爽| 久久久久国内视频| 亚洲一区高清亚洲精品| 亚洲精品乱码久久久v下载方式 | 看片在线看免费视频| 亚洲成人久久爱视频| 国产aⅴ精品一区二区三区波| 欧美日韩综合久久久久久 | 熟女少妇亚洲综合色aaa.| 欧美中文综合在线视频| 亚洲最大成人中文| 亚洲一区二区三区色噜噜| 国内揄拍国产精品人妻在线| 丰满人妻熟妇乱又伦精品不卡| a级毛片a级免费在线| 成人性生交大片免费视频hd| 日韩欧美在线二视频| 亚洲欧美日韩高清在线视频| 无限看片的www在线观看| 伊人久久大香线蕉亚洲五| 男人舔女人下体高潮全视频| 一区福利在线观看| 最新中文字幕久久久久 | 一a级毛片在线观看| 国产一区二区三区在线臀色熟女| 欧美性猛交╳xxx乱大交人| 九九久久精品国产亚洲av麻豆 | 2021天堂中文幕一二区在线观| 美女黄网站色视频| 中亚洲国语对白在线视频| 国内久久婷婷六月综合欲色啪| 国产精品一区二区精品视频观看| 人妻久久中文字幕网| 亚洲无线观看免费| 18禁黄网站禁片免费观看直播| 日本a在线网址| 午夜成年电影在线免费观看| 手机成人av网站| 99国产极品粉嫩在线观看| 亚洲成人久久性| 美女免费视频网站| 中文字幕熟女人妻在线| 久9热在线精品视频| 亚洲av熟女| 国产精品免费一区二区三区在线| 国产精品自产拍在线观看55亚洲| 欧美成狂野欧美在线观看| 欧美一级毛片孕妇| 国产精品1区2区在线观看.| 又爽又黄无遮挡网站| 国产黄片美女视频| 宅男免费午夜| 久久久久国内视频| 久9热在线精品视频| 久久久精品大字幕| 中文亚洲av片在线观看爽| 精品国产美女av久久久久小说| 99国产综合亚洲精品| 最近在线观看免费完整版| 国产真实乱freesex| 精品久久久久久久毛片微露脸| 精品熟女少妇八av免费久了| 精品一区二区三区av网在线观看| 99久久久亚洲精品蜜臀av| 老汉色av国产亚洲站长工具| 欧美日韩乱码在线| 色综合亚洲欧美另类图片| 母亲3免费完整高清在线观看| 成人亚洲精品av一区二区| 国产精品一区二区三区四区免费观看 | 亚洲国产欧美一区二区综合| 在线播放国产精品三级| 欧美xxxx黑人xx丫x性爽| 精品乱码久久久久久99久播| 国产精品女同一区二区软件 | 亚洲第一欧美日韩一区二区三区| 黑人操中国人逼视频| 男人舔奶头视频| 久久中文字幕一级| 又爽又黄无遮挡网站| 丁香六月欧美| 三级男女做爰猛烈吃奶摸视频| 国产精品永久免费网站| 老汉色∧v一级毛片| 床上黄色一级片| 久久久久国内视频| 亚洲在线自拍视频| 国产高潮美女av| 国产精品免费一区二区三区在线| 亚洲欧美日韩卡通动漫| 老司机福利观看| 欧美不卡视频在线免费观看| 欧美在线黄色| 亚洲av熟女| 国产又黄又爽又无遮挡在线| 在线观看美女被高潮喷水网站 | 久久人人精品亚洲av| 哪里可以看免费的av片| 两性夫妻黄色片| 国产亚洲精品久久久com| 19禁男女啪啪无遮挡网站| 亚洲精华国产精华精| 久久久久国产一级毛片高清牌| 欧美日本视频| 免费搜索国产男女视频| 婷婷六月久久综合丁香| cao死你这个sao货| 丰满人妻熟妇乱又伦精品不卡| 午夜两性在线视频| 成人特级av手机在线观看| 桃色一区二区三区在线观看| 国产精品乱码一区二三区的特点| 午夜a级毛片| 午夜福利在线观看吧| 午夜激情福利司机影院| 日本a在线网址| 国产成人啪精品午夜网站| 99久久久亚洲精品蜜臀av| 男女之事视频高清在线观看| 高清在线国产一区| 午夜福利高清视频| 久久中文字幕人妻熟女| 国产高清videossex| 国产一级毛片七仙女欲春2| 午夜免费观看网址| 午夜精品在线福利| 国产精品电影一区二区三区| 亚洲成人久久性| 国产黄a三级三级三级人| 99久久综合精品五月天人人| 亚洲国产高清在线一区二区三| 免费看美女性在线毛片视频| 国产精品一区二区免费欧美| 午夜影院日韩av| 99riav亚洲国产免费| 久久精品国产综合久久久| 在线观看一区二区三区| 免费观看的影片在线观看| 久久久久久久久中文| 少妇丰满av| 日韩有码中文字幕| 国产精品一区二区三区四区免费观看 | 欧美3d第一页| 一个人看的www免费观看视频| 岛国视频午夜一区免费看| 成人av在线播放网站| 9191精品国产免费久久| 亚洲成人久久性| 中文字幕熟女人妻在线| 国内久久婷婷六月综合欲色啪| 三级国产精品欧美在线观看 | 久久久精品大字幕| 亚洲成人久久性| 日本精品一区二区三区蜜桃| 日韩高清综合在线| 久久午夜综合久久蜜桃| 国产精品99久久99久久久不卡| 黄片小视频在线播放| 三级毛片av免费| 亚洲自偷自拍图片 自拍| 1024香蕉在线观看| 成年人黄色毛片网站| 国产真实乱freesex| 亚洲精品中文字幕一二三四区| 最近最新中文字幕大全电影3| 日韩欧美 国产精品| 国产黄色小视频在线观看| 国产在线精品亚洲第一网站| 97碰自拍视频| 99国产极品粉嫩在线观看| 日韩精品青青久久久久久| 51午夜福利影视在线观看| 成人18禁在线播放| 三级国产精品欧美在线观看 | 天堂影院成人在线观看| 日韩免费av在线播放| 日本 欧美在线| 亚洲欧美日韩高清在线视频| www日本黄色视频网| 可以在线观看的亚洲视频| 成人国产综合亚洲| 亚洲av电影在线进入| 久久久久精品国产欧美久久久| 国产男靠女视频免费网站| 丝袜人妻中文字幕| 黑人操中国人逼视频| 国产精品久久久久久精品电影| 最近在线观看免费完整版| 欧美精品啪啪一区二区三区| 天堂动漫精品| 成人鲁丝片一二三区免费| 亚洲专区国产一区二区| 久久中文字幕人妻熟女| 久久国产精品人妻蜜桃| 国产精品影院久久| 欧美黑人巨大hd| 美女被艹到高潮喷水动态| 亚洲欧美日韩无卡精品| 久久久久久久精品吃奶| 中出人妻视频一区二区| 91九色精品人成在线观看| 亚洲国产日韩欧美精品在线观看 | 桃红色精品国产亚洲av| 国产精品久久电影中文字幕| 亚洲中文日韩欧美视频| 午夜福利欧美成人| 99在线视频只有这里精品首页| 亚洲无线观看免费| 国产不卡一卡二| 男插女下体视频免费在线播放| 亚洲国产欧美网| 美女黄网站色视频| 91麻豆av在线| 三级国产精品欧美在线观看 | 看黄色毛片网站| 黄色女人牲交| 成年免费大片在线观看| 成人特级黄色片久久久久久久| 中文字幕熟女人妻在线| 淫妇啪啪啪对白视频| 在线观看一区二区三区| 国产69精品久久久久777片 | 国产午夜精品久久久久久| 国产亚洲精品久久久com| 岛国在线观看网站| 久久久精品欧美日韩精品| 一本精品99久久精品77| 亚洲乱码一区二区免费版| 一个人观看的视频www高清免费观看 | 成人三级做爰电影| 国产精品,欧美在线| 久久久久久久久中文| 久久香蕉国产精品| 男女视频在线观看网站免费| 欧美成狂野欧美在线观看| 国产97色在线日韩免费| 国产伦精品一区二区三区四那| 国产熟女xx| 三级国产精品欧美在线观看 | 18禁国产床啪视频网站| 午夜亚洲福利在线播放| 午夜激情欧美在线| 亚洲 欧美一区二区三区| 精品不卡国产一区二区三区| 国产不卡一卡二| 黄频高清免费视频| 日本撒尿小便嘘嘘汇集6| 中文在线观看免费www的网站| 国产毛片a区久久久久| 两个人的视频大全免费| 久久精品综合一区二区三区| 欧美丝袜亚洲另类 | 五月伊人婷婷丁香| 99精品久久久久人妻精品| 网址你懂的国产日韩在线| 噜噜噜噜噜久久久久久91| 观看免费一级毛片| 中文字幕人成人乱码亚洲影| av女优亚洲男人天堂 | 男女视频在线观看网站免费| 亚洲国产日韩欧美精品在线观看 | or卡值多少钱| 老司机福利观看| 97碰自拍视频| 久久性视频一级片| 欧美不卡视频在线免费观看| 长腿黑丝高跟| 色综合欧美亚洲国产小说| 国产麻豆成人av免费视频| 少妇熟女aⅴ在线视频| 欧美午夜高清在线| 女人被狂操c到高潮| 一进一出抽搐gif免费好疼| 天天躁日日操中文字幕| 成人精品一区二区免费| 日韩欧美一区二区三区在线观看| 美女高潮的动态| 色哟哟哟哟哟哟| 亚洲五月天丁香| 欧美日韩国产亚洲二区| 床上黄色一级片| 制服人妻中文乱码| 午夜福利欧美成人| 欧美色欧美亚洲另类二区| 日本黄大片高清| 最近最新免费中文字幕在线| 无限看片的www在线观看| 亚洲中文av在线| 午夜福利18| 男人和女人高潮做爰伦理| 一个人看的www免费观看视频| 久久久久久国产a免费观看| 国产成人av激情在线播放| 久久这里只有精品中国| 亚洲欧美激情综合另类| 午夜福利高清视频| 毛片女人毛片| 亚洲熟女毛片儿| 色播亚洲综合网| 国产午夜精品论理片| 亚洲中文av在线| 九色国产91popny在线| 一区福利在线观看| 午夜精品久久久久久毛片777| 亚洲午夜理论影院| 啦啦啦免费观看视频1| 18禁裸乳无遮挡免费网站照片| av女优亚洲男人天堂 | 久久亚洲精品不卡| 国产成人影院久久av| 免费观看精品视频网站| 人人妻,人人澡人人爽秒播| 亚洲九九香蕉| www国产在线视频色| 国产视频一区二区在线看| www.999成人在线观看| 女警被强在线播放| 90打野战视频偷拍视频| 女人被狂操c到高潮| 久久久久亚洲av毛片大全| tocl精华| 在线十欧美十亚洲十日本专区| 女生性感内裤真人,穿戴方法视频| 中文字幕最新亚洲高清| 亚洲成人精品中文字幕电影| 成人性生交大片免费视频hd| 精品一区二区三区视频在线观看免费| a级毛片在线看网站| 欧美日韩瑟瑟在线播放| 听说在线观看完整版免费高清| 啦啦啦韩国在线观看视频| 国产精品,欧美在线| www.熟女人妻精品国产| 少妇人妻一区二区三区视频| 熟女人妻精品中文字幕| 午夜两性在线视频| 亚洲aⅴ乱码一区二区在线播放| 欧美绝顶高潮抽搐喷水| 久久中文字幕人妻熟女| 中文资源天堂在线| 亚洲av日韩精品久久久久久密| 在线播放国产精品三级| 午夜影院日韩av| 国产成人一区二区三区免费视频网站| 身体一侧抽搐| 美女cb高潮喷水在线观看 | 日本 欧美在线| 国产成人精品久久二区二区免费| 亚洲精品美女久久久久99蜜臀| 最新美女视频免费是黄的| 国产主播在线观看一区二区| 国产97色在线日韩免费| 日韩免费av在线播放| 性色av乱码一区二区三区2| 久久久久久久久免费视频了| 日韩欧美精品v在线| 麻豆成人午夜福利视频| 国产黄片美女视频| 俄罗斯特黄特色一大片| 亚洲中文日韩欧美视频| 久久精品影院6| 久久热在线av| 精品乱码久久久久久99久播| 最新美女视频免费是黄的| 麻豆国产97在线/欧美| 国产午夜精品论理片| 久久久久九九精品影院| 久久久久亚洲av毛片大全| www.熟女人妻精品国产| 嫩草影院精品99| 又黄又爽又免费观看的视频| 1024手机看黄色片| 久久国产乱子伦精品免费另类| 男女那种视频在线观看| 搡老妇女老女人老熟妇| 精品国产亚洲在线| www国产在线视频色| 美女高潮喷水抽搐中文字幕| 香蕉久久夜色| 18禁裸乳无遮挡免费网站照片| 亚洲精品一区av在线观看| 日本三级黄在线观看| 亚洲无线在线观看| 小说图片视频综合网站| 国产亚洲欧美98| 欧美午夜高清在线| 国产午夜精品论理片| 美女cb高潮喷水在线观看 | www.熟女人妻精品国产| 1024香蕉在线观看| 日韩大尺度精品在线看网址| 国产一区二区三区在线臀色熟女| 老司机在亚洲福利影院| 一本久久中文字幕| 最好的美女福利视频网| 丰满人妻熟妇乱又伦精品不卡| 亚洲av日韩精品久久久久久密| 免费av不卡在线播放| 老司机午夜十八禁免费视频| 啦啦啦免费观看视频1| 免费在线观看日本一区| 岛国在线免费视频观看| 国产极品精品免费视频能看的| 很黄的视频免费| 99久久精品热视频| 国产aⅴ精品一区二区三区波| 天天添夜夜摸| 精品国产超薄肉色丝袜足j| 国产欧美日韩一区二区精品| 国产极品精品免费视频能看的| 在线观看日韩欧美| 人妻夜夜爽99麻豆av| 丝袜人妻中文字幕| 白带黄色成豆腐渣| 久久久久久国产a免费观看| 两性午夜刺激爽爽歪歪视频在线观看| 成年女人永久免费观看视频| 99久久精品一区二区三区| 亚洲av电影在线进入| 在线国产一区二区在线| 久久久国产精品麻豆| 国产精品亚洲美女久久久| 久久久久久大精品| 午夜两性在线视频| 中出人妻视频一区二区| 天堂√8在线中文| 国产精品一区二区三区四区免费观看 | 国产高潮美女av| 亚洲avbb在线观看| 欧美性猛交黑人性爽| 日本成人三级电影网站| e午夜精品久久久久久久| 97超视频在线观看视频| 不卡av一区二区三区| 一进一出抽搐动态| 在线永久观看黄色视频| 欧美日本视频| 两性夫妻黄色片| 99re在线观看精品视频| 看黄色毛片网站| 日韩高清综合在线| 2021天堂中文幕一二区在线观| 亚洲人成网站高清观看| а√天堂www在线а√下载| www.自偷自拍.com| 久久午夜综合久久蜜桃| 亚洲av熟女|