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

    元胞自動(dòng)機(jī)在模擬細(xì)胞行為方面的應(yīng)用

    2017-06-01 12:20:40李向陽崔甜甜
    關(guān)鍵詞:規(guī)則生物生長

    李向陽 孫 亮 崔甜甜

    (蘇州大學(xué)醫(yī)學(xué)部放射醫(yī)學(xué)與防護(hù)學(xué)院,江蘇 蘇州 215000)

    元胞自動(dòng)機(jī)在模擬細(xì)胞行為方面的應(yīng)用

    李向陽 孫 亮*崔甜甜

    (蘇州大學(xué)醫(yī)學(xué)部放射醫(yī)學(xué)與防護(hù)學(xué)院,江蘇 蘇州 215000)

    在生物細(xì)胞生長機(jī)理以及莫諾德微分方程模型的基礎(chǔ)上,推導(dǎo)出生物細(xì)胞生長以及底物消耗的微分方程,并基于此建立模擬生物細(xì)胞生長的元胞自動(dòng)機(jī)模型(CABCGM)。CABCGM采用二維細(xì)胞自動(dòng)機(jī)作為細(xì)胞的生長空間,采用Moore型鄰域作為細(xì)胞的鄰居規(guī)則,其演化規(guī)則根據(jù)細(xì)胞的生長機(jī)理和推導(dǎo)出的動(dòng)力學(xué)微分方程設(shè)計(jì)。在建立元胞自動(dòng)機(jī)模型的前提下,以Matlab為開發(fā)工具來實(shí)現(xiàn)模型編程,并進(jìn)行仿真實(shí)驗(yàn)。仿真結(jié)果表明,利用元胞自動(dòng)機(jī)建立的生物細(xì)胞的生長模型與實(shí)際微分方程細(xì)胞的生長曲線比較吻合(尤其在前期時(shí)間段),故元胞自動(dòng)機(jī)能較好地描述細(xì)胞生長的演化行為。

    元胞自動(dòng)機(jī)模型;生物細(xì)胞;細(xì)胞行為

    引言

    對(duì)于生物細(xì)胞生長動(dòng)態(tài)的描述模型,最常用的是純數(shù)學(xué)模型。純數(shù)學(xué)模型一般為非線性動(dòng)力學(xué)微分方程[1-2],通常是在理想情況下通過理論推導(dǎo)得出來的,是一類確定性的機(jī)理模型。非線性動(dòng)力學(xué)微分方程的這一特點(diǎn)有些與生物細(xì)胞生長的隨機(jī)性特點(diǎn)相悖,當(dāng)生物非線性系統(tǒng)的復(fù)雜性增加,以及求解微分方程或計(jì)算數(shù)值的過程變得越來越無望時(shí),元胞自動(dòng)機(jī)的優(yōu)點(diǎn)就突顯出來,它是非常靈活有效的建模工具。

    上世紀(jì)80年代末,Christopher Langton提出了人工生命,并指出生命系統(tǒng)的實(shí)現(xiàn)形式有多種,除已知的碳水化合物的形式,還可以利用別的物質(zhì)構(gòu)造出不同于自然生命的另類載體生命形式,賦予它們生命特征,使其具有進(jìn)化、變異等生命現(xiàn)象,得出生命的普遍行為。這類的形式可以是物理的、化學(xué)的以及程序的形式等[3-7]。元胞自動(dòng)機(jī),正是從人工生命衍生出的一個(gè)分支,其單位眾多、特定條件相互作用,具有量變到質(zhì)變的“涌現(xiàn)”現(xiàn)象等特性,是當(dāng)之無愧的復(fù)雜系統(tǒng),是自然生命系統(tǒng)的對(duì)偶系統(tǒng),為復(fù)雜系統(tǒng)建模、仿真和控制提供了一種理想的計(jì)算模型,其表現(xiàn)復(fù)雜系統(tǒng)的能力為細(xì)胞生長過程的建模仿真提供了新思路[8]。細(xì)胞的生長是一個(gè)復(fù)雜的生化過程,就單個(gè)細(xì)胞而言,其生物活性相對(duì)簡單,但其形成的整個(gè)生化過程卻表現(xiàn)出復(fù)雜的行為[9-10]。細(xì)胞群落與環(huán)境相互作用,發(fā)展并豐富自身結(jié)構(gòu)或形態(tài),也即自組織過程。這個(gè)過程機(jī)理復(fù)雜,動(dòng)態(tài)行為多變,其演化具有多樣性、隨意性、不確定性以及強(qiáng)非線性等,說明它符合復(fù)雜系統(tǒng)的特征。目前,研究細(xì)胞生長建立的模型不是以復(fù)雜系統(tǒng)的觀點(diǎn)建立的,所以建立的模型未能表現(xiàn)出復(fù)雜系統(tǒng)的特征,使得這類模型與細(xì)胞的生長實(shí)際存在較大的差距。

    本研究基于生物細(xì)胞生長的機(jī)理和推導(dǎo)出的動(dòng)力學(xué)微分方程,建立了模擬生物細(xì)胞生長的元胞自動(dòng)機(jī)模型(cellular automata model of biology cell growth,CABCGM),并根據(jù)設(shè)計(jì)的Matlab規(guī)則進(jìn)行了編程設(shè)計(jì),并進(jìn)行了仿真實(shí)驗(yàn)。仿真實(shí)驗(yàn)結(jié)果表明,CABCGM能較好地復(fù)現(xiàn)動(dòng)力學(xué)微分方程所描述的生物細(xì)胞生長過程,還可以使演化過程可視化。對(duì)于Matlab模擬生物細(xì)胞的研究鮮有報(bào)道,本次CABCGM的設(shè)計(jì),是用Matlab在生物細(xì)胞領(lǐng)域建模的新探索。

    1 實(shí)驗(yàn)方法

    1.1 細(xì)胞生長及基質(zhì)消耗的動(dòng)力學(xué)微分方程的推導(dǎo)

    現(xiàn)代細(xì)胞動(dòng)力學(xué)的奠基人Monod為代表的眾多科學(xué)工作者對(duì)各種細(xì)胞動(dòng)力學(xué)模型進(jìn)行了研究和建模,其中最為著名的是莫諾德方程[11-13],即

    (1)

    式中:μ為比生長速率;μm為最大比生長速率;Ks為半飽和常數(shù),其值為比生長速率等于最大比生長速率一半時(shí)的基質(zhì)濃度;S為生長的基質(zhì)濃度。

    按照比生長速率的定義,有

    (2)

    聯(lián)立式(1)、(2)便得到

    (3)

    若取足夠小,則

    (4)

    聯(lián)立式(3)、(4)便得到差分方程,即

    (5)

    同理,在基質(zhì)消耗方面,根據(jù)碳守恒原理,S底物=X細(xì)胞+P產(chǎn)物+m維持,即

    (6)

    由于本次模擬不考慮產(chǎn)物的生成,故不考慮P產(chǎn)物項(xiàng),得

    (7)

    同理,若取得足夠小,即

    (8)

    聯(lián)立式(3)、(7)、(8),得到關(guān)于基質(zhì)消耗的差分方程如下:

    (9)

    以上為細(xì)胞生長及基質(zhì)消耗的動(dòng)力學(xué)微分方程和差分方程的推導(dǎo)過程,可進(jìn)一步對(duì)其中的差分方程進(jìn)行理論特性統(tǒng)計(jì)分析,得到CABCGM與離散的生物細(xì)胞生長動(dòng)力學(xué)模型在數(shù)學(xué)期望意義下是一致的結(jié)論,這對(duì)本研究構(gòu)造的CABCGM具有指導(dǎo)意義?;诖耍€可以假設(shè)細(xì)胞生長過程中的濃度最大值為Xmax,細(xì)胞濃度在T0時(shí)間內(nèi)的最大變化量ΔXmax(t),根據(jù)細(xì)胞生長機(jī)理和實(shí)驗(yàn)實(shí)例,設(shè)置參數(shù)后,Xmax和ΔXmax(t)可由式(3)和式(5)求出,為后面演化概率的設(shè)計(jì)做準(zhǔn)備。

    1.2 CABCGM結(jié)構(gòu)設(shè)計(jì)

    本次模型采用二維元胞自動(dòng)機(jī)作為生物細(xì)胞生長模型,元胞自動(dòng)機(jī)中的每一個(gè)元胞代表一個(gè)生物細(xì)胞,每一個(gè)元胞具有不同的狀態(tài)。在元胞自動(dòng)機(jī)中,生物細(xì)胞的生長演化過程就是元胞自動(dòng)機(jī)各個(gè)元胞在不同狀態(tài)值之間不斷轉(zhuǎn)化的過程?;贑A的生物細(xì)胞生長模型被定義為一個(gè)五元組,有

    (10)

    式中:t為離散時(shí)間;Cells為CABCGM的基本元素,即元胞;Cellspace為細(xì)胞空間;Neighborhoods為CABCGM的鄰域,即演化規(guī)則的定義域;Rules為CABCGM的演化規(guī)則,規(guī)定著元胞狀態(tài)的演化,是CABCGM的核心。

    1.2.1 CABCGM的離散時(shí)間

    對(duì)生物細(xì)胞的整個(gè)生長過程進(jìn)行離散化,對(duì)于每一個(gè)離散時(shí)刻t=kT0,其中k={0,1,2,…,N}為離散時(shí)間序列,T0為離散時(shí)間間隔。

    1.2.2 CABCGM細(xì)胞和細(xì)胞空間

    元胞是CABCGM的基本元素,假設(shè)CABCGM的元胞是一個(gè)微小的正方形。本次模擬的是在理想情況下在直徑為d=10 cm的平皿上生長的細(xì)胞,最終長出5×106個(gè)細(xì)胞,據(jù)此把每個(gè)細(xì)胞由圓形折合成正方形,得到單個(gè)細(xì)胞的邊長為40 μm。整個(gè)細(xì)胞空間Cellspace={Cij|i,j∈{0,1,2,…,N}},其中N×N為元胞空間的元胞總數(shù)。

    1.2.3 CABCGM細(xì)胞鄰域

    CABCGM的元胞鄰域是其演化規(guī)則的演化作用范圍,也就是元胞演化規(guī)則的定義域。在此采用Moore型領(lǐng)域,定義如下:

    (11)

    1.2.4 CABCGM細(xì)胞狀態(tài)

    CA中的每一個(gè)元胞代表一個(gè)細(xì)胞,每個(gè)細(xì)胞具有不同的狀態(tài),令CA細(xì)胞Cell(i,j)在t時(shí)刻的狀態(tài)為Sij(t),取兩種狀態(tài)值(0和1)。其中,Sij(t)=0表示Cell(i,j)在t時(shí)刻為空或細(xì)胞死亡,Sij(t)=1表示新生細(xì)胞占據(jù)了該空格。如此,在CA模型中,細(xì)胞的生長過程就是CA各個(gè)細(xì)胞單元的狀態(tài)值在0、1之間不斷轉(zhuǎn)化的過程。

    1.2.5 CABCGM的演化規(guī)則

    CA的演化規(guī)則具有局部特性,CA的每個(gè)細(xì)胞演化均依賴于與之相鄰的細(xì)胞狀態(tài)。也就是說,細(xì)胞下一時(shí)刻的狀態(tài)只與它鄰域的細(xì)胞狀態(tài)有關(guān)。CA按照演化規(guī)則進(jìn)行并行演化,即在任意一個(gè)離散時(shí)刻,細(xì)胞空間中所有細(xì)胞同時(shí)依據(jù)規(guī)則演化。

    根據(jù)生物細(xì)胞生長機(jī)理以及推導(dǎo)出來的動(dòng)力學(xué)微分方程,設(shè)計(jì)了CABCGM的演化規(guī)則如下:

    1) IfSij(t)=0,then

    (12)

    式中,N(i,j,t)是t時(shí)刻CA元胞Cell(i,j)的Moore領(lǐng)域內(nèi)細(xì)胞的個(gè)數(shù),P0,1(i,j,t)為t時(shí)刻狀態(tài)為0的元胞在下一個(gè)時(shí)刻演化為狀態(tài)1的演化概率,P0,1,Thr(i,j,t)是t時(shí)刻狀態(tài)值為0的細(xì)胞下一個(gè)時(shí)刻演化為1的概率閾值。

    演化概率定義如下:

    (13)

    式中,ΔX(t)是CA模型中在t時(shí)刻的理想擴(kuò)增量,ΔXmax(t)是模型的最大擴(kuò)增量,L0,1是狀態(tài)0元胞向狀態(tài)1元胞演化的趨勢(shì)系數(shù),

    P0,1,Thr(i,j,t)、L0,1根據(jù)具體培養(yǎng)條件和機(jī)理確定,并在仿真實(shí)驗(yàn)中進(jìn)行必要調(diào)整。

    2) IfSij(t)=1,then

    (14)

    1.3 仿真實(shí)驗(yàn)

    利用Matlab仿真軟件,編制CA的演化程序進(jìn)行仿真實(shí)驗(yàn)。CA的演化程序包括下列功能:

    1)利用推導(dǎo)公式須計(jì)算Xmax、ΔXmax,此外還需計(jì)算每一kT0時(shí)刻的細(xì)胞濃度、kT0→(k+1)T0時(shí)間段內(nèi)的細(xì)胞濃度變化量ΔXmax(t);

    2)在二維狀態(tài)空間上,自動(dòng)描繪CA的元胞狀態(tài)演化圖案;

    3)在每一個(gè)離散時(shí)刻kT0對(duì)CA元胞進(jìn)行質(zhì)量統(tǒng)計(jì),并換算成CA的細(xì)胞濃度,描繪出各狀態(tài)元胞隨時(shí)間變化的曲線;

    4)對(duì)每一kT0離散時(shí)刻CA的元胞進(jìn)行質(zhì)量統(tǒng)計(jì),并換算成CA的細(xì)胞濃度,描繪其隨時(shí)間變化的曲線。

    CABCGM的活細(xì)胞質(zhì)量統(tǒng)計(jì)方法是:將狀態(tài)值1的細(xì)胞質(zhì)量均視為1,因此當(dāng)狀態(tài)值為1的元胞數(shù)量為N1時(shí),那么狀態(tài)值1的元胞質(zhì)量為N1,設(shè)CA的元胞質(zhì)量換算成濃度的換算系數(shù)為α,則CA的細(xì)胞濃度C=αN1,換算系數(shù)α可在實(shí)驗(yàn)仿真時(shí)進(jìn)行必要的調(diào)整。

    2 實(shí)驗(yàn)結(jié)果

    圖1 CABCGM模型演化過程中細(xì)胞濃度隨時(shí)間變化的曲線以及基質(zhì)濃度隨時(shí)間變化的曲線Fig.1 The cell concentration change over time in the evolution process of CABCGM model and the curve of the substrate concentration change over time

    由圖1可以看出,在開始的一段時(shí)間內(nèi),模擬的圖像與微分方程描述的曲線比較吻合,說明設(shè)計(jì)的元胞自動(dòng)機(jī)在規(guī)則和參數(shù)方面還較吻合,但隨著時(shí)間的進(jìn)行,圖像出現(xiàn)較大的分離,尤其是基質(zhì)消耗的曲線,說明當(dāng)細(xì)胞增至一定數(shù)量后,基質(zhì)瞬間被大量消耗,之前設(shè)計(jì)的元胞自動(dòng)機(jī)的規(guī)則和參數(shù)不再適應(yīng)。所以確切說,短期模擬細(xì)胞的生長,元胞自動(dòng)機(jī)較微分方程表現(xiàn)出諸多的優(yōu)勢(shì),但細(xì)胞數(shù)量到達(dá)一定值后,還需要調(diào)整元胞自動(dòng)機(jī)參數(shù)等,而微分方程就不需要調(diào)整。

    3 討論和結(jié)論

    本研究立足于生物細(xì)胞的生長機(jī)理,在莫諾德方程的基礎(chǔ)上,推導(dǎo)出描述生物細(xì)胞生長的非線性動(dòng)力學(xué)微分方程,通過考慮實(shí)際,對(duì)此微分方程設(shè)計(jì)參數(shù),建立了模擬生物細(xì)胞生長演化的元胞自動(dòng)機(jī)模型(CABCGM),通過對(duì)元胞自動(dòng)機(jī)進(jìn)行規(guī)則的設(shè)計(jì),該模型的仿真結(jié)果能較好地復(fù)現(xiàn)動(dòng)力學(xué)微分方程所描述的生物細(xì)胞的生長演化過程,且更能表現(xiàn)出復(fù)雜系統(tǒng)的復(fù)雜性(多樣性、隨機(jī)性、魯棒性、不確定性等)。本實(shí)驗(yàn)的新穎之處是:在莫諾德基礎(chǔ)上,推導(dǎo)出了細(xì)胞的生長曲線微分方程,這與其他同行利用現(xiàn)成的微分方程離散化后再模擬是不同的。此外,本次模擬討論點(diǎn)在基質(zhì)瞬間被大量消耗的細(xì)胞值,確定了這個(gè)點(diǎn)后,可以讓元胞自動(dòng)機(jī)的模擬過程分段進(jìn)行,該點(diǎn)之前的按照以前的規(guī)則和參數(shù)運(yùn)行,該點(diǎn)之后的根據(jù)實(shí)際設(shè)計(jì)的新規(guī)則及參數(shù)進(jìn)行。 本課題的下一個(gè)研究重點(diǎn)是:利用CA模型的優(yōu)點(diǎn),更具體、更深入地模擬生物細(xì)胞的生長演變,如具體到細(xì)胞的各個(gè)周期、細(xì)胞的成熟、成熟后的分裂分化等。

    [1] Bajpai RK, Reu BM. Evalution of feeding strategies in carbon-regulated secondary metabolite production through mathematical modeling. Biotechnology and Bioengineering. 1981, 23(4): 717-738.

    [2] Menezes J. Mathematical modeling of penicillin-G fed-batch fermentation[J]. Journal of Chemical Technology and Biotechnology.1994,61(2):123-138.

    [3] 李廣運(yùn),韓力群. 人工生命及其基本原理[J]. 計(jì)算機(jī)仿真,2004(5):121-125.

    [4] 喻海飛,汪定偉. 人工生命研究綜述[J]. 信息與控制,2004(4):434-439.

    [5] 盧欣華,孫吉貴,韓霄松,等. 人工生命:計(jì)算機(jī)與生物學(xué)交叉的前沿學(xué)科[J]. 計(jì)算機(jī)科學(xué),2008(5):9-13.

    [6] 黃文高,潘志庚. 人工生命在計(jì)算機(jī)圖形學(xué)中的應(yīng)用[J]. 計(jì)算機(jī)輔助設(shè)計(jì)與圖形學(xué)學(xué)報(bào),2005(7):1383-1388.

    [7] 張曉榮. 計(jì)算主義:從Cyborg走向人工生命[J]. 青海社會(huì)科學(xué),2007(1):133-135.

    [8] Ray TS. An evolutionary approach to synthetic biology: Zen and the art of creating life. Artificial Life Journal, 1994,1(1/2):179-209.

    [9] 于乃功,阮曉鋼. 模擬青霉素發(fā)酵過程中菌體生長動(dòng)態(tài)的細(xì)胞自動(dòng)機(jī)模型[J]. 生物物理學(xué)報(bào),2004(2):155-162.

    [10] 于乃功,阮曉鋼. 細(xì)胞自動(dòng)機(jī)及其在菌體生長建模仿真中的應(yīng)用[J]. 系統(tǒng)仿真學(xué)報(bào),2004(12):2651-2654.

    [11] 孔維寶,董妙音,李萬武,等. 從莫諾方程談起——生物類本科專業(yè)教學(xué)中3個(gè)相似方程的探討[J]. 微生物學(xué)通報(bào),2015(8):1599-1602.

    [12] 鄭學(xué)翔. 第二講 微生物生長動(dòng)力學(xué)[J]. 釀酒科技,1982(4):22-27.

    [13] 吳佩琮. 第四講 微生物生長代謝過程數(shù)學(xué)模型的建立[J]. 食品與發(fā)酵工業(yè),1982(4):73-81.

    [14] 王鵬,楊明,靳禎. 細(xì)胞自動(dòng)機(jī)及MATLAB仿真[J]. 中北大學(xué)學(xué)報(bào)(自然科學(xué)版),2006,27(5):388-391.

    Applications of Cellular Automata to Simulate Cell Behavior

    Li Xiangyang Sun Liang*Cui Tiantian

    (Schoolofradiationmedicineandprotection,MedicalDepartment,Suzhouuniversity,Jiangsu,Suzhou215000)

    cellular automata model;biological cells;cell behavior

    10.3969/j.issn.0258-8021. 2017. 02.018

    2016-01-18, 錄用日期:2016-09-30

    國家自然科學(xué)基金(K112805915)

    R318

    D

    0258-8021(2017) 02-0248-04

    *通信作者(Corresponding author), E-mail: slhmz@126.com

    猜你喜歡
    規(guī)則生物生長
    生物多樣性
    生物多樣性
    上上生物
    撐竿跳規(guī)則的制定
    數(shù)獨(dú)的規(guī)則和演變
    碗蓮生長記
    小讀者(2021年2期)2021-03-29 05:03:48
    第12話 完美生物
    航空世界(2020年10期)2020-01-19 14:36:20
    生長在哪里的啟示
    生長
    文苑(2018年22期)2018-11-19 02:54:14
    讓規(guī)則不規(guī)則
    Coco薇(2017年11期)2018-01-03 20:59:57
    美女午夜性视频免费| 亚洲中文日韩欧美视频| av超薄肉色丝袜交足视频| 一a级毛片在线观看| 午夜两性在线视频| 国产亚洲精品综合一区在线观看 | 黄色成人免费大全| 国产区一区二久久| 免费观看人在逋| 日韩欧美一区视频在线观看| 777久久人妻少妇嫩草av网站| 久久精品aⅴ一区二区三区四区| 美女福利国产在线| aaaaa片日本免费| 亚洲精品中文字幕一二三四区| 欧美黑人精品巨大| 夜夜看夜夜爽夜夜摸 | 又黄又粗又硬又大视频| 久久国产精品人妻蜜桃| 亚洲欧美激情综合另类| 久久性视频一级片| 免费高清在线观看日韩| 国产精品国产av在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 免费在线观看影片大全网站| 超碰成人久久| 怎么达到女性高潮| 成人三级做爰电影| 久99久视频精品免费| 国产精品九九99| 亚洲美女黄片视频| av电影中文网址| 午夜成年电影在线免费观看| 欧美人与性动交α欧美精品济南到| 久久久久久久午夜电影 | 91麻豆精品激情在线观看国产 | 久久亚洲真实| 日本黄色日本黄色录像| 天天躁狠狠躁夜夜躁狠狠躁| 久久欧美精品欧美久久欧美| 色老头精品视频在线观看| 午夜91福利影院| 免费观看精品视频网站| e午夜精品久久久久久久| 母亲3免费完整高清在线观看| 久久人妻福利社区极品人妻图片| www.999成人在线观看| 亚洲久久久国产精品| 成人三级黄色视频| 亚洲成人久久性| 久久久久国产精品人妻aⅴ院| 精品高清国产在线一区| 国产在线精品亚洲第一网站| 亚洲avbb在线观看| 国产亚洲欧美精品永久| 一级,二级,三级黄色视频| 嫩草影院精品99| 亚洲激情在线av| 亚洲激情在线av| 亚洲av成人av| 欧美日韩乱码在线| 久久99一区二区三区| 欧美黄色片欧美黄色片| 久久人人精品亚洲av| 久久精品91蜜桃| 老司机亚洲免费影院| 久久国产亚洲av麻豆专区| 日韩欧美在线二视频| 极品人妻少妇av视频| 欧美日韩乱码在线| 欧美日本亚洲视频在线播放| 成人特级黄色片久久久久久久| 中文字幕精品免费在线观看视频| 怎么达到女性高潮| 精品第一国产精品| av网站在线播放免费| 午夜老司机福利片| 69精品国产乱码久久久| 欧美激情高清一区二区三区| 日韩视频一区二区在线观看| 久久久久久大精品| 亚洲熟妇熟女久久| 国产成人精品在线电影| 亚洲中文av在线| 一进一出好大好爽视频| 国产又爽黄色视频| 日韩免费av在线播放| 激情在线观看视频在线高清| 国产成人av激情在线播放| 免费久久久久久久精品成人欧美视频| 99久久国产精品久久久| 欧美日韩福利视频一区二区| 色综合站精品国产| 国产高清激情床上av| 伊人久久大香线蕉亚洲五| 一区二区三区激情视频| 色老头精品视频在线观看| 日日夜夜操网爽| 国内久久婷婷六月综合欲色啪| 日韩一卡2卡3卡4卡2021年| 青草久久国产| 欧美日韩一级在线毛片| 国产精品偷伦视频观看了| 99香蕉大伊视频| 午夜老司机福利片| www.www免费av| 日韩人妻精品一区2区三区| 天天添夜夜摸| 黑丝袜美女国产一区| 黄色视频不卡| 亚洲熟妇中文字幕五十中出 | 亚洲黑人精品在线| 黄频高清免费视频| 国产精品av久久久久免费| 天堂中文最新版在线下载| 亚洲狠狠婷婷综合久久图片| 欧洲精品卡2卡3卡4卡5卡区| 一级片'在线观看视频| 女人被狂操c到高潮| 中出人妻视频一区二区| 欧美黄色淫秽网站| 欧美黄色淫秽网站| 一级,二级,三级黄色视频| а√天堂www在线а√下载| 精品国产乱码久久久久久男人| 国产有黄有色有爽视频| 国产黄a三级三级三级人| 欧美午夜高清在线| 国产在线观看jvid| 巨乳人妻的诱惑在线观看| 丰满的人妻完整版| 少妇被粗大的猛进出69影院| a级片在线免费高清观看视频| 午夜福利,免费看| 精品乱码久久久久久99久播| 在线免费观看的www视频| 宅男免费午夜| 久久天堂一区二区三区四区| 免费在线观看完整版高清| 亚洲精品一二三| 亚洲视频免费观看视频| 日本 av在线| 亚洲激情在线av| 日韩欧美三级三区| 宅男免费午夜| 久热这里只有精品99| 一a级毛片在线观看| 国产精品香港三级国产av潘金莲| 丁香欧美五月| 午夜免费观看网址| 日韩欧美一区视频在线观看| 亚洲国产毛片av蜜桃av| 欧美激情 高清一区二区三区| 老汉色av国产亚洲站长工具| 国产精品一区二区三区四区久久 | 黑丝袜美女国产一区| 亚洲九九香蕉| 母亲3免费完整高清在线观看| 久久午夜亚洲精品久久| 日本wwww免费看| 欧美另类亚洲清纯唯美| 高清在线国产一区| 久久久久精品国产欧美久久久| 久久国产精品人妻蜜桃| 高潮久久久久久久久久久不卡| 欧美不卡视频在线免费观看 | 亚洲国产欧美日韩在线播放| 成人三级黄色视频| 一进一出好大好爽视频| 91麻豆精品激情在线观看国产 | 狂野欧美激情性xxxx| 色综合欧美亚洲国产小说| 精品第一国产精品| 久久精品国产综合久久久| 国产精品久久电影中文字幕| 少妇的丰满在线观看| 后天国语完整版免费观看| 热99re8久久精品国产| 97超级碰碰碰精品色视频在线观看| 97超级碰碰碰精品色视频在线观看| 亚洲成a人片在线一区二区| 精品午夜福利视频在线观看一区| 人妻久久中文字幕网| 法律面前人人平等表现在哪些方面| 天堂影院成人在线观看| 亚洲欧美日韩高清在线视频| 一级毛片高清免费大全| 咕卡用的链子| 欧美日韩av久久| 激情视频va一区二区三区| 欧美一区二区精品小视频在线| 97碰自拍视频| 国产野战对白在线观看| 国产一区在线观看成人免费| 国产欧美日韩一区二区三| 亚洲熟妇中文字幕五十中出 | 欧美中文日本在线观看视频| 国产精品一区二区在线不卡| 99国产精品一区二区三区| 亚洲精品国产区一区二| 精品电影一区二区在线| 亚洲精品久久成人aⅴ小说| 一二三四社区在线视频社区8| 琪琪午夜伦伦电影理论片6080| 欧美日韩中文字幕国产精品一区二区三区 | 免费在线观看完整版高清| 成人免费观看视频高清| a级毛片在线看网站| 午夜免费激情av| 精品无人区乱码1区二区| 黄色丝袜av网址大全| 亚洲色图综合在线观看| 可以在线观看毛片的网站| 热99国产精品久久久久久7| 亚洲av五月六月丁香网| 757午夜福利合集在线观看| 少妇 在线观看| 99久久久亚洲精品蜜臀av| 亚洲国产毛片av蜜桃av| 国产主播在线观看一区二区| 亚洲精品美女久久久久99蜜臀| 十八禁网站免费在线| 丝袜美腿诱惑在线| 中文欧美无线码| 国产单亲对白刺激| 我的亚洲天堂| 亚洲五月色婷婷综合| 80岁老熟妇乱子伦牲交| 一夜夜www| 午夜精品久久久久久毛片777| 男女下面进入的视频免费午夜 | 免费av毛片视频| 久久香蕉国产精品| 久久精品成人免费网站| 国产极品粉嫩免费观看在线| www国产在线视频色| 黄片播放在线免费| 女人爽到高潮嗷嗷叫在线视频| 亚洲色图av天堂| 999久久久精品免费观看国产| 神马国产精品三级电影在线观看 | 国产99白浆流出| 涩涩av久久男人的天堂| 婷婷六月久久综合丁香| 十八禁网站免费在线| 51午夜福利影视在线观看| 午夜精品在线福利| 正在播放国产对白刺激| 97人妻天天添夜夜摸| 亚洲自拍偷在线| 午夜免费观看网址| 久久久精品欧美日韩精品| 免费看十八禁软件| 老司机午夜十八禁免费视频| 国产片内射在线| 日韩欧美在线乱码| 全区人妻精品视频| 日韩成人在线观看一区二区三区| 国产高清视频在线观看网站| 精品免费久久久久久久清纯| av女优亚洲男人天堂| 91久久精品国产一区二区成人| 简卡轻食公司| 级片在线观看| 亚洲国产日韩欧美精品在线观看| 在线观看免费视频日本深夜| 久久人人精品亚洲av| 亚洲av成人av| 女人十人毛片免费观看3o分钟| 国产欧美日韩一区二区精品| 色av中文字幕| 免费av观看视频| 久久精品91蜜桃| 五月伊人婷婷丁香| 国产一区二区激情短视频| 国产精品久久视频播放| 久久99热6这里只有精品| 又黄又爽又刺激的免费视频.| 成年版毛片免费区| av在线观看视频网站免费| 男女床上黄色一级片免费看| 亚洲黑人精品在线| 亚洲人与动物交配视频| 黄色丝袜av网址大全| 亚洲国产高清在线一区二区三| 搡女人真爽免费视频火全软件 | 三级毛片av免费| 夜夜躁狠狠躁天天躁| 国语自产精品视频在线第100页| 亚洲 国产 在线| 黄片小视频在线播放| 一边摸一边抽搐一进一小说| 性欧美人与动物交配| 怎么达到女性高潮| 变态另类成人亚洲欧美熟女| 中文字幕人成人乱码亚洲影| 90打野战视频偷拍视频| 美女 人体艺术 gogo| 成人欧美大片| 亚洲成人免费电影在线观看| 亚洲片人在线观看| 国产精品日韩av在线免费观看| 精品欧美国产一区二区三| 啦啦啦观看免费观看视频高清| 天天一区二区日本电影三级| 最好的美女福利视频网| 美女高潮喷水抽搐中文字幕| 日韩高清综合在线| 国产在线精品亚洲第一网站| 一区二区三区激情视频| 丰满的人妻完整版| www.www免费av| 亚洲第一电影网av| 婷婷亚洲欧美| 少妇丰满av| 天美传媒精品一区二区| 久久久久性生活片| 欧美乱色亚洲激情| 欧美+亚洲+日韩+国产| 国产黄a三级三级三级人| a在线观看视频网站| 丰满人妻熟妇乱又伦精品不卡| 黄色一级大片看看| 精品人妻视频免费看| 国产精品精品国产色婷婷| 中文字幕精品亚洲无线码一区| 亚洲三级黄色毛片| 日韩欧美 国产精品| 免费高清视频大片| 中文字幕精品亚洲无线码一区| 国产在线精品亚洲第一网站| 成人特级黄色片久久久久久久| 91麻豆精品激情在线观看国产| 美女cb高潮喷水在线观看| 中文字幕人妻熟人妻熟丝袜美| 在线观看66精品国产| 久久精品久久久久久噜噜老黄 | 精品国内亚洲2022精品成人| 精品国产亚洲在线| 久久精品国产亚洲av香蕉五月| 我要搜黄色片| 最新在线观看一区二区三区| 国产白丝娇喘喷水9色精品| 国产精品免费一区二区三区在线| 精品欧美国产一区二区三| 国产精品电影一区二区三区| 亚洲中文字幕一区二区三区有码在线看| 色吧在线观看| 欧美成狂野欧美在线观看| 最近最新中文字幕大全电影3| 亚洲国产高清在线一区二区三| 有码 亚洲区| 国内精品久久久久久久电影| 一本综合久久免费| av国产免费在线观看| 99国产精品一区二区蜜桃av| 校园春色视频在线观看| 黄色日韩在线| 久久精品国产亚洲av香蕉五月| 丝袜美腿在线中文| 色综合欧美亚洲国产小说| 欧美不卡视频在线免费观看| 黄色女人牲交| 欧美色视频一区免费| av专区在线播放| 男人舔奶头视频| 亚洲七黄色美女视频| 成熟少妇高潮喷水视频| 国产精品98久久久久久宅男小说| 美女高潮喷水抽搐中文字幕| 国产精品伦人一区二区| 亚洲专区国产一区二区| 欧美精品国产亚洲| 少妇人妻精品综合一区二区 | 啪啪无遮挡十八禁网站| 内地一区二区视频在线| 三级国产精品欧美在线观看| av在线天堂中文字幕| 国产大屁股一区二区在线视频| 国产白丝娇喘喷水9色精品| 99久久精品热视频| 最后的刺客免费高清国语| 51午夜福利影视在线观看| 色视频www国产| 国产精品乱码一区二三区的特点| 中文字幕人成人乱码亚洲影| 久久6这里有精品| 欧美极品一区二区三区四区| 伦理电影大哥的女人| 一进一出好大好爽视频| 久久性视频一级片| 可以在线观看毛片的网站| 嫩草影视91久久| 一级毛片久久久久久久久女| 91麻豆av在线| av福利片在线观看| 国产私拍福利视频在线观看| 成年人黄色毛片网站| 美女被艹到高潮喷水动态| 欧美又色又爽又黄视频| 深爱激情五月婷婷| 99久国产av精品| 国产乱人伦免费视频| 亚洲激情在线av| 日韩亚洲欧美综合| 小蜜桃在线观看免费完整版高清| 婷婷精品国产亚洲av在线| 亚洲av一区综合| 老司机福利观看| 在线观看66精品国产| 日本精品一区二区三区蜜桃| 简卡轻食公司| 韩国av一区二区三区四区| av国产免费在线观看| 嫩草影院入口| 精品久久久久久久久久免费视频| 成人国产一区最新在线观看| 精品免费久久久久久久清纯| 亚洲人成网站在线播| 在线天堂最新版资源| 久久久久久国产a免费观看| 蜜桃亚洲精品一区二区三区| 全区人妻精品视频| 舔av片在线| 久久久久精品国产欧美久久久| 老熟妇仑乱视频hdxx| 国产伦人伦偷精品视频| 高清在线国产一区| ponron亚洲| 91在线精品国自产拍蜜月| 久久国产精品人妻蜜桃| 日韩精品青青久久久久久| 两个人视频免费观看高清| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 99久久久亚洲精品蜜臀av| 欧美一区二区亚洲| 别揉我奶头 嗯啊视频| 九九在线视频观看精品| 国产一区二区亚洲精品在线观看| 精品日产1卡2卡| 亚洲 欧美 日韩 在线 免费| or卡值多少钱| 蜜桃久久精品国产亚洲av| 国产精品久久久久久久电影| 色综合站精品国产| 午夜福利成人在线免费观看| 亚洲国产高清在线一区二区三| 色综合站精品国产| 校园春色视频在线观看| 91麻豆av在线| 日韩国内少妇激情av| 久久久久久久久大av| 欧美xxxx性猛交bbbb| 成人国产综合亚洲| 国产av不卡久久| 欧美成人免费av一区二区三区| 国产精品久久久久久久电影| or卡值多少钱| 又爽又黄a免费视频| 亚洲av成人不卡在线观看播放网| 日韩av在线大香蕉| 成人一区二区视频在线观看| 色噜噜av男人的天堂激情| 少妇人妻精品综合一区二区 | 在线免费观看的www视频| 久久亚洲真实| 精品午夜福利视频在线观看一区| 一本综合久久免费| 午夜福利18| a级毛片a级免费在线| 简卡轻食公司| 精品久久久久久久久久免费视频| 在线a可以看的网站| 久久精品久久久久久噜噜老黄 | 五月伊人婷婷丁香| 岛国在线免费视频观看| 桃色一区二区三区在线观看| 午夜免费男女啪啪视频观看 | 人人妻,人人澡人人爽秒播| 变态另类丝袜制服| 男插女下体视频免费在线播放| 日本在线视频免费播放| 网址你懂的国产日韩在线| 午夜福利欧美成人| 51午夜福利影视在线观看| 婷婷丁香在线五月| 国产老妇女一区| 观看美女的网站| 狂野欧美白嫩少妇大欣赏| 亚洲专区国产一区二区| 我要看日韩黄色一级片| 十八禁国产超污无遮挡网站| 一级a爱片免费观看的视频| 久久久国产成人免费| 激情在线观看视频在线高清| 淫秽高清视频在线观看| 久久久久精品国产欧美久久久| 亚洲美女黄片视频| 亚洲中文字幕一区二区三区有码在线看| 成年女人永久免费观看视频| 蜜桃久久精品国产亚洲av| 好男人电影高清在线观看| 久久久久久久精品吃奶| 可以在线观看毛片的网站| 亚洲综合色惰| 成年免费大片在线观看| 久久香蕉精品热| 嫩草影院新地址| 免费人成在线观看视频色| 久久精品国产自在天天线| bbb黄色大片| 一个人看的www免费观看视频| 一个人免费在线观看的高清视频| 色尼玛亚洲综合影院| 免费看a级黄色片| 可以在线观看的亚洲视频| 欧美三级亚洲精品| 免费人成视频x8x8入口观看| 久久久久免费精品人妻一区二区| 天堂√8在线中文| 国产av一区在线观看免费| 亚洲成人精品中文字幕电影| 亚洲成人久久爱视频| 在线天堂最新版资源| 亚洲精品影视一区二区三区av| 亚洲无线观看免费| 国产成人福利小说| 乱码一卡2卡4卡精品| 少妇丰满av| 女人被狂操c到高潮| 99国产极品粉嫩在线观看| 99在线视频只有这里精品首页| 又紧又爽又黄一区二区| 欧美性猛交╳xxx乱大交人| 亚洲av成人精品一区久久| 久久精品久久久久久噜噜老黄 | 两性午夜刺激爽爽歪歪视频在线观看| 亚洲无线在线观看| 中文字幕av在线有码专区| 精品乱码久久久久久99久播| 一卡2卡三卡四卡精品乱码亚洲| 精品国产亚洲在线| 亚洲无线在线观看| 久久久久精品国产欧美久久久| av在线蜜桃| 欧美性猛交╳xxx乱大交人| 蜜桃亚洲精品一区二区三区| 两个人的视频大全免费| 极品教师在线免费播放| 级片在线观看| 波多野结衣高清无吗| 岛国在线免费视频观看| 国产大屁股一区二区在线视频| 色哟哟哟哟哟哟| 黄色视频,在线免费观看| aaaaa片日本免费| 变态另类成人亚洲欧美熟女| 亚洲精品粉嫩美女一区| 久久草成人影院| 全区人妻精品视频| 极品教师在线免费播放| 99热这里只有是精品50| 国产不卡一卡二| 欧美国产日韩亚洲一区| 久久精品国产清高在天天线| 色尼玛亚洲综合影院| 久久亚洲精品不卡| 人妻久久中文字幕网| 婷婷色综合大香蕉| 欧美成人免费av一区二区三区| 亚洲成人中文字幕在线播放| aaaaa片日本免费| 国产精品久久久久久亚洲av鲁大| 能在线免费观看的黄片| 两个人的视频大全免费| 欧美色视频一区免费| 国产免费一级a男人的天堂| 欧美成狂野欧美在线观看| 熟女电影av网| 日韩高清综合在线| 精品不卡国产一区二区三区| 97热精品久久久久久| 麻豆成人av在线观看| 99久久无色码亚洲精品果冻| 婷婷亚洲欧美| 亚洲欧美日韩卡通动漫| 成年版毛片免费区| 男女视频在线观看网站免费| 国产探花极品一区二区| 亚洲欧美日韩卡通动漫| 亚洲av五月六月丁香网| 久久天躁狠狠躁夜夜2o2o| 国产视频一区二区在线看| 99riav亚洲国产免费| 中文亚洲av片在线观看爽| 日本一二三区视频观看| 国产一区二区三区视频了| 别揉我奶头~嗯~啊~动态视频| 国产高潮美女av| 18禁黄网站禁片午夜丰满| 欧美+亚洲+日韩+国产| 亚洲成a人片在线一区二区| 日韩成人在线观看一区二区三区| 精品99又大又爽又粗少妇毛片 | 欧美+亚洲+日韩+国产| 午夜激情福利司机影院| 神马国产精品三级电影在线观看| 亚洲人与动物交配视频| 亚洲乱码一区二区免费版| 精品日产1卡2卡| 在线观看一区二区三区| 老司机福利观看| 婷婷六月久久综合丁香| 国产在线精品亚洲第一网站| 亚洲内射少妇av|