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

    Dynamic behavior of the cyanobacterial circadian clock with regulation of CikA?

    2021-10-28 07:18:30YingLi李瑩GuangKunZhang張廣鹍andYanMingGe葛焰明
    Chinese Physics B 2021年10期
    關(guān)鍵詞:李瑩

    Ying Li(李瑩), Guang-Kun Zhang(張廣鹍), and Yan-Ming Ge(葛焰明)

    College of Information Technology,Shanghai Ocean University,Shanghai 201306,China

    Keywords: cyanobacterial circadian clock,mathematical model,adaptability,sensitivity analysis

    1. Introduction

    Circadian clocks are the fundamental process of many organisms to keep the physical and behavior order and adapt to daily environmental variations, which are intrinsic selfsustained oscillators with a period of approximately 24 h. The great progress in molecular biology experiments has enabled researchers to explore the components and molecular mechanisms involved in circadian clocks.[1]Among a wide variety of organisms,cyanobacteria are the simplest organisms to have circadian clocks.

    In the cyanobacterial circadian clock, the gene clusterkaiABCand its product proteins KaiA, KaiB, and KaiC are essential for generating the rhythm.[2]KaiC plays a central role and undergoes rhythmic activities in transcription,translation, phosphorylation, and autodephosphorylation in the presence of both KaiA and KaiB.[3,4]The cyanobacterial circadian clock is composed of two coupled components: a post-translational oscillator (PTO) and a transcriptional/translational feedback loop(TTFL).These two components have a hierarchical relationship. The PTO is the core pacemaker to which the TTFL is attached. The TTFL feeds back to the PTO by the synthesis of clock proteins.[5]

    In the TTFL,the clock genes,kaiA,kaiB,andkaiC,form a gene cluster, wherekaiBandkaiCare co-transcribed askaiABCmRNA.[2]Biological experiment shows that KaiC over expression consistently reduces kaiBC promoter activity, which implies KaiC regulates kaiBC transcription negatively.[2]This self-repression of KaiC indicates that the negative feedback transcriptional regulation generates oscillation ofkaigene expression. Additionally, KaiA over expression shows different effects on kaiBC transcription depending on KaiC; the transcriptional regulation by KaiC also showskaiAdependency.[6]Based on the fact that KaiA enhances KaiC phosphorylation, this cooperative regulation by KaiA and KaiC is considered to be realized by KaiC phosphorylation.[7]

    The PTO is a self-sustained core pacemaker which consists of three core proteins: KaiA, KaiB, and KaiC.[8]KaiC can autophosphorylate and autodephosphorylate. The other clock proteins KaiA and KaiB modulate KaiC autophosphorylation. KaiA enhances the autophosphorylation of KaiC by binding to the A-loop, whereas KaiB inhibits this action of KaiA by sequestering KaiA from the A-loop,[6]which generates the oscillation of KaiC phosphorylation with a circadian period.[9]

    Although the Kai proteins alone produce a robust circadian rhythm in vitro,there are several other proteins that have a great influence on the rhythms.[10]CikA,a signal transduction protein, plays an important role in the three processes of cyanobacterial circadian clock: input, oscillator and output.In the input pathway, CikA senses environmental light information indirectly, which entrains the circadian oscillator.[11]CikA also functions as an output pathway component by acting as a phosphatase against RpaA,a transcription factor regulated by the circadian clock.[12]In the central oscillator,CikA affects the phosphorylation and dephosphorylation of KaiC indirectly. Structural analysis reveals that both CikA and KaiA bind competitively to the same site on KaiB that is used to sequester KaiA,[13]which changes the ability of KaiA to stimulate KaiC phosphorylation. This regulation of CikA affects both the amplitude and the period of the circadian oscillation.[14]Considering the important role of CikA on the circadian clock,many biological results have been found.However, the molecular mechanism is not fully understood.To our knowledge, there is no research work on it from the dynamical angle,which may help us more deeply understand the role of CikA in the cyanobacterial circadian clock.

    On the investigation of the regulatory mechanism of circadian clocks, mathematical models are useful and powerful,which have been widely used.[11,15–18]To explore the regulation mechanisms of CikA on the phosphorylation and dephosphorylation of KaiC,we bring forward a mathematical model with incorporation of additional regulation of CikA into the cyanobacterial circadian clock model. We will carry out a detailed study of the model to reveal the crucial role of CikA on the cyanobacterial circadian clock. The rest of this article is organized as follows. In Section 2, a detailed mathematical model is developed. Section 3 shows the numerical simulation and analysis. The conclusion and discussion are presented in Section 4.

    2. Description of the mathematical model

    When CikA regulation is incorporated into the PTO and TTFL circadian network,[8]the model of the cyanobacterial circadian oscillator is schematized in Fig.1. Cooperative regulation of KaiA and KaiC on kaiBC is realized by phosphorylated KaiC.KaiC can autophosphorylate and autodephosphorylate. KaiA promotes the phosphorylation of KaiC,and KaiB inhibits the phosphorylation in a competitive manner. However,CikA weakens the inhibition of KaiB.

    It has been observed that the amount of KaiA in a single cell remain constant at a low level in the cytosol, and the amount of KaiB is always kept to be proportional to that of KaiC.[19]The model is described by a three-dimensional system of ordinary differential equations which include the temporal dynamics of kaiBC mRNA and KaiC in two forms: unphosphorylated and phosphorylated. This system is described as follows:

    whereX1,X2, andX3are the concentration of kaiBC mRNA,unphosphorylated KaiC(NP-KaiC)and phosphorylated KaiC(P-KaiC),respectively;δ1,δ2,andδ3are the degradation rates of kaiBC,NP-KaiC,and P-KaiC,respectively;pis the translation rate of kaiBC.

    Fig.1. Schematic diagram for the cyanobacterial circadian clock with the regulation of CikA, which involves the processes of kaiBC transcription and KaiC phosphorylation.

    whereαandKare both nonnegative constants.

    Phosphorylation rate of KaiC is described by the Phos function. We suppose that KaiB interacts with KaiA which decreases KaiA amount that can stimulate the phosphorylation of KaiC,and CikA depresses the interaction. KaiA-mediated KaiC phosphorylation is described with the Michaelis–Menten function as follows:

    where [KaiA] represents the concentration of KaiA and it is a constant;kaandkmare the maximum phosphorylation rate per KaiA concentration activity and the Michaelis constant,respectively. The dephosphorylation rate is denoted askdX3,wherekdis the dephosphorylation rate per P-KaiC concentration. Without KaiB and CikA,the total change of phosphorylation and dephosphorylation is as follows:

    Considering the interaction of KaiB with KaiA decreases KaiA amount and [KaiB] is always proportional to the total amount of KaiC,we use the following equation to describe the total change of phosphorylation and dephosphorylation with KaiB:

    wheref,g,eandmare constants.

    In the following numerical simulations, we take the parameter values asδ1= 0.135 h?1,δ2= 0.1425 h?1,δ3=0.1425 h?1,p=36.1 h?1,α=5,K=15μM,ka=110 h?1,km= 1 μM,f= 16.48,g= 3.3,e= 1.6,m= 1,kd=0.105 h?1.

    3. Results

    3.1. The effect of CikA on the dynamics of the circadian clock

    For a mathematical model of the circadian clock, we should first validate whether its dynamic behaviors agree with experimental observations. In order to observe the effect of CikA on KaiC phosphorylation, we take various values of [CikA] in the circadian clock model with fixed [KaiA].Time evolution for the concentrations of P-KaiC is shown in Fig.2(a). The corresponding amplitude and period are shown in Fig. 2(b), from which we can see that CikA affects both the amplitude and the period of the circadian oscillation of KaiC phosphorylation. In particular, increasing amounts of CikA shorten the period of the KaiC phosphorylation oscillation,which is consistent with the experimental results.[14]

    Biological experiments showed that the amount of KaiA in the reaction mixture is an important factor for the circadian oscillator of KaiC phosphorylation.[21]Numerical simulation of the system shows that with the increase of [KaiA], the period decreases, and the amplitude increases first and then decreases (see Figs. 2(c)–2(d)). All these results are coincident with experimental findings.[14,21]

    Considering the competition of CikA with KaiA for the binding site on KaiB, to explore the effect of CikA, we decrease the concentration of KaiA in the reaction mixture while simultaneously increasing the concentration of CikA.The sum of[KaiA]and[CikA]is maintained. From Figs.2(e)–2(f),we can see that the phosphorylation state of KaiC maintains robust circadian rhythms with a stable circadian amplitude and a relatively stable period with the increase of [KaiA]. That is to say, in the presence of CikA, decreasing [KaiA] lengthens the period slightly. The ability of KaiA to stimulate KaiC phosphorylation is enhanced by CikA,which makes KaiA out of its sequestered state. Therefore,more CikA in the reaction mixture can release more free KaiA by occupying its binding site on KaiB,keeping the kinase activity of KaiC turned on by comparatively small amounts of KaiA.

    In order to see the dynamics more clearly, we give the bifurcation diagrams in Fig. 3, which show that only when[KaiA] is large enough the presented model produces selfsustained oscillations, no matter whether CikA presents or not. The bifurcation point in Fig. 3(b) is smaller than that in Fig.3(a),which means that the regulation of CikA expands the range of[KaiA]that makes the system oscillate.

    Fig. 2. Effects of varying the amounts of CikA and KaiA on KaiC phosphorylation. (a) Simulation of circadian oscillators of P-KaiC with various[CikA]and fixed[KaiA](1.2μM).(b)Amplitude and period of P-KaiC in(a). (c)Simulation of circadian oscillators of P-KaiC with various[KaiA]and without CikA.(d)Amplitude and period of P-KaiC in(c). (e)Simulation of circadian oscillators of P-KaiC with various[KaiA]and[CikA]. (f)Amplitude and period of P-KaiC in(e),with[CikA]+[KaiA]=1.8μM.

    Fig.3.Bifurcation diagrams with[KaiA]as the control parameter in the absence(a)and presence of CikA(b),respectively.

    3.2. Sensitivity analysis of the oscillation period to parameter values

    Robustness is an intrinsical property of a system, which describes the capability to keep performance facing up to the perturbation of system parameters. In order to study the robustness of the system with the regulation of CikA,we present the dynamic sensitivity analysis. The degree of robustness displays an inverse correlation with that of the corresponding sensitivity. The oscillator period is an important character to quantify the oscillator system,so we investigate the period sensitivity. The period sensitivity shows the variations of period coming from perturbations of the parameters in the model,which is described as follows:[22]

    whereTis the period of the circadian clock,andpiis the parameter with indexi. Considering that there are several parameters in the model,it is difficult to fully study the dynamic properties of the system in the whole parameter space.In order to evaluate the period’s sensitivity to changes of the parameters, we vary one parameter value at a time and keep the rest unchanged.

    The following system parameters are considered. The first kind is about the production of KaiC and dephosphorylation of KaiC, denoted by parameterspandkdrespectively,which reflect the roles of TTFL on the cyanobacterial circadian clock. The second kind is about the regulation of KaiA,KaiB and KaiC on the phosphorylation of KaiC, denoted by parameterka, which reflects the role of OPT on the circadian clock. The third kind is about the degradation rate of kaiBC,NP-KaiC,and P-KaiC,denoted byδ1,δ2,andδ3,respectively.

    Fig.4. Period sensitivity related to the parameters of the system. [(a),(b)] The conditions of +10% perturbation and ?10% perturbation to each parameter,respectively. The other parameters are fixed.

    Figure 4 shows the period sensitivity to the parameter perturbations with different concentrations of CikA. In Fig. 4(a)and 4(b), +10%perturbation and?10%perturbation to each parameter are introduced,respectively.From Fig.4,we obtain some interesting results. With the regulation of CikA,the system still oscillates when the parameters are set within a certain range. That is to say,the system is robust with the regulation of CikA. With or without CikA, the oscillator period is least sensitive to the perturbation of parameterspandkd, which are associated with TTFL.However,it is most sensitive to the perturbation of parameterka,which is associated to PTO.This finding indicates that the PTO as the core pacemaker plays the key role on deciding the oscillator period,while the TTFL is a slave oscillator. The presence of CikA decreases the sensitivity of the period to most of the parameters related to the TTFL and the sensitivity is inversely proportional to the concentration of CikA.On the contrary,the presence of CikA increases the sensitivity toka, which is proportional to the concentration of CikA.This result indicates that the regulation of CikA further strengthens the crucial role of the PTO. With the regulation of CikA,the system is more sensitive to the core PTO compared with that without CikA.That is to say,with the regulation of CikA,it is easier for the system to modulate its period against the environmental perturbations.

    3.3. CikA increases the adaptation to disrupted light-dark cycle

    Without the external time cues,the endogenous period of the cyanobacterial circadian clock is not exactly 24 h. The period becomes 24 h when it is exposed to natural environmental cycles. In other words, the oscillator is phase-locked or entrained in the external cycles.[23]The daily light-dark(LD) cycle is the strongest entraining signal. If entrainment is disrupted by a sudden shift in the phase of the LD cycle,the circadian oscillator will adjust its phase until the oscillator is reentrainmented or locked to the new LD cycle.[24]The time for the oscillator to get reentrainment reflects the ability of organisms to adapt to environmental changes,i.e.,adaptation.

    Experimental observations show that the light can increase KaiC phosphorylation rate,so we will simulate the effect of light by adding a term to the right-hand side of the differential equation for P-KaiC,i.e.,

    Fig.5.Reentrainment of the cyanobacterial circadian clock where the LD cycle is disturbed by 6 h delays.Simulation of circadian oscillators of P-KaiC in the absence of CikA(a)and in the presence of CikA(b),respectively. The solid and dotted lines indicate the oscillator in the normal LD cycle and the disturbed LD cycle respectively. (c)The corresponding circadian cycles in(a)–(b)are projected onto NP-KaiC-P-KaiC plane.(d)Phase shifts between the circadian oscillators in normal LD cycle and that in the disturbed LD cycle.

    Figures 5(a)and 5(b)show that the cyanobacterial circadian clock can be entrained by the LD cycle whenL0is set to appropriate values with or without CikA.Given the entrainment,we study the speed at which the system adjusts the phase to get reentrainment when there is a sudden shift in the phase of the LD cycle. The circadian oscillator runs for 20 days in a normal LD cycle. Then 6 h delay for light are added on the 21-th day. The dotted lines in Figs. 5(a) and 5(b) show the process of adjusting the phase of the system to reentrainment with and without CikA, respectively. Six-hour delay of light leads to phase delay of the oscillator. Comparing Figs. 5(a)and 5(b),we find that the phase shift is larger in the situation with the regulation of CikA than that without CikA. In order to see the adjustment speed more clearly, taking the peak of P-KaiC concentration as the phase of the circadian oscillator,we calculate the phase shifts between the circadian oscillator under the new LD cycle and the normal LD cycle, as shown in Fig.5(c). We can see that the presence of CikA speeds up the reentrainment of the circadian system, i.e., the regulation of CikA strengthens the adaptation of the system.

    In order to explore the dynamical mechanisms of the effect of CikA,we present the phase plane diagram in Fig.5(c).The presence of CikA weakens the inhibition of KaiB on KaiA which increases the activation of KaiA on promoting the phosphorylation of KaiC; however, the phosphorylated KaiC inhibits the transcription of KaiBC. Consequently, the regulation of CikA brings about a lower concentration of the total KaiC compared with the situation without CikA. Thus, the limit cycle is smaller in the situation with CikA than that without CikA, as shown in Fig. 5(c). On the other hand, the LD cycle directly affects the phosphorylation of KaiC.The phase changing of the LD cycle leads to the phase adjustments until the system resynchronizes with the new LD cycle. However,the amount of total KaiC that can be changed during a period is limited. Therefore, it will take longer time for the larger cycle to adjust the same amount of phase shift. This analysis explains why it is easier for the system to get reentrainment with CikA than that without CikA.

    Instead of light delay in the LD cycle, we also studied the reentrainment where the light is advanced by 6 h. We get the same conclusion that CikA increases the speed of adjusting phase shifts,which increase the system’s adaptation to the disrupted LD cycle.

    3.4. Amplitude response curve

    In addition to period, amplitude also plays an important role in the circadian clock. The degree of circadian rhythm was related to the amplitude of circadian oscillator,which can be response to external stimuli.[25]Researchers experimentally found that higher amplitude of circadian oscillator enhances fitness in cyanobacteriain silico.[26]For mammals, it was found that people with schizophrenia,depression or cancer have no circadian clock or very weak circadian clock(low amplitude). Thus, it is important to explore at what time the external stimulation added helps to increase amplitude and when stimulation inhibits amplitude which should be avoided.In clinical medicine,external stimulation is an important measure to treat related diseases.[27,28]

    Fig.6. (a)The amplitude response curve(ARC)of the circadian clock. (b)–(d)The partial enlarged details of(a)during different intervals of CT.The horizontal axis represents the circadian time when light pulse is added. Negative(positive)amplitude changes correspond to the increase(decrease)of amplitudes after stimuli.The black line indicates the situation without the regulation of CikA,and the red dotted line indicates that with[CikA]=0.6μM.The pulse duration is 1 h.

    Making amplitude response curves(ARCs)is an effective way to study the amplitude changes induced by environmental signals added at different circadian times(CTs).In order to see the difference between the amplitude changes clearly,we draw the ARCs corresponding the situations with and without CikA shown in Fig. 6, where 1 h light pulse is added at different CTs. A complete circadian rhythm cycle includes two parts,the phosphorylation process and the dephosphorylation process,which correspond to the increase and decrease processes of the P-KaiC concentration respectively.The time at which PKaiC reaches its peak level is defined as CT0. The periods of the circadian system are 22.77 h and 21.47 h in the situations with the regulation of CikA and without CikA, respectively.Thus, the interval from CT0 to CT11 corresponds roughly to the dephosphorylation process;and the interval from CT11 to CT22 corresponds roughly to the phosphorylation process.We measured the amplitude of P-KaiC before and after the light pulse. From Fig. 6(a), we can see that the system receiving light stimulation exhibits increase of amplitude at most CTs.However,the system receiving light stimulation at CTs(11 h–13 h) around the trough level of P-KaiC shows significant amplitude reduction (see Fig. 6(c)). This behavior should be avoided because reduction of amplitude decreases the fitness of cyanobacteria. Comparing the black line (without CikA)and the red dotted line(with CikA)in Fig.6(a),the peak and trough of the red dotted curve are slightly larger than those of the black curve. That is to say, with the regulation of CikA,the amplitude increases slightly more and decreases slightly less, which means that CikA slightly increases the fitness of cyanobacteria. Combining Figs. 6(b) and 6(d), we can see that the amplitude increases more when the light stimulation is added during the dephosphorylation process than that when the light stimulation is added during the phosphorylation process.

    4. Conclusion and discussion

    As an input pathway component, CikA plays essential roles in the cyanobacterial circadian clock by regulating the phosphorylation state of KaiC. However, to our knowledge,the current mathematical models of the cyanobacterial circadian clock have been made only considering interactions between the Kai proteins. In the present paper,to investigate the dynamical potential of CikA in regulating the cyanobacterial circadian clock, a more accurate model of introducing CikA into the cyanobacterial circadian clock is built based in part on the model proposed by Takigawaimamuraet al.[7]A detailed study of the model is carried out and the uncovered dynamic mechanism fits the biological observations on the roles of CikA in regulating circadian clock. Therefore, the model captures the main features of CikA regulation, which can be used to study other biological problems associated with CikA as an input component.

    With the model, we get some interesting results about CikA, which affects both the period and amplitude of the cyanobacterial circadian clock. The presence of CikA expands the range of parameters that make the system produce oscillators, which indicates stronger robustness of circadian oscillator. Sensitivity analysis shows that CikA also strengthens the central role of the PTO as a core loop. At the same time, CikA improves the adaptability of the system. In other words,with the regulation of CikA,it is easier for the system to adapt the change of environment. The ARCs tell us that CikA enhances the fitness of cyanobacterial. These findings may help us understand the mechanism of CikA-regulated circadian clock more deeply, which also provides a theoretical clue for biological researchers.

    Besides input functionality of CikA as an input pathway component, it also acts to deipnosophist the output pathway transcription factor RpaA.[12]The technique used is also helpful to analyze the output functionality of CikA. In the near future, with the help of mathematical model and dynamical analysis, we will further explore the dynamical mechanisms of CikA as an output pathway component in the cyanobacterial circadian clock.

    猜你喜歡
    李瑩
    李瑩:為雷達(dá)焊“心臟”的“軍工繡娘”
    李瑩:為雷達(dá)焊“心臟”的“軍工繡娘”
    基于投資者保護(hù)的公平與效率關(guān)系
    智富時代(2019年3期)2019-04-30 11:11:14
    “盈之寶”創(chuàng)業(yè)之道:用藝術(shù)尋找客戶
    大學(xué)生(2018年2期)2018-02-02 00:41:19
    習(xí)近平文藝思想淺析
    和“母親”一起懷孕,房東租客演繹大愛親情
    “小丈夫” 屢屢玩失蹤, 浪漫姐弟戀一地狼煙
    分憂(2016年7期)2016-07-14 02:52:14
    一位公益律師的反家暴情結(jié)
    民主與法制(2016年7期)2016-05-14 10:25:11
    “走教”也精彩
    湖南教育(2016年8期)2016-04-16 03:20:28
    Enoch Emery and His‘Wise Blood’
    综合色av麻豆| 两个人的视频大全免费| 久久久成人免费电影| av在线天堂中文字幕| 97热精品久久久久久| 国产视频内射| 99久久精品一区二区三区| 在线观看免费视频日本深夜| avwww免费| 观看美女的网站| 国产成人aa在线观看| 亚洲av成人精品一区久久| 免费无遮挡裸体视频| av在线观看视频网站免费| 国产黄a三级三级三级人| 高清毛片免费观看视频网站| 九色国产91popny在线| 亚洲色图av天堂| 一级a爱片免费观看的视频| 亚洲内射少妇av| 国产亚洲精品久久久久久毛片| 天美传媒精品一区二区| 色哟哟·www| 亚洲性夜色夜夜综合| 两个人视频免费观看高清| videossex国产| 欧美色视频一区免费| 97超级碰碰碰精品色视频在线观看| 欧美另类亚洲清纯唯美| 又爽又黄a免费视频| 国产探花在线观看一区二区| 小蜜桃在线观看免费完整版高清| 国产精品国产高清国产av| 999久久久精品免费观看国产| 草草在线视频免费看| 69av精品久久久久久| 成年女人毛片免费观看观看9| 欧美日本视频| 精华霜和精华液先用哪个| 校园人妻丝袜中文字幕| 成人二区视频| 老熟妇乱子伦视频在线观看| 国产精品,欧美在线| 高清毛片免费观看视频网站| 精品国内亚洲2022精品成人| 此物有八面人人有两片| 狂野欧美激情性xxxx在线观看| 久99久视频精品免费| 中文字幕久久专区| 婷婷精品国产亚洲av在线| 色综合站精品国产| 一级av片app| 国产精品亚洲美女久久久| 国产精品无大码| 国产伦人伦偷精品视频| 欧美日韩亚洲国产一区二区在线观看| 国产黄片美女视频| 久久精品夜夜夜夜夜久久蜜豆| 一区二区三区四区激情视频 | www日本黄色视频网| 精品一区二区三区人妻视频| 亚洲最大成人手机在线| 乱系列少妇在线播放| 欧美极品一区二区三区四区| 亚洲av成人精品一区久久| 久久精品国产亚洲av涩爱 | 很黄的视频免费| 人妻久久中文字幕网| 啦啦啦观看免费观看视频高清| 国产精品乱码一区二三区的特点| 99久久久亚洲精品蜜臀av| 婷婷色综合大香蕉| 国产国拍精品亚洲av在线观看| 女人十人毛片免费观看3o分钟| 国产淫片久久久久久久久| 蜜桃久久精品国产亚洲av| 亚洲av中文字字幕乱码综合| 亚洲av五月六月丁香网| 日韩欧美 国产精品| 久久亚洲精品不卡| 日韩强制内射视频| 一本久久中文字幕| 女人被狂操c到高潮| 亚洲在线观看片| 老熟妇乱子伦视频在线观看| 亚洲av中文av极速乱 | 老女人水多毛片| 老司机深夜福利视频在线观看| 人妻夜夜爽99麻豆av| 在线观看一区二区三区| 亚洲无线在线观看| 国产精华一区二区三区| 蜜桃亚洲精品一区二区三区| 99精品在免费线老司机午夜| 国产黄片美女视频| 亚洲美女黄片视频| 少妇的逼好多水| 亚洲五月天丁香| 桃色一区二区三区在线观看| 亚洲三级黄色毛片| 亚洲专区国产一区二区| 日本与韩国留学比较| 五月玫瑰六月丁香| 91久久精品国产一区二区成人| 很黄的视频免费| 日韩大尺度精品在线看网址| 精品福利观看| 最好的美女福利视频网| 国产男人的电影天堂91| 日韩高清综合在线| 一卡2卡三卡四卡精品乱码亚洲| av专区在线播放| 变态另类成人亚洲欧美熟女| 色在线成人网| 国产一区二区激情短视频| 在线观看午夜福利视频| 婷婷丁香在线五月| 看黄色毛片网站| av.在线天堂| 久久亚洲精品不卡| 天美传媒精品一区二区| 99热6这里只有精品| 亚洲欧美日韩高清在线视频| 国产激情偷乱视频一区二区| 久久精品国产亚洲av涩爱 | 麻豆国产av国片精品| 欧美日韩亚洲国产一区二区在线观看| 亚洲国产日韩欧美精品在线观看| 两个人的视频大全免费| 国产精品国产三级国产av玫瑰| 狠狠狠狠99中文字幕| 99精品久久久久人妻精品| 尤物成人国产欧美一区二区三区| 天堂影院成人在线观看| 真人做人爱边吃奶动态| 在线天堂最新版资源| eeuss影院久久| 亚洲性夜色夜夜综合| 能在线免费观看的黄片| 亚洲av成人精品一区久久| av天堂中文字幕网| 搡老岳熟女国产| 欧美性感艳星| 在线免费观看不下载黄p国产 | 长腿黑丝高跟| 日韩高清综合在线| 亚洲不卡免费看| 老司机深夜福利视频在线观看| 久久久午夜欧美精品| 99在线人妻在线中文字幕| 亚洲精品久久国产高清桃花| 欧美另类亚洲清纯唯美| 不卡视频在线观看欧美| 国产三级在线视频| 男女边吃奶边做爰视频| 午夜免费激情av| 日韩欧美在线二视频| 亚洲最大成人中文| 村上凉子中文字幕在线| 露出奶头的视频| 亚洲国产色片| 国产色婷婷99| 国产午夜福利久久久久久| 丰满的人妻完整版| 亚洲国产精品成人综合色| 免费看日本二区| 久久久久久久久久成人| 男人舔奶头视频| 少妇的逼好多水| 国产淫片久久久久久久久| 亚洲av成人精品一区久久| 99视频精品全部免费 在线| 成人国产麻豆网| 69人妻影院| 亚洲一区高清亚洲精品| 亚洲性久久影院| 赤兔流量卡办理| 999久久久精品免费观看国产| 免费电影在线观看免费观看| 日本 欧美在线| 色av中文字幕| 天天一区二区日本电影三级| 欧美日本亚洲视频在线播放| 99九九线精品视频在线观看视频| 一进一出抽搐动态| 窝窝影院91人妻| 国内精品美女久久久久久| 91麻豆精品激情在线观看国产| 精品久久久久久久末码| 国产精品久久电影中文字幕| 一个人看视频在线观看www免费| 亚洲欧美日韩无卡精品| 深爱激情五月婷婷| 97超级碰碰碰精品色视频在线观看| 日韩欧美在线乱码| 成人国产一区最新在线观看| 久久精品综合一区二区三区| 国产真实乱freesex| 成人欧美大片| 91午夜精品亚洲一区二区三区 | 国产精品人妻久久久久久| 亚洲 国产 在线| 91久久精品国产一区二区三区| 日韩 亚洲 欧美在线| 色综合婷婷激情| 欧美+亚洲+日韩+国产| 亚洲无线在线观看| 免费在线观看成人毛片| 韩国av一区二区三区四区| 亚洲精品一区av在线观看| 噜噜噜噜噜久久久久久91| 麻豆成人午夜福利视频| a级一级毛片免费在线观看| 国产三级中文精品| 91久久精品国产一区二区三区| 国内精品久久久久久久电影| 看黄色毛片网站| 最好的美女福利视频网| 18禁在线播放成人免费| 男人狂女人下面高潮的视频| 亚洲午夜理论影院| 午夜福利在线在线| 最近在线观看免费完整版| 亚洲国产精品合色在线| av在线天堂中文字幕| 国产 一区 欧美 日韩| 午夜影院日韩av| 一本久久中文字幕| 国产精品1区2区在线观看.| 欧美日韩中文字幕国产精品一区二区三区| 一本精品99久久精品77| 自拍偷自拍亚洲精品老妇| 久久九九热精品免费| 亚洲国产精品合色在线| 琪琪午夜伦伦电影理论片6080| 亚洲第一电影网av| 美女大奶头视频| 在现免费观看毛片| 三级男女做爰猛烈吃奶摸视频| 露出奶头的视频| 国产又黄又爽又无遮挡在线| 国内揄拍国产精品人妻在线| 欧美3d第一页| 日韩精品青青久久久久久| 亚洲五月天丁香| 久久久久免费精品人妻一区二区| 国产成年人精品一区二区| 一级av片app| 久久久久久九九精品二区国产| 免费人成视频x8x8入口观看| 久久精品国产亚洲av涩爱 | 日韩欧美三级三区| 别揉我奶头 嗯啊视频| 人人妻,人人澡人人爽秒播| 国产亚洲精品久久久com| 老熟妇乱子伦视频在线观看| 狠狠狠狠99中文字幕| 午夜视频国产福利| 欧美精品国产亚洲| 老师上课跳d突然被开到最大视频| 国产 一区精品| 成人性生交大片免费视频hd| 国产极品精品免费视频能看的| 亚洲av第一区精品v没综合| 免费人成在线观看视频色| 热99在线观看视频| 精品午夜福利视频在线观看一区| 午夜影院日韩av| 午夜福利视频1000在线观看| av福利片在线观看| 欧美黑人巨大hd| 中文资源天堂在线| 偷拍熟女少妇极品色| 国产av一区在线观看免费| 亚洲精品一区av在线观看| 嫁个100分男人电影在线观看| 伦精品一区二区三区| 国产亚洲91精品色在线| 国产高清三级在线| 香蕉av资源在线| 国产精品久久久久久久久免| 日本与韩国留学比较| 午夜爱爱视频在线播放| 免费在线观看影片大全网站| 一边摸一边抽搐一进一小说| 成人国产综合亚洲| 一区二区三区四区激情视频 | 美女被艹到高潮喷水动态| 99久久精品热视频| 日本撒尿小便嘘嘘汇集6| 国产免费一级a男人的天堂| 午夜老司机福利剧场| 日本熟妇午夜| 一级毛片久久久久久久久女| 欧美潮喷喷水| 十八禁网站免费在线| 国内久久婷婷六月综合欲色啪| 村上凉子中文字幕在线| 夜夜夜夜夜久久久久| 国产爱豆传媒在线观看| 国产欧美日韩一区二区精品| 亚洲第一区二区三区不卡| 久久天躁狠狠躁夜夜2o2o| 搞女人的毛片| 超碰av人人做人人爽久久| 亚洲中文日韩欧美视频| 淫秽高清视频在线观看| 午夜免费男女啪啪视频观看 | bbb黄色大片| 日本黄色视频三级网站网址| 久久久久九九精品影院| 夜夜看夜夜爽夜夜摸| 亚洲五月天丁香| 欧美精品啪啪一区二区三区| 亚洲av熟女| 韩国av一区二区三区四区| 国产精品日韩av在线免费观看| 亚洲av成人精品一区久久| 亚洲综合色惰| av黄色大香蕉| 久久中文看片网| 日本在线视频免费播放| 久久精品国产自在天天线| 两人在一起打扑克的视频| 久久热精品热| 国产中年淑女户外野战色| 91久久精品电影网| 国产精品一区www在线观看 | 99久国产av精品| 午夜福利在线在线| 桃色一区二区三区在线观看| 国产午夜精品久久久久久一区二区三区 | 97人妻精品一区二区三区麻豆| 99国产精品一区二区蜜桃av| 在线播放国产精品三级| 国产精品一区二区性色av| 美女高潮喷水抽搐中文字幕| 在线看三级毛片| 欧美激情久久久久久爽电影| 国产精品免费一区二区三区在线| 午夜激情福利司机影院| 老熟妇乱子伦视频在线观看| 长腿黑丝高跟| 欧美中文日本在线观看视频| 欧美日韩黄片免| 丰满乱子伦码专区| 亚洲第一电影网av| 亚洲精品在线观看二区| 欧美成人a在线观看| 欧美潮喷喷水| 男女下面进入的视频免费午夜| 免费观看的影片在线观看| 免费av观看视频| 一进一出好大好爽视频| 成人av在线播放网站| 国产毛片a区久久久久| 国产成人福利小说| 国内久久婷婷六月综合欲色啪| 日本与韩国留学比较| 午夜福利视频1000在线观看| 赤兔流量卡办理| 自拍偷自拍亚洲精品老妇| videossex国产| 在线国产一区二区在线| avwww免费| 国产麻豆成人av免费视频| 国产精品一区二区三区四区免费观看 | 欧美三级亚洲精品| av在线天堂中文字幕| 老女人水多毛片| 欧美日本亚洲视频在线播放| 一本久久中文字幕| 欧美激情国产日韩精品一区| 国产高清视频在线观看网站| 中文字幕免费在线视频6| 一区福利在线观看| 成人亚洲精品av一区二区| 日韩人妻高清精品专区| 亚洲18禁久久av| 亚洲18禁久久av| 国产亚洲欧美98| 我的老师免费观看完整版| 久久亚洲真实| 在线国产一区二区在线| 国产精品国产三级国产av玫瑰| 蜜桃亚洲精品一区二区三区| 久久热精品热| 男人的好看免费观看在线视频| 亚洲,欧美,日韩| 欧美精品国产亚洲| 亚洲图色成人| 国产午夜精品久久久久久一区二区三区 | 国产一级毛片七仙女欲春2| 久久6这里有精品| 亚洲精华国产精华液的使用体验 | 国产精品98久久久久久宅男小说| 国产精品精品国产色婷婷| 中文字幕久久专区| www.色视频.com| 成年女人毛片免费观看观看9| 久久久午夜欧美精品| 女人十人毛片免费观看3o分钟| 色综合婷婷激情| 五月玫瑰六月丁香| 又爽又黄a免费视频| 又粗又爽又猛毛片免费看| 久久6这里有精品| 亚洲欧美日韩高清在线视频| 男人舔女人下体高潮全视频| 久久午夜亚洲精品久久| 可以在线观看毛片的网站| 在线天堂最新版资源| 免费在线观看日本一区| 国产精品久久久久久亚洲av鲁大| 午夜免费激情av| 国产精品一区二区三区四区免费观看 | 精品久久久久久久人妻蜜臀av| 黄片wwwwww| 狂野欧美白嫩少妇大欣赏| av福利片在线观看| 黄色视频,在线免费观看| 欧美一级a爱片免费观看看| 精品久久久久久久久亚洲 | 精品人妻1区二区| 男人的好看免费观看在线视频| 亚洲最大成人中文| 成年女人永久免费观看视频| 久久久久久久久久久丰满 | 欧美日韩乱码在线| 真人做人爱边吃奶动态| 久久久成人免费电影| 女人被狂操c到高潮| 成年女人看的毛片在线观看| 99九九线精品视频在线观看视频| 成人综合一区亚洲| 亚洲综合色惰| 成年女人毛片免费观看观看9| 老熟妇仑乱视频hdxx| 亚洲va在线va天堂va国产| 免费看光身美女| 亚洲精品国产成人久久av| 黄色丝袜av网址大全| 97人妻精品一区二区三区麻豆| bbb黄色大片| 欧美三级亚洲精品| 中亚洲国语对白在线视频| 熟妇人妻久久中文字幕3abv| 国产成人福利小说| 国产高清三级在线| 他把我摸到了高潮在线观看| 老司机福利观看| 久久99热这里只有精品18| 真实男女啪啪啪动态图| 久久精品综合一区二区三区| 欧美精品啪啪一区二区三区| 美女cb高潮喷水在线观看| 国产乱人视频| 亚洲av电影不卡..在线观看| 亚洲自偷自拍三级| 搡老熟女国产l中国老女人| 狠狠狠狠99中文字幕| 欧美区成人在线视频| 中国美女看黄片| 国产毛片a区久久久久| 国产精品精品国产色婷婷| 亚洲图色成人| 国产三级在线视频| 韩国av在线不卡| 日韩精品青青久久久久久| 国产美女午夜福利| 国产爱豆传媒在线观看| 精品一区二区三区视频在线| www.色视频.com| 成人二区视频| 精品久久久久久,| 99热6这里只有精品| 老熟妇仑乱视频hdxx| 简卡轻食公司| 人妻夜夜爽99麻豆av| 久久午夜亚洲精品久久| 999久久久精品免费观看国产| 男插女下体视频免费在线播放| 欧美日韩瑟瑟在线播放| 一个人免费在线观看电影| netflix在线观看网站| 国产国拍精品亚洲av在线观看| 真实男女啪啪啪动态图| 精品一区二区免费观看| 婷婷丁香在线五月| 国产午夜精品久久久久久一区二区三区 | 精华霜和精华液先用哪个| 麻豆av噜噜一区二区三区| 久久久午夜欧美精品| 三级毛片av免费| 国产精品综合久久久久久久免费| 久久国内精品自在自线图片| 国产老妇女一区| 欧美xxxx黑人xx丫x性爽| 亚洲人成网站在线播放欧美日韩| 少妇人妻精品综合一区二区 | 成人国产综合亚洲| 最近在线观看免费完整版| 日韩欧美精品v在线| 乱码一卡2卡4卡精品| 国国产精品蜜臀av免费| 亚洲色图av天堂| 国产真实乱freesex| 18禁黄网站禁片免费观看直播| 观看免费一级毛片| 一级毛片久久久久久久久女| 亚洲国产高清在线一区二区三| aaaaa片日本免费| 久久精品国产亚洲av香蕉五月| 欧美潮喷喷水| 免费搜索国产男女视频| 日韩人妻高清精品专区| 亚洲美女搞黄在线观看 | 51国产日韩欧美| www日本黄色视频网| 亚洲欧美日韩高清在线视频| 日本a在线网址| 身体一侧抽搐| 女人被狂操c到高潮| 99热只有精品国产| 欧美人与善性xxx| 18禁裸乳无遮挡免费网站照片| 91午夜精品亚洲一区二区三区 | 亚洲熟妇熟女久久| 国产精品一区二区三区四区免费观看 | aaaaa片日本免费| 丰满的人妻完整版| 搡老熟女国产l中国老女人| 18禁裸乳无遮挡免费网站照片| av中文乱码字幕在线| 亚洲天堂国产精品一区在线| 欧美一区二区亚洲| 999久久久精品免费观看国产| 国产精品伦人一区二区| 久久人人爽人人爽人人片va| 国产精品国产高清国产av| 乱系列少妇在线播放| 日韩精品有码人妻一区| 亚洲中文日韩欧美视频| 久久久久国内视频| 村上凉子中文字幕在线| 亚洲欧美日韩高清在线视频| 久久6这里有精品| 22中文网久久字幕| 日韩欧美在线二视频| 亚洲无线在线观看| 久久香蕉精品热| 干丝袜人妻中文字幕| 亚洲精品色激情综合| 国产精品乱码一区二三区的特点| av女优亚洲男人天堂| 一区福利在线观看| 午夜日韩欧美国产| 嫁个100分男人电影在线观看| 国产精品无大码| 黄色欧美视频在线观看| 色播亚洲综合网| 亚洲av中文av极速乱 | 日韩强制内射视频| 在线天堂最新版资源| 一级av片app| 床上黄色一级片| 黄色日韩在线| 很黄的视频免费| 99热只有精品国产| 欧美日韩瑟瑟在线播放| 51国产日韩欧美| 国产精品一及| 国产探花极品一区二区| 乱系列少妇在线播放| 国产精品久久久久久精品电影| 91在线观看av| 桃红色精品国产亚洲av| 五月伊人婷婷丁香| 赤兔流量卡办理| 一本久久中文字幕| 亚洲精品成人久久久久久| 国产激情偷乱视频一区二区| 久久久久国内视频| 91av网一区二区| 看片在线看免费视频| 久久久久国内视频| 97人妻精品一区二区三区麻豆| av专区在线播放| 一区二区三区免费毛片| 色在线成人网| 成人欧美大片| 18禁裸乳无遮挡免费网站照片| 久久国产乱子免费精品| 麻豆一二三区av精品| 亚洲欧美清纯卡通| 美女大奶头视频| 成人精品一区二区免费| 在线观看午夜福利视频| 亚洲精品日韩av片在线观看| 国产激情偷乱视频一区二区| a在线观看视频网站| 国产一区二区三区av在线 | 精品一区二区三区av网在线观看| 99国产精品一区二区蜜桃av| 午夜福利在线在线| 黄色欧美视频在线观看| 久久午夜福利片| 日本色播在线视频| 看黄色毛片网站| 成年人黄色毛片网站| 久久精品国产亚洲av香蕉五月| 国产精品99久久久久久久久| 国产精品一区www在线观看 | 国产精品不卡视频一区二区| 毛片女人毛片|