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

    Analysis on the Network Node of Konjac Glucomannan Molecular Chain①

    2014-12-17 03:04:08WENChengRongSUNZhongQiPANGJieMAZhenSHENBenShuXIEBingQingJINGPu
    結(jié)構(gòu)化學(xué) 2014年8期

    WENCheng-RongSUNZhong-QiPANGJieMAZhenSHENBen-ShuXIEBing-QingJINGPu

    ?

    Analysis on the Network Node of Konjac Glucomannan Molecular Chain①

    WEN Cheng-RongaSUN Zhong-QiaPANG JieaMA ZhenaSHEN Ben-Shua②XIE Bing-QingaJING Pub

    a(350002)b(200240)

    The dynamic changes of the complex network and the material form and function were actuated by the molecular chains. The interaction behavior between molecular chains was difficult to illuminate because the dynamic changes of macromolecules were observed difficultly by normal spectrum method and the methods to test and evaluate the complex network evolution prediction and intervention are rare. The mathematic model of domino offect of molecular chains was established based on the topological structure of molecular chain aggregation of Konjac glucomannan, and the molecular entanglement mechanism of Konjac glucamannan blends was studied through molecular simulation and knot theory analysis combined with experimental verification. The results suggested that two network models (topological entanglement and solid knot) of Konjac glucomannon blends were formed through hydrogen bond nodes. The topological entanglement was strengthened with the increase of concentration and the form of molecular chains was Gaussian chain which could not allow traverse moving owing to the intermolecular cross and entanglement and the shield of intramolecular interaction. Besides, the structures of Konjac glucomannon blends became more stable due to the solid knot. Both of them were verified by the experimental results. This experimental method simplifies the microscopic description of Konjac glucomannon, and there is important guiding significance of the experimental results for the prediction and control of polysaccharides’ structure and function.

    knot, Konjac glucomannon, topological entanglement, network node

    1 INTRODUCTION

    As a typical mathematical method, topology has important applications on the aspect of chain ring formed by polymer, such as nucleic acids, proteins, polysaccharides,[1-3]. The knot theory, one of topologies, is a new research field developed rapidly in recent years. It helps explain how interaction causes change of polymer chain spatial structure and form complex network, then it helps forecast and regulate the polymer’s function[4-5]. Molecular chain is a new research subject in knot theory. The forming of polymer molecular chain model can help form prediction mechanism of the path of chain, which is expected to provide a simple, unified platform andfairer comparison for the evolution of complex net-works. What’s more, it will greatly promote the theory of evolution of complex networks model. The polymer’s physicochemical properties can be pre- dicted by its structure characteristic parameters[6].

    How to drive the evolution of a complex network of macromolecules at the molecular level, to explore the interactions between molecules, to find be- havioral effects of macromolecular molecular chains and to provide prediction and regulation of poly- mer’s structure and function have become a hot and difficult issue of studies reference. Now the researches on polymer’s chain are almost focused on the experimental results and test methods, and those on molecular simulation almost focus on the si- mulation system, single mechanism and simulation methods. Owing to the complexity of biological macromolecules, the existing methods are not good to clarify the intermolecular interaction behaviors, or to achieve the target of regulating its structure and function in dose-effect.

    Polysaccharide is a kind of macromolecular sub- stances with complex structure and functions[7,8], whose physical and chemical properties and bi- ological activity are affected by their network structure[9,10]. How to characterize and design the molecular structure of polysaccharide and regulate its function has become difficult and hot issue. In some cases, the polysaccharide will show a ring structure of end to end, which is not a simple planar ring structure, but topological structure in different environments[11]. Some of its properties and be- haviors are independent on their chemical com- position and other details, but only on its own topological structure.These properties and behavior of polysaccharide are between micro and macro, which need to be studied from mesoscale. So, it can be studied by approximated research of molecular simulation, without considering the details[12,13]. Although the basis of research is experiment, we need to explore the nature of the phenomenon from theory and calculations, and then calculate the unknown quantity.

    Konjac glucomannan (KGM) is one of the natural polysaccharides in which glucose and mannose at the molar ratio of 1:1.6~1.7 are connected via-1,4 glycosidic bonds[14]. KGM has been widely applied in the fields of food, biomedical, materials,. as water conservation materials, cartilage, drug and nutrient carriers, and so on[6,14]. Owing to the limi- tations of experimental studies, the conformation of KGM molecular chain and its stable formation mechanism, especially the interaction of other poly- mer and inorganic small molecule compounds, have not been studied in depth[15]. Based on the topology of KGM, combined with molecular simulation and mathematical model establishment, the molecular links of KGM, KGM compounds, and KGM blends are analyzed. The processes of changing from trivial knot to the topology entanglement and topological knot and building molecular chain network mathe- matical model to predict the formation of a link mechanism are studied. And the mathematical model of molecular chain network and prediction me- chanism of chain path are established as well. In order to achieve predicting results, the verification of experimental research is adopted, too. This study fuses so many multidisciplinary advantages that it provides a new thought of the functional regulation of poly- saccharides and other macromolecules.

    2 MATERIALS AND METHODS

    2.1 Materials

    Konjac refined powder was obtained from San Ai Organic Konjac Development Co., Ltd., Zhaotong, Yunnan, China; Curdlan (CUD) was purchased from Japanese Takeda kirin food co., LTD., Japan; Other chemical reagents (analytical grade) were obtained from Sinopharm Chemical Reagent Co. Ltd.

    2.2 Main instruments and equipments

    EquipmentTypeOrigin Circular dichroism, CDBioLogic MOS 450Thermo Electron Corporation, USA Magnetic stirring85-2Changzhou sino instrument co., LTD Digital thermostat water bathHH-2Jiangyin Poly Scientific Instruments Co. Ltd. Electronic analytical balancePL402-CMETTLER TOLEDO Instruments Co., Ltd.

    2.3 Experimental methods

    2.3.1 Preparation of KGM/borax blend

    1.0 g KGM was dispersed in 100 mL distilled deionized water to obtain 0.01 g/mL solutions, and stirred by a magnetic stirrer for 2 h at 45 ℃, then vacuum filtrated by a 200 mesh filter cloth. 5 mL 0.6% borax solution was added to 95 mL KGM solution which was stored at 4 ℃ for 1 d, then stirred for 20 min and kept at room temperature for 40 min[16].

    2.3.2 Circular dichroism (CD) spectra

    The same volume of 80mol/L congo dye solution and KGM/ borax blend solution were mixed. After being stored for 24 h, bioLogic MOS 450 circular dichroism spectrometer was used. Deter- mination conditions: temperature: 25 ℃, radius of color dish: 0.2 cm, wavelength: 180~400 nm, distinguishability: 0.5 nm, scanning speed: 30 nm/ min.

    2.3.3 Molecular model optimization

    The molecular dynamic simulation of the solution of KGM and its blends in water (4 ℃) in addition to individual aqueous solutions was carried out on O2working station (American HP) using Hyper- Chem7.0 release soft ware (American Hypercube Corporation). Force-field function included bond, angle and non-bonding; the Langevin temperature regulation method was used to adjust the tempera- ture; the ensemble was canonical ensemble; the step size of calculation was 0.001 ps and the relaxation time 0.1 ps; the system was run for 100 ns; and the interval of data collection was 2 ps. In order to understand the law of molecular conformation change in optimizing procession intuitively, the H atoms were concealed.

    2.3.4 Molecular chain model building of KGM

    The kink theory was analyzed by simplifying the complex structure of KGM to flat topological structure combined with molecular chain folding, based on the topology analysis of P. C. Burns[17]. Then the mathematic model was established and the knot was found. The results were verified by experi- mental research at last.

    We assumed that the zero chain knot coincided with the negative half shaft of theaxis in the three-dimensional rectangular coordinate system, whose terminal point is the origin of coordinate. So, the terminal point of the ithchain knot was marked as Mi= (xi, yi, zi), i = 0, 1, 2, …, n, which means M0= (0, 0, 0). And we also assumed that QiMi= { cosφi, cosψi, cosθi}, because the following chain knot is relative to the prior’s free rotatability, of which free rotation angle Φ is a uniformly distributed random number in the range of [0, 2π)[18].

    (1) When i = 1, the uniformly distributed random number in the range of [0, 2π) is Φ1, andcosφ1is confirmed.

    Then we get a formula as follows: cosψ1= ±√ (1-cos2φ1- cos2θ1)= ±√ (1-cosφ1),(θ1=π/2). It is “+” as 0 ≤ φ1≤ π, while it is “-” as π ≤ φ1≤ 2π.

    And we get the following formulas, x1= l*sin70°32*cosφ1, y1= l*sin70°32*cosψ1, z1= l*sin70°32*cosθ1, which means the point of M1is confirmed.

    Popularly, if the point of Mi= (xi, yi, zi) is confirmed, we can count the results of cosines in the direction of Mi-1Mias below: cosα1= (xi-xi-1)/l, cosβ1= (yi-yi-1)/l, cosγ1= (zi-zi-1)/l.

    Then the coordinate values of planar circle (Qi+1) are as follows:

    ui+1= xi+1*cosα1, vi+1= yi+1*cosβ1, wi+1= zi+1*cosγ1.

    And the plane equation is perpendicular to the direction of Mi-1Miand crosses the point of Qi+1as below: cosαi*(x-ui+1)+cosβi*(y-vi+1)+cosγi*(z-wi+1)= 0.

    For instance, if i = 1 and φ1= π/4, we get the following results: cosφi=√2/2, cosψi=√2/2, cosθi=0;

    x1= 0.6717* l, y1= 0.6717* l, z1= 0;

    cosα1= 0.6717, cosβ1= 0.6717, cosγ1= 0;

    u2= 0.67172* l, v2= 0.67172* l, w2= 0;

    0.6717*(x- 0.67172* l) + 0.6717*(y- 0.67172* l) = 0, which can be simplified to x+y = 0.67172* l.

    (2) If we get the uniformly distributed random number (Φi+1) in the range of [0, 2π), we can work out

    cosθi+1= (-cosαi* cosφi+1* cosγi±| cosβi|*√ (|cos2βi+ cos2γi- cos2φi+1|))/( cos2βi+ cos2γi) (The probability of "+" and "-" is the same) and cosψi+1= ±√ (1-cos2φi+1- cos2θi+1), (It is “+” as 0 ≤ φi+1≤ π, but “–” as π ≤ φi+1≤ 2π) from the simultaneous equations as follows: cosαi* cosφi+1+ cosβi* cosψi+1+ cosγi* cosθi+1= 0 and cos2φi+1+ cos2ψi+1+cos2θi+1= 0.

    Then the point values of Mi+1= (xi+1, yi+1, zi+1) are worked out as below:

    xi+1= ui+1+r* cosφi+1, yi+1= vi+1+r* cosψi+1, zi+1= wi+1+r* cosθi+1.

    3 RESULTS AND ANALYSIS

    3.1 Establishment of the KGM molecularchain model

    The force of intermolecular chains is hydrogen bonds of KGM-KGM and KGM-water when KGM is dissolved in water. The stable conformation and topological structure of KGM in water are shown in Figs. 1 and 2, respectively.

    The knot is a closed curve which never crosses itself in three-dimensional space. When a knot is equivalent with the planar circumference, we call it a trivial knot which plays a role of zero[19]. KGM is composed of-(1→4) linked D-mannose and D-glucose in a molar ratio of 1.6:1 or 1.4:1, with about 1 in 19 units being acetylated. KGM may contain short side branches at the C-3 position of mannoses, and acetyl groups are randomly present at the C-6 position of sugar units[20]. The C and O atoms are represented by salient points and end- points, respectively. The bonds between C and C or O are represented by straight line, while those between H and O or C are ignored. The molecular chain of KGM in water solution is approximately circular from overlooking (shown as Fig. 1), and its helix chain is a trivial knot (shown as Fig. 2).

    3.2 Establishment of the KGM/CUD molecular model

    3.2.1 Entanglement of KGM/CUD blends

    The mesoscopic molecular chain of KGM/CUD blends in water solution is shown in Fig. 3, from which we can observe the mesoscopic morphology of KGM/CUD blend. Fig. 3 tells us that the KGM chains are entangled and intertwined with the CUD chains tightly, existing in a form of random coil state. The overlapping effect of the molecular chain has formed a strong interaction entanglement between the molecular chains that are conducive to enhancing the interaction of intermolecular chain and the pro- perties of KGM/CUD blend.

    Fig. 1. Stable conformation of KGM in water solution

    Fig. 2. Knot structure of KGM in water solution

    Fig. 3. Mesoscopic molecular chain form of KGM-CUD blends (Red is CUD; pink, blue and green are KGM)

    3.2.2 Topology entanglement model of the KGM/CUD blends

    Entanglement is one of the most important charac- ters of polymer’s condensed state. It is a network structure constituted with physical junction points between the polymer chains, and the movement of the molecular chain is fettered and restricted by the surrounding molecules. According to Flory’s deduc- tion[21], the intermolecular chain’s entanglement is more stable due to its orientation when the entangled molecular chains form a network structure. But the orientation will take time, and the hindering of entanglement between the molecular chains makes it longer. If the entanglement of molecular chains has a synergistic effect, it is beneficial to the stability of structure; otherwise it is harmful.

    When KGM and CUD are blended in the external force, the entanglement processes of molecular chains are hindered with each other, and the time of forming a network structure is prolonged. Otherwise, the entanglement of molecular chains is affected by the concentration. When in low concentration, the entanglement tends to be folded by itself; when the concentration is increased, the polymer chains are throughout entangled with each other to form an entanglement network, as shown in Fig. 3; when in extreme high concentration, the polymer chains penetrate fully with each other, which are Gaussian chains, so that intramolecular interaction is replaced by intermolecular interaction and completely shiel- ded. The entanglement and penetration of molecular chains are respectively shown as the box and circle parts in Fig. 4, because of which the traverse move- ment of molecular chains is restricted. Therefore, the molecular structure is stable.

    3.3 Establishment of the KGM/borax molecular model

    3.3.1 Molecular simulation of the KGM/borax molecule

    The conformation and hydrogen bonding sites of KGM/borax blends are shown in Fig. 5. The mole- cular dynamics simulations of KGM/boric blends were taken out under periodic boundary conditions in different temperature. It shows that the molecular chains were helical structures under the effect of hydrogen bonds in low and normal temperature. But, as the temperature increased, the hydrogen bonding sites between the KGM molecule and boric acid ions changed slightly. Although the hydrogen bonds of boric acid ions with acetyl were broken, those of the boric acid ions with -OH of Man and Glu glycone remained stable. It suggests that the helical struc- tures of KGM/borax blends were stable with less influence of temperature.

    3.3.2 Solid knot of KGM/borax blends

    There is some difference between the mechanism of molecular chain entanglement between KGM/ borax blends with KGM/other polymer blends. The KGM molecular chain forms helix rings through self winding, knots through weak bonds and supra- molecular force such as hydrogen bond, and solid knots through coordination reaction between boric acid ions and glycosyl, in KGM/borax blends. The molecular conformation is stable because of the acetyl on the helix structure, which owns interspace in the structure. The mathematic modelof KGM/ borax blends is a torus knot ring, as shown in Fig. 6.

    3.3.3 Circular dichroism spectrum of KGM/borax blends

    The circular dichroism spectrum is an important method to research the conformation of polymer molecule including the discipline of conformation change. As shown in Fig. 7, the variation tendencies of the circular dichroism spectrum of KGM/borax blends are almost the same at 30 and 60 ℃ from 180 to 400 nm wavelength. Both of them have positive cotton domino offect near 190, 210, 250, and 320 nm, and negative cotton domino offect around 200, 220, 300, and 345 nm. The difference lies in that the positive cotton domino offect of KGM/borax blends at 60℃ is weaker than that at 30℃, while it is on the contrary for the negative, especially in the 200 nm wavelength. Overall, the effect of temperature on the conformation of KGM/borax blends is slight, suggesting that the helical structures of KGM/borax blends were stable and the analysis of solid knot was feasible.

    Fig. 4. Topological entanglement of the KGM/CUD blends (Red is CUD, black is KGM; circle part is penetration, and box part is entanglement)

    Fig. 5. Conformation and hydrogen bonding sites of the KGM/borax blends(dotted line is hydrogen bond)

    Fig. 6. Torus knot ring of the KGM/borax blends

    Fig. 7. Effect of temperature on the circular dichroismspectrum of KGM/borax blends

    4 DISCUSSION AND CONCLUTION

    Nowadays, there is a further deep research on the prediction of complex networks through topological entanglement and path of chain, as a promotion of molecular dynamics method[22,23]. There are a lot of advantages on using computer simulation to research the network of polymers as well as using the concept of system or entirety to study some complex phenomena, which make the prediction results more accurate because more factors could be considered. It is absolutely a good method to research the complex network structure combined with computer simulation with knot analysis, which has a certain guiding role in predicting and controlling the structure and function of polysaccharides. This means not only makes up for some deficiency of the experimental results which are difficult to observe with the existing experimental methods, but also makes up for the difficulties of complex systems which can not be solved analytically in theory. However, it helps us observe polymer structures and functions more in-depth[24-26].

    The complex networks of KGM and CUD were divided into solid molecule through conbining topological entanglement and path of chain together. It can disassemble polysaccharide movement to a lot of single movements and then superposed, describe the structure of natural polysaccharide as a dynamic and elastic network structure, and simplify a polysaccharide to a line and a micromolecule to a dot. Harmonic potential with a same elastic coefficient was used in the interaction between the molecular chains of polysaccharide. We have found that KGM’s helix chain is a trivial knot and proposed and validated the mathematical model of KGM blends molecule chain network (topology entanglement and topological knot model). This is benefit on understanding the complex network of molecular chains in theory, and this also provides network evolution mechanism for the same class or network. The conformation and molecular chain of KGM blends are representated and designed very well and reasonable models of this chain are esta- blished to guide experiment and control its function.

    Knot theory and mathematic model simplify the microscopic description of KGM molecular chain, and the predicted results have guiding effect on the structure and function prediction and regulation of polysaccharides. It can also achieve rational reconstruction in this process and the optimization of molecular chain network. Finally, it makes the achievement of full dose regulation possible. It provides a theoretical reference of regulating and controlling the fold, entanglement and properties of polysaccharide molecular chain by using the net- work evolution to simulate the molecular chain of optimal physical and chemical state, but it also provides reference effect of methods on the research of interaction between polysaccharides and rational reconstruction and inorganic micromolecule. The topological entanglement is proved objectively to be existence by experimental facts. But, there should be a new breakthrough in the observations and mole- cular dynamics simulations of microscopic images of polymer topological entanglements. And there should be more databases to solve the problem of complex network control through the prediction theory of path of chain so as to predict the in- formation of structure, which will be the focus of future research directions.

    (1) Duan, J. W.; Zheng, Z.; Zhou, P. P.; Qiu, W. Y. The architecture of Sierpinski knots.. 2012, 68, 595-610.

    (2) Redner, S. Teasing out the missing links.2008, 453, 47.

    (3) Li, H.; Helling, R.; Tang, C.; Wingreen, N. S. Emergence of preferred structures in a simple model of protein folding.1996, 273, 666-669.

    (4) Darcy, I. K. Biological distances on DNA knots and links, applications to XER recombination.2001, 10, 269-294.

    (5) Vazquez, M.; Colloms, S. D.; Sumners, D. W. Tangle analysis of Xer recombination reveals only three solutions, all consistent with a single three-dimensional topological pathway.2005, 346, 493-504.

    (6) Haliloglu, T.; Bahar, I.; Erman, B. Gaussian dynamics of folded proteins.1997, 79, 3090-3093.

    (7) Dumitriu, S. Eds.. CRC Press 2010.

    (8) Stephen, A. M.; Phillips, G. O.; Williams, P. A. Eds.CRC Press 2010.

    (9) Zhang, B. B.; Su, H. Q.; Gu, X. J.; Huang, X. L.; Wang, H. M. Effect of structure and charge of polysaccharide flocculants on their flocculation performance for bentonite suspensions.2013, 436, 443-449.

    (10) Tzianabos, A. O. Polysaccharide immunomodulators as therapeutic agents, structural aspects and biologic 33.

    (11) Orellana, A. Biosynthesis of non-cellulosic polysaccharides in the Golgi apparatus. Topological considerations.2005, 139, 42-45.

    (12) Lyman, E.; Cui, H. S.; Voth, G. A. Reconstructing protein remodeled membranes in molecular detail from mesoscopic models.2011, 13, 10430-10436.

    (13) Draviam, V. M.; Funahashi, A.; Hiroi, N.; Kimura, A.; Kobayashi, T. J.Meeting report of the international workshop on quantitative biology 2012, mesoscopic and microscopic worlds meet.2012, 3.

    (14) Jian, W. J.; Zeng, Y.; Xiong, H. J.; Pang, J. Molecular simulation of the complex of konjac glucomannan–borate in water.2011, 85, 452-456.

    (15) Melinda, C.; Timothy, C.; Baldwin, T. J.; Hocking, K. C. Traditional uses and potential health benefits of Amorphophallus konjac K. Koch ex N.E.Br.2010, 128, 268-278.

    (16) Forget, A.; Christensen, J.; Lüdeke, S.; Kohler, E.; Tobias, S.; Matloubi, M.; Shastri, V. P. Polysaccharide hydrogels with tunable stiffness and provasculogenic properties via-helix to-sheet switch in secondary structure.2013, 110, 12887-12892.

    (17) Jian, W. J.Fujian Agriculture and Forestry University2011.

    (18) Burns, P. C.; Deely, K. M.; Skanthakumar, S. Neptunium incorporation into uranyl compounds that form as alteration products of spent nuclear fuel, implications for geologic repository performance.2004, 92, 151-160.

    (20) Kinoshita, S.; Terasaka, H. On unions of knots.1957, 9, 131-153.

    (21) Katsuraya, K.; Okuyama, K.; Hatanaka, K.; Oshima, R.; Sato, T.; Matsuzaki, K. Constitution of konjac glucomannan, chemical analysis and13C NMR spectroscopy. Carbohydrate Polymers2003, 53, 183-189.

    (22) Flory, P. J.; Rehner, J. Statistical mechanics of cross-linked polymer networks II. welling.2004, 11, 521-526.

    (23) Walker, D. M.; Tordesillas, A. Topological evolution in dense granular materials, a complex networks perspective.2010, 47, 624-639.

    (24) Tuszynska, I.; Bujnicki, J. M. Predicting atomic details of the unfolding pathway for YibK, a knotted protein from the SPOUT superfamily.2010, 27, 511-520.

    (25) Lee, J. Y.; Chun, M. S.; Jung, H. W.; Hyun, J. C. Conformational dynamics of sub-micron sized wormlike polyelectrolyte polymer in flow fields.2012, 20, 1163-1172.

    (26) Perkyns, J. S.; Lynch, G. C.; Howard, J. J.; Pettitt, B. M. Protein solvation from theory and simulation, exact treatment of Coulomb interactions in three-dimensional theories.2010, 132, 064106.

    21 April 2014;

    16 July 2014

    ① This project was supported by the National Natural Science Foundation of China (31271837), Specialized Research Fund for the Doctoral Program of Higher Education jointly funded by Ministry of Education (20113515110010), Special Research Funds from Ministry of Science and Technology (2012GA7200022), Major projects of industries, universities and research in Fujian Province (2013N5003), and Natural Science Foundation of Fujian Province (2011J0101)

    . Pang Jie, professor, majoring in food chemistry and nutrition. E-mail: pang3721941@163.com

    亚洲av成人精品一区久久| 日本免费一区二区三区高清不卡| 日韩精品青青久久久久久| 此物有八面人人有两片| 91麻豆精品激情在线观看国产| 搡老岳熟女国产| 国产一区二区三区在线臀色熟女| 国产精品久久久久久人妻精品电影| 中国美女看黄片| 狂野欧美白嫩少妇大欣赏| 麻豆成人午夜福利视频| 国产精品电影一区二区三区| 91在线观看av| 动漫黄色视频在线观看| 国产一区在线观看成人免费| 免费看美女性在线毛片视频| 久久午夜亚洲精品久久| 欧美极品一区二区三区四区| 亚洲人成网站高清观看| 久久精品91蜜桃| 国产精品影院久久| 久久精品国产清高在天天线| 最近最新免费中文字幕在线| 久久久久久亚洲精品国产蜜桃av| 黄色a级毛片大全视频| 国产激情欧美一区二区| 亚洲精品中文字幕在线视频| 午夜a级毛片| 长腿黑丝高跟| 大型av网站在线播放| 午夜福利视频1000在线观看| 两人在一起打扑克的视频| 亚洲精品国产一区二区精华液| 亚洲av成人精品一区久久| 国产激情偷乱视频一区二区| 神马国产精品三级电影在线观看 | 丝袜美腿诱惑在线| 欧美成人性av电影在线观看| 老司机午夜福利在线观看视频| 国产又黄又爽又无遮挡在线| 亚洲最大成人中文| 伊人久久大香线蕉亚洲五| 天堂√8在线中文| 在线观看美女被高潮喷水网站 | 国产免费av片在线观看野外av| 亚洲欧美一区二区三区黑人| 免费在线观看亚洲国产| 中文字幕久久专区| 国产精品1区2区在线观看.| 亚洲成人精品中文字幕电影| 久久久久久久午夜电影| 桃红色精品国产亚洲av| 成熟少妇高潮喷水视频| 天堂av国产一区二区熟女人妻 | 国产私拍福利视频在线观看| 亚洲av第一区精品v没综合| 欧美日韩一级在线毛片| 国产久久久一区二区三区| 免费一级毛片在线播放高清视频| 午夜两性在线视频| 久久国产乱子伦精品免费另类| 亚洲精品粉嫩美女一区| 黄片大片在线免费观看| 亚洲精品国产一区二区精华液| 国产1区2区3区精品| 日日爽夜夜爽网站| 好男人电影高清在线观看| 久久亚洲精品不卡| 国产亚洲精品久久久久5区| 狠狠狠狠99中文字幕| 国产av不卡久久| 欧美黑人巨大hd| 狠狠狠狠99中文字幕| 别揉我奶头~嗯~啊~动态视频| 色av中文字幕| 中文资源天堂在线| 三级国产精品欧美在线观看 | 宅男免费午夜| 男女之事视频高清在线观看| 亚洲欧美激情综合另类| 黄色丝袜av网址大全| 国产成+人综合+亚洲专区| 两个人的视频大全免费| av中文乱码字幕在线| 嫩草影院精品99| 高清毛片免费观看视频网站| 少妇的丰满在线观看| 日韩大尺度精品在线看网址| 国产三级中文精品| 亚洲熟妇熟女久久| 久久九九热精品免费| 亚洲一区高清亚洲精品| 亚洲国产精品合色在线| 2021天堂中文幕一二区在线观| 国产精品久久久久久精品电影| 色综合站精品国产| 亚洲七黄色美女视频| 狂野欧美白嫩少妇大欣赏| 最近最新中文字幕大全电影3| 一本大道久久a久久精品| 国产黄片美女视频| 制服诱惑二区| 日韩欧美在线二视频| 变态另类成人亚洲欧美熟女| 国产又色又爽无遮挡免费看| 国产黄a三级三级三级人| 欧美日韩福利视频一区二区| 亚洲人成网站高清观看| 久久久久久免费高清国产稀缺| 久久人人精品亚洲av| 午夜亚洲福利在线播放| 久久久久久亚洲精品国产蜜桃av| 亚洲国产欧美网| 婷婷精品国产亚洲av| 999精品在线视频| av有码第一页| 国产熟女午夜一区二区三区| 黑人欧美特级aaaaaa片| 亚洲va日本ⅴa欧美va伊人久久| 97碰自拍视频| 国产精品九九99| 久久久久国内视频| 少妇的丰满在线观看| 亚洲成av人片在线播放无| 两个人的视频大全免费| 99re在线观看精品视频| 国产av又大| 级片在线观看| 亚洲色图av天堂| 国产蜜桃级精品一区二区三区| 亚洲九九香蕉| 亚洲人与动物交配视频| 欧美一级a爱片免费观看看 | 国产成人aa在线观看| 成人三级做爰电影| 成人av一区二区三区在线看| 国产欧美日韩精品亚洲av| 午夜福利在线观看吧| 午夜成年电影在线免费观看| 在线观看66精品国产| 无遮挡黄片免费观看| 丰满人妻一区二区三区视频av | 熟妇人妻久久中文字幕3abv| 国产爱豆传媒在线观看 | 97超级碰碰碰精品色视频在线观看| 特大巨黑吊av在线直播| 国产一区二区三区视频了| 1024香蕉在线观看| 日本 欧美在线| 无人区码免费观看不卡| 精品福利观看| 国产精品av视频在线免费观看| 狠狠狠狠99中文字幕| av片东京热男人的天堂| 成人永久免费在线观看视频| 国产亚洲精品第一综合不卡| 久久久精品欧美日韩精品| 可以免费在线观看a视频的电影网站| 在线国产一区二区在线| av天堂在线播放| 91麻豆精品激情在线观看国产| 久久性视频一级片| 黄色a级毛片大全视频| 国产亚洲精品久久久久久毛片| 国产三级在线视频| 一级片免费观看大全| 久久草成人影院| 91在线观看av| 午夜福利18| 亚洲国产欧美一区二区综合| 国产高清videossex| 欧美不卡视频在线免费观看 | 日本一二三区视频观看| 国产精品久久电影中文字幕| 国产午夜福利久久久久久| 亚洲男人天堂网一区| 18美女黄网站色大片免费观看| 亚洲人成伊人成综合网2020| 亚洲午夜精品一区,二区,三区| 亚洲国产看品久久| 麻豆av在线久日| 熟女少妇亚洲综合色aaa.| 欧美日韩黄片免| 国产精品免费一区二区三区在线| 午夜精品在线福利| 国产又黄又爽又无遮挡在线| 美女免费视频网站| 老鸭窝网址在线观看| 婷婷丁香在线五月| 国内少妇人妻偷人精品xxx网站 | 久久中文字幕人妻熟女| 欧美黄色淫秽网站| 91九色精品人成在线观看| 男女下面进入的视频免费午夜| 国产成人精品久久二区二区91| 少妇的丰满在线观看| 免费一级毛片在线播放高清视频| 亚洲国产精品999在线| 久久人妻福利社区极品人妻图片| 久久久久久亚洲精品国产蜜桃av| 亚洲午夜精品一区,二区,三区| 国产主播在线观看一区二区| 啪啪无遮挡十八禁网站| 中文字幕av在线有码专区| 日本免费一区二区三区高清不卡| 欧美性猛交黑人性爽| 男人的好看免费观看在线视频 | 国产成人欧美在线观看| 9191精品国产免费久久| 中文字幕人成人乱码亚洲影| 色播亚洲综合网| 丰满人妻熟妇乱又伦精品不卡| 精品福利观看| 国产精品一区二区三区四区久久| 国产成人精品无人区| 欧美另类亚洲清纯唯美| 国产黄片美女视频| 国产一区在线观看成人免费| 麻豆成人午夜福利视频| 亚洲一区二区三区不卡视频| 男女下面进入的视频免费午夜| 国产av又大| 可以在线观看毛片的网站| 男人舔女人的私密视频| 国产主播在线观看一区二区| 91国产中文字幕| 亚洲av成人精品一区久久| 在线观看舔阴道视频| 一区福利在线观看| av片东京热男人的天堂| 精品国产亚洲在线| 国产黄a三级三级三级人| 国模一区二区三区四区视频 | 亚洲av电影不卡..在线观看| 曰老女人黄片| 欧美在线黄色| 日本免费a在线| 久久草成人影院| 亚洲欧美精品综合一区二区三区| 长腿黑丝高跟| 观看免费一级毛片| 黄色成人免费大全| 亚洲精品美女久久av网站| 成人av在线播放网站| 国产精品野战在线观看| 中文字幕人成人乱码亚洲影| 国产精华一区二区三区| 人人妻人人看人人澡| 两性午夜刺激爽爽歪歪视频在线观看 | 黑人巨大精品欧美一区二区mp4| 国内少妇人妻偷人精品xxx网站 | 亚洲精品色激情综合| tocl精华| 一级毛片精品| 亚洲精品一区av在线观看| 俺也久久电影网| 国产探花在线观看一区二区| www.自偷自拍.com| 亚洲狠狠婷婷综合久久图片| www.熟女人妻精品国产| 久久久久久人人人人人| 久久久国产精品麻豆| 91大片在线观看| 一进一出抽搐gif免费好疼| 久久久国产精品麻豆| 久久久久性生活片| 日韩欧美 国产精品| 免费看美女性在线毛片视频| 日本在线视频免费播放| 国产视频一区二区在线看| 国产一区二区三区在线臀色熟女| 色av中文字幕| 国产99白浆流出| 日日干狠狠操夜夜爽| 日韩欧美在线二视频| 亚洲精品在线美女| 日韩欧美免费精品| 看免费av毛片| 国产一区二区在线观看日韩 | 国产精品野战在线观看| 丁香欧美五月| 欧美黑人巨大hd| 一本大道久久a久久精品| 久久精品国产亚洲av香蕉五月| 中国美女看黄片| 国产精品野战在线观看| 亚洲自拍偷在线| 久久久水蜜桃国产精品网| 欧美中文日本在线观看视频| 欧美三级亚洲精品| 免费看十八禁软件| 亚洲成av人片在线播放无| 免费在线观看影片大全网站| 男女视频在线观看网站免费 | 国产亚洲精品av在线| 国产私拍福利视频在线观看| 99久久国产精品久久久| 婷婷亚洲欧美| 婷婷六月久久综合丁香| 欧美黄色淫秽网站| 精品一区二区三区视频在线观看免费| 变态另类丝袜制服| 亚洲国产精品久久男人天堂| 男人舔奶头视频| 欧美性猛交╳xxx乱大交人| 日韩欧美在线二视频| 亚洲一区二区三区色噜噜| 日韩成人在线观看一区二区三区| 99在线视频只有这里精品首页| 99国产精品99久久久久| 村上凉子中文字幕在线| 欧美日韩一级在线毛片| 在线a可以看的网站| 18禁美女被吸乳视频| 国产aⅴ精品一区二区三区波| 老熟妇仑乱视频hdxx| 男女之事视频高清在线观看| 一级毛片精品| 久久精品国产亚洲av香蕉五月| 午夜精品在线福利| 老汉色∧v一级毛片| 黄色视频,在线免费观看| 色尼玛亚洲综合影院| 在线观看免费视频日本深夜| av视频在线观看入口| 91大片在线观看| 午夜福利在线在线| 国产亚洲精品久久久久5区| 日日夜夜操网爽| 男女下面进入的视频免费午夜| 一级片免费观看大全| 国产午夜精品论理片| 国产av一区二区精品久久| 久久九九热精品免费| 国产亚洲av高清不卡| 国产精品日韩av在线免费观看| 国产欧美日韩一区二区精品| 一级毛片女人18水好多| 丰满人妻一区二区三区视频av | 久久久久久久久免费视频了| ponron亚洲| 一夜夜www| 精品国内亚洲2022精品成人| 亚洲国产精品999在线| 最新在线观看一区二区三区| 日韩中文字幕欧美一区二区| 最近视频中文字幕2019在线8| 精品高清国产在线一区| 他把我摸到了高潮在线观看| 国产97色在线日韩免费| 午夜影院日韩av| 久久九九热精品免费| 欧美乱色亚洲激情| 少妇被粗大的猛进出69影院| 中出人妻视频一区二区| 亚洲全国av大片| 一本久久中文字幕| 老熟妇乱子伦视频在线观看| 一a级毛片在线观看| 色综合欧美亚洲国产小说| 亚洲精品在线美女| 亚洲国产欧美一区二区综合| 国产精品一区二区精品视频观看| 欧美性猛交╳xxx乱大交人| 亚洲天堂国产精品一区在线| 久久精品国产亚洲av香蕉五月| 国产成人av教育| av免费在线观看网站| 51午夜福利影视在线观看| 一级毛片女人18水好多| 精华霜和精华液先用哪个| 一级毛片高清免费大全| 国产午夜精品久久久久久| 女警被强在线播放| 97碰自拍视频| 久久午夜亚洲精品久久| 免费在线观看影片大全网站| 亚洲黑人精品在线| 国产熟女xx| 国产精品九九99| 色尼玛亚洲综合影院| 久久亚洲精品不卡| 亚洲第一欧美日韩一区二区三区| 国产高清视频在线播放一区| 免费看美女性在线毛片视频| 可以在线观看毛片的网站| 久久这里只有精品中国| 99热6这里只有精品| 亚洲人成网站在线播放欧美日韩| 99国产综合亚洲精品| 亚洲国产日韩欧美精品在线观看 | 90打野战视频偷拍视频| 淫妇啪啪啪对白视频| 免费看a级黄色片| 国产精品美女特级片免费视频播放器 | 狂野欧美激情性xxxx| 国产一区二区三区视频了| 高清毛片免费观看视频网站| www日本在线高清视频| 一夜夜www| 看黄色毛片网站| 国产欧美日韩一区二区三| 成年免费大片在线观看| 成人国产一区最新在线观看| 黄片大片在线免费观看| 人妻夜夜爽99麻豆av| 欧美成人一区二区免费高清观看 | 国产av麻豆久久久久久久| 亚洲精品在线观看二区| 欧美色欧美亚洲另类二区| 国产一级毛片七仙女欲春2| 亚洲欧美一区二区三区黑人| 777久久人妻少妇嫩草av网站| 极品教师在线免费播放| 一级毛片女人18水好多| 巨乳人妻的诱惑在线观看| 美女 人体艺术 gogo| 国产视频一区二区在线看| 9191精品国产免费久久| 日本五十路高清| 后天国语完整版免费观看| 中出人妻视频一区二区| 精品人妻1区二区| 亚洲免费av在线视频| 999精品在线视频| 亚洲国产精品成人综合色| 久热爱精品视频在线9| 一个人免费在线观看电影 | 国产人伦9x9x在线观看| av天堂在线播放| 老司机午夜福利在线观看视频| av欧美777| 国产又黄又爽又无遮挡在线| 女生性感内裤真人,穿戴方法视频| 天堂av国产一区二区熟女人妻 | 九色国产91popny在线| АⅤ资源中文在线天堂| 国产熟女xx| 国产成人精品无人区| 精品少妇一区二区三区视频日本电影| 久久久国产成人免费| 国产精品久久久久久亚洲av鲁大| 欧美一级毛片孕妇| 他把我摸到了高潮在线观看| 成人一区二区视频在线观看| 我要搜黄色片| 午夜激情av网站| 毛片女人毛片| 99国产极品粉嫩在线观看| 又爽又黄无遮挡网站| 一a级毛片在线观看| 中文在线观看免费www的网站 | 欧美绝顶高潮抽搐喷水| 69av精品久久久久久| 国产免费男女视频| 精品久久久久久久末码| 欧美日韩瑟瑟在线播放| 亚洲狠狠婷婷综合久久图片| 757午夜福利合集在线观看| 91成年电影在线观看| 观看免费一级毛片| 国产精品日韩av在线免费观看| 久久久精品大字幕| 男女做爰动态图高潮gif福利片| 91大片在线观看| 男女做爰动态图高潮gif福利片| 两人在一起打扑克的视频| av欧美777| 特大巨黑吊av在线直播| 2021天堂中文幕一二区在线观| 亚洲精品国产一区二区精华液| 亚洲欧美日韩无卡精品| 99在线视频只有这里精品首页| ponron亚洲| 99在线视频只有这里精品首页| 波多野结衣高清无吗| 黄色片一级片一级黄色片| 精品一区二区三区视频在线观看免费| bbb黄色大片| 两个人的视频大全免费| 91国产中文字幕| 神马国产精品三级电影在线观看 | 在线观看舔阴道视频| 亚洲av熟女| 不卡av一区二区三区| 舔av片在线| 免费在线观看视频国产中文字幕亚洲| 天天添夜夜摸| 欧美激情久久久久久爽电影| 日日干狠狠操夜夜爽| 香蕉久久夜色| 免费搜索国产男女视频| 亚洲色图av天堂| 成年人黄色毛片网站| 亚洲中文字幕一区二区三区有码在线看 | 日本精品一区二区三区蜜桃| 级片在线观看| 丰满的人妻完整版| 国产激情欧美一区二区| 国产一区二区激情短视频| 亚洲自偷自拍图片 自拍| 男女之事视频高清在线观看| 女人爽到高潮嗷嗷叫在线视频| 免费人成视频x8x8入口观看| 在线十欧美十亚洲十日本专区| 亚洲精品美女久久av网站| 色综合婷婷激情| av有码第一页| 亚洲欧美日韩东京热| 搡老妇女老女人老熟妇| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲国产精品成人综合色| 欧美 亚洲 国产 日韩一| 久久久久久久久久黄片| 成人永久免费在线观看视频| 少妇裸体淫交视频免费看高清 | 18禁国产床啪视频网站| 欧洲精品卡2卡3卡4卡5卡区| av天堂在线播放| 亚洲av成人精品一区久久| 宅男免费午夜| 亚洲人成网站高清观看| 欧美性猛交╳xxx乱大交人| 成人永久免费在线观看视频| a级毛片a级免费在线| 国产精品1区2区在线观看.| 色精品久久人妻99蜜桃| 亚洲午夜理论影院| 日韩大码丰满熟妇| 欧美三级亚洲精品| 国产成人欧美在线观看| www日本黄色视频网| 又大又爽又粗| 成人特级黄色片久久久久久久| 国产精华一区二区三区| 午夜免费成人在线视频| 免费看日本二区| 最好的美女福利视频网| 欧洲精品卡2卡3卡4卡5卡区| 免费av毛片视频| 中文字幕av在线有码专区| 叶爱在线成人免费视频播放| 亚洲国产欧美网| 这个男人来自地球电影免费观看| 麻豆成人av在线观看| 国产日本99.免费观看| 欧美日本亚洲视频在线播放| 观看免费一级毛片| 国产午夜精品久久久久久| 日日夜夜操网爽| 啪啪无遮挡十八禁网站| 女同久久另类99精品国产91| 国产区一区二久久| 五月伊人婷婷丁香| 久久久精品大字幕| videosex国产| 免费观看精品视频网站| 叶爱在线成人免费视频播放| 亚洲专区字幕在线| 麻豆国产av国片精品| 大型黄色视频在线免费观看| 亚洲国产中文字幕在线视频| 香蕉av资源在线| 国产成人啪精品午夜网站| 90打野战视频偷拍视频| 欧美乱色亚洲激情| 1024手机看黄色片| 伊人久久大香线蕉亚洲五| 日本三级黄在线观看| 美女高潮喷水抽搐中文字幕| 黄色a级毛片大全视频| 久久久精品欧美日韩精品| 亚洲中文日韩欧美视频| www.www免费av| 久久精品91蜜桃| 成年女人毛片免费观看观看9| 亚洲美女视频黄频| 一卡2卡三卡四卡精品乱码亚洲| 麻豆国产97在线/欧美 | 悠悠久久av| 精品午夜福利视频在线观看一区| 男女下面进入的视频免费午夜| 精品欧美一区二区三区在线| 国产欧美日韩一区二区三| 国产aⅴ精品一区二区三区波| 国产黄色小视频在线观看| 国内精品一区二区在线观看| 免费在线观看黄色视频的| 精品久久久久久,| 老司机午夜十八禁免费视频| 黄色女人牲交| 亚洲欧美日韩高清在线视频| 丝袜美腿诱惑在线| 不卡av一区二区三区| 国产精品 国内视频| 亚洲熟女毛片儿| 亚洲人与动物交配视频| 丁香欧美五月| 成人18禁在线播放| 欧美日韩乱码在线| 日韩精品中文字幕看吧| 很黄的视频免费| 亚洲国产精品999在线| 免费在线观看日本一区| 久久中文字幕一级| 黄色 视频免费看| 亚洲免费av在线视频| 日本五十路高清| 国产aⅴ精品一区二区三区波| 亚洲熟妇中文字幕五十中出| 在线国产一区二区在线| 亚洲欧美一区二区三区黑人| 亚洲熟女毛片儿| 国产成人一区二区三区免费视频网站| 熟女电影av网|