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

    Zn(Ⅱ)金屬有機(jī)骨架作為基質(zhì)金屬蛋白酶模擬物催化水解微囊藻毒素-LR

    2023-12-21 01:02:10方艷芬賀玉婷牛慧斌陳春城黃應(yīng)平
    無機(jī)化學(xué)學(xué)報 2023年12期
    關(guān)鍵詞:微囊三峽庫區(qū)宜昌

    方艷芬 賀玉婷 趙 麗 ?;郾?陳春城 李 悅*, 黃應(yīng)平*,

    (1三峽大學(xué)材料與化工學(xué)院,宜昌 443002)

    (2三峽大學(xué)三峽庫區(qū)生態(tài)環(huán)境教育部工程研究中心,宜昌 443002)

    (3南開大學(xué)化學(xué)學(xué)院,天津 300071)

    (4中國科學(xué)院化學(xué)研究所,北京 100190)

    0 Introduction

    Microcystins (MCs) are a large type of cyclic heptapeptide toxins in eutrophic waters, which are commonly considered hepatotoxic and carcinogenic because of the strong inhibition of protein phosphatases and interference with cell signaling pathways. In some bodies of water that have become eutrophicated,the levels of MCs have increased to levels as high as dozens of parts per billion, which greatly exceeds the World Health Organization′s water quality standard of 1.0 parts per billion[1-4]. Advanced oxidation processes(AOPs)are capable of efficiently oxidizing MCs by producing reactive oxygen species (ROS). Nevertheless,due to the high presence of natural organic matter(NOM) in natural water bodies, AOPs are not selective to MCs because most ROS are consumed by NOM[5-7].Hydrolysis, on the other hand, is a gentler and safer way of treating pollutants,making it a preferable option over AOPs[8-11].Based on our earlier research,it was discovered that siderite, a mineral that contains iron, can break down microcystin-LR (MC-LR) peptide bonds,even in the presence of humic acid (HA), which is 10 times its amount[12]. Although siderite′s effectiveness is not very high, the results indicate that catalytic hydrolysis is a feasible approach to achieve the desired degradation of MCs when NOM is present.

    Matrix metalloproteinases (MMPs), a group of proteinases dependent on Zn(Ⅱ), play a significant role in breaking down polypeptide bonds in connective tissues[13-14]. This breakdown can potentially aid in the hydrolysis of peptide bonds in MCs. Metal-organic frameworks(MOFs)hold great potential as a foundation for creating enzyme mimics or nanozymes. This is due to their remarkable versatility in composition, which enables the modification of metal ion identity and coordination environment to replicate enzymatic active sites. The functions of a series of metalloenzymes, such as tyrosinase[15], CO-dehydrogenase[16], and phosphotriesterase[17], have been performed by MOF mimics.Therefore, MOFs having functional groups can be organized in a structured manner, resulting in synergistic action, similar to MMPs. However, there are currently no artificial MMP mimics reported for MC hydrolysis,despite the widespread use of MMPs.

    In this work, the physical and chemical properties of the Zn (Ⅱ)-based MOF nanomaterials (Zn-MOF-1-NSs) were meticulously analyzed and successfully used to hydroxylate MC-LR.The Zn-MOF-1-NS nanomaterial, acting as an MMP mimic, showed remarkable efficiency in catalyzing MC-LR hydroxylation with high activity and substrate specificity. Further investigation revealed that both the Zn(Ⅱ)ion and the carboxyl group of the MOF were involved in activating the peptide bond.

    1 Experimental

    1.1 Chemical and materials

    MC-LR (Fig.S1, Supporting information) standard(1 mg solid, purity not less than 95% by HPLC) was purchased from Express Technology Co., Ltd. and stored at -25 ℃. MC-LR stock solution was prepared by adding H2O to dissolve the MC-LR standard and then stored at 4 ℃. The concentration of the MC-LR solution used inin-situATR-FTIR experiments was 250 mg·L-1and that in the degradation experiments was 2 mg·L-1. Deuterium oxide (D2O, 99.9%) was purchased from Cambridge Isotope Laboratories, Inc.Glutamic acid (D-Glu), arginine (L-Arg), alanine (DAla) and leucine (L-Leu) were obtained from Sigma-Aldrich. River humic acid standard (IR105H) was purchased from the International Humic Substances Society (IHSS). Methanol and acetonitrile were HPLC grade, and water was purified by reverse osmosis and deionization. All chemicals were obtained commercially and used as received without further purification.

    1.2 Preparation of Zn-MOF-1-NS

    Zn-MOF-1 precursor was synthesized by the following procedures[18]: a mixture of ZnSO4·6H2O (143.8 mg, 0.5 mmol), 3-amino-1,2,4-triazole-5-carboxylic acid (38.5 mg, 0.3 mmol), and 1,3,5-benzenetricarboxylic acid (42 mg, 0.2 mmol) was dissolved in water(10.0 mL), and placed in a 23 mL Teflon-lined stainless steel vessel.The pH of the solution was adjusted to 6.0 by triethylamine.This mixture was sealed and heated at 150 ℃for 48 h.After cooling to 25 ℃,the precipitate of the crude product was collected, washed with H2O and methanol, and dried in a vacuum at 60 ℃for 6 h. Anal. Calcd. for Zn-MOF-1 (C22H18N8O15Zn4)(%):C, 29.49; H, 2.02; N, 12.51. Found(%): C, 29.17; H,2.08;N,12.65.

    The exfoliation of Zn-MOF-1 was carried out by sonicating the ethanol suspension of its precursor (8 mg in 7 mL)for 3 h.After the centrifugation at 4 000 r·min-1for 10 min to remove unexfoliated particles,Zn-MOF-1-NSs were separated from the ethanol supernatant by rotatory evaporation.

    1.3 Characterization

    Powder X-ray diffraction (PXRD) patterns were recorded on a Rigaku D/max-2500 diffractometer with CuKαradiation (λ=0.154 06 nm) at 40 kV and 100 mA with a scanning range of 5°-50°.Scanning electron microscopy (SEM) images were taken using a JEOL JSM-7500F scanning electron microscope at an accelerating voltage of 30 kV.Elemental analysis (C, H, and N) was performed on a Perkin-Elmer 240C analyzer.Thermogravimetric analysis (TGA) was operated with a Rigaku standard TG-DTA analyzer from ambient temperature to 700 ℃with a heating rate of 10 ℃·min-1in the air, and an empty Al2O3crucible was used as the reference. Atomic force microscopy (AFM) was carried out on a Bruker Icon probe microscope in tapping mode in the air under ambient conditions, and the delaminated MOF samples were deposited on aluminum-coated silicon cantilevers. Theζpotential of the catalyst was measured using a Malvern Zeta Sizer. N2sorption isotherms were measured with a Micromeritics ASAP 2020 analyzer at 77 K. X-ray photoelectron spectroscopy (XPS) was conducted on an AXIS Supra spectrometer equipped with 300 W AlKαradiation.The hydrocarbon C1sline at 284.8 eV from adventitious carbon was used for energy referencing.A Nicolet iS50 Fourier transform infrared spectrometer collected infrared spectra(FTIR).

    1.4 Hydrolysis experiments

    Reactions were carried out in a 10 mL vial placed in a dark box with magnetic stirring at (25±0.1) ℃. In each experiment, 3 mL aqueous solution of MC-LR with the initial concentration of 2 mg·L-1was placed in the vial, and 5 mg catalyst was dispersed to it. The pH of the suspension was not adjusted. At specified reaction times, 400 μL samples of the reaction mixture were collected, centrifuged to remove the catalyst, and analyzed using HPLC.

    The substrate specificity of Zn-MOF-1-NS toward MC-LR in the presence of NOM was tested with the addition of HA with 1, 2, 10, and 20 mg·L-1. Other reaction and analytical conditions were unchanged.

    1.5 Analytical methods

    MC-LR. A Waters 2690 HPLC determined the concentration of MC-LR with a PDA detector and a Kromasil C18 column (4.6 mm×250 mm, 10 μm particle size). The mobile phase was water and methanol(60∶40,V/V) containing 0.05% trifluoroacetic acid(TFA). The flow rate, injection volume, and column temperature were 0.8 mL·min-1, 20 μL, and 35 ℃,respectively. MC-LR was quantified by its absorption at 238 nm. The retention time of MC-LR was 8.3 min.The quantification range was from 0 to 10 μmol·L-1for MC-LR with the correlation coefficient(R2)for the standard of 0.998(Fig.S2a).

    Amino acids. Before measurement, 0.6 mL of 4∶1(V/V) water-ethanol mixture was added to wash the adsorbed amino acids from the MOF catalyst (Fig.S3)[19].The details about desorption experiments are provided in the Supporting information. The hydrolysis products(D-Glu,L-Arg,D-Ala, andL-Leu) of MC-LR were determined by a Waters 2690 HPLC with an FLR detector and a Nova-PakTMC18 column (3.9 mm×150 mm, 4 μm particle size). The retention times ofD-Glu,L-Arg,D-Ala, andL-Leu were 8.3, 14.3, 19.8, 21.3,and 31.8 min, respectively. The hydrolysis products were quantified using fluorescent detection withλex=330 nm andλem=445 nm. These compounds were identified and quantified by the external standard curves of the peak area. The quantification ranges for four amino acids were from 0.5 to 20 μmol·L-1with theR2for the standards above 0.99(Fig.S2b).

    1.6 In-situ ATR-FTIR analysis

    In-situinfrared spectroscopy was used to monitor the surface interaction between Zn-MOF-1-NS film and MC-LR in different solvents (H2O,D2O,and C6H6).The ATR cell was placed in the internal compartment of a Nicolet iS50 FTIR spectrometer (Thermo, USA)equipped with an MCT-A detector.Thin Zn-MOF-1-NS film was prepared by drying Zn-MOF-1-NS suspension(0.4 mL, 10 g·L-1) on a ZnSe crystal strip (80 mm×10 mm×4 mm). After adding 0.8 mL of solvent to the flow cell, ATR-FTIR spectra were collected continuously for 2 or 4 h until the establishment of the sorption equilibrium between solvent and Zn-MOF-1-NS film. After MC-LR (100 μL, 250 mg·L-1) was added, the timeresolved ATR-FTIR spectra were recorded at the interval of 0.283 s.

    1.7 Theoretical calculation

    The geometry optimization, frequency calculation,and natural bond orbital (NBO) analysis were all performed at the M06/def2-TZVP level[20-21]. To make the calculation affordable, besides the central Zn2+ion,only the ligands coordinating directly with it were considered in treating Zn3 and Zn4. For Zn2, because of the possibility of H2O/reactant forming an H bond with nearby carboxyl groups, two additional btc3-ligands were taken into the model. The initial configurations of Zn2+and ligands were picked up from the crystalline structure of Zn-MOF-1. Since the peripheral carboxyl groups of btc3-ligands are far from the reaction center of Zn2+, they were removed with H atoms supplemented to the vacancies. During optimization, the C atoms in the benzene rings and the N atoms coordinating to peripheral Zn2+ions (not included in the model) were fixed to replicate the bulk behavior, while the other atoms were relaxed (Fig.1). An implicit solvent model of SMD (solvation model density) modulated the solvent effects.

    Fig.1 Coordination environments of(a)Zn1,(b)Zn2,(c)Zn3,and(d)Zn4 in Zn-MOF-1-NS

    2 Results and discussion

    2.1 Selection of MOF catalyst

    Due to the large molecular size of MC-LR (2.4 nm), it is difficult to design channels large enough to accommodate MC-LR molecules without affecting the stability of the MOF. The results reported by Tissot et al.[22]and our previous works[23]have shown the inefficiency of internal active sites towards large-molecular reactants. Therefore, we chose 2D nanosheets, which expose the maximum quantity of active sites on the external surface, as the morphology for the MOF catalyst.Combining with the structure and catalytic mechanism of MMP,the selection criteria of the MOF catalyst were: (1) the backbone is constructed by the connection among Zn2+ions and five-membered N-heterocyclic bridging ligands, and carboxyl groups are contained for the activation of peptide bond;(2)in addition to heterocyclic N atoms, Zn2+ion is coordinated with at least one solvent molecule, which would be easily replaced by substrate and does not destroy the integrity of the MOF; (3) the framework is 2D connected, which could be delaminated to nanosheets. After the screening of the 424 MOFs containing Zn2+and N-pentacyclic ligand containing MOFs summarized in the CCDC database, [Zn4(atz)2(btc)2(H2O)3]n(Zn-MOF-1), a MOF assembled by the ligands of 3-amino-1,2,4-triazolate(atz-) and 1,3,5-benzenetricarboxylate (btc3-), was chosen as a potential catalyst for the following investigation (Fig.2)[18]. This MOF contains four crystallographically independent Zn(Ⅱ)centers (Fig.1). The benzene ring and triazole ring projections on the ab plane are staggered, but there isn′t any effective overlap. Therefore,theπ-πstacking doesn′t play a significant role in the crystal stability of this MOF.All of these Zn(Ⅱ)centers are surrounded by triazole N atoms and carboxyl O atoms, and three of them (Zn2, Zn3, and Zn4) are coordinated by one H2O molecule. The interconnection among Zn2+ions and two kinds of ligands form a double-layer structure, which expands into the supramolecular framework of the crystal only by the weak connection of H bonds. The 2D connected structure of Zn-MOF-1 was drawn with Mercury 4.2.0.

    Fig.2 Two-dimensional connected structure of Zn-MOF-1

    2.2 Preparation and characterization of MOF nanosheets

    Zn-MOF-1 was synthesized generally according to the reported procedure with minor modifications to minimize the particle size[18]. PXRD was applied to characterize the obtained sample (Fig. 3a). All the recorded diffraction peaks could be readily indexed to the simulated pattern, showing that our sample presents the same framework structure and has good phase purity. This backbone was further supported by the results of elemental analysis, thermogravimetric (Fig.S4),and FTIR analyses(Fig.S5).Because of its 2D connected framework, the sample of this MOF presented a layer-stacking morphology under SEM observation(Fig.3b). In consistent with its nonporous structure,Zn-MOF-1 presented a small BET surface area of 8 m2·g-1(Fig.S6a and Table S1). This result implied that Zn-MOF-1 has limited contact with the reactants,so its delamination is necessary.

    Fig.3 (a)XRD patterns and(b)SEM image of Zn-MOF-1

    The Zn-MOF-1 precursor was exfoliated simply by the sonication in ethanol (the product is denoted as Zn-MOF-1-NS).After the removal of sediments,a transparent supernatant was obtained. The exfoliation process was evidenced by an apparent Tyndall effect under the irradiation of a laser beam (Fig.4a).The morphology of the exfoliated sample was examined by AFM.As shown in Fig.4b and 4c,this sample was composed of irregularly shaped nanosheets with several hundred nanometers lateral dimensions.The thicknesses of these nanosheets were generally about 4 nm,indicating that each nanosheet comprises only six two-layer structures of Zn-MOF-1. Zn-MOF-1-NS presented a negativeζpotential of -9.2 mV at pH=7.4 (Fig.S7),which ensured its good stability in the aqueous suspension[24].

    Fig.4 (a)Photograph of the Tyndall effect of the ethanol suspension of Zn-MOF-1-NS;(b)AFM image of Zn-MOF-1-NS;(c)Height profile of the corresponding section of figure b

    2.3 Catalytic degradation of MC-LR

    With 5 mg of Zn-MOF-1-NS, 82.6% of MC-LR was degraded in 7.5 h(Fig.5a).MC-LR decay generally followed first-order kinetics withR2=0.94 (Fig. 5b).Within 7.5 h, only 20.9% of MC-LR was degraded by siderite, indicating that the activity of Zn-MOF-1-NS was much higher than those of siderite that showed the highest activity reported for the hydrolysis of MC-LR[12].To further examine the reuse potential of Zn-MOF-1-NS, cyclic degradation experiments were carried out(Fig.5d).It was observed that Zn-MOF-1-NS continued to display high catalytic activity in five successive cycles. However, due to the high dispersity of nanosheets,part of the catalyst was lost during the recycling process, which induced the slowdown of the reaction(by 45% in five cycles).Furthermore,as evidenced by the PXRD pattern (Fig.S8),the Zn-MOF-1-NS catalyst maintained good stability during the reaction.

    Fig.5 (a)Degradation curves and(b)reaction kinetics of MC-LR using various catalysts;(c)Hydrolysis yield production of amino acids using Zn-MOF-1-NS(Inset:MC-LR degradation curve);(d)Recycling of Zn-MOF-1-NS;(e)Michaelis-Menten kinetic curve of the initial reaction rate against MC-LR concentration;(f)Degradation curves of MC-LR(2 mg·L-1)at different concentrations of HA in Zn-MOF-1-NS(Inset:effect of HA concentration on the reaction rate constant of MC-LR hydrolysis)

    To illustrate the reaction mechanism, another 2D Zn(Ⅱ)based MOF,[Zn7(OH)2(Hbta)6(H2O)6]n(H3bta=N,Nbis-(1(2)H-tetrazol-5-yl)-amine) (named Zn-MOF-2)was also tested for comparison[25]. The characterization of Zn-MOF-2 is also shown in Fig.S9-S11. Under Zn-MOF-2-NS catalysis, the degradation efficiency of MC-LR in 7.5 h could only reach 54.0%,obviously lower than the result of Zn-MOF-1-NS. The Zn2+density could not explain the superior activity of Zn-MOF-1-NS because the surficial Zn2+density of Zn-MOF-1-NS(1.67×10-5nm-2)was even lower than that of Zn-MOF-2-NS (2.31×10-5nm-2). Considering that the carboxyl group was only contained in Zn-MOF-1-NS, we speculated that this group participates in the hydrolysis of MC-LR. However, the addition of acetic acid to the Zn-MOF-2-NS catalytic system could not remarkably accelerate MC - LR consumption. This phenomenon indicated that the relative positions of Zn2+and the carboxyl group are crucial.

    The intermediates produced during the reaction were tested with an initial concentration of MC-LR of 10 mg·L-1.D-Glu,L-Arg,D-Ala, andL-Leu were detected during the decomposition of MC-LR (Fig.5c).At 10 h, the concentrations were 0.653, 0.466, 1.20 and 0.0183 μmol·L-1, respectively. The quantity ofD-Ala could account for 50.0% of MC-LR consumption, confirming that the hydrolysis of peptide bonds was the main pathway for MC-LR degradation.

    A steady-state kinetic assay was performed by varying the initial concentration of MC-LR. The relationship between the initial reaction rate and substrate concentration could be well described by the Michaelis-Menten equation withR2=0.97, showing the enzyme mimic feature of Zn-MOF-1-NS (Fig.5e). The fitting gave the Michaelis constant (Km) of 6.2 μmol·L-1and the maximum reaction rate (Vmax) of 0.61 μmol·L-1·min-1. SinceKmis the substrate concentration at which half the enzymes are bound and its values for MMPs range from tenths to tens of millimoles[26-29], this result indicated that the affinity of Zn-MOF-1-NS for MC-LR is comparable with those of MMPs towards proteins.

    Various NOMs are contained in environmental water samples and could potentially inhibit MC-LR degradation. HA was added to the reaction system to test the ability to resist NOM interference of Zn-MOF-1-NS. Even when the HA concentration was increased to 20 mg·L-1, the rate change of MC-LR (2 mg·L-1) was still negligible (Fig.5f), showing that Zn-MOF-1-NS can specifically degrade MC-LR in the presence of NOMs. In contrast, with Zn-MOF-2-NS as a catalyst,the same concentration of HA led to a significant slowdown of the reaction (the percentage of degradation decreased from 46.7% to 29.5%, Fig.S12). Based on the structure of Zn-MOF-1-NS, this catalyst was expected to catalyze the hydroxylation of the peptide bond by the Lewis acidity of Zn2+ions. However, our previous research illustrated that a single Lewis acid site of mineral catalysts did not display specificity for MC-LR[12].Thus, the phenomenon in the Zn-MOF-1-NS catalytic system indicates the existence of multiple binding processes in the reaction.

    2.4 Hydrolysis of MC-LR on Zn-MOF-1-NS surface

    In the ATR-FTIR spectrum of Zn-MOF-1-NS(Fig.S13a),the wide band at 3 260 cm-1includes O—H and N—H stretching vibrations[30]. The 1 615 and 1 431 cm-1peaks belong to the C=O and C—O stretching vibrations of COO-that coordinate with Zn2+, respectively[12,24]. The vibration at 1 554 cm-1is attributed to N—H[31]. The band at 1 369 cm-1belongs to the C—C bonds of the benzene ring[32]. The absorptions at 1 211 and 1 062 cm-1correspond to the C—N bonds on and in the triazole rings, respectively[33]. The stretching vibration of Zn—O appeared at 734 cm-1[34]. The spectrum of MC-LR is given in Fig.S13b. In this spectrum,the absorption bands at 3 260 and 1 515 cm-1refer to the stretching and bending vibrations of N—H[12]. The bands at 2 925 and 2 851 cm-1are assigned to the stretching vibrations of —CH2—[35]. The C=O stretching vibration of COOH was centered at 1 731 cm-1[36].The peak at 1 642 cm-1(Ⅰband) and 1 450 cm-1(Ⅱband)was ascribed to the amide C=O stretching vibration[37]. The peak at 1 642 cm-1is ascribed to the amide C=O stretching vibration (Ⅰband). The vibrations of C—H in-plane bending, C—O—C stretching, and C—N stretching correspond to the peaks at 1 409, 1 178,and 1 082 cm-1,respectively[25].

    Fig.6a and 6b show thein-situATR-FTIR spectra of the Zn-MOF-1-NS/H2O equilibrium process, which were recorded at intervals after MC-LR addition (100 μL, 250 mg·L-1). The positive bands at 3 501, 3 409,1 648, 1 537, 1 487, 1 401, and 1 087 cm-1showed that MC-LR is adsorbed on the surface of Zn-MOF-1-NS.A redshift of the peak of the amide C=O of MC-LR from 1 648 to 1 637 cm-1was observed,which could be explained by its coordination with the Zn2+of Zn-MOF-1-NS. Because Zn2+works as a Lewis acid site, such coordination would reduce the π electron density of the amide group and thus weaken the C=O bond[38-39]. In addition,the negative peak at 1 615 cm-1indicated that the—COOH of Zn-MOF-1-NS is involved in the interaction with MC-LR in addition to Zn2+. Meanwhile, the new band at 3 183 cm-1illustrates the formation of a hydrogen bond between the amide N—H of MC-LR and the carboxyl oxygen of Zn-MOF-1-NS[27]. With the proceeding of the reaction, two new peaks of —COOH positively occurred at 1 601 and 1 664 cm-1[40-41], while two negative bands appeared at 1 216 and 1 065 cm-1,which are assigned to the cleavage of the amide C—N bond of MC-LR. These changes indicated that the peptide bond of MC-LR is broken into—COOH and N—H under the catalysis of Zn-MOF-1-NS.

    To verify the hydrolysis mechanism,analogous experiments in D2O were performed.In the spectra of Zn-MOF-1-NS, the bands of O—H at 3 297 and 1 644 cm-1were red-shifted to 2 470 and 1 201cm-1for the O—D bond (Fig.S14), which conformed to the Hooke′s law[42]. Upon the addition of MC-LR, the same trend as that in the Zn-MOF-1-NS/H2O system was observed(Fig.6c). A crest appeared at 1 664 cm-1, showing the formation of —COOH by the cleavage of amide C—N[43-44]. The peak of N—H bending vibration at 1 330 cm-1and the valley of the break of the C—N bond at 1 047 cm-1confirm the cleavage of the peptide bond of MC-LR (Fig.6d). The redshifts of the —COO-peak of Zn-MOF-1-NS from 1 431 cm-1(in the pure Zn-MOF-1-NS/D2O system, Fig.S13a) to 1 450 cm-1and the N—H peak of MC-LR from 1 554 to 1 330 cm-1confirm the formation of hydrogen bond between them. Together with the results of the H2O system,these changes manifest that both Zn2+and—COO-of Zn-MOF-1-NS participate in the interaction with MC-LR.The former acts as a Lewis acid site to coordinate with amide C=O and the latter binds to N—H through a hydrogen bond.

    Fig.6 In-situ ATR-FTIR spectra of(a,b)MC-LR(100 μL,250 mg·L-1)in H2O system and(c,d)MC-LR(100 μL,250 mg·L-1)in D2O system

    In addition, the infrared spectra of Zn-MOF-1-NS separately immersed in H2O and D2O were recorded for comparison (Fig.S15). Relative to the system of H2O,two new peaks were observed at 2 488 and 1 246 cm-1in D2O, indicating that D2O replaced the coordinating H2O molecule. Therefore, it is necessary to study the fate of MC-LR after binding with Zn2+sites in the D2O system.After the addition of MC-LR to the D2O system,thein-situATR-FTIR spectra presented a valley at 1 205 cm-1(O—D, for D2O desorption) and the crests at 3 400 cm-1(N—H of adsorbed MC-LR) and 1 450 cm-1(C=O of amide Ⅱ) from 0.2 to 119.5 min(Fig.6d), indicating that the D2O coordinating to Zn2+is replaced by MC-LR.

    Similar phenomena were also observed in the benzene system. Zn-MOF-1-NS equilibrated at C6H6for 2 h (Fig.7a). The loss of the coordinating H2O molecule in Zn-MOF-1-NS gave negative bands of H—O stretching at 3 386 and 1 632 cm-1(Fig.7b).At the same time,it was observed that the Zn—O band presented a blue shift from 774 to 804 cm-1. This change confirms the replacement of the H2O molecule by MC-LR on the Zn2+site.

    Fig.7 In-situ ATR-FTIR spectra of(a)Zn-MOF-1-NS equilibrated for 2 h and(b)MC-LR(100 μL,250 mg·L-1)/Zn-MOF-1-NS in C6H6 system(Solvent-Zn-MOF-1-NS spectrum as background);XPS spectra of Zn-MOF-1-NS before and after the reaction with MC-LR for 10 h:(c)O1s;(d)Zn2p

    We further used XPS spectroscopy to investigate the interaction between MC-LR and Zn-MOF-1-NS.For pristine Zn-MOF-1-NS, three peaks of O1swere recorded at 531.44, 532.08, and 532.52 eV, corresponding to Zn—O[45], O—H[46], and H2O groups[47],respectively (Fig.7c). After the addition of MC-LR, the peak of H2O shifted to 533.49 eV and its area decreased(from 36.1% to 16.0%), verifying the replacement of H2O molecule by MC-LR. Because of the formation of hydrogen bonds with MC-LR,the peak of the O atom of Zn—O presented a blue shift to 532.28 eV. Considering that Zn2+works as the Lewis acid site for MC-LR,we also examined the Zn2pspectrum(Fig.7d).The pristine MOF gave two peaks at 1 022.4 and 1 045.5 eV,which could be assigned to the Zn2p3/2and Zn2p1/2spinorbit components of O/N coordinated Zn atom, respectively[48-49]. Since replacing the coordinating aqueous O atom with the amide O atom of MC-LR does not change the first coordination sphere of Zn2+, the addition of MC-LR did not significantly alter the positions of these two peaks.

    The experimental observations indicate: (1) the much higher activity of Zn-MOF-1-NS than that of Zn-MOF-2-NS proved the participation of the carboxyl group; (2) the substrate specificity of Zn-MOF-1-NS towards MC-LR illustrates the formation of multiple bonds between them; (3) FTIR analysis shows that the amide O atom coordinates to Zn2+ion, and N—H is hydrogen bounded; (4)in-situATR-FTIR and XPS spectra indicate the coordination of amide O atom of MC-LR to Zn2+with the replacement of H2O molecule.So, we speculate that the catalytic mechanism of Zn-MOF-1-NS is similar to that of MMP and involves the participation of both Zn2+ions and carboxyl groups.Thus, we resorted to the theoretical calculation to give an atomic insight into the activation mechanism of MC-LR by Zn-MOF-1-NS. Since the reaction occurred mainly around the peptide bond, N-methylacetamide was taken as a simplified model for MC-LR. As introduced above, among the four crystallographically independent Zn(Ⅱ)centers in Zn-MOF-1-NS,three of them,Zn2, Zn3, and Zn4, are coordinated by H2O molecules and can provide coordination sites for amide O atom.Therefore, we first calculated the Gibbs free energy change (ΔG) values for the substitution of H2O molecule byN-methylacetamide on these metal sites, and the results are listed in Table 1. Due to the similar coordination environments of Zn3 and Zn4, they gave close ΔGvalues of -10.4 and -11.6 kJ·mol-1, respectively.The negative values of these free energy changes show that the coordination of amide to these Zn sites is thermodynamically feasible. The scrutinization of the geometries of these complexes shows that besides the coordination of the amide O atom to Zn2+(Zn—O:0.208 nm for Zn3 and 0.204 nm for Zn4), the N—H also forms a strong hydrogen bond with the carboxyl O atom in the ligand (O…H: 0.198 and 0.202 nm for Zn3 and Zn4, respectively, and ∠N—H…O=163.4° and 157.9°,Fig.S16).

    Table 1 Results of the calculation about amide adsorption to Zn sites of Zn-MOF-1-NS

    In contrast, Zn2 provided a positive ΔGfor the binding with amide (28.3 kJ·mol-1), which may result from the position of its coordinating H2O molecule.Unlike the coordinating H2O molecules of Zn3 and Zn4, which point to the outside of the double-layer structure, the coordinating H2O molecule of Zn2 is located on the opposite layer (Fig.1 and 2). This position makes the substituting amide group crowd with nearly btc3-ligands and thus weakens the formed Zn—O bond (0.226 nm). The steric hindrance also hinders the H bond formation with the carboxyl group coordinating directly with Zn2+. The amide N—H can only form a weak H bond with a carboxyl group in the second coordination sphere with longer O…H distance(0.214 nm) and narrower ∠N—H…O (141.8°). Therefore, the main active sites in Zn-MOF-1-NS are the Zn3 and Zn4 on the surface of the nanosheet, and their carboxyl ligands contribute to the binding with the amide.

    Since the first step of amide hydroxylation is the nucleophilic addition of H2O to the amide C atom, the bond order and polarity of C=O[50]and the LUMO level of the amide group[51]would be parameters to predict the reactivity of the substrate. It is observed that the C=O bonds in three amide-Zn complexes are all longer than those in free N-methylacetamide, and their total NBO polarities (Δt, the difference in the NBO charge between C and O atoms)also show a remarkable increase. The LUMO level of theN-methylacetamide presents aca. 0.5 eV decrease upon the coordination with Zn sites.These results prove that,after the adsorption to the catalyst,the amide C=O bond is weakened,its C atom becomes more positive-charged and electrophilic, and the antibonding orbital of this group becomes more feasible to accommodate the lone pair electrons of nucleophilic agent (H2O), and thus the reactivity is enhanced.

    3 Conclusions

    In summary, simulating the structures of MMPs, a Zn (Ⅱ)- based MOF with nanosheet morphology was developed to catalyze MC-LR hydrolysis. This catalyst can rapidly remove MC-LR even in the presence of a high concentration of HA. The comparison with Zn-MOF-2-NS showed that the carboxyl groups in Zn-MOF-1-NS contribute much to the superior activity and substrate specificity.In-situATR-FTIR,XPS,and theoretical calculation proved that multiple binding processes exist between the MC-LR and Zn-MOF-1-NS and contribute to the peptide bond activation. This work demonstrates the potential of mimic enzyme catalysis in environmental treatment and that MOF is an ideal platform to simulate the multi-group synergism in enzymes.

    Supporting information is available at http://www.wjhxxb.cn

    猜你喜歡
    微囊三峽庫區(qū)宜昌
    宜昌“清單之外無事項”等
    放歌宜昌(女聲獨唱)
    湖北宜昌卷
    三峽庫區(qū)萬家壩滑坡變形區(qū)穩(wěn)定性復(fù)核研究
    奇姝雙生 畫滿宜昌
    中華奇石(2016年9期)2016-12-15 14:20:46
    三峽庫區(qū)產(chǎn)業(yè)培育及結(jié)構(gòu)調(diào)整的思考
    微囊懸浮-懸浮劑和微囊懸浮劑不是同種劑型
    微囊藻毒素-LR對秀麗線蟲精子形成的毒性作用
    馬錢子生物堿微囊的制備及評價
    ELISA試劑盒法測定水中LR型微囊藻毒素
    日韩制服骚丝袜av| 妹子高潮喷水视频| 国产真人三级小视频在线观看| 国产亚洲欧美在线一区二区| av福利片在线| 91精品国产国语对白视频| 秋霞在线观看毛片| 人妻一区二区av| 欧美精品人与动牲交sv欧美| 欧美少妇被猛烈插入视频| 欧美性长视频在线观看| 日本五十路高清| 啦啦啦 在线观看视频| 国产精品久久久久久人妻精品电影 | 丝袜喷水一区| 国产高清不卡午夜福利| 少妇人妻 视频| xxxhd国产人妻xxx| 久久久久久久大尺度免费视频| 亚洲男人天堂网一区| 日韩人妻精品一区2区三区| 久久久久国产一级毛片高清牌| 丝袜美足系列| 波多野结衣av一区二区av| 国产男女超爽视频在线观看| 最新在线观看一区二区三区 | 午夜福利,免费看| 母亲3免费完整高清在线观看| 少妇的丰满在线观看| 午夜91福利影院| www.自偷自拍.com| 丝袜美腿诱惑在线| 又粗又硬又长又爽又黄的视频| 我要看黄色一级片免费的| 亚洲国产精品国产精品| 精品久久久久久久毛片微露脸 | 成年女人毛片免费观看观看9 | 国产精品.久久久| 如日韩欧美国产精品一区二区三区| 中文字幕最新亚洲高清| 又紧又爽又黄一区二区| 在线精品无人区一区二区三| 精品国产乱码久久久久久小说| 人人妻人人添人人爽欧美一区卜| 三上悠亚av全集在线观看| 久久鲁丝午夜福利片| 波多野结衣一区麻豆| 黄色a级毛片大全视频| 亚洲精品成人av观看孕妇| 午夜免费男女啪啪视频观看| 欧美成狂野欧美在线观看| 久久久久久久久免费视频了| 亚洲av综合色区一区| 欧美黄色片欧美黄色片| 18禁黄网站禁片午夜丰满| 国精品久久久久久国模美| 免费看av在线观看网站| www.av在线官网国产| 中文字幕人妻熟女乱码| 亚洲色图综合在线观看| 久久久久久人人人人人| 免费观看人在逋| 国产精品亚洲av一区麻豆| 中文字幕精品免费在线观看视频| 久久精品熟女亚洲av麻豆精品| 男女免费视频国产| 国产高清videossex| 曰老女人黄片| 欧美激情 高清一区二区三区| 又大又爽又粗| 国产av国产精品国产| 999久久久国产精品视频| 成人免费观看视频高清| 18禁裸乳无遮挡动漫免费视频| 中文字幕精品免费在线观看视频| 国产av一区二区精品久久| 亚洲激情五月婷婷啪啪| av在线播放精品| 亚洲 国产 在线| 人人妻人人澡人人爽人人夜夜| 又大又黄又爽视频免费| 老司机深夜福利视频在线观看 | 最近最新中文字幕大全免费视频 | 久久久久精品人妻al黑| 桃花免费在线播放| 亚洲欧美精品综合一区二区三区| 一本一本久久a久久精品综合妖精| 青草久久国产| 女人爽到高潮嗷嗷叫在线视频| 新久久久久国产一级毛片| 久久中文字幕一级| 成人影院久久| 国产精品久久久久久精品古装| 国产精品国产av在线观看| 亚洲人成电影观看| 国产成人免费无遮挡视频| 国产深夜福利视频在线观看| 欧美激情高清一区二区三区| 久久久久精品国产欧美久久久 | 在线亚洲精品国产二区图片欧美| 精品国产乱码久久久久久男人| 91老司机精品| 中文字幕制服av| 久久精品人人爽人人爽视色| bbb黄色大片| 亚洲国产欧美网| 免费黄频网站在线观看国产| 亚洲精品日韩在线中文字幕| 久久久久久久精品精品| av天堂在线播放| 欧美日韩亚洲综合一区二区三区_| 爱豆传媒免费全集在线观看| 中文乱码字字幕精品一区二区三区| 欧美在线黄色| 香蕉国产在线看| 国产1区2区3区精品| 女人高潮潮喷娇喘18禁视频| 男女之事视频高清在线观看 | 乱人伦中国视频| 少妇裸体淫交视频免费看高清 | 高清不卡的av网站| 免费少妇av软件| 麻豆av在线久日| 久久综合国产亚洲精品| 十八禁高潮呻吟视频| 亚洲av日韩在线播放| svipshipincom国产片| 亚洲欧美激情在线| 久久精品人人爽人人爽视色| 夫妻午夜视频| 纯流量卡能插随身wifi吗| 99精国产麻豆久久婷婷| 久久影院123| 每晚都被弄得嗷嗷叫到高潮| 性色av乱码一区二区三区2| 亚洲国产中文字幕在线视频| 999久久久国产精品视频| 一级毛片黄色毛片免费观看视频| 男女国产视频网站| 成年av动漫网址| 国产亚洲精品久久久久5区| 色综合欧美亚洲国产小说| 久久国产精品影院| 国产在线视频一区二区| 高清av免费在线| 黄网站色视频无遮挡免费观看| 五月开心婷婷网| 精品久久蜜臀av无| 亚洲av综合色区一区| 看十八女毛片水多多多| www.精华液| 日韩中文字幕欧美一区二区 | 2021少妇久久久久久久久久久| av欧美777| 欧美精品一区二区免费开放| 在现免费观看毛片| 日韩,欧美,国产一区二区三区| 国产成人免费观看mmmm| 亚洲精品在线美女| 免费观看人在逋| 狂野欧美激情性xxxx| 婷婷色综合www| 国产深夜福利视频在线观看| 五月开心婷婷网| 日本wwww免费看| 美女国产高潮福利片在线看| 亚洲精品美女久久av网站| 亚洲av男天堂| 免费看av在线观看网站| 免费观看人在逋| 男女高潮啪啪啪动态图| 在线观看人妻少妇| 99精品久久久久人妻精品| 一本久久精品| 成人午夜精彩视频在线观看| 久久久亚洲精品成人影院| 午夜老司机福利片| 多毛熟女@视频| 欧美性长视频在线观看| 国产亚洲av片在线观看秒播厂| 国产成人精品在线电影| 视频区图区小说| 校园人妻丝袜中文字幕| 另类精品久久| 久久久欧美国产精品| 日韩人妻精品一区2区三区| 午夜久久久在线观看| 成年人午夜在线观看视频| 美女午夜性视频免费| 天堂中文最新版在线下载| 成人国语在线视频| 欧美乱码精品一区二区三区| 久久久国产欧美日韩av| 一级毛片黄色毛片免费观看视频| 成人国产av品久久久| 久久国产精品大桥未久av| 欧美精品一区二区免费开放| 精品一区在线观看国产| 一级黄色大片毛片| 日韩欧美一区视频在线观看| 大型av网站在线播放| 夫妻午夜视频| 91精品国产国语对白视频| 在线亚洲精品国产二区图片欧美| 国产又色又爽无遮挡免| 人成视频在线观看免费观看| av又黄又爽大尺度在线免费看| 国产免费一区二区三区四区乱码| www.自偷自拍.com| 成人亚洲精品一区在线观看| av福利片在线| 少妇猛男粗大的猛烈进出视频| 狠狠精品人妻久久久久久综合| 亚洲精品美女久久av网站| 成人黄色视频免费在线看| 一级毛片女人18水好多 | 亚洲五月婷婷丁香| 精品少妇内射三级| 国产91精品成人一区二区三区 | 婷婷色综合大香蕉| 久久人人爽人人片av| 天天添夜夜摸| av国产精品久久久久影院| 国产精品香港三级国产av潘金莲 | 男女午夜视频在线观看| 成年人黄色毛片网站| 欧美国产精品va在线观看不卡| 亚洲精品美女久久av网站| 美女扒开内裤让男人捅视频| 日本91视频免费播放| 两性夫妻黄色片| 免费日韩欧美在线观看| 成人亚洲欧美一区二区av| 久久国产精品人妻蜜桃| 国产精品麻豆人妻色哟哟久久| 国产成人欧美| 纵有疾风起免费观看全集完整版| 午夜日韩欧美国产| 欧美人与善性xxx| 国产男人的电影天堂91| av在线app专区| 国产精品国产三级专区第一集| 日韩免费高清中文字幕av| 亚洲国产精品999| 精品久久蜜臀av无| 精品亚洲成a人片在线观看| 久久99精品国语久久久| 欧美日韩av久久| 天天添夜夜摸| 高清不卡的av网站| 亚洲久久久国产精品| 免费久久久久久久精品成人欧美视频| 国产亚洲av片在线观看秒播厂| 如日韩欧美国产精品一区二区三区| 久久久久网色| 九草在线视频观看| 日本av手机在线免费观看| 欧美97在线视频| 日本一区二区免费在线视频| 亚洲精品中文字幕在线视频| 欧美亚洲 丝袜 人妻 在线| 2021少妇久久久久久久久久久| 国产在线免费精品| 国产成人av激情在线播放| 精品一品国产午夜福利视频| 久久久精品国产亚洲av高清涩受| 中文字幕亚洲精品专区| 大香蕉久久成人网| 老司机在亚洲福利影院| 少妇裸体淫交视频免费看高清 | 日韩 亚洲 欧美在线| 欧美日韩亚洲综合一区二区三区_| 亚洲欧美日韩高清在线视频 | 无遮挡黄片免费观看| 中文字幕人妻熟女乱码| 欧美日韩黄片免| 亚洲黑人精品在线| 国产片特级美女逼逼视频| 亚洲人成77777在线视频| 日本av手机在线免费观看| 赤兔流量卡办理| 一本色道久久久久久精品综合| 狠狠精品人妻久久久久久综合| 欧美日韩成人在线一区二区| 秋霞在线观看毛片| 国产黄色免费在线视频| 午夜激情av网站| 免费观看人在逋| 国产1区2区3区精品| 亚洲国产精品999| √禁漫天堂资源中文www| 丝袜喷水一区| 十分钟在线观看高清视频www| 午夜视频精品福利| 香蕉国产在线看| 久久久久国产精品人妻一区二区| 永久免费av网站大全| 美国免费a级毛片| 女性生殖器流出的白浆| 久久久久久久精品精品| 一二三四社区在线视频社区8| 男女高潮啪啪啪动态图| 赤兔流量卡办理| 亚洲综合色网址| 国产亚洲精品第一综合不卡| 国产成人免费观看mmmm| 女人高潮潮喷娇喘18禁视频| 色婷婷av一区二区三区视频| 亚洲精品一区蜜桃| 精品国产乱码久久久久久小说| 国产日韩欧美视频二区| 亚洲av国产av综合av卡| 亚洲专区国产一区二区| av网站免费在线观看视频| 中文字幕精品免费在线观看视频| 国产有黄有色有爽视频| 99九九在线精品视频| 成年人免费黄色播放视频| 亚洲国产精品一区三区| 国产av一区二区精品久久| 国产精品人妻久久久影院| 黑人巨大精品欧美一区二区蜜桃| 日本一区二区免费在线视频| 18禁观看日本| 日本wwww免费看| 狠狠精品人妻久久久久久综合| 亚洲色图 男人天堂 中文字幕| 亚洲熟女毛片儿| 久久久亚洲精品成人影院| 美女国产高潮福利片在线看| 欧美xxⅹ黑人| a级毛片黄视频| 欧美xxⅹ黑人| 久久久久网色| 母亲3免费完整高清在线观看| 老司机亚洲免费影院| 国产熟女欧美一区二区| 老司机深夜福利视频在线观看 | 国产精品亚洲av一区麻豆| 精品亚洲乱码少妇综合久久| 国产成人免费观看mmmm| 国产精品久久久av美女十八| 亚洲综合色网址| 国产日韩欧美在线精品| 高清av免费在线| 国产真人三级小视频在线观看| 少妇的丰满在线观看| 丁香六月天网| xxx大片免费视频| 成年av动漫网址| 婷婷色麻豆天堂久久| 国产成人a∨麻豆精品| 亚洲精品一区蜜桃| 亚洲av欧美aⅴ国产| 国产日韩一区二区三区精品不卡| 亚洲少妇的诱惑av| 丝袜人妻中文字幕| 99久久综合免费| 赤兔流量卡办理| 老司机午夜十八禁免费视频| 九色亚洲精品在线播放| 90打野战视频偷拍视频| 欧美精品av麻豆av| 涩涩av久久男人的天堂| 只有这里有精品99| 婷婷色综合www| 十八禁网站网址无遮挡| 自拍欧美九色日韩亚洲蝌蚪91| 老司机靠b影院| 丝袜美腿诱惑在线| 青青草视频在线视频观看| 国产精品麻豆人妻色哟哟久久| 欧美亚洲 丝袜 人妻 在线| 亚洲欧美一区二区三区久久| 国产欧美日韩一区二区三 | 国产免费现黄频在线看| svipshipincom国产片| 亚洲伊人色综图| 丝袜人妻中文字幕| 无限看片的www在线观看| 欧美黑人欧美精品刺激| 18禁裸乳无遮挡动漫免费视频| 热re99久久精品国产66热6| 免费观看人在逋| 亚洲色图综合在线观看| av欧美777| 老司机午夜十八禁免费视频| 久久久精品区二区三区| 男女免费视频国产| 真人做人爱边吃奶动态| 黄片小视频在线播放| 亚洲国产中文字幕在线视频| 亚洲少妇的诱惑av| 别揉我奶头~嗯~啊~动态视频 | 一边摸一边抽搐一进一出视频| 在线观看免费视频网站a站| 亚洲欧美成人综合另类久久久| 亚洲一区中文字幕在线| 久久久国产精品麻豆| 制服人妻中文乱码| 亚洲国产毛片av蜜桃av| 精品福利观看| 国产精品九九99| 欧美大码av| 免费少妇av软件| 午夜老司机福利片| 亚洲三区欧美一区| 色网站视频免费| 久久ye,这里只有精品| 亚洲国产看品久久| 国产成人91sexporn| 免费在线观看影片大全网站 | 成年女人毛片免费观看观看9 | 狂野欧美激情性bbbbbb| 黄色视频不卡| 国产高清国产精品国产三级| 欧美av亚洲av综合av国产av| 一区在线观看完整版| 国产精品免费视频内射| 你懂的网址亚洲精品在线观看| 啦啦啦在线免费观看视频4| 亚洲精品一二三| 精品欧美一区二区三区在线| 久久久久久久大尺度免费视频| 又大又爽又粗| 手机成人av网站| 最近中文字幕2019免费版| xxxhd国产人妻xxx| 国产91精品成人一区二区三区 | 亚洲人成77777在线视频| xxxhd国产人妻xxx| 亚洲国产精品999| 久久青草综合色| 国产精品一区二区在线观看99| 夜夜骑夜夜射夜夜干| 纵有疾风起免费观看全集完整版| 天堂俺去俺来也www色官网| av福利片在线| 国产一区有黄有色的免费视频| 午夜久久久在线观看| 亚洲第一青青草原| 啦啦啦视频在线资源免费观看| 老司机影院毛片| 超碰成人久久| 欧美日韩成人在线一区二区| 成年av动漫网址| 国产精品亚洲av一区麻豆| 一区福利在线观看| 久久人人爽av亚洲精品天堂| 国产一区二区激情短视频 | 国产又爽黄色视频| 久久影院123| 美女高潮到喷水免费观看| 天天躁日日躁夜夜躁夜夜| 男女下面插进去视频免费观看| 久热爱精品视频在线9| 后天国语完整版免费观看| 久久精品国产a三级三级三级| 精品欧美一区二区三区在线| 你懂的网址亚洲精品在线观看| 最近手机中文字幕大全| 亚洲一区二区三区欧美精品| 国产欧美日韩一区二区三区在线| 国产成人欧美在线观看 | 国产成人a∨麻豆精品| 免费观看人在逋| 日韩 亚洲 欧美在线| 欧美日韩视频精品一区| 国产野战对白在线观看| 男女边吃奶边做爰视频| 亚洲精品美女久久av网站| 人妻人人澡人人爽人人| 久久精品成人免费网站| 成年av动漫网址| 精品欧美一区二区三区在线| 午夜福利,免费看| 色播在线永久视频| av网站免费在线观看视频| 久久久久久免费高清国产稀缺| 99热国产这里只有精品6| 国产精品麻豆人妻色哟哟久久| 精品一区在线观看国产| 性高湖久久久久久久久免费观看| 欧美人与性动交α欧美软件| 免费看十八禁软件| 丝瓜视频免费看黄片| 久久久国产一区二区| 深夜精品福利| 婷婷色综合www| 又大又爽又粗| 国产在线一区二区三区精| 2018国产大陆天天弄谢| 99re6热这里在线精品视频| 天天躁日日躁夜夜躁夜夜| 久久天躁狠狠躁夜夜2o2o | 国产亚洲欧美精品永久| 欧美人与善性xxx| 精品久久久精品久久久| 18在线观看网站| 黑人巨大精品欧美一区二区蜜桃| www.精华液| 看免费成人av毛片| 亚洲精品日本国产第一区| 欧美成狂野欧美在线观看| 丁香六月天网| 两个人看的免费小视频| 久久人妻福利社区极品人妻图片 | 黄色片一级片一级黄色片| 亚洲国产看品久久| 脱女人内裤的视频| 国产福利在线免费观看视频| 亚洲国产欧美日韩在线播放| 91麻豆av在线| 日韩,欧美,国产一区二区三区| 超碰成人久久| 波多野结衣一区麻豆| 国产成人系列免费观看| 50天的宝宝边吃奶边哭怎么回事| 你懂的网址亚洲精品在线观看| 亚洲熟女毛片儿| 两人在一起打扑克的视频| 免费看不卡的av| 色播在线永久视频| 国产黄色免费在线视频| 日本五十路高清| 国产视频首页在线观看| 欧美日韩综合久久久久久| 国产视频一区二区在线看| 老司机午夜十八禁免费视频| 国产不卡av网站在线观看| av国产久精品久网站免费入址| 色精品久久人妻99蜜桃| 精品亚洲乱码少妇综合久久| 欧美大码av| 叶爱在线成人免费视频播放| 一区二区三区激情视频| 国产成人欧美在线观看 | 最近手机中文字幕大全| 国产色视频综合| 欧美av亚洲av综合av国产av| 黄色怎么调成土黄色| 国产免费福利视频在线观看| 欧美成人精品欧美一级黄| av网站在线播放免费| 国产在线免费精品| 天天躁日日躁夜夜躁夜夜| 男女免费视频国产| 人人澡人人妻人| 在线观看免费高清a一片| 国产欧美日韩一区二区三 | 又大又爽又粗| 激情视频va一区二区三区| 丝袜人妻中文字幕| 亚洲一区中文字幕在线| 麻豆av在线久日| 午夜视频精品福利| 国产成人精品在线电影| 一边摸一边抽搐一进一出视频| 亚洲五月婷婷丁香| 精品亚洲乱码少妇综合久久| 2018国产大陆天天弄谢| 亚洲av美国av| 99久久99久久久精品蜜桃| 高清黄色对白视频在线免费看| 亚洲精品第二区| 亚洲一卡2卡3卡4卡5卡精品中文| 日韩制服丝袜自拍偷拍| 中国美女看黄片| 18禁观看日本| av又黄又爽大尺度在线免费看| 一区二区日韩欧美中文字幕| 国产精品亚洲av一区麻豆| 国产午夜精品一二区理论片| 亚洲欧洲精品一区二区精品久久久| 亚洲成人国产一区在线观看 | 亚洲国产精品国产精品| 在线av久久热| 美女午夜性视频免费| 一级毛片电影观看| 欧美老熟妇乱子伦牲交| 亚洲欧洲精品一区二区精品久久久| 国产av一区二区精品久久| 久久99热这里只频精品6学生| 性高湖久久久久久久久免费观看| 亚洲自偷自拍图片 自拍| 亚洲国产精品一区二区三区在线| 9191精品国产免费久久| 免费观看av网站的网址| 黄网站色视频无遮挡免费观看| 久久这里只有精品19| 久久人人爽人人片av| 亚洲精品自拍成人| 飞空精品影院首页| 一区二区av电影网| 精品一品国产午夜福利视频| 久久精品久久精品一区二区三区| 男女之事视频高清在线观看 | 九草在线视频观看| 又黄又粗又硬又大视频| a级毛片在线看网站| 狠狠精品人妻久久久久久综合| 精品欧美一区二区三区在线| 国产精品二区激情视频| 人妻人人澡人人爽人人| 亚洲精品久久久久久婷婷小说| 欧美+亚洲+日韩+国产| 美女扒开内裤让男人捅视频| 看免费成人av毛片| a级毛片黄视频| 99香蕉大伊视频| 欧美97在线视频| 脱女人内裤的视频| 中文字幕人妻丝袜制服| 国产精品国产三级国产专区5o| 老司机午夜十八禁免费视频|