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

    Characterization of a Novel Esterase Belonging to Family V from Marinobacter flavimaris

    2024-03-12 11:14:24HEJingjingZHANGYunhuiWULeileiWANGYaruZHANGHeLIUZhengangandSHIXiaochong
    Journal of Ocean University of China 2024年1期

    HE Jingjing, ZHANG Yunhui , WU Leilei, WANG Yaru, ZHANG He,LIU Zhengang, and SHI Xiaochong, ,

    1) Frontiers Science Center for Deep Ocean Multispheres and Earth System, and College of Marine Life Sciences,Ocean University of China, Qingdao 266003,China

    2) Laboratory for Marine Ecology and Environmental Science, Laoshan Laboratory, Qingdao 266237, China

    3) Institute of Evolution & Marine Biodiversity, Ocean University of China, Qingdao 266003, China

    Abstract Lipolytic enzymes have attracted enormous attentions because of their ability in ester hydrolysis, ester synthesis, transesterification and other biochemical reactions. Bacteria are important sources of lipolytic enzymes applied in industry. Here, a novel lipolytic enzyme encoded by esterase gene est1347 was identified in Marinobacter flavimaris WLL162, and was purified and characterized. The lipolytic enzyme Est1347 consisted of 312 amino acid residues and a 21-amino-acids N-terminal signal peptide with a predicted molecular weight of 34.2 kDa. It belongs to family V of bacterial lipolytic enzymes based on the amino acid sequence homology analysis. Est1347 is a mesophilic and alkali-resistant enzyme with the highest activity at 45℃ and pH 8.5; it is stable at temperatures below 50℃ and pH 7.5 – 11.0. Est1347 showed a preference for middle-length chain substrate p-NPC10 and a wide range of other substrates. The Km, Vmax, Kcat and Kcat/Km values of Est1347 for p-NPC10 in pH 8.5 at 45℃ were 0.9411 mmol L?1, 1285 μmol min?1 mg?1, 698.91 s?1 and 743.65 s?1 (mmol L?1)?1, respectively. It is also tolerant to the metal ions, organic solvents and detergents. In conclusion, the esterase Est1347 laid a foundation for further study of bacterial lipolytic enzyme family V.

    Key words esterase; Marinobacter flavimaris; enzymatic properties; lipolytic enzyme family V

    1 Introduction

    Lipolytic enzymes can catalyze the hydrolysis and synthesis of the ester bond, and have a function in transesterification and other reactions. It distributes extensively in bacteria, archaea, viruses and eukaryote organisms, while bacteria are the main sources of industry-applied lipolytic enzymes. Bacterial lipolytic enzymes include carboxylesterases (EC 3.1.1.1, esterase) and true lipases (EC 3.1.1.3, lipase), and have putative physiological functions with virulence, metabolism and growth promotion (Kovacicet al.,2019). Esterases prefer water-soluble short-chain to medium-chain esters, whereas lipases hydrolyze the water-insoluble long-chain triglycerides (Jaegeret al., 1994). The catalysis of lipolytic enzymes is characterized by mild reaction conditions, high efficiency and no need for coenzyme, which can be widely applied to food, pharmaceutical and cosmetic industries (Bornscheuer, 2002).

    Bacterial lipolytic enzymes are classified into eight families (Family I – VIII) based on the amino acid sequences similarities and physiological properties (Arpigny and Jaeger, 1999). Recently, an increasing number of lipolytic enzymes have been characterized and classified, while 35 esterase families and 11 true lipase families have been newly introduced (Hitch and Clavel, 2019). The catalytic domain of lipolytic enzymes is composed of a classic α/βhydrolase fold, with eight β-sheets connected by α-helices in the central position and α-helices flanked on both sides(M?lgaardet al., 2000). The catalytic triad Ser, Asp/Glu and His are located in the catalytic domain and the active site Ser is in the conserved GXSXG pentapeptide (Olliset al.,1992). According to the structural homology, the ESTHER database reclassifies proteins containing the α/β-hydrolase fold into four blocks: C, H, L and X. Family V is classified into block X and it is divided into three subfamilies: Family V.1, Family V.2, Family V.3 (Lenfantet al., 2012).

    Microorganisms from extreme environments are found to possess lipolytic enzymes with special characteristics and these lipolytic enzymes can be used as uniquely valuable biocatalysts (Vermaet al., 2021). A thermostable lipolytic enzyme PLP from the metagenomic library of deep-sea hydrothermal vent demonstrates optimal activity at 70℃ and this property implies it can be used in the high-temperature food processing industry (Fuet al., 2015). The hyperthermostable carboxylesterase EstD9 from thermophilicAnoxybacillus geothermalisD9 presents excellent stability at high-temperature range (70 – 100℃) and high tolerance to various surfactants and metal ions; these enzymatic features suggest it is a highly suitable catalyst for industrial bioprocesses (Johanet al., 2022). The cold-adapted esterase EstDR4 from extremophileDeinococcus radioduransdisplays regioselectivity, enantioselectivity, and degradation effects on four insecticides, which may have the potential to be used in the pharmaceutical industry and the treatment for pollution caused by insecticides (Zhanget al.,2021). Another psychrophilic esterase Est19 from the Antarctic bacteriumPseudomonassp. E2-15 shows the stereo-specificity that can identify the methyl-L lactate and methyl-D lactate, and this characteristic suggests Est19 may be useful for cold and chiral catalyses (Liuet al., 2021).Additionally, the halophilic esterase EstA1 isolated from the marine bacteriumAlteromonassp. 39-G1 has excellent salttolerance characteristic and its activity is continuously increased with increasing NaCl concentration, and this unique property enhances its academic and industrial values(Wonet al., 2020). Another novel high salt-tolerant esterase EstSHJ2 remains high activity in the solution with 1.0 – 3.5 mol L?1NaCl while the activity is very low without NaCl(Wanget al., 2020). It is previously reported some properties of halophilic lipase can improve the flavor of fish sauce(Kanlayakrit and Boonpan, 2007), so the study of the halophilic lipolytic enzyme may have applicable values in the food industry. The lipolytic enzymes with diverse enzymatic characteristics from different environments have important applications in the food, dairy, detergent and pharmaceutical industries (Fuci?oset al., 2012). Therefore, mining unique esterase resources has raised enormous attention.

    Marinobacterbelongs toAlteromonadaceaein the class ofγ-Proteobacteria, which can degrade hydrocarbon and utilize the lipid (Duran, 2010; Boninet al., 2015). It is previously reported thatMarinobacterplays an indispensable role in the remediation of hydrocarbon pollution in highsalt environments (Lyuet al., 2022; Zhanget al., 2022).Meanwhile, mostMarinobacterisolates are found to have halophilic lipolytic ability, such asMarinobacter lipolyticusCEES 33 (Oveset al., 2017),Marinobacter litoralisSW-45 (Musaet al., 2019) andMarinobacter lipolyticusSM19 (Martínet al., 2003). Therefore, exploring the potential lipolytic enzymes inMarinobacteris of great significance. The halophilic lipase, LipBL, isolated from theMarinobacter lipolyticusSM19, belongs to family VIII and is highly efficient in producing eicosapentaenoic acid (EPA),whose properties are of great interest in the food industry(Pérezet al., 2011). The active sites of LipBL are investigated through site-directed mutagenesis, the critical S-XX-K motif is found to have a crucial effect on the catalysis and the conserved GXSXG pentapeptide presumably represents a substrate-binding site, which is different from the classical catalytic residues (Pérezet al., 2012).

    Some reported lipolytic enzymes have poor stability,which is an important factor limiting their industrial application. Therefore, exploring novel lipolytic enzymes with great stability from microorganisms is a feasible method.Recently,Marinobacter flavimarisWLL162 with strong lipolytic enzyme activity was isolated from the Western Pacific Cobalt-rich crust (Wu, 2020) and 12 putative lipolytic enzyme genes were identified in its genome sequence.In this study, a novel esterase geneest1347fromM. flavimarisWLL162 was cloned and the esterase Est1347 was expressed, purified and characterized. These results propose a novel esterase from the marine source and provide a supplement for family V of the bacterial lipolytic enzyme.

    2 Materials and Methods

    2.1 Bacterial Strains, Plasmids, Media and Growth Conditions

    Marinobacter flavimarisWLL162 (=Marinobacter flavimarisSW-145T) (Yoonet al., 2004) was isolated from the sediment (5443 m water depth) of Western Pacific Cobalt-rich crust area at 157.2144?E, 19.60048?N and cultured using 2216E (MA, marine agar) at 28℃. The cloning and expression vectors were pUCm-T (Sangon Biotech,Shanghai, China) and pET-24a (+) (Novagen, Beijing, China), respectively.Escherichia coliDH5α was used for gene clone andE. coliBL21 (DE3) was used for the protein expression, while both of them were cultured at 37℃ in Luria-Bertani (LB) agar.

    2.2 Sequence and Phylogenetic Analysis of Esterase Gene

    The genomic DNA from WLL162 was extracted using a DNA extraction kit (TIANGEN BIOTECH, Beijing, China). The sequencing and assembling of total genomic DNA were based on the methods of Linet al. (2018). The prediction and annotation of the Open reading frame (ORF)were carried out by the RAST server (Overbeeket al., 2014).The amino acid sequence was analyzed by BLASTP against the protein data bank (PDB) database and UniProt Knowledgebase (UniProtKB) (Batemanet al., 2021). The multiple sequence alignment between Est1347 and reference sequences was performed using ClustalX 2.1 (Larkinet al.,2007) and the result was rendered using the ESPript 3.0 online tool (ESPript 3.x / ENDscript 2.x (ibcp.fr)). The phylogenetic tree based on the amino acid sequence was constructed by MEGA X (Kumaret al., 2018), with 1000 bootstraps, and sequences of other lipolytic enzymes were downloaded from GenBank. The three-dimensional model was predicted using Swiss-Model online tool and analyzed by PyMOL. The molecular weight and pI of esterase were predicted by the ExPASy database (Gasteigeret al., 2003)and the signal peptide was predicted by the SignalP 5.0 server (Almagro Armenteroset al., 2019).

    2.3 Cloning, Expression and Purification of Est1347

    To obtain the recombinant proteins without signal peptide, theest1347gene (GenBank accession number OQ29 5858) was amplifiedviapolymerase chain reaction (PCR)with the primer pairs,est1347-F(CGCATGTG CTCTCGTCAG) andest1347-R(CCCCTCGA CCGGCTGCC), which containedBamHI andHindⅢ digestion sites (underlined), respectively. The PCR products were cloned to the pUCm-T vector and the resulting plasmid was transformed intoE. coliDH5α competent cells.The recombinant plasmid pUCm-T/est1347was extracted and digested byBamHI andHindⅢ restriction endonucleases after sequencing confirmation; the products were purified and inserted into the linearized plasmid pET-24a (+).The correct expression vector pET-24a (+)/est1347was introduced intoE. coliBL21 (DE3) competent cells for protein expression. The recombinant cells were grown in LB medium supplemented with 50 μg mL?1kanamycin at 37℃.At the mid-exponential growth phase (OD600nm= 0.4 – 0.6),isopropyl-β-D-thiogalactopyranoside (IPTG) was added to a final concentration of 0.05 mmol L?1, and the cultivation was continued at 16℃, with shaking at 150 r min?1for another 12 – 14 h. The induced cells were harvested by centrifugation at 8000 r min?1, 4℃ for 10 min, washed with prechilled PBS buffer, and resuspended with prechilled binding buffer (20 mmol L?1Tris-HCl, 0.5 mol L?1NaCl, pH 8.0).After ultrasonic disruption on ice, the cell debris was removed by centrifugation at 12000 r min?1, 4℃ for 10 min.The supernatant was filtered through a 0.22-μm-pore-size filter and loaded onto the Ni-NTA affinity column (Qiagen,Hilden, Germany). Proteins were purified according to the manufacturer’s recommendations. The eluate containing the target protein was dialyzed with dialysate (20 mmol L?1Tris-HCl, 0.85% NaCl, 10% glycerol, pH 8.0) for 24 h at 4℃. The expression and purification of protein were verified by 12% sulfate-polyacrylamide gel electrophoresis(SDS-PAGE) (Laemmli, 1970).

    2.4 Esterase Activity Detection of Est1347

    The enzymatic activity of purified recombinant Est1347 was tested by the reported method with some modification(Wanget al., 2016). Esterase activity againstp-nitrophenyl (p-NP) esters was determined by measuring the amount ofp-nitrophenol released by catalytic hydrolysis. The production ofp-nitrophenol was measured by the absorbance at 405 nm. The standard reaction system contained 100 m mol L?1ofp-nitrophenyl butyrate (p-NPC4) in acetonitrile,ethanol, and 50 mmol L?1of Tris-HCl buffer (pH 8.5) at a ratio of 1:4:94 (V/V/V). The reaction system was preincubated at 45℃ for 5 min and initiated by adding 10 μL of purified enzyme solution to 1 mL of substrate mixture. After 10 min, the absorbance of the reaction system at 405 nm was measured. One unit of esterase activity was defined as the amount of esterase capable of releasing 1 μmol ofp-nitrophenol per minute. Protein concentration was measured by the Bradford method (Bradford, 1976) using bovine serum albumin (BSA) as the standard.

    2.5 Biochemical Characterization of Est1347

    The optimal temperature of Est1347 was determined in optimal pH at 0, 5, 15, 20, 30, 40, 45, 50, 55, 60, and 70℃.For detecting the thermostability, the enzyme samples were incubated at 0, 5, 15, 20, 28, 40, 45, 50, 55, 60, 70, 80,and 90℃ for 1 h, and then the residual activities were measured. The optimal pH of Est1347 was measured with pH range of 3.0 – 11.0 using four kinds of buffer systems: 50 mmol L?1citric acid/sodium citrate buffer (pH 3.0 – 6.0), 50 mmol L?1Na2HPO4-NaH2PO4buffer (pH 6.0 – 8.0), 50 mmol L?1Tris-HCl (pH 7.5 – 9.0) and 50 mmol L?1Na2HPO4-NaOH (pH 10.0 and 11.0). To determine the pH stability of esterase, the purified enzyme samples were preincubated in the above-mentioned buffers at 4℃ for 12 h,and the residual activities were measured at 45℃. Substrate specificity was assayed with the followingp-nitrophenyl esters:p-nitrophenyl acetate (p-NPC2),p-nitrophenyl butyrate (p-NPC4),p-nitrophenyl caproate (p-NPC6),p-nitrophenyl caprylate (p-NPC8),p-nitrophenyl decanoate (p-NPC10),p-nitrophenyl laurate (p-NPC12) andp-nitrophenyl palmitate (p-NPC16).

    The kinetic parameter of esterase was determined by measuring enzymatic activity usingp-NPC10 as substrate at 45℃ in pH 8.5 for 10 min with substrate concentrations of 0.01, 0.05, 0.1, 0.15, and 0.2 – 1.2 mmol L?1(at an interval of 0.2). TheKmandVmaxvalues were calculated using the Michaelis-Menten equation in GraphPad Prism 9.0. And theKcatandKcat/Kmwere calculated based on the value ofVmax.

    The effect of metal ions (Fe3+, Co2+, Mg2+, Cu2+, Ni2+,Mn2+, Na+, K+, Ca2+, Zn2+) and metal ion chelator EDTA on esterase activity was determined by adding different ions and EDTA at a final concentration of 1 mmol L?1and 10 mmol L?1into the reaction system. The effect of salt concentration on Est1347 activity was measured with the concentration of NaCl of 0 – 5 mol L?1in the reaction system.The organic solvents (acetonitrile, DMF, isopropanol, methanol, acetone, DMSO, glycerol) were added respectively to the standard reaction system (10%, V/V) to investigate the effect of organic solvents on Est1347. The detergents (SDS, Triton X-100, Tween-20, Tween-80) were added respectively to the standard reaction system (1%, V/V)to study the effect of detergents on the enzyme.

    3 Results

    3.1 Protein Sequence and Phylogenetic Analysis of Est1347

    Based on the result of genome annotation of the putative esterase geneest1347, a 939-bp open reading frame(ORF) was identified. The ORF encodes a protein containing 312 amino acid residues with a predicted molecular weight of 34.2 kDa and a predicted pI of 4.51. The N-terminal signal peptide was composed of 1 – 21 amino acids.Sequence alignment revealed that this protein shared the highest identity of 41.42% with the membrane phospholipase A (UniProtKB: Q9KJG6) fromPseudomonas aeruginosaPAO1. The phylogenetic analysis of Est1347 and known bacterial lipolytic enzyme families showed that Est1347 belongs to family V of bacterial lipolytic enzymes(Fig.1). Multiple sequence alignment of Est1347 with highly related homologs showed that the canonical catalytic triad containing Ser139, Asp261and His289, and Ser139is located in the GXSXGG motif (Fig.2). A three-dimensional structure of Est1347 was predicted based on the template 6I8W (the crystal structure of the protein Q9KJG6) and the model structure revealed that Est1347 had a typical α/β-hydrolase fold (Fig.3a). The multiple sequence alignment and phylogenetic analysis results confirmed that Est1347 is a novel lipolytic enzyme belonging to bacterial lipolytic enzyme family V.

    Fig.1 Phylogenetic tree of lipolytic enzymes based on the homology of protein sequences, constructed using the neighborjoining method by MEGA X. The sequences used in the analysis were obtained from the GenBank database. Bootstrap values are based on 1000 replicates, and only values > 50% are shown. The scale bar indicates the number of amino acid substitutions per site.

    Fig.2 Multiple sequence alignments of amino acid sequences of Est1347 and homologs from different bacteria. Accession numbers of proteins in the UniProtKB are given for Q9KJG6 (Pseudomonas aeruginosa PAO1), P24640 (Moraxella sp.TA144), Q02104 (Psychrobacter immobilis B10), Q1QEU6 (Psychrobacter cryohalolentis K5), B0V9K7 (Acinetobacter baumannii AYE). The active-site residues are indicated by circles below the sequences.

    Fig.3 Predicted three-dimensional structure of Est1347 (a); SDS-PAGE analysis of recombinant esterase Est1347 (b); substrate specificity of Est1347 (c); and enzyme kinetics of recombinant Est1347 (d). (a) The three-dimensional structure of Est1347 was modeled against template 6I8W (the crystal structure of the protein Q9KJG6). The catalytic triad Ser, Asp and His is shown in yellow sticks. (b) SDS-PAGE of recombinant esterase Est1347. Lanes: M, molecular mass markers;lane 1, cell-free extracts of E. coli BL21 (DE3) harboring plasmid pET-24a (+)/est1347; lane 2, purified Est1347; lane 3,cell-free extracts of E. coli BL21 (DE3) harboring plasmid pET-24a (+). (c) Substrate specific of Est1347 toward p-NP esters. Substrate specific was determined using p-NP esters with various chain lengths (C2 – C16) at a final concentration of 1 mmol L?1 at 45℃ in 50 mmol L?1 Tris-HCl buffer (pH 8.5). (d) Determination of enzyme kinetics of Est1347 toward p-NPC10. All data are shown as mean ± SD (n = 3).

    3.2 Expression and Activity of Esterase Est1347

    The vector pET-24a (+) was used to express the esterase in theE. coliBL21(DE3) cells. The purified recombinant proteins were evaluated by SDS-PAGE (Fig.3b),which displayed a single band with an apparent molecular weight of approximately 32 kDa, accordant to the predicted molecular mass (32.02 kDa). About 7 mg of the recombinant protein was obtained after His-tag affinity purification from 300 mL LB broth. Est1347 could efficiently hydrolyze short- to middle-lengthp-NP esters (C2 – C10) with the highest activity towardp-NPC10 (6991 U mg?1) followed byp-NPC4 (4927.3 U mg?1) andp-NPC6 (2036.0 U mg?1)among thep-NP esters examined (Fig.3c). The specific activity of Est1347 was observed usingp-NPC10 as substrate in 50 mmol L?1Tris-HCl (pH 8.5) at 45℃, and theKmandVmaxfor esterase were 0.9411 mmol L?1and 1285 μmol min?1mg?1(Fig.3d), respectively. TheKcatandKcat/Kmvalues of Est1347 were calculated as 698.91 s?1and 743.65 s?1(mmol L?1)?1, respectively. In addition, Est1347 also exhibited hydrolytic activity against tributyrin and Tween-20, 40 and 80 (Fig.4).

    3.3 Effects of Temperature and pH on Activity and Stability of Est1347

    The purified esterase was evaluated for its activity and stability at different temperatures and pH values. The optimum temperature of Est1347 was 45℃ and it was active in a wide range of 0 – 60℃ (Fig.5a). For the stability of Est1347, it could retain more than 80% of its highest activity after 1 h incubation at temperatures 0 – 50℃, but lost almost all activity after 1 h incubation at 90℃ (Fig.5b).

    Fig.5 Responses of temperature of recombinant esterase Est1347. (a) Effect of temperature on Est1347 activity. The highest activity of Est1347 at 45℃ was taken as 100%. (b) Thermostability of Est1347. The enzyme was incubated at 0 – 90℃for 1 h. The remaining activity was measured under optimal conditions. The highest activity at 15℃ was taken as 100%.The data are shown as mean ± SD (n = 3).

    Est1347 showed the highest activity at pH 8.5 and exhibited activity in a narrow pH range (6.0 – 9.0, Fig.6a). Moreover, it could maintain most of its initial activity after 12 h incubation at pH 7.5 – 11.0 (Fig.6b), which indicated that Est1347 had relatively high pH stability and it preferred alkaline environment.

    Fig.6 Responses of pH of recombinant esterase Est1347. (a) Effect of pH on the activity of Est1347. The absorbance of the reaction system was measured at 405 nm. (b) pH-stability of Est1347. The enzyme was incubated at pH from 3.0 to 11.0 for 12 h at 4℃. The remaining activity was measured under optimal conditions. The absorbance of the reaction system with preincubated enzyme was measured at 405 nm. The data are shown as mean ± SD (n = 3).

    3.4 Effects of Metal Ions, EDTA and Salinity on Activity of Est1347

    The effects of metal ions and ion chelator EDTA on the activity of Est1347 were also investigated (Fig.7a). The results showed that only K+(1 mmol L?1) could improve the enzyme activity. Mg2+, Mn2+, Na+, Ca2+(1 mmol L?1and 10 mmol L?1) and K+(10 mmol L?1) had no obvious effect on Est1347. Fe3+, Co2+, Cu2+, Ni2+and Zn2+(1 mmol L?1and 10 mmol L?1) could inhibit the activity of Est1347 significantly. The ion chelator EDTA had no significant effect on Est1347 at 1 mmol L?1while inhibited enzymatic activity at 10 mmol L?1.

    Fig.7 Response of Est1347 to metal ions and EDTA (a) and salinity (b). (a) Effects of 1 mmol L?1 and 10 mmol L?1 metal ions and EDTA on the activity of Est1347; ** P < 0.01, * P < 0.05; (b) Effects of different NaCl concentrations on the activity of Est1347. The data are shown as mean ± SD (n = 3).

    The activity of Est1347 decreased (Fig.7b) with the increase of salinity. Est1347 could remainca.83% andca.69% of the initial activity with 0.5 mol L?1and 1 mol L?1NaCl in the reaction system respectively, which indicates that Est1347 is a moderate salt-tolerant esterase.

    3.5 Effects of Organic Solvents and Detergents on the Activity of Est1347

    The effects of different organic solvents on the activity of Est1347 were measured (Fig.8). The methanol, DMSO and glycerol had no obvious effect on the enzymatic activity of Est1347. While acetonitrile, DMF, isopropanol, ethanol and acetone inhibited the enzymatic activity significantly, and acetonitrile showed the strongest inhibition effect which could reduce the activity of Est1347 to onlyca.16%.The detergents also showed significant effects on the activity of Est1347 (Table 1). The ionic detergent SDS and nonionic detergents Tween-20 and Tween-80 inhibited the activity of Est1347, wherein SDS and Tween-80 could inhibit the activity completely. However, the non-ionic detergent Triton X-100 could increase the enzyme activity to 106.9%.

    Table 1 Effects of various detergents on the activity of Est1347

    Fig.8 Response of Est1347 to various organic solvents. The value observed without organic solvent was taken as control (100%). ** P < 0.01, * P < 0.05.

    4 Discussion

    In this study, a novel esterase Est1347 was identified from the marine bacteriumMarinobacter flavimarisWLL162 and classified into the bacterial lipolytic enzyme family V based on the conserved motif and homologs of amino acid sequences. Among the characterized lipolytic enzymes,Est1347 had the highest identity (41.24%) to the α/β fold hydrolase ofPseudomonas aeruginosaPAO1 (Bleffertet al.,2022), which was a novel bacterial virulence factor and the esterase activity of which was characterized using 4-methylumbelliferyl palmitate andp-nitrophenyl butyrate as substrates. Est1347 also shared similarities with esterase EstPc fromPsychrobacter cryohalolentisK5 (accession number ABE73807.1), lipase Lip1 fromPsychrobacter immobilisB10 (accession number CAA47949.1) and lipase Lip3 fromPsychotroph moraxellaTA144 (accession number CAA37863.1), with 36.88%, 35.21% and 34.86% amino acid sequence identity, respectively. These three lipolytic enzymes were characterized by psychrophilic activeties (Felleret al., 1991; Arpignyet al., 1993; Novototskaya-Vlasovaet al., 2012). However, Est1347 showed the highest activity at moderate temperature (Table 2), which indicates that Est1347 probably represents the moderate temperature type esterase in family V.

    Table 2 Comparison of properties of known esterases from different bacterial lipolytic enzyme families

    Lipolytic enzymes grouped in family V were originated from various organisms in different environments and showed diverse characteristics, they remained to be further explored (Ruizet al., 2007). The novel esterase Est1347 in this study displayed great thermostability at its optimal temperature range (retaining over 80% activity after incubation at 0 – 50℃ for 1 h). The esterase Est_p1, a mesophilic enzyme and a member of family V, displayed the highest activity at 40℃ and retained about 50% of activity after incubated at 40℃ for 1 h while the residual activity at 45℃or 50℃ dropped rapidly within 15 min (Penget al., 2011).Another esterase EstIM1 showed its highest activity at 40℃ but it was sensitive to 40℃, whose activity could retain no more than 20% after incubation at 40℃ for 20 min(Koet al., 2012). These results indicate Est1347 is more suitable for medium-temperature catalysis and is capable of maintaining catalytic activity at the optimum reaction temperature for a long time.

    For the optimal pH and pH stability, the maximum activity of Est1347 was observed at pH 8.5, similar to other alkaline esterases belonging to family V, which include EstPc fromPsychrobacter cryohalolentisK5 (Novototskaya-Vlasovaet al., 2012), Est16 from a metagenomic library of a microbial consortium specialized for diesel oil degradation (Pereiraet al., 2015) and EstKa fromKlebsiella aerogenesGF-0526 (Gaoet al., 2021). Interestingly, Est1347 retained most of the enzyme activity at pH 7.5 – 11.0 after incubation at relevant buffers for 12 h, which indicated it had obvious advantages in the preparation of enzyme and industrial application. Alkali-stable enzymes usually possess more arginine and neutral hydrophilic amino acid (histidine, asparagine and glutamine) residues, and the arginine residues can form hydrogen bonds or ion pairs to adapt to the alkaline conditions (Shiraiet al., 1997). Est1347 had a proportion of 15.4% of arginine and neutral hydrophilic amino acid residues, which may be related to its alkali stability.

    Most lipolytic enzymes in family V preferred the shortor middle-length (C < 10) chainsp-NP esters as the substrate. For example, EstV showed a preference for shortlength (C2 – C6)p-NP esters (Ruizet al., 2007), est-OKK had the maximal activity towardp-NPC4 (Yanget al., 2018),BroH displayed the highest activity towardp-NPC6 (Chenet al., 2013) and FCLip1 preferentially hydrolyzedp-NPC8(Caiet al., 2011). However, Est1347 had the highest activity againstp-NPC10, which was different from other members and provided new insight into the lipolytic enzyme family V. According to former reports, the substrate-binding pocket of theE. colithioesterase TesA was reshaped by structure-guided mutational screening to improve the C8 substrate selectivity (Denget al., 2020). Considering the preference of Est1347 to substratep-NPC10, we deduced the substrate-binding pocket of Est1347 probably prefer-red binding the C10 substrate over the C4 substrate. This deduction requires further confirmation. Furthermore, Est1347 prefers short- to middle-length chain acyl substrates and could degrade tributyrin, which confirms the function of Est1347 as an esterase. In addition, Est1347 also could degrade non-ionic detergents Tween-20, Tween-40 and Tween-80, indicating Est1347 might be employed in detergent pollution treatment.

    The specific activities of enzymes from different sources were summarized and the kinetic paraments of Est1347 were compared with other esterases. The best substrate of Est1347 wasp-NPC10 and theVmaxvalue (1285 μmol min?1mg?1) of Est1347 againstp-NPC10 was at a moderate level (Table 2). It had been reported that medium chain length(C6 – C10) fatty acids and their derivatives had unique application values, such as antibacterial, lubricant, substitutes for types of gasoline, and intermediates for herbicide synthesis (Sarriaet al., 2017). The property of Est1347 which preferred hydrolyzingp-NPC10 implied it might have the potential in the production of middle-chain fat acid. Besides,the cleaning experiment was also performed according to the method of Zafaret al.(2022). And esterase Est1347 combined with detergent could removeca. 14.12% more stained oil than detergent alone. The oil in this study was with relatively long chain length. Est1347 possibly has higher cleaning efficiency on the middle-length oil stain. The potential of Est1347 applied to commercial laundry detergents still needs to be further investigated.

    5 Conclusions

    In this study, a novel esterase geneest1347was identified and cloned fromMarinobacter flavimarisWLL162 and the encoding protein Est1347 was expressed successfully usingE. coli. The esterase Est1347 belongs to the bacterial lipolytic enzyme family V and is a mesophilic and alkali-resistant enzyme. Est1347 showed a preference for middle-length chain substratep-NPC10 and exhibited a large range of substrates. Furthermore, Est1347 showed great thermostability and pH stability, which could maintain more than 80% of its maximal activity at 0 – 50℃ and pH 7.5 –11.0. These characteristics of Est1347 prompted it to be a candidate for biotechnological applications in industry.

    Acknowledgements

    This work was supported by the projects from the National Natural Science Foundation of China (No. 42230411),and the China Ocean Mineral Resources R and D Association (COMRA) Special Foundation (No. DY135-B2-10).

    国产成人一区二区三区免费视频网站| 最近在线观看免费完整版| 日日摸夜夜添夜夜添小说| 欧美乱妇无乱码| 欧美黑人欧美精品刺激| 成人av一区二区三区在线看| 变态另类成人亚洲欧美熟女| 国产精品亚洲美女久久久| 国产熟女xx| 国产精品美女特级片免费视频播放器 | 一夜夜www| 12—13女人毛片做爰片一| 怎么达到女性高潮| 色av中文字幕| 色视频www国产| 成年人黄色毛片网站| 亚洲国产精品成人综合色| 岛国视频午夜一区免费看| 国产精品,欧美在线| 精品国产美女av久久久久小说| 国产一区二区在线观看日韩 | 一卡2卡三卡四卡精品乱码亚洲| 两性夫妻黄色片| 黄片小视频在线播放| 男人舔女人下体高潮全视频| 国产精品久久久久久亚洲av鲁大| 国产高潮美女av| 国产免费男女视频| 午夜成年电影在线免费观看| 真人一进一出gif抽搐免费| 1024手机看黄色片| 久久亚洲真实| 久久天躁狠狠躁夜夜2o2o| 亚洲精华国产精华精| 毛片女人毛片| 哪里可以看免费的av片| 1024香蕉在线观看| 国产蜜桃级精品一区二区三区| 欧美黄色片欧美黄色片| av国产免费在线观看| 男人舔女人的私密视频| 麻豆成人av在线观看| 99国产精品99久久久久| 看黄色毛片网站| 看片在线看免费视频| 国产综合懂色| 国产高清有码在线观看视频| 欧美黄色片欧美黄色片| www.999成人在线观看| 91麻豆精品激情在线观看国产| 午夜a级毛片| 五月玫瑰六月丁香| 国产亚洲欧美在线一区二区| 1024香蕉在线观看| 18禁国产床啪视频网站| 亚洲av免费在线观看| 日韩欧美三级三区| 18禁国产床啪视频网站| 欧美日韩中文字幕国产精品一区二区三区| 国产精品精品国产色婷婷| 久久久久久久精品吃奶| 亚洲精华国产精华精| 偷拍熟女少妇极品色| 久久亚洲真实| 午夜福利视频1000在线观看| 精品一区二区三区视频在线 | 在线观看日韩欧美| 老司机午夜十八禁免费视频| 精品久久蜜臀av无| 亚洲色图 男人天堂 中文字幕| 欧美成人性av电影在线观看| 最近最新中文字幕大全免费视频| 亚洲国产精品成人综合色| 亚洲aⅴ乱码一区二区在线播放| 国产高清有码在线观看视频| 国内少妇人妻偷人精品xxx网站 | av欧美777| 不卡一级毛片| 国产一区二区三区视频了| 夜夜看夜夜爽夜夜摸| 99视频精品全部免费 在线 | 免费看美女性在线毛片视频| 一夜夜www| 两个人视频免费观看高清| 男人舔女人的私密视频| 成人精品一区二区免费| 噜噜噜噜噜久久久久久91| www日本黄色视频网| 国产精华一区二区三区| 欧美黄色片欧美黄色片| 国产精品九九99| 亚洲性夜色夜夜综合| 俄罗斯特黄特色一大片| 色播亚洲综合网| 欧美日韩乱码在线| 一边摸一边抽搐一进一小说| 国产真实乱freesex| 国产精品亚洲av一区麻豆| 国产精品久久电影中文字幕| 久久伊人香网站| 脱女人内裤的视频| 麻豆一二三区av精品| 国产精品亚洲美女久久久| 在线观看免费视频日本深夜| 亚洲一区二区三区色噜噜| 午夜福利高清视频| 又粗又爽又猛毛片免费看| av在线蜜桃| 亚洲专区中文字幕在线| 国产真人三级小视频在线观看| 宅男免费午夜| 欧美性猛交黑人性爽| 男女之事视频高清在线观看| 最近最新中文字幕大全免费视频| 国产精华一区二区三区| 夜夜夜夜夜久久久久| xxx96com| 国产成+人综合+亚洲专区| 不卡一级毛片| 搡老岳熟女国产| 亚洲熟妇中文字幕五十中出| 欧美日韩中文字幕国产精品一区二区三区| 日本a在线网址| 亚洲五月天丁香| 午夜激情福利司机影院| 成人鲁丝片一二三区免费| 亚洲在线自拍视频| 国产精品乱码一区二三区的特点| av在线蜜桃| 中文字幕人成人乱码亚洲影| 90打野战视频偷拍视频| 最近最新免费中文字幕在线| 欧美日韩一级在线毛片| 国产精品久久电影中文字幕| 小蜜桃在线观看免费完整版高清| 国内精品一区二区在线观看| 亚洲 国产 在线| 他把我摸到了高潮在线观看| 一个人免费在线观看电影 | 亚洲专区国产一区二区| 又粗又爽又猛毛片免费看| 级片在线观看| 久久伊人香网站| 亚洲午夜理论影院| 亚洲中文字幕日韩| 99re在线观看精品视频| 成人三级黄色视频| 在线观看66精品国产| 日本黄色视频三级网站网址| 亚洲国产欧洲综合997久久,| 欧美成人性av电影在线观看| 国产一级毛片七仙女欲春2| 一本一本综合久久| 国产成人福利小说| 久久国产乱子伦精品免费另类| 一级毛片精品| 欧美另类亚洲清纯唯美| 午夜激情欧美在线| 国产精品久久久久久人妻精品电影| 国产精品日韩av在线免费观看| 亚洲自拍偷在线| 成年女人看的毛片在线观看| 精品不卡国产一区二区三区| 亚洲色图av天堂| 国内精品一区二区在线观看| 一区二区三区国产精品乱码| 怎么达到女性高潮| 亚洲欧美日韩高清在线视频| 国产精品一区二区免费欧美| 99久久成人亚洲精品观看| 在线播放国产精品三级| www国产在线视频色| 精品午夜福利视频在线观看一区| 日本撒尿小便嘘嘘汇集6| 欧美色欧美亚洲另类二区| 国产麻豆成人av免费视频| 精品无人区乱码1区二区| 99久久国产精品久久久| 欧美成狂野欧美在线观看| 国产精华一区二区三区| 欧美在线一区亚洲| 男女视频在线观看网站免费| 热99re8久久精品国产| 热99在线观看视频| 19禁男女啪啪无遮挡网站| 黄色视频,在线免费观看| av视频在线观看入口| 麻豆成人午夜福利视频| 真人一进一出gif抽搐免费| 91麻豆av在线| 午夜精品一区二区三区免费看| 麻豆成人av在线观看| 亚洲乱码一区二区免费版| 亚洲18禁久久av| 午夜精品在线福利| 欧美日韩亚洲国产一区二区在线观看| 搡老熟女国产l中国老女人| 麻豆成人午夜福利视频| 两人在一起打扑克的视频| 亚洲欧美日韩高清在线视频| 亚洲五月天丁香| 男女之事视频高清在线观看| 国产不卡一卡二| 久久久国产欧美日韩av| 在线a可以看的网站| 久久伊人香网站| 国产69精品久久久久777片 | 亚洲午夜理论影院| 757午夜福利合集在线观看| 亚洲国产看品久久| 国产三级黄色录像| 在线观看免费午夜福利视频| 国产久久久一区二区三区| 欧美日韩瑟瑟在线播放| 国产蜜桃级精品一区二区三区| 午夜免费激情av| 九色国产91popny在线| 成人一区二区视频在线观看| 国产成人av激情在线播放| 日日摸夜夜添夜夜添小说| 99热精品在线国产| 一区二区三区高清视频在线| 午夜a级毛片| 国产精华一区二区三区| 高潮久久久久久久久久久不卡| 欧美日韩黄片免| 国产精品免费一区二区三区在线| 给我免费播放毛片高清在线观看| 高潮久久久久久久久久久不卡| 国产亚洲精品久久久久久毛片| 丰满的人妻完整版| 国产精华一区二区三区| 高潮久久久久久久久久久不卡| 动漫黄色视频在线观看| 亚洲精品国产精品久久久不卡| 毛片女人毛片| 久久久久久久久中文| 亚洲真实伦在线观看| 岛国在线免费视频观看| 很黄的视频免费| 午夜福利成人在线免费观看| 欧美乱码精品一区二区三区| 午夜福利在线观看免费完整高清在 | 色在线成人网| 俄罗斯特黄特色一大片| 免费看日本二区| 看免费av毛片| 女人高潮潮喷娇喘18禁视频| 亚洲精品美女久久av网站| 久久久久国产一级毛片高清牌| 淫妇啪啪啪对白视频| 亚洲精品在线观看二区| 叶爱在线成人免费视频播放| 亚洲午夜理论影院| 国产伦人伦偷精品视频| 99国产精品一区二区三区| av在线天堂中文字幕| 欧美日韩国产亚洲二区| 性色av乱码一区二区三区2| 欧美激情久久久久久爽电影| a级毛片在线看网站| 麻豆成人av在线观看| 后天国语完整版免费观看| 亚洲av免费在线观看| 淫秽高清视频在线观看| 久久久精品欧美日韩精品| 香蕉丝袜av| 高清毛片免费观看视频网站| 青草久久国产| 黄色片一级片一级黄色片| 国内精品久久久久久久电影| 99久久综合精品五月天人人| 女人被狂操c到高潮| 美女 人体艺术 gogo| 亚洲av日韩精品久久久久久密| 亚洲欧美精品综合一区二区三区| 日韩欧美在线二视频| 老汉色av国产亚洲站长工具| 亚洲性夜色夜夜综合| 高清在线国产一区| 国产黄a三级三级三级人| 午夜福利高清视频| 两人在一起打扑克的视频| 一区二区三区国产精品乱码| 欧美一区二区精品小视频在线| АⅤ资源中文在线天堂| 亚洲色图 男人天堂 中文字幕| 色av中文字幕| 久久香蕉精品热| 久久亚洲精品不卡| 日本一二三区视频观看| 久久久久久九九精品二区国产| 欧美国产日韩亚洲一区| 人妻久久中文字幕网| 精品一区二区三区av网在线观看| 国产成+人综合+亚洲专区| 网址你懂的国产日韩在线| 夜夜看夜夜爽夜夜摸| 国产精品久久久久久久电影 | 中文字幕久久专区| 亚洲天堂国产精品一区在线| 日韩欧美在线二视频| 观看免费一级毛片| 不卡av一区二区三区| 久久久久久久久免费视频了| 可以在线观看的亚洲视频| 91av网站免费观看| 在线观看一区二区三区| 欧美一级a爱片免费观看看| 亚洲成人久久爱视频| 免费观看人在逋| 丰满人妻一区二区三区视频av | 国产男靠女视频免费网站| 日韩人妻高清精品专区| 亚洲欧美精品综合一区二区三区| 久久国产精品人妻蜜桃| 色av中文字幕| 一个人免费在线观看电影 | 啦啦啦韩国在线观看视频| 人妻丰满熟妇av一区二区三区| 欧洲精品卡2卡3卡4卡5卡区| 亚洲精品456在线播放app | 久久久精品大字幕| 欧洲精品卡2卡3卡4卡5卡区| 亚洲avbb在线观看| 国产成人一区二区三区免费视频网站| 亚洲精品久久国产高清桃花| 岛国在线免费视频观看| cao死你这个sao货| 欧美成人性av电影在线观看| 极品教师在线免费播放| 亚洲九九香蕉| 国产精品自产拍在线观看55亚洲| 亚洲av免费在线观看| 88av欧美| 午夜精品久久久久久毛片777| 亚洲avbb在线观看| 天天躁日日操中文字幕| 床上黄色一级片| 五月伊人婷婷丁香| 欧美3d第一页| 色精品久久人妻99蜜桃| 欧美三级亚洲精品| 国产免费男女视频| 日本三级黄在线观看| 一级作爱视频免费观看| 国产黄色小视频在线观看| 深夜精品福利| av国产免费在线观看| 1000部很黄的大片| 国产免费男女视频| 欧美一级a爱片免费观看看| 亚洲黑人精品在线| 国内精品一区二区在线观看| 国产视频内射| 欧美成人性av电影在线观看| 十八禁网站免费在线| 麻豆久久精品国产亚洲av| 免费无遮挡裸体视频| 久99久视频精品免费| 啦啦啦免费观看视频1| 欧美不卡视频在线免费观看| 夜夜躁狠狠躁天天躁| 精品一区二区三区视频在线 | 亚洲美女视频黄频| 午夜福利高清视频| 亚洲人成网站在线播放欧美日韩| 真实男女啪啪啪动态图| 久久草成人影院| 国产欧美日韩一区二区精品| 99国产精品一区二区蜜桃av| 成人特级av手机在线观看| 国产一区二区在线观看日韩 | 最好的美女福利视频网| 国产一区二区三区视频了| 手机成人av网站| netflix在线观看网站| 亚洲精品乱码久久久v下载方式 | 免费看光身美女| 国产精品电影一区二区三区| 高清在线国产一区| 欧美日韩一级在线毛片| 男人的好看免费观看在线视频| 国产精品日韩av在线免费观看| 久久久成人免费电影| 美女午夜性视频免费| 久久人人精品亚洲av| 午夜视频精品福利| 不卡一级毛片| 一个人看的www免费观看视频| 香蕉久久夜色| 一级毛片精品| 丰满人妻熟妇乱又伦精品不卡| 色综合欧美亚洲国产小说| 午夜激情福利司机影院| 1024手机看黄色片| 99riav亚洲国产免费| 国产成年人精品一区二区| 午夜福利18| 亚洲国产精品合色在线| 国产综合懂色| 国产精品 欧美亚洲| 欧美3d第一页| 午夜精品在线福利| 99久久精品热视频| 国产精品 国内视频| 九九热线精品视视频播放| 国产又色又爽无遮挡免费看| 三级国产精品欧美在线观看 | 男人舔奶头视频| 身体一侧抽搐| 99热这里只有是精品50| www.熟女人妻精品国产| 亚洲成人久久爱视频| 日日干狠狠操夜夜爽| 欧美性猛交╳xxx乱大交人| АⅤ资源中文在线天堂| 国产精品一区二区三区四区久久| 香蕉av资源在线| 亚洲精品美女久久久久99蜜臀| 男人的好看免费观看在线视频| 亚洲第一欧美日韩一区二区三区| 久久人妻av系列| 不卡av一区二区三区| 亚洲一区二区三区色噜噜| 色噜噜av男人的天堂激情| 亚洲性夜色夜夜综合| 日本免费a在线| 欧美zozozo另类| 桃色一区二区三区在线观看| 中文字幕久久专区| 韩国av一区二区三区四区| 国产极品精品免费视频能看的| 国产精品香港三级国产av潘金莲| a级毛片a级免费在线| 18美女黄网站色大片免费观看| 国产成人欧美在线观看| 成人欧美大片| 久久精品91蜜桃| 黄色女人牲交| 国产精品日韩av在线免费观看| 午夜福利视频1000在线观看| 99国产精品一区二区三区| 国产成人福利小说| 制服人妻中文乱码| 亚洲精品一区av在线观看| 国产高清有码在线观看视频| 最新美女视频免费是黄的| 一级黄色大片毛片| 国内揄拍国产精品人妻在线| 一级a爱片免费观看的视频| 国产黄a三级三级三级人| 日本 欧美在线| 精品午夜福利视频在线观看一区| 国产1区2区3区精品| 亚洲精品在线美女| 国产主播在线观看一区二区| 国产精品自产拍在线观看55亚洲| 别揉我奶头~嗯~啊~动态视频| 国产精品 欧美亚洲| 久久国产精品人妻蜜桃| 久久国产精品影院| 久久这里只有精品中国| 成年女人永久免费观看视频| 国产精品久久视频播放| 亚洲国产中文字幕在线视频| 久久精品亚洲精品国产色婷小说| 欧美又色又爽又黄视频| 日韩三级视频一区二区三区| 亚洲国产精品成人综合色| 欧美成狂野欧美在线观看| 国产精品 欧美亚洲| 激情在线观看视频在线高清| 丁香欧美五月| 国产精品自产拍在线观看55亚洲| 成年女人永久免费观看视频| 在线观看日韩欧美| 啦啦啦观看免费观看视频高清| 一区二区三区激情视频| 国产又色又爽无遮挡免费看| 国产高清激情床上av| 母亲3免费完整高清在线观看| 在线免费观看不下载黄p国产 | 性色av乱码一区二区三区2| 日本撒尿小便嘘嘘汇集6| 99久久精品一区二区三区| 我的老师免费观看完整版| www.自偷自拍.com| 久久久久久大精品| 免费在线观看成人毛片| 色播亚洲综合网| 美女扒开内裤让男人捅视频| 网址你懂的国产日韩在线| 九九久久精品国产亚洲av麻豆 | 欧美一级毛片孕妇| 久久人人精品亚洲av| 国产单亲对白刺激| 2021天堂中文幕一二区在线观| 一进一出抽搐动态| 好看av亚洲va欧美ⅴa在| 久久国产乱子伦精品免费另类| av在线天堂中文字幕| 国内精品久久久久久久电影| 国产熟女xx| 美女大奶头视频| 啪啪无遮挡十八禁网站| 国产高清videossex| 国产极品精品免费视频能看的| 亚洲熟妇中文字幕五十中出| 18美女黄网站色大片免费观看| 国产高清videossex| 国产一区二区三区视频了| 亚洲精品国产精品久久久不卡| 国产精品美女特级片免费视频播放器 | 熟女少妇亚洲综合色aaa.| 精品福利观看| 国产免费av片在线观看野外av| 久久久久国内视频| 淫妇啪啪啪对白视频| 男人舔奶头视频| 欧美成人性av电影在线观看| 国产精品 欧美亚洲| 亚洲电影在线观看av| 欧美日韩中文字幕国产精品一区二区三区| 国产真实乱freesex| 欧美高清成人免费视频www| 成人三级黄色视频| 色播亚洲综合网| 亚洲在线自拍视频| 精品人妻1区二区| 岛国在线免费视频观看| 国产综合懂色| 欧美+亚洲+日韩+国产| 天天添夜夜摸| 欧洲精品卡2卡3卡4卡5卡区| 999久久久精品免费观看国产| 长腿黑丝高跟| 欧美3d第一页| 日本成人三级电影网站| 天堂√8在线中文| 久久久国产成人精品二区| 国产乱人伦免费视频| 国产欧美日韩精品亚洲av| 国产97色在线日韩免费| 国产毛片a区久久久久| 精品乱码久久久久久99久播| 精品无人区乱码1区二区| 曰老女人黄片| 国产午夜精品久久久久久| 欧美丝袜亚洲另类 | 国产亚洲精品久久久com| 国产成人精品久久二区二区免费| 亚洲精品一区av在线观看| 日本a在线网址| 天堂影院成人在线观看| 不卡一级毛片| 日日干狠狠操夜夜爽| 高潮久久久久久久久久久不卡| 精品久久久久久,| 青草久久国产| 亚洲va日本ⅴa欧美va伊人久久| 亚洲,欧美精品.| 亚洲欧美激情综合另类| 国产精品亚洲一级av第二区| 亚洲狠狠婷婷综合久久图片| 欧美国产日韩亚洲一区| 神马国产精品三级电影在线观看| 两个人视频免费观看高清| 久久久久久久久免费视频了| 2021天堂中文幕一二区在线观| 18美女黄网站色大片免费观看| 国产亚洲精品综合一区在线观看| 国产精品av视频在线免费观看| 国产真实乱freesex| 国产精品98久久久久久宅男小说| 成人特级av手机在线观看| 欧美中文综合在线视频| 此物有八面人人有两片| 欧美激情在线99| 99久久综合精品五月天人人| 久久精品91蜜桃| 国产精品日韩av在线免费观看| 精品99又大又爽又粗少妇毛片 | 99热6这里只有精品| 日本成人三级电影网站| 日本免费a在线| 亚洲av中文字字幕乱码综合| 久久精品国产清高在天天线| 狂野欧美激情性xxxx| 国产精品一区二区精品视频观看| 久久久国产精品麻豆| 国产 一区 欧美 日韩| 啦啦啦免费观看视频1| 美女 人体艺术 gogo| 精品一区二区三区视频在线观看免费| 两个人视频免费观看高清| 国产成人精品久久二区二区免费| 亚洲,欧美精品.| 午夜成年电影在线免费观看| 亚洲乱码一区二区免费版| 国产精品自产拍在线观看55亚洲| av黄色大香蕉| 久久人人精品亚洲av| 午夜福利在线观看吧| 又黄又粗又硬又大视频| 亚洲精品美女久久久久99蜜臀| 色尼玛亚洲综合影院| 国产精品久久久av美女十八| 久99久视频精品免费| 亚洲精品在线美女| 成人无遮挡网站| 91九色精品人成在线观看| 18禁美女被吸乳视频| 黄色成人免费大全| www.www免费av|