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

    Regulating Electron Transport Band Gaps of Bovine Serum Albumin by Binding Hemin

    2021-04-02 02:23:54ChuanliWuWenhuiLiangJingjingFanYuxianCaoPingWuChenxinCai
    物理化學學報 2021年3期

    Chuanli Wu , Wenhui Liang , Jingjing Fan , Yuxian Cao , Ping Wu , Chenxin Cai ,*

    1 Jiangsu Key Laboratory of New Power Batteries, Jiangsu Collaborative Innovation Center of Biomedical Functional Materials,College of Chemistry and Materials Science, Nanjing Normal University, Nanjing 210023, China.

    2 College of Education, Nanjing Normal University, Nanjing 210023, China.

    Abstract: The small size (nanoscale) of proteins and their favorable electron transport (ETp) properties make them suitable for various types of bioelectronic devices and offer a solution for miniaturizing these devices to nanoscale dimensions.The performance of protein-based devices is predominantly affected by the ETp property of the proteins,which is largely determined by the band gaps of the proteins, i.e., the energy difference between the conduction band (CB) and valence band (VB).Regulating the protein ETp band gaps to appropriate values is experimentally demanding and hence remains a significant challenge.This study reports a facile method for modulating the ETp band gaps of bovine serum albumin (BSA), via its binding with a foreign molecule, hemin.The formation of the hemin-BSA complex was initially confirmed by theoretical simulation (molecular docking) and experimental characterization(fluorescence and absorption spectra), which indicated that the hemin is positioned inside a hydrophobic cavity formed by hydrophobic amino acid residues and near Trp213, at subdomain IIA of BSA, with no significant effects on the structure of BSA.Circular dichroism (CD) spectra indicated that the BSA conformation remains essentially unaltered following the formation of the hemin-BSA complex, as the helicities of the free BSA (non-binding) and the hemin-BSA complex were estimated to be 66% and 65%, respectively.Moreover, this structural conformation remains preserved after the hemin-BSA complex is immobilized on the Au substrate surface.The hemin-BSA complex is immobilized onto the Au substrate surface along a single orientation, via the -SH group of Cys34 on the protein surface.Atomic force microscopy (AFM) images indicate that hemin-BSA forms a dense layer on the surface of the Au substrate with a lateral size of ~3.2-3.7 nm, which is equivalent to the actual size of BSA, ~4.0 nm × 4.0 nm × 14.0 nm.The current-voltage (I-V) responses were measured using eutectic gallium-indium (EGaIn) as the top electrode and an Au film as the substrate electrode, revealing that the ETp processes of BSA and hemin-BSA on the Au surface have distinct semiconducting characteristics.The CB and VB were estimated by analysis of the differential conductance spectra, and for the free BSA, they were ~0.75 ± 0.04 and ~-0.75 ± 0.08 eV, respectively, being equally distributed around the Fermi level (0 eV), with a band gap of ~1.50 ± 0.05 eV.Following hemin binding, the CB (~0.64 ± 0.06 eV) and VB (~ -0.29 ± 0.07 eV) of the protein were closer to the Fermi level, resulting in a band gap of ~0.93 ± 0.05 eV.These results demonstrated that hemin molecules can effectively regulate ETp characteristics and the transport band gap of BSA.This methodology may provide a general approach for tuning protein ETp band gaps, enabling broad variability by the preselection of binding molecules.The protein and foreign molecule complex may further serve as a suitable material for configuring nanoscale solid-state bioelectronic devices.

    Key Words: Protein electron transport; Protein-based electronic device; Bovine serum albumin; Hemin

    1 Introduction

    The use of protein molecules as electronic components offers a solution to miniaturize electronic devices to nanoscale dimensions, which is essential for the future of the electronics industry because these types of circuits would have the benefits of reduced energy consumption and enhanced device speed and functionality1-5.The small size (nano-scale) of proteins and their effective electron transport (ETp) make them suitable for many types of bioelectronic devices1,6-8, including biosensors for biological molecule monitoring8-15, biomolecular transistors for data storage16, and biofuel cells for microscale implantable17,18.The properties of these devices are mainly affected by the ETp performance of the protein, which is largely determined by the protein’s band gaps,i.e., the energy difference between the conduction band (CB) and valence band (VB).From a practical and theoretical point of view, it is necessary to adjust the protein ETp band gaps to appropriate values to improve the performance of protein-related electronic devices and expand their application fields.

    Although each protein has an inherent ETp and characteristic band gap, they can still be regulated to some extent by changing parameters that modulate ETp of the proteins, such as the conformational structure and core compositions.Sheves and Cahenetal.reported that the ETp of human serum albumin(HSA) is improved by noncovalently binding of retinoate19and hemin20.Richter and Cohenetal.21,22showed that the addition of a C60molecule can change the gate voltage of BSA-transistors.However, these groups did not explore how the ETp band gap of proteins changes after binding with foreign molecules.We studied the regulation of protein (enzyme) conformation to regulate the ETp characteristics of glucose oxidase and myoglobin23,24.However, the band gap and ETp features of these denatured proteins are not reflective of the intrinsic properties of the proteins.Mukhopadhyay and coworkers regulated the band gaps of the ferritin protein from ~0 to 1.17 V by replacing natural iron cores with other metal cores such as Mn, Cu, Ni, Co, and Zn25,26.Although core reconstitution methods can preserve the native conformation of proteins, they require complicated procedures for operation.Therefore, it is necessary to study a method for regulating the characteristics of ETp and to adjust the band gap of protein easily and conveniently.

    Our previous work demonstrated that the ETp band gap of BSA was reduced by noncovalently binding it with cyanocobalamin (Vb12)27.However, the Vb12-bound BSA was immobilized on the surface of a silicon (111) substrate through the reaction of the -COOH groups of BSA and the -NH2group of the 3-aminopropyl trimethoxysilane (3-APTMS)modified on the Si surface.Because there is a high number of-COOH groups on the BSA surface, located at different positions28, the orientation of the immobilized BSA molecules on the Si surface could not be controlled well, which caused the obtainedI–Vcharacteristics of BSA to be an average of the ETp from different directions.Because protein orientation on the electrode surface plays an important role in ETp4,8,29,30, this effort sought to control the immobilized BSA on the substrate surface with the same orientation through the formation of an Au-S bond between the -SH groups of the cysteines (Cys34)on the BSA molecule and the Au substrate.BSA molecule contains 35 cysteines, however, 34 of them are linked to form 17 S-S bonds embedded in the inner cavity of the protein molecule31and are not available for the formation of Au-S bonds.Cys34 is the only cysteine residue that is located at the surface of the protein and does not participate in a disulfide bridge32.The surface-located Cys34 is favorable for forming an Au-S bond,leading to the protein film formation with controlled orientation.Moreover, we replaced the Hg top electrode with eutectic gallium-indium (EGaIn) to avoid Hg electrode amalgamating with Au film33.The work function of EGaIn (~4.2 eV) is close to that of Hg (~4.5 eV), but EGaIn is nontoxic and does not alloy with Au34.Therefore, the obtainedI–Vfeatures from the EGaIn top electrode and Au substrate electrode can reflect the ETp of BSA at one particular orientation.

    2 Experimental

    2.1 Preparation and characterization of the hemin-BSA

    We chose hemin as the binding molecule.Hemin is the prosthetic group of many porphyrin-containing metalloproteins and serves as an electron sink during electron transfer or redox chemistry mediated by those proteins.Moreover, it also binds noncovalently to HSA35, and the structure of the hemin-HSA complex has been determined by X-ray diffraction36.Based on the structural similarity between HSA and BSA, hemin is also expected to bind BSA.For preparation of the hemin-BSA, BSA(5 μmol·L–1, Sigma-Aldrich) was allowed to interact with hemin(5 μmol·L–1) in phosphate buffer solution (PBS, 10 mmol·L–1,pH 7.4) for more than 4 h.The unbound hemin molecules were removed by dialyzing with 10 kDa cutoff against PBS buffer for at least 48 h to remove these unbound hemin molecules; the PBS buffer was changed every 12 h.The formation of the hemin-BSA was characterized with fluorescence spectra and absorption spectra, which were recorded on a Cary Eclipse fluorescence spectrophotometer (Varian) and a Cary 5000 UV-Vis-NIR spectrometer (Varian), respectively.BSA was excited at 295 nm in order to minimize the contribution from tyrosine.To evaluate the structural and conformational changes of BSA caused by binding with hemin, circular dichroism (CD) measurements were performed on a circular dichroism spectropolarimeter(Chriascan, Applied Photophysics).The Raman spectra of the hemin-BSA complex were measured on a Labram HR 800 microspectrometer (Jobin Yvon, France) with an excitation source of 785 nm and the laser power was 50 mW.All spectra were baseline corrected.

    2.2 Immobilization of the hemin-BSA onto gold substrate

    To measure theI–Vresponse, the hemin-BSA was immobilized on an Au film substrate surface.For preparation of the Au film substrate, the indium tin oxide (ITO) glass plate (20 mm × 20 mm, Changzhou Best Photoelectronics Technology Co.Ltd., China) were first cleaned by dipping it into a Piranha solution (a mixture of the 30% (volume percentage) H2O2and 70% H2SO4) at 80 ℃ for 30 min (Warning: the Piranha solution reacts strongly with organic compounds and should not be stored in closed containers.It must be handled with extreme caution!) and followed by rinsing with deionized water thoroughly and absolute ethanol in an ultrasonic bath and drying under a steam of nitrogen.Before sputtering Au, a layer of chromium (Cr) was first sputtered on glass slide to promote adhesion of Au (step a, Fig.1).Au film was then sputtered onto the ITO glass surface (step b, Fig.1).The Au and the Cr films had a thickness of 50 and 10 nm, respectively.The freshly prepared Au film surface was cleaned sequentially with deionized water, acetone and then annealed to obtain a relatively flat Au surface.Finally, the Au film substrate was immersed into a hemin-BSA solution (5 μmol·L–1, in a sealed vial) for 6 h for the hemin-BSA immobilization (step c, Fig.1), which proceedsviathe reaction of the -SH of the Cys34 on the surface of BSA and the Au film to form Au-S bonds, yielding a protein layer with a thickness of ~3.2 to 3.7 nm.After that, the hemin-BSA-coated Au substrate was washed thoroughly with distilled water to remove those unbound hemin-BSA molecules.The immobilization processes and the morphology of the protein on the Au surface were observed by atomic force microscopy(AFM) with a Nanoscope IIIa scanning probe microscope(Digital Instruments) in tapping mode (step d, Fig.1).

    Fig.1 Illustration of the formation of the Au film on the ITO surface,the immobilization of hemin-BSA on the Au film surface, and configuration of the fabricated device for the I–V measurements.

    2.3 I–V response measurements

    I–Vresponses were recorded with EGaIn severed as a topcontact electrode using an Autolab PGSTAT302N electrochemical station (Metrohm) at an ambient temperature by employing a two-electrode system in the bias range of -2.5 to+2.5 V (step e, f, Fig.1).The EGaIn tip (the preparation of the EGaIn tip has been depicted in Supporting Information) was putted in contact with the hemin-BSA layer using the micromanipulator.The EGaIn tip converges with its reflection on the surface of the gold substrate indicated their contacts.

    3 Results and discussion

    Before recording theI–Vresponses, the formation of the hemin-BSA was confirmed.Molecular docking was first performed to theoretically predict the possibility of BSA interacting with hemin and the binding sites of hemin in BSA(the details for simulation are depicted in Supporting Information).Analysis of the docking structures revealed that hemin binding produces a convergent structure only in the IIA subdomain (Fig.2a).This binding site of hemin in BSA differs from the HSA molecule on which the hemin is believed to bind within subdomain IB36.After binding with BSA, the hemin is located inside a hydrophobic cavity (Fig.2b), formed by hydrophobic amino acid residues, such as Trp213 (0.23 nm, the nearest distance of hemin to amino residue), Glu152 (0.26 nm),Ser191 (0.24 nm), Arg194, 198, 217 (0.23, 0.29, 0.22 nm),His287 (0.24 nm), Cys288 (0.30 nm), Val292 (0.28 nm), and Lys294 (0.19 nm; Fig.2c), indicating that van der Waals interactions dominate the binding.Even though the hemin in the hydrophobic cavity is positioned near to the Trp213 (Fig.2d),the hemin binding does not significantly change the basic environment of Trp213.

    Fig.2 (a) The most stable docking structure of the hemin molecule at subdomain IIA of the BSA molecule.The molecular structure of hemin is also depicted.(b) Bound hemin positioning inside the hydrophobic cavity composed of the hydrophobic amino acid residues in subdomain IIA of the BSA.Red represents negative charge, and blue represents positive charge.The surface potential ranges from -0.1e to 0.1e.(c) Illustration of the amino acid residues in the vicinity of the bound hemin.The hydrogen bonds are shown by dashed green lines.(d) Illustration of the relative position of the bound hemin in subdomain IIA of the BSA to the Trp213 residue.

    The formation of hemin-BSA was then verifiedviaabsorption spectra, which showed that both hemin and BSA have significant absorptions in the UV-Vis region (Fig.3a).BSA exhibits two absorption peaks, located at ~215 and 278 nm (a1, Fig.3a).The peak at ~215 nm is the characteristic absorption of the protein framework conformation and corresponds to the peptide bond,and the peak at ~278 nm is from the absorption of the aromatic amino acids, such as Trp, Tyr or Phe, in the protein molecule37.Hemin displayed a broad characteristic absorption with a center at ~385 nm (a2, Fig.3a), corresponded to the porphyrin Soret band, as well as the Q-bands38.Hemin binding caused to an increase in the intensity of the peak at ~278 nm that is slightly blue shifted (~3 nm; a3, Fig.3a).Moreover, a new strong absorption peak at ~406 nm was observed.This new peak should be resulted from the Soret band of the hemin molecule bound with BSA, and has a significant red shift comparing with that for the free hemin (not binding with BSA).These findings suggest the formation of hemin-BSA, which was consistent with the theoretical predictions, and also imply that the protein molecule in hemin-BSA complex essentially retains its native structure.

    Fig.3 (a) Absorption spectra of free BSA (a1, 5 μmol·L-1), hemin (a2, 5 μmol·L-1), and hemin-BSA (a3, 5 μmol·L-1) in PBS, and the hemin-BSA immobilized on the Au substrate surface (a4).(b) Fluorescence emission spectra of BSA (b1, 5 μmol·L–1), hemin-BSA (b2,5 μmol·L-1), and hemin alone (b3, 5 μmol·L-1) in phosphate buffer solution (PBS, 10 mmol·L-1, pH 7.4).The spectra were recorded under excitation of 295 nm.(c) CD spectra of free BSA (c1, 1 μmol·L-1), hemin-BSA (c2, 1 μmol·L-1), and hemin (c3, 1 μmol·L-1) in PBS, and the hemin-BSA immobilized on the Au substrate surface (c4).

    The formation of hemin-BSA was further supported by fluorescence spectroscopy.Excitation of the unbound BSA (free BSA) at 295 nm resulted in a strong fluorescence peak at ~350 nm (b1, Fig.3b), which is the characteristic fluorescence peak of BSA, corresponding to the Trp in a hydrophobic environment in the protein.The fluorescence quenching of BSA results from binding with hemin (b2, Fig.3b) at a 1 : 1 ratio of BSA to hemin.In the studied region, we did not observe any fluorescence signal of the hemin molecule (b3, Fig.3b).These results indicate the formation of hemin-BSA and were consistent with the results predicted by molecular docking and absorption spectrum results again.We further conjecture that, as expected from the simulation by molecular docking (Fig.2d), hemin is in the vicinity of Trp213 in the protein because the fluorescence of the BSA is mainly attributed to the Trp213.In addition, the fluorescence peak position of the protein almost remains unchanged after the formation of the hemin-BSA (b1, b2, Fig.3b), indicating that BSA in hemin-BSA complex essentially retains its native structure.

    After theoretically and experimentally confirming the formation of hemin-BSA, we carried out CD spectral analysis to gain a better understanding of the structural and conformational changes in BSA induced by hemin binding.The CD spectra of BSA in the far-UV region showed two bands at ~208 and 222 nm (c1, Fig.3c), which are the characteristic of the α-helix structure39,40, corresponding to theπ-π* andn-π* transitions,respectively, implying that the native BSA is mainly characterized by an α-helix conformation, consistent with the previous reported results (the α-helix structure of the BSA is~67%)41.Compared with free BSA, the intensities of the two negative minima remained to be mostly the same for hemin-BSA(c2, Fig.3c), indicating that the helical structure of the protein remained stable after hemin was bound into the protein (in this far-UV region, hemin itself did not showed any CD response.See c3, Fig.3c).CD spectral analysis using the CDNN program(version 2.0) showed that the helicity of BSA in its native and hemin-BSA states was estimated to be 66 (1.0 μmol·L-1, 25 ℃)and 65%, respectively, confirming that BSA conformation remained essentially invariant after forming the hemin-BSA.

    Next, immobilization of hemin-BSA on an Au film substrate and its ETp characteristics were investigated.Due to the presence of the -SH group of Cys34 on the surface of BSA,hemin-BSA can be immobilized on the Au substrate by an Au-S bond.The formation of Au-S bonds between the Au substrate and hemin-BSA are proved by the band of 270 cm-142,which is the characteristic Raman peak of the Au-S bond.Moreover, the “fingerprint region” between 200 and 2000 cm-1of the Raman spectrum provides secondary and tertiary structural information or the microenvironment of each functional group: the amide I band at 1646 cm-1, the CH2stretching band at 1448 cm-1, and the bands of hydrophobic residues such as Phe at 1006 cm-1and Tyr at 850 and 830 cm-1(curve b, Fig.4).These Raman also observed in the spectrum of the hemin-BSA in solution (curve a, Fig.4) except the peak of the Au-S bond, indicating the immobilization of the hemin-BSA on the Au substrateviathe Au-S bond.

    Fig.4 The Raman spectra of hemin-BSA in solution (a) and immobilized on Au substrate (b).

    The immobilization of the hemin-BSA on the Au substrate can also observed from the atomic force microscopy (AFM) images,which showed hemin-BSA could form a dense layer on the surface of the Au substrate (Fig.5a,b).These features were not observed on the surface of a pure Au substrate (without immobilization of proteins, Fig.5c,d), indicating that the interaction of the protein with the substrate led to the absorption of the molecule.The lateral dimension analysis based on the AFM image revealed that on the top Au electrode-sample convolution, the lateral size the protein was ~3.2-3.7 nm, which is equivalent to the actual size of ~4.0 nm × 4.0 nm × 14.0 nm7.These results demonstrate that hemin-BSA has been immobilized on the Au surface successfully.

    Fig.5 (a) AFM topography image of the hemin-BSA immobilized on the Au substrate surface and (b) its corresponding 3D image.(c) AFM topography image of Au substrate and (d) its 3D image.

    Before performing the ETp experiments on hemin-BSA, we investigated whether the conformation structures of BSA immobilized on an Au substrate remained by measuring its absorption and CD spectra and comparing them with values recorded in solution.The resulting absorption spectrum showed peaks at ~276 and 406 nm (a4, Fig.3a), similar to the absorption of hemin-BSA in solution (a3, Fig.3a).Two negative bands appeared in the CD spectrum of the immobilized hemin-BSA (c4,Fig.3c) was also similar to those for the hemin-BSA in solution(c2, Fig.3c).Based on the absorption and CD spectral results, we conclude that the native conformation of BSA was retained after hemin-BSA was immobilized on the Au substrate surface.

    The ETp characteristics of the BSA and hemin-BSA were investigated by recording theirI–Vresponses.According to theI–Vresponses (Fig.6a), the ETp processes of BSA and the hemin-BSA on the Au substrate surface have obvious semiconductor characteristics.I–Vresponses were derived from the protein ETp rather than its redox reaction because BSA (with no redox center) is electrochemically inactive.The current corresponding to pure hemin immobilized on the surface of the Au substrate (Fig.6b) was higher than the measured ETp current with hemin-BSA (Fig.6a).This means that hemin-BSA formed a less electronically conductive film on the Au substrate surface,consistent with our AFM observations.Compared to free BSA,the ETp current of hemin-BSA was higher by nearly an order of magnitude (Fig.6c), indicating the bind of hemin in BSA can significantly improve the ETp.

    Fig.6 (a) I–V responses of free BSA, hemin-BSA, and apo-hemin-BSA immobilized on the Au substrate surface.(b) I–V responses of the hemin immobilized on the Au substrate surface.I–V responses presented here are the average of 15 different measurements on three separate samples (five measurements at different points on each sample).(c) Semi-logarithmic plot of the I–V response curves for free BSA and hemin-BSA immobilized on the Au substrate surface.(d and e) Differential conductance spectra of free BSA (d) and hemin-BSA (e).The vertical lines in (d) and (e) represent Fermi energies, which are aligned to 0 V.(f) Schematic illustration of the CB and VB band levels and the heights of the band gaps of BSA and hemin-BSA.The values of the CB and VB presented here are average values of the 15 measurements.

    Hemin has a Fe ion redox center (Fig.2a).To check whether the enhancement of the ETp current of the hemin-BSA is caused by the Fe ion redox center, we employed the apo-hemin (i.e.,without Fe ion redox center) to bind the BSA, forming apohemin-BSA.The recordedI–Vresponse of the apo-hemin-BSA is depicted in Fig.6a, which indicates the ETp current of the apohemin-BSA is almost in the same as that of hemin-BSA,implying that the Fe ion center of the hemin hardly affects the ETp current of the BSA.This result suggests that the enhancement of the ETp current of bound BSA (Fig.6c) is not caused by the Fe ion redox center of the hemin, and further verify that theI–Vresponses in Fig.6a result from the ETp of protein,not from the electrochemical redox reaction.

    Finally, we estimated the changes of ETp band gaps of BSA after binding with hemin.The differential conductance spectra,i.e., differential conductivity as functions of the sample bias voltage (dI/dVvsV), were plotted (Fig.6d,e).In these plots, the bias voltage of zero (V= 0) corresponds to the Fermi energy,negative (V< 0) and positive (V> 0) voltages indicate energies below and above the Fermi energy, respectively.According to WKB (Wentzel Kramers Brillouin) approximation43-45, dI/dVis a convolution of the surface DOS,i.e., is proportional to the density of states (DOS) of the semiconductor, and thus the DOS can be obtained from the dI/dVvsVcurves.The dI/dVbelow the Fermi energy indicates the occupied DOS (corresponding to VB)and the dI/dVabove the Fermi energy indicates the unoccupied DOS (corresponding to CB).Therefore, these plots can also be used to predict the CB and VB levels of BSA in hemin-BSA because the plots reflect the characteristics of each protein.The peaks in the negative voltage region of the DOS spectra can be expressed as the VB of the semiconductor, similarly, the position of CB is provided by the peak at the positive voltages.In this way, band gap values are calculated utilizing the position of CB and VB edges, and the influence of binding molecules on the protein can be determined based on the difference in band gaps.Fig.6d shows that the CB and VB of the BSA is ~0.75 ± 0.04(average of 15 measurements) and ~ -0.75 ± 0.08 eV (Fig.6d),respectively, and is equally distributed around 0 eV.The band gap is ~1.50 ± 0.05 eV.After binding with hemin, the protein CB(~0.64 ± 0.06 eV) and VB (~ -0.29 ± 0.07 eV) are closer to the Fermi levels (Fig.6e), resulting in a band gap of ~0.93 ± 0.05 eV (Fig.6f).The decrease in the band gap resulted to an increase in ETp, demonstrating that hemin molecules can effectively regulate the ETp characteristics and the transport band gap of BSA.In addition, the protein ETp band gap generally fell within the range of 0.8-2.6 eV25, which is effective for common semiconductor technologies (for example,p–njunction-based transistors), suggesting that controlling the electronic band gap of the protein may be further applied to suitable materials for constructing novel bioelectronic devices at nanoscale levels.

    4 Conclusions

    In conclusion, we developed a simple method for tuning the ETp band gap of BSA by binding with hemin.Both molecular docking and absorption spectroscopy results demonstrated that conformational changes caused by the binding of hemin to the IIA subdomain of BSA were negligible.I–Vcurves indicated that binding of BSA with hemin reduced the protein band gaps and significantly increased the ETp.These results could easily be extended to tune the ETp band gaps of other proteins by selecting the foreign binding molecules.More importantly, these protein bound with foreign molecules have possibilities to be used as ETp mediators in future nanoscale bioelectronic devices.

    Supporting Information:available free of chargeviathe internet at http://www.whxb.pku.edu.cn.

    一本大道久久a久久精品| 久久99热这里只频精品6学生| 制服诱惑二区| 亚洲九九香蕉| 99re6热这里在线精品视频| 男女午夜视频在线观看| 国产精品久久久久久精品古装| 一本久久精品| 一级毛片精品| 久久久国产一区二区| 日韩人妻精品一区2区三区| 久久av网站| 80岁老熟妇乱子伦牲交| 99国产精品一区二区蜜桃av | 王馨瑶露胸无遮挡在线观看| 人妻久久中文字幕网| 三级毛片av免费| 人成视频在线观看免费观看| 90打野战视频偷拍视频| kizo精华| 亚洲精品久久成人aⅴ小说| 久久人人97超碰香蕉20202| 老汉色∧v一级毛片| 午夜两性在线视频| 欧美精品啪啪一区二区三区| 亚洲av片天天在线观看| 久久 成人 亚洲| 我的亚洲天堂| 正在播放国产对白刺激| 成人三级做爰电影| 高清av免费在线| 国产亚洲欧美精品永久| 欧美另类亚洲清纯唯美| 午夜成年电影在线免费观看| 成人国语在线视频| 国产成人一区二区三区免费视频网站| 别揉我奶头~嗯~啊~动态视频| 成年人黄色毛片网站| 久久精品亚洲av国产电影网| 一级片免费观看大全| 精品福利永久在线观看| 欧美精品一区二区大全| 久久99一区二区三区| 亚洲av第一区精品v没综合| 性色av乱码一区二区三区2| 精品国产一区二区三区久久久樱花| 亚洲精品中文字幕一二三四区 | 国产精品久久电影中文字幕 | 18在线观看网站| tocl精华| 人妻久久中文字幕网| 美女福利国产在线| 在线观看66精品国产| 18禁观看日本| av欧美777| 女人久久www免费人成看片| 肉色欧美久久久久久久蜜桃| 午夜福利在线免费观看网站| 美女高潮到喷水免费观看| 国产一卡二卡三卡精品| 久久亚洲精品不卡| 精品人妻在线不人妻| 精品国产一区二区三区四区第35| 国产亚洲精品第一综合不卡| 男人舔女人的私密视频| 少妇粗大呻吟视频| 精品一区二区三区视频在线观看免费 | 成在线人永久免费视频| 丝袜喷水一区| 飞空精品影院首页| 97人妻天天添夜夜摸| 亚洲一区中文字幕在线| 国产在视频线精品| 国产精品国产高清国产av | 国产片内射在线| 亚洲欧洲日产国产| 99精品在免费线老司机午夜| 热re99久久国产66热| 一二三四在线观看免费中文在| 美女国产高潮福利片在线看| 在线十欧美十亚洲十日本专区| 夜夜骑夜夜射夜夜干| 女性被躁到高潮视频| 热99国产精品久久久久久7| 757午夜福利合集在线观看| 中文字幕精品免费在线观看视频| 咕卡用的链子| 怎么达到女性高潮| 久久热在线av| 淫妇啪啪啪对白视频| 欧美久久黑人一区二区| 欧美精品av麻豆av| 亚洲欧洲日产国产| 50天的宝宝边吃奶边哭怎么回事| 亚洲精品久久午夜乱码| 亚洲情色 制服丝袜| 9色porny在线观看| 成人18禁高潮啪啪吃奶动态图| 两个人看的免费小视频| 一级毛片精品| 色综合婷婷激情| 久久精品亚洲精品国产色婷小说| 国产在视频线精品| 欧美激情久久久久久爽电影 | 欧美国产精品va在线观看不卡| 热re99久久国产66热| 女人爽到高潮嗷嗷叫在线视频| 免费在线观看视频国产中文字幕亚洲| 亚洲午夜理论影院| 国产成+人综合+亚洲专区| 老司机在亚洲福利影院| 黑人巨大精品欧美一区二区mp4| 亚洲熟女毛片儿| 国产无遮挡羞羞视频在线观看| 欧美 日韩 精品 国产| videos熟女内射| 天天影视国产精品| 免费少妇av软件| 18禁国产床啪视频网站| 国产在线免费精品| 欧美激情高清一区二区三区| 一级片免费观看大全| 黄色片一级片一级黄色片| 国产一区二区在线观看av| 久久精品国产99精品国产亚洲性色 | 久久婷婷成人综合色麻豆| 欧美精品高潮呻吟av久久| 美女高潮喷水抽搐中文字幕| videosex国产| 亚洲中文av在线| 亚洲九九香蕉| 最新的欧美精品一区二区| 久久午夜亚洲精品久久| 国产精品美女特级片免费视频播放器 | 欧美日韩精品网址| 97人妻天天添夜夜摸| 久久婷婷成人综合色麻豆| 亚洲av成人不卡在线观看播放网| 日本五十路高清| 国产精品免费一区二区三区在线 | 五月天丁香电影| 午夜激情av网站| 国产视频一区二区在线看| 久久人人97超碰香蕉20202| 日本五十路高清| 午夜福利欧美成人| 91精品三级在线观看| 亚洲九九香蕉| 欧美日韩视频精品一区| 久久久久久亚洲精品国产蜜桃av| 色综合欧美亚洲国产小说| 看免费av毛片| 亚洲国产欧美网| 亚洲va日本ⅴa欧美va伊人久久| 超碰97精品在线观看| 亚洲av成人不卡在线观看播放网| 97在线人人人人妻| 在线观看一区二区三区激情| 少妇被粗大的猛进出69影院| 久久天堂一区二区三区四区| 91成人精品电影| 精品一区二区三卡| 久久国产精品影院| 亚洲精品中文字幕在线视频| 在线亚洲精品国产二区图片欧美| 日韩制服丝袜自拍偷拍| 女人高潮潮喷娇喘18禁视频| 男女无遮挡免费网站观看| 欧美午夜高清在线| 色94色欧美一区二区| 成在线人永久免费视频| 又大又爽又粗| 一进一出抽搐动态| 亚洲免费av在线视频| 久久ye,这里只有精品| 99国产精品99久久久久| 老司机深夜福利视频在线观看| 一级黄色大片毛片| 欧美+亚洲+日韩+国产| 一区二区日韩欧美中文字幕| 91成年电影在线观看| 黄频高清免费视频| 免费人妻精品一区二区三区视频| 欧美日韩一级在线毛片| 国产成人免费观看mmmm| 国产aⅴ精品一区二区三区波| 大片免费播放器 马上看| a级毛片黄视频| 国产一区二区三区视频了| 国产亚洲精品第一综合不卡| 国产又爽黄色视频| avwww免费| 国产欧美日韩精品亚洲av| 亚洲欧美一区二区三区黑人| 高清毛片免费观看视频网站 | 多毛熟女@视频| 我的亚洲天堂| 国产精品熟女久久久久浪| 久久婷婷成人综合色麻豆| 一个人免费看片子| av视频免费观看在线观看| 男人操女人黄网站| 一级片免费观看大全| 两人在一起打扑克的视频| 久久99热这里只频精品6学生| av免费在线观看网站| 亚洲伊人色综图| 亚洲男人天堂网一区| 高清黄色对白视频在线免费看| 极品人妻少妇av视频| 黄色a级毛片大全视频| 激情在线观看视频在线高清 | 久久婷婷成人综合色麻豆| 国产欧美日韩一区二区三| 国产黄色免费在线视频| 99热网站在线观看| 免费在线观看日本一区| 我要看黄色一级片免费的| 三上悠亚av全集在线观看| 高清视频免费观看一区二区| 久久人妻熟女aⅴ| 国产成人一区二区三区免费视频网站| 亚洲情色 制服丝袜| 美女国产高潮福利片在线看| 日韩欧美国产一区二区入口| 欧美变态另类bdsm刘玥| 老司机影院毛片| 久久精品亚洲av国产电影网| 高清av免费在线| 在线观看舔阴道视频| 亚洲三区欧美一区| 男女无遮挡免费网站观看| 黄色视频,在线免费观看| 国产av又大| 国产人伦9x9x在线观看| 老司机亚洲免费影院| 国产高清视频在线播放一区| 一本一本久久a久久精品综合妖精| 交换朋友夫妻互换小说| 精品一区二区三区四区五区乱码| 正在播放国产对白刺激| 色婷婷久久久亚洲欧美| 亚洲伊人色综图| 亚洲精品美女久久av网站| 69av精品久久久久久 | 妹子高潮喷水视频| 美女午夜性视频免费| 叶爱在线成人免费视频播放| bbb黄色大片| 99精品在免费线老司机午夜| 麻豆乱淫一区二区| 欧美黑人精品巨大| 俄罗斯特黄特色一大片| 成人手机av| www.熟女人妻精品国产| 成年版毛片免费区| 欧美在线黄色| 日韩欧美免费精品| 国产精品秋霞免费鲁丝片| 成人18禁高潮啪啪吃奶动态图| 国产精品1区2区在线观看. | 国产精品电影一区二区三区 | 国产视频一区二区在线看| 久久亚洲精品不卡| 高清黄色对白视频在线免费看| 女人精品久久久久毛片| 精品国内亚洲2022精品成人 | 中文字幕最新亚洲高清| 久久国产精品影院| 国产一区二区三区视频了| 日本欧美视频一区| 五月开心婷婷网| 69av精品久久久久久 | 日韩欧美一区二区三区在线观看 | 91麻豆精品激情在线观看国产 | 黑人巨大精品欧美一区二区mp4| av在线播放免费不卡| 成在线人永久免费视频| 69精品国产乱码久久久| 精品一区二区三区av网在线观看 | 男女无遮挡免费网站观看| 国产亚洲欧美在线一区二区| 久久亚洲真实| 亚洲欧洲日产国产| a级毛片黄视频| 美女午夜性视频免费| 久久99一区二区三区| 精品欧美一区二区三区在线| 亚洲午夜理论影院| 欧美+亚洲+日韩+国产| 肉色欧美久久久久久久蜜桃| 制服人妻中文乱码| 国产精品亚洲一级av第二区| 精品一品国产午夜福利视频| 精品国产一区二区三区久久久樱花| 精品熟女少妇八av免费久了| 淫妇啪啪啪对白视频| 黑人巨大精品欧美一区二区mp4| 国产精品香港三级国产av潘金莲| 亚洲国产精品一区二区三区在线| 国产精品亚洲av一区麻豆| 黑人操中国人逼视频| 久久久久久人人人人人| 黄片小视频在线播放| 一区二区三区精品91| 两个人看的免费小视频| 另类精品久久| 十八禁网站网址无遮挡| 99精品久久久久人妻精品| 丁香六月欧美| 老司机靠b影院| 美女高潮到喷水免费观看| 精品国产乱码久久久久久小说| 黄色怎么调成土黄色| 黄色 视频免费看| 久久久水蜜桃国产精品网| 久久人人爽av亚洲精品天堂| 侵犯人妻中文字幕一二三四区| 亚洲欧美色中文字幕在线| 国产有黄有色有爽视频| 18禁裸乳无遮挡动漫免费视频| 日日夜夜操网爽| 日韩制服丝袜自拍偷拍| 欧美日韩亚洲综合一区二区三区_| 18禁裸乳无遮挡动漫免费视频| 欧美精品啪啪一区二区三区| 啦啦啦免费观看视频1| 久久精品国产亚洲av高清一级| 自拍欧美九色日韩亚洲蝌蚪91| 建设人人有责人人尽责人人享有的| 老汉色av国产亚洲站长工具| 成人黄色视频免费在线看| 国产成人精品久久二区二区免费| 精品亚洲成a人片在线观看| 肉色欧美久久久久久久蜜桃| 菩萨蛮人人尽说江南好唐韦庄| 黑人操中国人逼视频| 国产成人系列免费观看| 国产亚洲精品久久久久5区| 精品久久蜜臀av无| 又黄又粗又硬又大视频| 国产精品久久久av美女十八| 男女下面插进去视频免费观看| 2018国产大陆天天弄谢| 亚洲色图av天堂| 热99久久久久精品小说推荐| 丝袜人妻中文字幕| 国产日韩一区二区三区精品不卡| 在线观看免费视频日本深夜| 欧美黑人欧美精品刺激| 久久久精品国产亚洲av高清涩受| 国产成人精品久久二区二区免费| 极品少妇高潮喷水抽搐| 满18在线观看网站| 国产高清videossex| 女人被躁到高潮嗷嗷叫费观| 91成年电影在线观看| 成人三级做爰电影| 亚洲欧美日韩另类电影网站| 热99久久久久精品小说推荐| 两性夫妻黄色片| 99国产精品99久久久久| 免费黄频网站在线观看国产| 天堂中文最新版在线下载| 高清在线国产一区| 热re99久久精品国产66热6| 日韩三级视频一区二区三区| 国产成人av教育| 欧美国产精品va在线观看不卡| 热99久久久久精品小说推荐| 亚洲人成电影观看| 久久青草综合色| 97在线人人人人妻| 国产精品久久久久久人妻精品电影 | 久久国产精品人妻蜜桃| 黄色 视频免费看| 国产有黄有色有爽视频| 在线播放国产精品三级| 国产无遮挡羞羞视频在线观看| 国产一区二区三区在线臀色熟女 | 乱人伦中国视频| 亚洲精品中文字幕一二三四区 | 美女午夜性视频免费| 中国美女看黄片| 国产av国产精品国产| 涩涩av久久男人的天堂| 男女边摸边吃奶| 一区二区三区国产精品乱码| 亚洲国产毛片av蜜桃av| 午夜福利在线观看吧| 啦啦啦在线免费观看视频4| 999久久久国产精品视频| 亚洲精品中文字幕一二三四区 | 啦啦啦免费观看视频1| 欧美日韩成人在线一区二区| 国产1区2区3区精品| 十八禁高潮呻吟视频| 99精品欧美一区二区三区四区| 操出白浆在线播放| 色播在线永久视频| 亚洲精品成人av观看孕妇| 亚洲精品粉嫩美女一区| 18禁裸乳无遮挡动漫免费视频| 一级a爱视频在线免费观看| 日韩欧美三级三区| 亚洲免费av在线视频| 侵犯人妻中文字幕一二三四区| videosex国产| 黑人猛操日本美女一级片| 12—13女人毛片做爰片一| 国产精品av久久久久免费| 一级毛片电影观看| 男人舔女人的私密视频| 男女下面插进去视频免费观看| 久久午夜综合久久蜜桃| 精品亚洲乱码少妇综合久久| 久久精品成人免费网站| 一区二区三区国产精品乱码| 国产高清激情床上av| 日韩有码中文字幕| 老司机深夜福利视频在线观看| 久热这里只有精品99| 性少妇av在线| 精品国产亚洲在线| 91av网站免费观看| 国产精品一区二区免费欧美| 精品亚洲成a人片在线观看| 欧美亚洲 丝袜 人妻 在线| 精品卡一卡二卡四卡免费| 18禁国产床啪视频网站| 极品人妻少妇av视频| 一边摸一边做爽爽视频免费| 国产有黄有色有爽视频| 久久久久久久久免费视频了| 在线观看一区二区三区激情| 男人操女人黄网站| 国产精品熟女久久久久浪| 国产欧美日韩精品亚洲av| 久久久久国产一级毛片高清牌| 亚洲精品国产区一区二| 99国产精品一区二区蜜桃av | 久久人妻av系列| 久久精品国产a三级三级三级| 欧美精品一区二区大全| 91精品国产国语对白视频| 国产精品香港三级国产av潘金莲| 欧美日韩精品网址| 天堂俺去俺来也www色官网| 天堂中文最新版在线下载| 黄片播放在线免费| 亚洲国产毛片av蜜桃av| 亚洲色图 男人天堂 中文字幕| 99国产精品一区二区三区| 人人妻人人澡人人爽人人夜夜| 国产aⅴ精品一区二区三区波| 亚洲av成人不卡在线观看播放网| 又大又爽又粗| 成人永久免费在线观看视频 | 免费久久久久久久精品成人欧美视频| 丝袜在线中文字幕| cao死你这个sao货| 久久国产亚洲av麻豆专区| 考比视频在线观看| 成人国产av品久久久| 欧美黑人欧美精品刺激| 777米奇影视久久| 久久久久久久大尺度免费视频| 母亲3免费完整高清在线观看| 久久 成人 亚洲| 热99久久久久精品小说推荐| 日韩大片免费观看网站| 12—13女人毛片做爰片一| 麻豆av在线久日| 美国免费a级毛片| 一本一本久久a久久精品综合妖精| 大片电影免费在线观看免费| 国产精品久久久av美女十八| 成人国产av品久久久| 国产成人精品久久二区二区91| 亚洲av片天天在线观看| 国产伦理片在线播放av一区| 国产激情久久老熟女| 精品人妻1区二区| 久久午夜亚洲精品久久| 曰老女人黄片| 叶爱在线成人免费视频播放| 日日夜夜操网爽| a级片在线免费高清观看视频| 狂野欧美激情性xxxx| 国产欧美日韩一区二区三区在线| 一级,二级,三级黄色视频| 在线十欧美十亚洲十日本专区| 一个人免费在线观看的高清视频| 国产片内射在线| 久久久国产欧美日韩av| 黄色片一级片一级黄色片| 国产亚洲一区二区精品| 电影成人av| 亚洲自偷自拍图片 自拍| 国产欧美亚洲国产| 窝窝影院91人妻| 大香蕉久久网| 99九九在线精品视频| 午夜福利在线免费观看网站| 亚洲精品国产色婷婷电影| 国产免费av片在线观看野外av| 午夜成年电影在线免费观看| 婷婷成人精品国产| 国产黄色免费在线视频| 欧美激情久久久久久爽电影 | 免费高清在线观看日韩| 啦啦啦在线免费观看视频4| 欧美黄色片欧美黄色片| 亚洲欧洲日产国产| 波多野结衣av一区二区av| 国产麻豆69| 大型黄色视频在线免费观看| 中文欧美无线码| 一边摸一边做爽爽视频免费| 18禁美女被吸乳视频| 一本久久精品| 亚洲人成77777在线视频| 日日夜夜操网爽| 美女高潮喷水抽搐中文字幕| 高清在线国产一区| avwww免费| 婷婷丁香在线五月| 精品亚洲成国产av| 国产欧美日韩综合在线一区二区| 妹子高潮喷水视频| 亚洲人成电影观看| 国产精品免费视频内射| 在线看a的网站| 无限看片的www在线观看| 黄色视频,在线免费观看| 中文字幕高清在线视频| 99在线人妻在线中文字幕 | 国产一区二区 视频在线| 亚洲国产看品久久| 少妇的丰满在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 国产区一区二久久| 国产野战对白在线观看| 无人区码免费观看不卡 | 大香蕉久久成人网| 亚洲中文字幕日韩| 18在线观看网站| 黑人巨大精品欧美一区二区mp4| 一级片免费观看大全| 成年版毛片免费区| 最近最新中文字幕大全免费视频| 久久精品亚洲熟妇少妇任你| 久久人妻熟女aⅴ| 久久精品亚洲精品国产色婷小说| 99精国产麻豆久久婷婷| 香蕉久久夜色| 中文字幕人妻丝袜制服| 美国免费a级毛片| 极品教师在线免费播放| 大片免费播放器 马上看| 成人免费观看视频高清| 日日爽夜夜爽网站| 久久久水蜜桃国产精品网| 啦啦啦中文免费视频观看日本| 99国产精品免费福利视频| 国产成人av激情在线播放| 19禁男女啪啪无遮挡网站| 91精品国产国语对白视频| videosex国产| 天天躁日日躁夜夜躁夜夜| 国产av精品麻豆| 夜夜骑夜夜射夜夜干| 国产免费福利视频在线观看| 免费在线观看视频国产中文字幕亚洲| 天堂中文最新版在线下载| 久久久久网色| 久久久久久久精品吃奶| 老司机在亚洲福利影院| 18在线观看网站| 激情视频va一区二区三区| 精品熟女少妇八av免费久了| 久久久久久人人人人人| 日韩有码中文字幕| 亚洲av成人不卡在线观看播放网| 极品人妻少妇av视频| 久久久国产欧美日韩av| 一级片'在线观看视频| 国产淫语在线视频| 色尼玛亚洲综合影院| 久久中文看片网| 啦啦啦中文免费视频观看日本| 一个人免费在线观看的高清视频| 搡老熟女国产l中国老女人| 肉色欧美久久久久久久蜜桃| 亚洲va日本ⅴa欧美va伊人久久| 男女免费视频国产| 黄色成人免费大全| 18在线观看网站| 国产精品久久久久久精品古装| 亚洲色图 男人天堂 中文字幕| 每晚都被弄得嗷嗷叫到高潮| 亚洲伊人久久精品综合| 国产日韩欧美亚洲二区| 久久精品国产99精品国产亚洲性色 | 菩萨蛮人人尽说江南好唐韦庄| 国产野战对白在线观看| 欧美成人免费av一区二区三区 | 下体分泌物呈黄色| 中国美女看黄片| 欧美亚洲日本最大视频资源| 女警被强在线播放| 国产日韩欧美亚洲二区| 免费久久久久久久精品成人欧美视频| 国产精品国产av在线观看| 在线播放国产精品三级|