• <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高清涩受| 国产在线免费精品| 国产精品久久久av美女十八| 国产精品一二三区在线看| a 毛片基地| 美女福利国产在线| 人人澡人人妻人| 国产男女内射视频| 亚洲av男天堂| 国产欧美日韩综合在线一区二区| 一级毛片黄色毛片免费观看视频| kizo精华| tube8黄色片| 久久精品亚洲熟妇少妇任你| 日本欧美国产在线视频| 国产精品一区二区在线不卡| 男女下面插进去视频免费观看| 国产一区二区三区综合在线观看| 精品亚洲成a人片在线观看| av在线播放精品| 成人三级做爰电影| 69精品国产乱码久久久| 天天躁夜夜躁狠狠久久av| 亚洲成色77777| 久久久精品免费免费高清| 国产欧美日韩一区二区三区在线| 亚洲精品av麻豆狂野| 操出白浆在线播放| 电影成人av| 超色免费av| 性高湖久久久久久久久免费观看| 国产xxxxx性猛交| 91成人精品电影| 国产精品国产三级国产专区5o| 国产一区二区 视频在线| 美女扒开内裤让男人捅视频| 麻豆国产av国片精品| 亚洲精品自拍成人| 久久精品国产亚洲av高清一级| 99re6热这里在线精品视频| av欧美777| 无遮挡黄片免费观看| 亚洲精品久久午夜乱码| 18禁裸乳无遮挡动漫免费视频| 亚洲精品自拍成人| 母亲3免费完整高清在线观看| 19禁男女啪啪无遮挡网站| 精品少妇黑人巨大在线播放| 深夜精品福利| 日韩 欧美 亚洲 中文字幕| 啦啦啦啦在线视频资源| 亚洲欧美日韩高清在线视频 | 国产男女内射视频| 一级毛片 在线播放| 久久精品久久久久久噜噜老黄| 交换朋友夫妻互换小说| 精品少妇黑人巨大在线播放| 国产野战对白在线观看| 女性被躁到高潮视频| 国产精品人妻久久久影院| 国产欧美日韩综合在线一区二区| 成年美女黄网站色视频大全免费| 久久亚洲国产成人精品v| 男女午夜视频在线观看| 国产精品99久久99久久久不卡| 精品国产国语对白av| 我的亚洲天堂| 免费av中文字幕在线| 午夜老司机福利片| 又大又爽又粗| 久久国产精品大桥未久av| 亚洲精品一二三| 丁香六月天网| 美女视频免费永久观看网站| a级毛片黄视频| 91精品伊人久久大香线蕉| av天堂久久9| 爱豆传媒免费全集在线观看| 日本午夜av视频| 亚洲一码二码三码区别大吗| videosex国产| 在线精品无人区一区二区三| avwww免费| 狠狠精品人妻久久久久久综合| 韩国高清视频一区二区三区| 飞空精品影院首页| 欧美日韩亚洲高清精品| 老司机靠b影院| 亚洲 国产 在线| www.自偷自拍.com| 国产在线观看jvid| 日韩视频在线欧美| 一区二区日韩欧美中文字幕| 考比视频在线观看| 热99久久久久精品小说推荐| 在线精品无人区一区二区三| 亚洲中文av在线| 美女脱内裤让男人舔精品视频| 国产成人精品无人区| 国产成人精品久久久久久| 日韩av不卡免费在线播放| 老司机影院毛片| 精品久久久久久久毛片微露脸 | a 毛片基地| 精品福利观看| 精品国产国语对白av| 久久精品国产综合久久久| 国产精品熟女久久久久浪| 午夜福利影视在线免费观看| 高清欧美精品videossex| 亚洲国产欧美网| 人妻人人澡人人爽人人| 美女福利国产在线| 又黄又粗又硬又大视频| 亚洲欧美一区二区三区久久| 欧美性长视频在线观看| 欧美日韩国产mv在线观看视频| 大片免费播放器 马上看| 99九九在线精品视频| 在线观看免费高清a一片| 午夜免费鲁丝| 午夜91福利影院| 国产精品一区二区精品视频观看| 一边摸一边做爽爽视频免费| 侵犯人妻中文字幕一二三四区| 侵犯人妻中文字幕一二三四区| 精品人妻一区二区三区麻豆| 纯流量卡能插随身wifi吗| 久久精品国产a三级三级三级| 丁香六月天网| 欧美日韩av久久| 女人爽到高潮嗷嗷叫在线视频| 99久久综合免费| 美国免费a级毛片| 亚洲精品日本国产第一区| 一级a爱视频在线免费观看| 国产黄色视频一区二区在线观看| 超碰97精品在线观看| 蜜桃国产av成人99| 亚洲熟女毛片儿| 婷婷成人精品国产| 国产成人欧美| 中文字幕亚洲精品专区| 悠悠久久av| 在线观看www视频免费| 亚洲专区国产一区二区| 国产日韩欧美亚洲二区| 曰老女人黄片| 成年人午夜在线观看视频| 无限看片的www在线观看| 免费在线观看视频国产中文字幕亚洲 | 夫妻午夜视频| 91精品国产国语对白视频| 婷婷色综合www| av不卡在线播放| 国产高清不卡午夜福利| 好男人电影高清在线观看| 少妇粗大呻吟视频| 亚洲欧美一区二区三区久久| 久久午夜综合久久蜜桃| 爱豆传媒免费全集在线观看| 免费久久久久久久精品成人欧美视频| 中文字幕最新亚洲高清| 99国产精品99久久久久| 精品国产国语对白av| 777久久人妻少妇嫩草av网站| 亚洲国产欧美日韩在线播放| 日本欧美视频一区| 午夜免费男女啪啪视频观看| 午夜福利在线免费观看网站| 波野结衣二区三区在线| 一区二区三区精品91| 国产欧美日韩综合在线一区二区| 久久这里只有精品19| 99久久精品国产亚洲精品| 午夜久久久在线观看| www.熟女人妻精品国产| 在线观看国产h片| 脱女人内裤的视频| www日本在线高清视频| 国产精品一区二区精品视频观看| 亚洲av综合色区一区| 极品少妇高潮喷水抽搐| 日韩电影二区| 国产精品av久久久久免费| 亚洲天堂av无毛| 国产成人一区二区三区免费视频网站 | 午夜两性在线视频| 国产精品熟女久久久久浪| 777久久人妻少妇嫩草av网站| 少妇粗大呻吟视频| 日本91视频免费播放| 午夜免费男女啪啪视频观看| 视频区欧美日本亚洲| 欧美激情极品国产一区二区三区| 免费观看人在逋| 国产亚洲精品久久久久5区| 日本a在线网址| 亚洲美女黄色视频免费看| 欧美日韩福利视频一区二区| 2018国产大陆天天弄谢| 丝袜美足系列| 好男人电影高清在线观看| 精品人妻熟女毛片av久久网站| 永久免费av网站大全| 99热网站在线观看| 亚洲 欧美一区二区三区| 丰满少妇做爰视频| 婷婷色综合www| 久久久久久人人人人人| 欧美黑人精品巨大| 黄频高清免费视频| 99热国产这里只有精品6| 亚洲国产精品一区三区| 手机成人av网站| 精品亚洲乱码少妇综合久久| 欧美老熟妇乱子伦牲交| 五月开心婷婷网| 亚洲av欧美aⅴ国产| 免费高清在线观看视频在线观看| 欧美精品一区二区大全| 一区二区三区精品91| e午夜精品久久久久久久| 精品一区在线观看国产| 亚洲欧美一区二区三区久久| 水蜜桃什么品种好| 嫩草影视91久久| 久久亚洲精品不卡| 成在线人永久免费视频| 欧美人与性动交α欧美精品济南到| 高清av免费在线| 女人爽到高潮嗷嗷叫在线视频| 亚洲欧美精品综合一区二区三区| 国产在线一区二区三区精| 蜜桃在线观看..| 90打野战视频偷拍视频| 色精品久久人妻99蜜桃| 在线观看免费午夜福利视频| 大码成人一级视频| 黄色a级毛片大全视频| 午夜福利影视在线免费观看| 久久国产精品影院| 国产深夜福利视频在线观看| 久热这里只有精品99| 亚洲av日韩精品久久久久久密 | 纵有疾风起免费观看全集完整版| 9191精品国产免费久久| 91九色精品人成在线观看| 欧美日韩精品网址| 91精品伊人久久大香线蕉| 女性被躁到高潮视频| 婷婷成人精品国产| 午夜免费观看性视频| 久久国产精品大桥未久av| 久9热在线精品视频| 国产片内射在线| 亚洲av片天天在线观看| 亚洲国产欧美日韩在线播放| 宅男免费午夜| 日本一区二区免费在线视频| 欧美精品一区二区免费开放| 国产精品人妻久久久影院| 午夜日韩欧美国产| 美女视频免费永久观看网站| 亚洲视频免费观看视频| 99精国产麻豆久久婷婷| 波多野结衣一区麻豆| 午夜免费鲁丝| 亚洲中文字幕日韩| 精品一品国产午夜福利视频| 久久女婷五月综合色啪小说| 欧美精品一区二区免费开放| 国产高清videossex| 男人爽女人下面视频在线观看| 男人舔女人的私密视频| 色94色欧美一区二区| 日韩制服骚丝袜av| 一级毛片电影观看| 91九色精品人成在线观看| 纯流量卡能插随身wifi吗| 精品一品国产午夜福利视频| 99国产综合亚洲精品| 日本午夜av视频| www.精华液| 我要看黄色一级片免费的| 韩国精品一区二区三区| 午夜福利,免费看| 亚洲五月婷婷丁香| 亚洲精品国产一区二区精华液| 91九色精品人成在线观看| 亚洲国产成人一精品久久久| 啦啦啦 在线观看视频| 久9热在线精品视频| 亚洲 欧美一区二区三区| 久久久久久久大尺度免费视频| av有码第一页| 成年人黄色毛片网站| 亚洲 欧美一区二区三区| 免费在线观看影片大全网站 | 波野结衣二区三区在线| 一级,二级,三级黄色视频| 亚洲国产欧美在线一区| 一级毛片我不卡| av线在线观看网站| 国产精品九九99| 只有这里有精品99| 国产精品 欧美亚洲| 国产成人a∨麻豆精品| 国产日韩欧美视频二区| 国产xxxxx性猛交| 99国产综合亚洲精品| 欧美日韩一级在线毛片| 中国国产av一级| www日本在线高清视频| av在线app专区| 成人免费观看视频高清| 日本午夜av视频| 亚洲情色 制服丝袜| 精品福利观看| 美女午夜性视频免费| 精品国产一区二区三区久久久樱花| 国产熟女午夜一区二区三区| 制服人妻中文乱码| 一二三四在线观看免费中文在| 少妇精品久久久久久久| 啦啦啦 在线观看视频| 大香蕉久久成人网| 多毛熟女@视频| 老司机深夜福利视频在线观看 | 亚洲专区国产一区二区| 国产真人三级小视频在线观看| 99精国产麻豆久久婷婷| 亚洲av成人精品一二三区| 人人妻人人添人人爽欧美一区卜| 一级黄片播放器| 80岁老熟妇乱子伦牲交| 亚洲七黄色美女视频| 欧美黄色淫秽网站| 又黄又粗又硬又大视频| 午夜福利视频在线观看免费| 色网站视频免费| 亚洲av成人精品一二三区| 男男h啪啪无遮挡| 日韩制服丝袜自拍偷拍| 99国产精品一区二区蜜桃av | 狂野欧美激情性xxxx| 老司机亚洲免费影院| 天天影视国产精品| 亚洲欧美精品自产自拍| 精品欧美一区二区三区在线| 在现免费观看毛片| 美女午夜性视频免费| 欧美精品一区二区大全| 国产亚洲午夜精品一区二区久久| 国产精品免费视频内射| 国产精品国产av在线观看| 亚洲精品日韩在线中文字幕| 伦理电影免费视频| 亚洲精品成人av观看孕妇| 日本黄色日本黄色录像| 日本一区二区免费在线视频| 精品少妇内射三级| 国产在线一区二区三区精| av欧美777| 美女高潮到喷水免费观看| bbb黄色大片| 好男人视频免费观看在线| 少妇 在线观看| 免费黄频网站在线观看国产| 欧美激情极品国产一区二区三区| 亚洲欧美一区二区三区黑人| 高清av免费在线| 少妇 在线观看| 人妻一区二区av| 99热国产这里只有精品6| 99久久综合免费| 精品亚洲成a人片在线观看| 亚洲av日韩在线播放| 精品少妇久久久久久888优播| 极品少妇高潮喷水抽搐| 黄频高清免费视频| 多毛熟女@视频| 免费看av在线观看网站| 大香蕉久久网| 亚洲伊人久久精品综合| 日本欧美国产在线视频| 久久精品aⅴ一区二区三区四区| 18在线观看网站| 高潮久久久久久久久久久不卡| 国产爽快片一区二区三区| 99国产综合亚洲精品| 亚洲国产欧美在线一区| 久久久久国产精品人妻一区二区| 人人妻人人爽人人添夜夜欢视频| 国产精品一国产av| 国精品久久久久久国模美| 女人高潮潮喷娇喘18禁视频| 亚洲精品自拍成人| 丝袜美足系列| 少妇精品久久久久久久| 国产成人精品久久二区二区91| √禁漫天堂资源中文www| 精品少妇一区二区三区视频日本电影| 男女床上黄色一级片免费看| 日本一区二区免费在线视频| 91九色精品人成在线观看| 免费少妇av软件| 日韩人妻精品一区2区三区| 人人妻人人添人人爽欧美一区卜| www.自偷自拍.com| 亚洲av日韩在线播放| 国产片内射在线| 老司机靠b影院| 香蕉丝袜av| 2018国产大陆天天弄谢| 午夜两性在线视频| 性高湖久久久久久久久免费观看| 中文字幕色久视频| 日韩 亚洲 欧美在线| 日本黄色日本黄色录像| 成人三级做爰电影| 热re99久久精品国产66热6| 美女脱内裤让男人舔精品视频| 大片电影免费在线观看免费| 在线观看www视频免费| 免费高清在线观看日韩| 国产视频首页在线观看| 乱人伦中国视频| 中文欧美无线码| 无限看片的www在线观看| 国产午夜精品一二区理论片| 如日韩欧美国产精品一区二区三区| 亚洲成人手机| 国产精品香港三级国产av潘金莲 | 人人妻人人添人人爽欧美一区卜| 亚洲国产成人一精品久久久| 母亲3免费完整高清在线观看| 99久久综合免费| 岛国毛片在线播放| 国产高清videossex| 韩国精品一区二区三区| 亚洲精品一区蜜桃| 亚洲人成电影观看| 成年人午夜在线观看视频| 大话2 男鬼变身卡| 亚洲自偷自拍图片 自拍| 亚洲av电影在线进入| av在线播放精品| 国产伦理片在线播放av一区| 中文字幕最新亚洲高清| 欧美97在线视频| 色婷婷av一区二区三区视频| 欧美 日韩 精品 国产| 欧美日韩视频精品一区| 天堂俺去俺来也www色官网| 亚洲五月婷婷丁香| 建设人人有责人人尽责人人享有的| 汤姆久久久久久久影院中文字幕| 久热爱精品视频在线9| 亚洲黑人精品在线| 久久精品国产综合久久久| 欧美激情极品国产一区二区三区| 黄片播放在线免费| 欧美97在线视频| 一区在线观看完整版| 国产成人欧美| 国产高清视频在线播放一区 | 国产亚洲欧美在线一区二区| 国产精品 国内视频| 丝袜美足系列| 亚洲第一青青草原| 亚洲av男天堂| 午夜两性在线视频| 免费女性裸体啪啪无遮挡网站| 欧美久久黑人一区二区| 亚洲国产毛片av蜜桃av| 别揉我奶头~嗯~啊~动态视频 | 手机成人av网站| 亚洲七黄色美女视频| 午夜福利影视在线免费观看| 欧美日本中文国产一区发布| 欧美日韩亚洲国产一区二区在线观看 | 亚洲久久久国产精品| 久久久精品国产亚洲av高清涩受| 在线 av 中文字幕| 国产欧美亚洲国产| 欧美精品av麻豆av| 90打野战视频偷拍视频| 午夜激情久久久久久久| 999久久久国产精品视频| 国产黄色免费在线视频| 国产片特级美女逼逼视频| 视频区欧美日本亚洲| 久久综合国产亚洲精品| 国产主播在线观看一区二区 | 国产视频一区二区在线看| 久久人妻熟女aⅴ| 亚洲中文日韩欧美视频| 亚洲av美国av| 熟女少妇亚洲综合色aaa.| 人人妻,人人澡人人爽秒播 | 精品一区在线观看国产| 日本欧美视频一区| 国产一区二区在线观看av| 建设人人有责人人尽责人人享有的| 满18在线观看网站| 99久久99久久久精品蜜桃| 欧美日本中文国产一区发布| 在线观看国产h片| 久久狼人影院| 国产成人欧美在线观看 | 少妇粗大呻吟视频| 亚洲黑人精品在线| 一级毛片我不卡| 日韩,欧美,国产一区二区三区| 国产高清videossex| 国产又爽黄色视频| 亚洲av男天堂| 狠狠精品人妻久久久久久综合| 丝袜人妻中文字幕| www日本在线高清视频| 美女脱内裤让男人舔精品视频| 人人妻人人添人人爽欧美一区卜| 考比视频在线观看| 国产精品亚洲av一区麻豆| 18禁观看日本| 观看av在线不卡| 在线观看www视频免费| 男女边摸边吃奶| 亚洲一区中文字幕在线| 国产有黄有色有爽视频| 在线观看免费高清a一片| 人妻一区二区av| 国产日韩一区二区三区精品不卡| 18禁国产床啪视频网站| 成人18禁高潮啪啪吃奶动态图| 午夜福利乱码中文字幕| 亚洲国产欧美在线一区| 国产在线一区二区三区精| 婷婷色综合大香蕉| 天堂俺去俺来也www色官网| 亚洲av欧美aⅴ国产| a 毛片基地| 1024视频免费在线观看| 男人爽女人下面视频在线观看| 成年美女黄网站色视频大全免费| 精品国产国语对白av| 叶爱在线成人免费视频播放| 女人爽到高潮嗷嗷叫在线视频| 国产一区二区三区综合在线观看| 亚洲情色 制服丝袜| 国产av一区二区精品久久| 色综合欧美亚洲国产小说| 国产又爽黄色视频| 午夜激情久久久久久久| 在线观看免费视频网站a站| 高清av免费在线| 精品视频人人做人人爽| 女警被强在线播放| 亚洲精品在线美女| 亚洲激情五月婷婷啪啪| 亚洲欧美一区二区三区国产| 亚洲av欧美aⅴ国产| 水蜜桃什么品种好| 丝袜人妻中文字幕| 欧美日韩黄片免| 香蕉国产在线看| 老汉色∧v一级毛片| 日本午夜av视频| 欧美 日韩 精品 国产| 后天国语完整版免费观看| 婷婷成人精品国产| 国产一区二区三区av在线| 妹子高潮喷水视频| 亚洲色图 男人天堂 中文字幕| 在线观看www视频免费| 亚洲美女黄色视频免费看| 在线观看www视频免费| 国产精品成人在线| 国产欧美日韩一区二区三区在线| 亚洲欧洲国产日韩| 亚洲av成人不卡在线观看播放网 | 视频在线观看一区二区三区| 亚洲欧美激情在线| 大香蕉久久成人网| 亚洲国产成人一精品久久久| 中文字幕人妻熟女乱码| 成人国产一区最新在线观看 | 久久久久久免费高清国产稀缺| 久久九九热精品免费| 三上悠亚av全集在线观看| 午夜老司机福利片| 韩国精品一区二区三区| 亚洲成国产人片在线观看| tube8黄色片| 成年女人毛片免费观看观看9 | 欧美日韩av久久| 日韩av在线免费看完整版不卡| 秋霞在线观看毛片| 男女边摸边吃奶| 欧美日韩黄片免| 99热国产这里只有精品6| 成人国产一区最新在线观看 | 婷婷色av中文字幕| 婷婷色综合大香蕉| 日韩制服丝袜自拍偷拍| 十分钟在线观看高清视频www| 免费少妇av软件| 一级毛片电影观看| 欧美黄色片欧美黄色片| 日本av手机在线免费观看| 国产成人精品久久二区二区免费| 精品国产乱码久久久久久男人| 亚洲精品日本国产第一区|