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

    納/微秒復(fù)合脈沖提高細(xì)胞電融合效率的仿真與實(shí)驗(yàn)研究

    2022-04-08 05:41:58李成祥李新皓姚陳果
    電工技術(shù)學(xué)報(bào) 2022年6期
    關(guān)鍵詞:融合實(shí)驗(yàn)

    李成祥 李新皓 柯 強(qiáng) 姚 成 姚陳果

    納/微秒復(fù)合脈沖提高細(xì)胞電融合效率的仿真與實(shí)驗(yàn)研究

    李成祥 李新皓 柯 強(qiáng) 姚 成 姚陳果

    (輸配電裝備及系統(tǒng)安全與新技術(shù)國(guó)家重點(diǎn)實(shí)驗(yàn)室(重慶大學(xué)) 重慶 400030)

    傳統(tǒng)的細(xì)胞電融合技術(shù)通常采用持續(xù)時(shí)間為微秒級(jí)的短時(shí)高壓電脈沖誘導(dǎo)細(xì)胞電穿孔來(lái)實(shí)現(xiàn)。然而,高幅值微秒脈沖在促進(jìn)細(xì)胞電融合的同時(shí)易導(dǎo)致細(xì)胞過(guò)度穿孔而死亡,從而限制了細(xì)胞的融合效率。為進(jìn)一步提高細(xì)胞電融合率,該文基于納秒脈沖電場(chǎng)所誘導(dǎo)的“超穿孔現(xiàn)象”,結(jié)合常規(guī)的微秒脈沖電融合技術(shù),提出納/微秒復(fù)合脈沖協(xié)同誘導(dǎo)細(xì)胞電融合的方法。以SP2/0小鼠骨髓瘤細(xì)胞為實(shí)驗(yàn)對(duì)象,利用COMSOL軟件分析復(fù)合脈沖作用下細(xì)胞膜上跨膜電壓與穿孔密度的分布,并進(jìn)行初步的電融合實(shí)驗(yàn),探究該融合方法對(duì)細(xì)胞融合率和死亡率的影響。仿真結(jié)果顯示,傳統(tǒng)微秒脈沖作用下,細(xì)胞膜的兩極與接觸區(qū)域的跨膜電壓均超過(guò)穿孔閾值,顯著穿孔面積超過(guò)10%。而等劑量的納/微秒復(fù)合脈沖作用下,跨膜電壓僅在細(xì)胞接觸區(qū)域達(dá)到較高值,電穿孔也集中于該區(qū)域。實(shí)驗(yàn)結(jié)果同樣證明,該方法能夠保證較低細(xì)胞死亡率(10%)的同時(shí)達(dá)到較高的細(xì)胞融合效率(75%)。該研究為細(xì)胞電融合技術(shù)的進(jìn)一步發(fā)展提供了新的思路。

    細(xì)胞電融合 納/微秒復(fù)合脈沖 融合效率 跨膜電壓 穿孔密度

    0 引言

    細(xì)胞融合作為細(xì)胞工程的一項(xiàng)核心技術(shù),在生產(chǎn)生活中發(fā)揮著越來(lái)越重要的作用[1]。采用人工干預(yù)的細(xì)胞融合為某些疾病提供了全新的治療方式,如糖尿病[2]、中樞神經(jīng)系統(tǒng)軸突再生[3]、自身免疫疾病[4]以及惡性腫瘤[5]等疾病均可通過(guò)細(xì)胞融合產(chǎn)物得到有效治療。根據(jù)所采用干預(yù)方式的不同,細(xì)胞融合技術(shù)可分為化學(xué)融合法[6](聚乙二醇)、生物融合法[7](仙臺(tái)病毒、日本血凝病毒)和物理融合法[8](離心、振動(dòng)、電融合)等。其中,細(xì)胞電融合作為物理融合方法的一種,相較化學(xué)融合法、生物融合法以及其他物理因子所誘導(dǎo)的細(xì)胞融合方法,具有可重復(fù)性高、操作簡(jiǎn)便、對(duì)細(xì)胞無(wú)毒、便于實(shí)時(shí)觀(guān)察等突出優(yōu)勢(shì),已經(jīng)逐漸發(fā)展成為最有效、最具前景的細(xì)胞融合技術(shù)之一[1]。

    事實(shí)上,細(xì)胞電融合技術(shù)是由人們所熟知的微秒脈沖電場(chǎng)(microsecond Pulsed Electric Fields,msPEF)誘導(dǎo)細(xì)胞電穿孔現(xiàn)象發(fā)展而來(lái)。維爾茲堡大學(xué)的U. Zimmermann[8]利用msPEF進(jìn)行電穿孔實(shí)驗(yàn)時(shí),發(fā)現(xiàn)經(jīng)過(guò)脈沖處理的部分細(xì)胞發(fā)生了融合,由此便創(chuàng)立了一種全新的細(xì)胞融合方法——電融合法?;趍sPEF的細(xì)胞電融合技術(shù)一經(jīng)提出便得到了國(guó)內(nèi)外學(xué)者的廣泛關(guān)注和持續(xù)研究,日本學(xué)者M(jìn). Senda等[9]利用msPEF成功實(shí)現(xiàn)了植物原生質(zhì)體的電融合,并迅速將該技術(shù)推廣到動(dòng)物、微生物的原生質(zhì)體融合中,推動(dòng)了電融合技術(shù)的跨越式發(fā)展;南開(kāi)大學(xué)汪和睦團(tuán)隊(duì)[10-11]實(shí)現(xiàn)了鼠源B淋巴細(xì)胞和骨髓瘤細(xì)胞的電融合,并在此基礎(chǔ)上制備了多種不同抗性的雜交瘤細(xì)胞株;盧布雅爾納大學(xué)的D. Miklavcic等[12]采用電場(chǎng)強(qiáng)度1.2kV/cm、脈寬100ms、頻率1Hz的8個(gè)脈沖分別單獨(dú)處理B16-F1細(xì)胞和CHO細(xì)胞,發(fā)現(xiàn)兩種細(xì)胞都能夠進(jìn)行融合,但融合效率有所差別;瓦倫西亞理工大學(xué)的J. Cardona- Costa等[13]采用不同電場(chǎng)強(qiáng)度、脈寬的msPEF處理太平洋牡蠣受精卵,探討受精卵融合效率與微秒脈沖參數(shù)的依賴(lài)關(guān)系;上海瑞金醫(yī)院李雁等[14]研究發(fā)現(xiàn),200V、25ms的直流脈沖是人體細(xì)胞核移植的最佳電融合參數(shù),融合率達(dá)到53.6%。除了上述理論研究成果之外,基于msPEF的細(xì)胞電融合在設(shè)備研制方面也取得了積極的進(jìn)展。重慶大學(xué)Hu Ning等[15]為提高融合效率,研制出高通量細(xì)胞電融合芯片;東京大學(xué)K. Sugahara等[16]開(kāi)發(fā)了有利于大細(xì)胞高效配對(duì)的動(dòng)態(tài)微陣列電融合設(shè)備。

    盡管基于msPEF的電融合技術(shù)已經(jīng)取得了長(zhǎng)足的進(jìn)展,但是在融合過(guò)程中仍存在較高的細(xì)胞死亡率,極大地限制了融合效率的進(jìn)一步提升。為此,M. Kanduser等[17]提出通過(guò)添加脂質(zhì)抗氧化劑對(duì)細(xì)胞進(jìn)行保護(hù)以有效降低細(xì)胞死亡率;Qu Pengxiang等[18]則嘗試通過(guò)添加褪黑素來(lái)達(dá)到降低死亡率的效果。然而,化學(xué)試劑的添加無(wú)疑會(huì)導(dǎo)致步驟的增加、生產(chǎn)工序的繁雜以及經(jīng)濟(jì)成本的提高。M.Kanduser等[19]研究發(fā)現(xiàn),融合細(xì)胞的死亡率之所以居高不下,其關(guān)鍵在于msPEF施加過(guò)程中導(dǎo)致細(xì)胞膜上出現(xiàn)大面積穿孔。因此,欲從根本上提升細(xì)胞電融合效率,應(yīng)在保證產(chǎn)生滿(mǎn)足電融合需求的穿孔的同時(shí),盡可能地避免細(xì)胞膜過(guò)度穿孔情況的出現(xiàn)。

    近年來(lái),脈寬相對(duì)更窄的納秒脈沖電場(chǎng)(nano- second Pulsed Electric Fields, nsPEF)以其獨(dú)特的“細(xì)胞內(nèi)電處理”效應(yīng)在生物醫(yī)學(xué)領(lǐng)域得到了越來(lái)越多的關(guān)注[20-21]。同時(shí),也有大量的仿真和實(shí)驗(yàn)研究證實(shí),nsPEF同樣可以作用在細(xì)胞外膜,并使其產(chǎn)生大量納米級(jí)尺寸的微孔,即所謂的“超穿孔”(supra- electroporation)現(xiàn)象。Hu Qingfang等[22]基于分子動(dòng)力學(xué)仿真發(fā)現(xiàn),電場(chǎng)強(qiáng)度100kV/cm、脈寬10ns的nsPEF能在細(xì)胞外膜上快速產(chǎn)生納米級(jí)微孔;L. M. Mir等[23]采用電場(chǎng)強(qiáng)度40kV/cm、脈寬10ns的單個(gè)nsPEF處理DC-3F細(xì)胞,發(fā)現(xiàn)nsPEF同樣可以使細(xì)胞膜穿孔,但穿孔尺寸較小,常用的大分子染色劑并不能通過(guò)細(xì)胞膜進(jìn)入細(xì)胞內(nèi)部;G. Saulis 等[24]也同樣證實(shí)nsPEF作用下細(xì)胞膜上產(chǎn)生孔洞數(shù)量遠(yuǎn)大于msPEF,但尺寸較小,恢復(fù)時(shí)間較短,不易造成細(xì)胞死亡。

    鑒于此,本文嘗試將nsPEF引入細(xì)胞電融合,使其首先在細(xì)胞膜上產(chǎn)生尺寸相對(duì)較小的微孔,繼而再施加低幅值msPEF對(duì)孔洞進(jìn)行維持和擴(kuò)張,以避免直接施加高幅值msPEF造成細(xì)胞膜大面積穿孔,從而降低死亡率,最終達(dá)到提高細(xì)胞融合效率的目的。為此,本文選擇小鼠骨髓瘤細(xì)胞SP2/0為研究對(duì)象,建立融合細(xì)胞動(dòng)態(tài)電穿孔模型,對(duì)比分析msPEF及納/微秒復(fù)合脈沖作用下,細(xì)胞膜的跨膜電壓(Transmembrane Voltage, TMV)分布及穿孔變化情況,并初步開(kāi)展細(xì)胞融合實(shí)驗(yàn),利用熒光顯微鏡觀(guān)察染色細(xì)胞的融合情況,探究納/微秒復(fù)合脈沖提高細(xì)胞電融合效率的可行性。

    1 模型建立與實(shí)驗(yàn)設(shè)計(jì)

    1.1 數(shù)值仿真模型的建立

    在COMSOL Multiphysics中建立相互接觸的雙細(xì)胞模型,采用有限元方法,對(duì)脈沖作用下細(xì)胞跨膜電壓的分布及細(xì)胞膜穿孔情況進(jìn)行仿真分析。

    1.1.1 雙細(xì)胞幾何模型

    由于細(xì)胞融合實(shí)驗(yàn)是在低滲融合緩沖液中進(jìn) 行[25],待融合細(xì)胞會(huì)略微膨脹接近球體,因此,本文將細(xì)胞簡(jiǎn)化為標(biāo)準(zhǔn)球體[19]。標(biāo)準(zhǔn)球體具有空間軸對(duì)稱(chēng)性,可以建立二維幾何模型來(lái)減少計(jì)算量,以縮短仿真時(shí)間。細(xì)胞融合仿真模型示意圖如圖1所示,用一個(gè)長(zhǎng)200mm、寬200mm的空間來(lái)模擬細(xì)胞融合時(shí)的外環(huán)境,兩側(cè)的條狀矩形用于模擬銅制極板,其中,左側(cè)極板施加高電位,右側(cè)極板接地。細(xì)胞的物理尺寸按SP2/0骨髓瘤細(xì)胞的實(shí)際尺寸進(jìn)行設(shè)置,其中,細(xì)胞半徑為5.6mm,細(xì)胞核半徑為4.723mm,細(xì)胞核質(zhì)比約為70%。令兩細(xì)胞相互接觸并沿電場(chǎng)方向排列,接觸區(qū)域垂直于電場(chǎng),長(zhǎng)度設(shè)置為2mm。將細(xì)胞各部分簡(jiǎn)化為各向同性的線(xiàn)性均勻介質(zhì),分別用電導(dǎo)率和介電常數(shù)表示細(xì)胞模型各部分的電特性,各區(qū)域邊界設(shè)置為電絕緣,計(jì)算時(shí)用三角形進(jìn)行剖分,并對(duì)需要重點(diǎn)求解的細(xì)胞膜部分進(jìn)行二次細(xì)化,剖分單元為12 368個(gè)域單元(999個(gè)邊界元)。

    1.1.2 動(dòng)態(tài)電穿孔模型

    圖1 細(xì)胞融合仿真模型示意圖

    式中,DY為跨膜電壓;i為膜的內(nèi)電位;o為膜的外電位。

    在外加電場(chǎng)的作用下,細(xì)胞內(nèi)外離子運(yùn)動(dòng)并聚集在膜兩側(cè),細(xì)胞膜和核膜的穿孔可以描述為磷脂雙分子層受電位變化影響形成大量親水微孔。穿孔的形成為跨膜電流提供新的通道,大量的穿孔將顯著提高膜的導(dǎo)電性。此時(shí),整體的跨膜電流密度可以表示[26-27]為

    式中,為電流密度;m0為膜的初始電導(dǎo)率;m為細(xì)胞外溶液介電常數(shù);m為細(xì)胞膜厚度;ep為跨膜電流密度增量,其計(jì)算公式由K. A. De Bruin和W. Krassowska提出[28],有

    其中

    式中,ep為流過(guò)單個(gè)微孔的電流;為孔密度;p為孔隙電導(dǎo)率;為微孔半徑;0為孔隙阻礙能量系數(shù);為孔隙相對(duì)密度;為法拉第常數(shù);為氣體常數(shù);為熱力學(xué)溫度??酌芏瘸蕜?dòng)態(tài)變化,可通過(guò)Smoluchowsky方程[28]求解,有

    式中,0為初始孔密度;ep為電穿孔特征電壓;為穿孔有效形成率系數(shù);為穿孔系數(shù)。ep決定了細(xì)胞膜的穿孔閾值T,文獻(xiàn)[27]給出了二者之間的關(guān)系,生物質(zhì)膜的T通常取值為500~1 000mV,本文選取T=1 000mV(ep=258mV)[19, 29]。

    將式(5)~式(7)代入式(4),可得到電流增量為

    膜電導(dǎo)率m與孔密度呈正相關(guān),膜動(dòng)態(tài)電導(dǎo)率由初始電導(dǎo)率及穿孔后的電導(dǎo)率增量相加得到[30-31],有

    電場(chǎng)作用下穿孔半徑的變化情況[30-31]為

    式中,為孔半徑擴(kuò)散系數(shù);為玻耳茲曼常數(shù);為空間排斥能;為邊緣能量;h和t為平流速度常數(shù);max為T(mén)MV=1V時(shí)的最大電動(dòng)力;*為親水孔最小半徑;eff為膜有效張力系數(shù),可求得eff為

    仿真研究中細(xì)胞所有幾何參數(shù)和電氣參數(shù)的取值及其來(lái)源見(jiàn)表1。

    表1 仿真模型幾何參數(shù)與電氣參數(shù)

    Tab.1 The geometric and electrical parameters used in the simulation

    1.1.3 復(fù)合脈沖電穿孔效應(yīng)仿真

    利用COMSOL軟件中的AC-DC模塊和PDE模塊對(duì)兩個(gè)接觸細(xì)胞在脈沖作用下的電穿孔情況進(jìn)行瞬態(tài)仿真。為了準(zhǔn)確地對(duì)比單msPEF和納/微秒復(fù)合脈沖的作用效果,本文采用等劑量的脈沖施加方式,脈沖劑量[32, 35]為

    式中,V為第個(gè)脈沖的幅值;T為第個(gè)脈沖的脈寬;為脈沖個(gè)數(shù)。

    脈沖波形示意圖如圖2所示。參考細(xì)胞融合儀(ECM2001,美國(guó)BTX Corporation)參數(shù)推薦,選用電場(chǎng)強(qiáng)度2.21kV/cm,脈寬40ms的單個(gè)msPEF作為對(duì)照組的脈沖施加方案(見(jiàn)圖2a)。鑒于單個(gè)nsPEF對(duì)細(xì)胞膜作用有限[19],本文經(jīng)過(guò)初步實(shí)驗(yàn)探究,選取5個(gè)電場(chǎng)強(qiáng)度6kV/cm,脈寬200ns的nsPEF,結(jié)合一個(gè)電場(chǎng)強(qiáng)度2kV/cm,脈寬40ms的msPEF作為復(fù)合脈沖的組合方式(見(jiàn)圖2b)。nsPEF的施加間隔為1ms,兩種脈沖轉(zhuǎn)換間隔設(shè)置為5ms。對(duì)整個(gè)模型區(qū)域進(jìn)行有限元仿真計(jì)算,由于雙細(xì)胞模型的對(duì)稱(chēng)性,為避免重復(fù),在進(jìn)行細(xì)胞膜分析時(shí)僅對(duì)部分膜區(qū)域進(jìn)行展示,具體區(qū)域如圖3加粗部分所示。

    圖2 脈沖波形示意圖

    圖3 仿真分析區(qū)域示意圖

    1.2 實(shí)驗(yàn)材料與方法

    本文通過(guò)開(kāi)展細(xì)胞電融合實(shí)驗(yàn),進(jìn)一步從實(shí)驗(yàn)的角度驗(yàn)證納/微秒復(fù)合脈沖作用細(xì)胞電融合的可行性。為更好地對(duì)比納/微秒復(fù)合脈沖及msPEF的作用效果,采用三種不同的脈沖劑量分別進(jìn)行實(shí)驗(yàn)并統(tǒng)計(jì)細(xì)胞融合效率和死亡率,所涉及的實(shí)驗(yàn)方法和材料如下。

    1.2.1 融合緩沖液配置

    低滲融合緩沖液通常用于脈沖處理細(xì)胞的實(shí)驗(yàn)過(guò)程,該緩沖液(90mOsm/L)包含Mg2+(0.1mmol/L)、Ca2+(0.1mmol/L)、牛血清蛋白(1mg/mL)以及D-甘露醇(所有試劑均購(gòu)自美國(guó)Sigma Corporation)。等滲融合緩沖液含120mmol/L D-甘露醇,該緩沖液僅用于洗滌細(xì)胞。用電導(dǎo)率儀測(cè)試兩種緩沖液的電導(dǎo)率,電導(dǎo)率均為0.01S/m。

    1.2.2 細(xì)胞處理

    實(shí)驗(yàn)所用的小鼠骨髓瘤細(xì)胞(SP2/0)購(gòu)自上海雅潔生物科技有限公司,培養(yǎng)在37℃、5% CO2加濕的環(huán)境中,含10%胎牛血清的RPMI1640培養(yǎng)基上。在脈沖處理前,先將清洗、消化、潤(rùn)洗后的細(xì)胞通過(guò)血細(xì)胞計(jì)數(shù)板進(jìn)行計(jì)數(shù),并配制成細(xì)胞密度為1×106/mL的細(xì)胞懸液。將0.5mL的細(xì)胞懸液和0.5mL的Hoechst 33342(10mg/mL,美國(guó)Solarbio Corporation)加入培養(yǎng)皿中充分混合,并置于培養(yǎng)箱中培養(yǎng)20min充分染色。染色完成后使用5mL磷酸緩沖鹽溶液(Phosphate Buffered Saline, PBS)洗滌兩次,以去除剩余染料。再用5mL等滲融合緩沖液洗滌一至兩次,即可將細(xì)胞置于500mL低滲融合緩沖液中,進(jìn)行脈沖處理。

    1.2.3 細(xì)胞電融合實(shí)驗(yàn)驗(yàn)證

    取Hoechst染色完成的細(xì)胞接種于平行電極(45-0105,美國(guó)BTX Corporation)進(jìn)行細(xì)胞融合實(shí)驗(yàn),該電極由帶有矩形極板的玻璃微載玻片制成,電極間距為3.2mm,體積為650mL。細(xì)胞電融合實(shí)驗(yàn)平臺(tái)為課題組自行設(shè)計(jì),如圖4所示。數(shù)字信號(hào)發(fā)生器和功率放大器用于生成正弦波形,并輸入脈沖發(fā)生器進(jìn)行統(tǒng)一控制。脈沖發(fā)生器中有納微秒脈沖發(fā)生模塊及切換模塊,可以輸出不同脈寬的方波脈沖或正弦電壓波形。不同波形之間的切換由繼電器(JPK43A234,中國(guó) WPVAC)實(shí)現(xiàn),并通過(guò)現(xiàn)場(chǎng)可編程門(mén)陣列(Field Programmable Gate Array, FPGA)(EP4CE6,中國(guó) ALINX)進(jìn)行控制。示波器(DPO3024,美國(guó)Tektronix)用于顯示電流傳感器(2100,美國(guó)Pearson)和高壓探頭(PPE5kV,美國(guó)LeCroy)采集的電流和電壓信號(hào)。融合實(shí)驗(yàn)所施加的脈沖參數(shù)見(jiàn)表2,脈沖施加結(jié)束繼續(xù)培養(yǎng)10min,在細(xì)胞懸液中加入20mL碘化丙啶(PI;1mg/mL,美國(guó)Life Technologies Corporation)并靜置8min,用熒光顯微鏡(美國(guó)Life Technologies Corporation)觀(guān)察Hoechst 33342以及PI對(duì)細(xì)胞的染色情況。Hoechst 33342用于細(xì)胞核染色,融合實(shí)驗(yàn)后通過(guò)觀(guān)察熒光場(chǎng)和明場(chǎng)圖像可以判斷兩個(gè)或多個(gè)細(xì)胞是否融合。PI用于死亡細(xì)胞的染色,若細(xì)胞被PI染成紅色,則認(rèn)為細(xì)胞已經(jīng)死亡。

    圖4 細(xì)胞融合實(shí)驗(yàn)平臺(tái)

    表2 實(shí)驗(yàn)脈沖參數(shù)

    Tab.2 Experimental pulse parameter

    1.2.4 細(xì)胞電融合效率統(tǒng)計(jì)分析

    細(xì)胞融合率通常被定義為參與融合的細(xì)胞數(shù)與細(xì)胞總數(shù)之比[19, 32],細(xì)胞死亡率為被PI染色的細(xì)胞數(shù)與細(xì)胞總數(shù)之比。若參與融合的細(xì)胞同時(shí)被PI染色,則計(jì)入死亡細(xì)胞,不計(jì)入融合細(xì)胞。

    所有數(shù)據(jù)均以四組獨(dú)立實(shí)驗(yàn)的均數(shù)±標(biāo)準(zhǔn)差表示,并對(duì)不同脈沖組合的實(shí)驗(yàn)結(jié)果進(jìn)行顯著性檢驗(yàn),用Origin 2019b進(jìn)行繪制。

    2 結(jié)果與分析

    2.1 跨膜電壓仿真分析

    在外加脈沖電場(chǎng)的作用下,細(xì)胞膜由于電容效應(yīng)的存在表現(xiàn)出充放電過(guò)程,跨膜電壓隨之發(fā)生變化,本節(jié)對(duì)單msPEF及納/微秒復(fù)合脈沖作用下,雙細(xì)胞模型跨膜電壓分布情況進(jìn)行分析。圖5所示為兩種脈沖作用后細(xì)胞膜表面跨膜電壓沿弧長(zhǎng)的分布情況(分析區(qū)域見(jiàn)圖3)。根據(jù)仿真結(jié)果可知,單msPEF作用下細(xì)胞膜跨膜電壓分布(虛線(xiàn)部分)呈“W”型,細(xì)胞模型的左右兩極及接觸區(qū)域(黑色點(diǎn)劃線(xiàn)范圍)的跨膜電壓均達(dá)到了較高值(>1V),說(shuō)明這些區(qū)域均存在較大可能被誘導(dǎo)產(chǎn)生顯著穿孔。在納/微秒復(fù)合脈沖作用時(shí)(實(shí)線(xiàn)部分),細(xì)胞接觸區(qū)域產(chǎn)生的跨膜電壓略大于msPEF作用,而在細(xì)胞兩極的跨膜電壓則僅為0.3~0.5V,遠(yuǎn)小于細(xì)胞膜的穿孔閾值。

    圖5 細(xì)胞膜跨膜電壓分布情況

    2.2 穿孔密度仿真分析

    當(dāng)細(xì)胞跨膜電壓達(dá)到穿孔閾值時(shí),會(huì)在膜上誘導(dǎo)產(chǎn)生大量的親水孔,這些孔洞的產(chǎn)生與分布和所施加的脈沖參數(shù)密切相關(guān)。本節(jié)分別選取細(xì)胞模型的兩個(gè)極點(diǎn)以及接觸區(qū)域的中點(diǎn)(見(jiàn)圖3)來(lái)研究脈沖作用時(shí)膜上孔密度隨時(shí)間的變化情況,并分析脈沖作用結(jié)束后發(fā)生顯著穿孔(孔密度>1×1013/m2)的區(qū)域[19,32]。

    圖6為msPEF及納/微秒復(fù)合脈沖分別誘導(dǎo)細(xì)胞融合產(chǎn)生的孔密度隨時(shí)間的變化情況。由圖6a可見(jiàn),msPEF單獨(dú)作用時(shí),前3ms孔密度并沒(méi)有發(fā)生明顯變化,隨后B點(diǎn)(見(jiàn)圖3)率先出現(xiàn)穿孔,并于7ms左右率先達(dá)到了顯著穿孔臨界點(diǎn),并最終維持在2.5×1014/m2左右。隨著脈沖的持續(xù)作用,兩級(jí)位置的孔密度分別于8ms(A點(diǎn))和10ms(C點(diǎn))開(kāi)始增加,并先后超過(guò)了1×1013/m2,但上升趨勢(shì)較緩,最終穩(wěn)定的孔密度在1013/m2量級(jí),略小于接觸區(qū)域孔密度。

    納/微秒復(fù)合脈沖作用下的細(xì)胞膜穿孔情況如圖6b所示。由圖可見(jiàn),在nsPEF作用時(shí)間段(前5ms),細(xì)胞膜接觸區(qū)域的孔密度迅速上升,并最終穩(wěn)定在1.17×1015/m2,在此期間,兩極位置的孔密度幾乎沒(méi)有變化。隨著msPEF的施加,A點(diǎn)和C點(diǎn)的孔密度開(kāi)始緩慢增加,增加的速率明顯小于單msPEF作用結(jié)果,直至脈沖施加結(jié)束,兩極的孔密度均沒(méi)有超過(guò)1×1013/m2,而B(niǎo)點(diǎn)的孔密度在msPEF作用下一直保持到脈沖作用結(jié)束。

    圖6 細(xì)胞膜孔密度隨時(shí)間變化情況

    顯著穿孔分布如圖7所示。為了更加直觀(guān)地展示細(xì)胞膜上顯著穿孔的分布,在圖7中將膜上孔密度大于1×1013/m2的區(qū)域用粗實(shí)線(xiàn)標(biāo)出。由圖7a可知,msPEF作用下,細(xì)胞膜表面約有10%的膜區(qū)域發(fā)生了顯著穿孔,除了接觸區(qū)域,細(xì)胞的左右兩極也發(fā)生了大面積的顯著穿孔,且靠近正電極的左側(cè)極點(diǎn)區(qū)域的顯著穿孔面積要大于靠近接地電極的右側(cè)極點(diǎn)區(qū)域。而納/微秒復(fù)合脈沖作用下(見(jiàn)圖7b),細(xì)胞膜顯著穿孔都集中在細(xì)胞膜接觸區(qū)域,其余區(qū)域并沒(méi)有發(fā)生顯著穿孔。

    圖7 顯著穿孔分布

    仿真結(jié)果表明,納/微秒復(fù)合脈沖這種脈沖組合方式能有效地將穿孔集中于細(xì)胞接觸區(qū)域并可以在較長(zhǎng)時(shí)間內(nèi)有效地維持,而對(duì)其余膜區(qū)域的影響較小。從理論上分析,該現(xiàn)象有利于降低融合細(xì)胞的死亡率,提高細(xì)胞的融合效率。

    2.3 細(xì)胞電融合實(shí)驗(yàn)結(jié)果分析

    使用已經(jīng)過(guò)Hoechst 33342染料對(duì)細(xì)胞核進(jìn)行染色的SP2/0小鼠骨髓瘤細(xì)胞進(jìn)行電融合實(shí)驗(yàn)。在進(jìn)行脈沖處理前需先施加20s的高頻正弦電壓(100V/cm,1MHz)進(jìn)行細(xì)胞排隊(duì),正弦電壓作用結(jié)束后在較短時(shí)間內(nèi)切換為脈沖電壓誘導(dǎo)細(xì)胞融合。細(xì)胞在電處理結(jié)束后繼續(xù)培養(yǎng)10min,再經(jīng)過(guò)PI染色,即可分別獲取經(jīng)Hoechst和PI染色的細(xì)胞明場(chǎng)圖像及紅色和藍(lán)色的細(xì)胞熒光圖像。每種脈沖施加方案分別進(jìn)行四組獨(dú)立實(shí)驗(yàn),每組實(shí)驗(yàn)中隨機(jī)選取3個(gè)視野進(jìn)行融合細(xì)胞、死亡細(xì)胞及細(xì)胞總量的計(jì)數(shù),取平均數(shù)計(jì)算各組細(xì)胞融合率和死亡率。

    實(shí)驗(yàn)中獲得的細(xì)胞排隊(duì)圖像及熒光成像結(jié)果如圖8所示,其中,圖8a為細(xì)胞排隊(duì)圖像,在交變電壓導(dǎo)致的介質(zhì)電泳效應(yīng)作用下,細(xì)胞沿電場(chǎng)方向排列成珍珠串狀。圖8b為假處理對(duì)照組,該組不進(jìn)行脈沖處理,其余步驟及培養(yǎng)時(shí)間均與實(shí)驗(yàn)組相同。圖8c和圖8d分別為總劑量為286kV2·ms/cm2的msPEF和復(fù)合脈沖作用電融合的實(shí)驗(yàn)結(jié)果。

    圖8 熒光染色實(shí)驗(yàn)成像

    圖9 細(xì)胞融合率、死亡率統(tǒng)計(jì)結(jié)果

    msPEF和納/微秒復(fù)合脈沖作用細(xì)胞電融合時(shí),細(xì)胞死亡率的變化情況如圖9b所示。由對(duì)照組死亡率可知,細(xì)胞排隊(duì)所用的高頻正弦電壓也會(huì)對(duì)細(xì)胞造成一定傷害,死亡率約為6%;在單msPEF作用時(shí),隨著脈沖劑量的上升,細(xì)胞死亡率從19.12%上升到了27.83%、53.16%,死亡率增長(zhǎng)十分顯著;而納/微秒復(fù)合脈沖作用下,隨著脈沖劑量的提升,細(xì)胞的死亡率分別為13.58%、9.76%、39.35%。由此可見(jiàn),相較于單個(gè)較高幅值的msPEF,高幅值nsPEF結(jié)合較低幅值的msPEF這種脈沖施加方式能夠顯著降低細(xì)胞死亡率,并且在一定的脈沖劑量范圍內(nèi),復(fù)合脈沖作用下的細(xì)胞死亡率不會(huì)隨著劑量的增加發(fā)生明顯變化。

    3 討論

    外加脈沖電場(chǎng)會(huì)改變細(xì)胞膜內(nèi)外原有的離子分布狀態(tài),造成跨膜電壓的變化。而跨膜電壓的提高會(huì)降低膜脂重排所需的能量,從而導(dǎo)致親水孔的形成[36]。通過(guò)分析細(xì)胞膜上跨膜電壓分布(見(jiàn)圖5)可知,持續(xù)時(shí)間足夠短的脈沖(納秒級(jí)別)可以在細(xì)胞膜接觸區(qū)域有選擇性地進(jìn)行大量穿孔,通過(guò)計(jì)算脈沖作用引起的跨細(xì)胞膜的孔密度變化進(jìn)一步證明了這一點(diǎn)(見(jiàn)圖6)。在脈寬40ms的msPEF作用下,細(xì)胞膜接觸區(qū)域的電穿孔同時(shí)伴隨著細(xì)胞兩極區(qū)域的較大范圍穿孔(見(jiàn)圖7a),靠近正極板的細(xì)胞左極由于靜息電位的影響比右極附近區(qū)域穿孔更加嚴(yán)重[37]。而在脈沖劑量相同的情況下,先施加5個(gè)高幅值nsPEF,再使用低幅值msPEF作用只會(huì)導(dǎo)致接觸區(qū)域發(fā)生顯著電穿孔,對(duì)細(xì)胞膜上其他區(qū)域則影響較?。ㄒ?jiàn)圖7b)。過(guò)度的細(xì)胞電穿孔容易導(dǎo)致細(xì)胞死亡,本文的數(shù)值仿真結(jié)果與實(shí)驗(yàn)觀(guān)察到,單msPEF作用時(shí),SP2/0骨髓瘤細(xì)胞的高死亡率相一致。

    有研究表明[24, 38],細(xì)胞膜上的電穿孔數(shù)量由施加脈沖的幅值決定,而孔洞的尺寸則更多受脈沖施加時(shí)間的影響。在較高幅值的納秒級(jí)別脈沖作用下,孔隙的生成速率會(huì)遠(yuǎn)遠(yuǎn)大于孔擴(kuò)張的速率(超穿孔現(xiàn)象),但較短的脈沖持續(xù)時(shí)間不足以支撐孔的擴(kuò)張,生成的微孔半徑僅為0.8nm左右,且極容易重新封閉。根據(jù)電融合基本條件,細(xì)胞融合過(guò)程的持續(xù)時(shí)間約為數(shù)分鐘[39],在這個(gè)過(guò)程中需要微孔一直保持“打開(kāi)”的狀態(tài)才能保證細(xì)胞的成功融合,若通過(guò)提高nsPEF的頻率和數(shù)量來(lái)延長(zhǎng)孔隙的維持時(shí)間,則可能因?yàn)閚sPEF的高頻透膜性誘導(dǎo)其他生物效應(yīng)的產(chǎn)生而對(duì)細(xì)胞活性產(chǎn)生影響[40-41]。因此,本文利用msPEF僅靶向作用細(xì)胞膜的特點(diǎn)[42],對(duì)已經(jīng)生成的孔隙進(jìn)行維持和擴(kuò)張。低幅值msPEF作用產(chǎn)生的跨膜電壓很難達(dá)到穿孔閾值,因此不易在其他區(qū)域產(chǎn)生大面積穿孔,造成細(xì)胞死亡。

    本文對(duì)SP2/0小鼠骨髓瘤細(xì)胞進(jìn)行了電融合實(shí)驗(yàn),從實(shí)驗(yàn)角度驗(yàn)證nsPEF協(xié)同msPEF進(jìn)行細(xì)胞電融合的可行性。圖9表明,在SP2/0細(xì)胞電融合實(shí)驗(yàn)中,復(fù)合脈沖作用的細(xì)胞融合效率普遍大于傳統(tǒng)msPEF作用結(jié)果,而細(xì)胞死亡率則普遍小于msPEF作用結(jié)果,該現(xiàn)象證實(shí)本文提出的納/微秒復(fù)合脈沖協(xié)同作用的細(xì)胞電融合效率高于傳統(tǒng)微秒脈沖。與msPEF作用相似,隨著脈沖劑量的增加,復(fù)合脈沖作用下細(xì)胞融合效率也經(jīng)歷了先增大后減小的過(guò)程,但通過(guò)細(xì)胞死亡率統(tǒng)計(jì)結(jié)果可知,復(fù)合脈沖作用時(shí),在一定的電場(chǎng)強(qiáng)度范圍內(nèi)細(xì)胞死亡率不會(huì)隨著電場(chǎng)強(qiáng)度的提升而發(fā)生顯著改變。此現(xiàn)象說(shuō)明,納/微秒復(fù)合脈沖可以在提高細(xì)胞融合效率的同時(shí)不增加細(xì)胞的死亡率。

    若將本文的研究?jī)?nèi)容擴(kuò)展到實(shí)際應(yīng)用中,可為脈沖組合的選取提供依據(jù)和理論支撐。本實(shí)驗(yàn)中所采用的5個(gè)脈寬200ns,電場(chǎng)強(qiáng)度6kV/cm的nsPEF結(jié)合1個(gè)脈寬40ms,電場(chǎng)強(qiáng)度2.5kV/cm的msPEF的復(fù)合脈沖組合方式雖然已經(jīng)超過(guò)了傳統(tǒng)msPEF的電融合效率,但是并沒(méi)有考慮到脈沖參數(shù)的優(yōu)化問(wèn)題,40ms的微秒脈沖也僅為細(xì)胞融合儀推薦使用的脈沖參數(shù),并不一定為復(fù)合脈沖作用細(xì)胞電融合的最優(yōu)脈沖參數(shù),最優(yōu)脈沖組合問(wèn)題仍需后續(xù)進(jìn)一步開(kāi)展實(shí)驗(yàn)進(jìn)行探究。另一方面,本文僅通過(guò)細(xì)胞融合率來(lái)表征復(fù)合脈沖作用細(xì)胞電融合的結(jié)果,而在實(shí)際雜交瘤制備過(guò)程中,融合的細(xì)胞不一定都能有效存活并發(fā)展成為雜交瘤細(xì)胞,因此,在后續(xù)研究中需要引入雜交瘤培養(yǎng)技術(shù),來(lái)探究納/微秒復(fù)合脈沖電融合產(chǎn)生雜交瘤細(xì)胞的實(shí)際生產(chǎn)效率。

    4 結(jié)論

    本文從數(shù)值仿真和實(shí)驗(yàn)分析兩個(gè)方面,開(kāi)展了納/微秒復(fù)合脈沖作用細(xì)胞電融合融合效率及細(xì)胞死亡率的研究,證明了該融合方法提高細(xì)胞融合效率、降低細(xì)胞死亡率的可行性。

    對(duì)脈沖作用下細(xì)胞跨膜電壓分布及穿孔動(dòng)態(tài)過(guò)程的仿真研究中,發(fā)現(xiàn)納/微秒復(fù)合脈沖具有明顯的細(xì)胞接觸區(qū)域靶向穿孔的特性,能夠?qū)⒋┛准杏诩?xì)胞接觸位置并維持較長(zhǎng)時(shí)間,而對(duì)細(xì)胞膜其他區(qū)域的影響較小。該研究結(jié)果證明,納/微秒復(fù)合脈沖能夠結(jié)合兩種脈沖各自的優(yōu)勢(shì),避免了高幅值微秒脈沖作用下細(xì)胞膜過(guò)度穿孔情況的出現(xiàn)。通過(guò)SP2/0小鼠骨髓瘤細(xì)胞的電融合實(shí)驗(yàn)研究,證實(shí)納/微秒復(fù)合脈沖與微秒脈沖相比,可以降低細(xì)胞死亡率,提高細(xì)胞融合效率。同時(shí),該融合方法的作用效果同樣對(duì)脈沖劑量具有一定的依賴(lài)性,最優(yōu)的電融合脈沖組合仍需進(jìn)一步探索。

    [1] Bajaj Y. Plant protoplasts and genetic engineering II[M]. Berlin: Springer-Verlag, 1989.

    [2] Mccluskey J T, Hamid M, Guo-Parke H, et al. Development and functional characterization of insulin-releasing human pancreatic beta cell lines produced by electrofusion[J]. Journal of Biological Chemistry, 2011, 286(25): 21982-21992.

    [3] Sretavan D W, Chang W, Hawkes E, et al. Microscale surgery on single axons[J]. Neurosurgery, 2005, 57(4): 635-646.

    [4] Avigan D, Rosenblatt J, Kufe D. Dendritic/tumor fusion cells as cancer vaccines[J]. Seminars in Oncology, 2012, 39(3): 287-295.

    [5] Trontelj K, Reber?ek M, Kandu?er M, et al. Optimi- zation of bulk cell electrofusion in vitro for pro- duction of human-mouse heterohybridoma cells[J]. Bioelectrochemistry, 2008, 74(1): 124-129.

    [6] Golestani R, Pourfathollah A A, Moazzeni S M. Cephalin as an efficient fusogen in hybridoma technology: can it replace poly ethylene glycol?[J]. Hybridoma, 2007, 26(5): 296-301.

    [7] Okada Y, Koseki I, Kim J, et al. Modification of cell membranes with viral envelopes during fusion of cells with HVJ (Sendai virus): I. Interaction between cell membranes and virus in the early stage[J]. Experimental Cell Research, 1975, 93(2): 368-378.

    [8] Zimmermann U, Neil G A. Electromanipulation of Cells[M]. New York: CRC Press, 1996.

    [9] Senda M, Takeda J, Abe S, et al. Induction of cell fusion of plant protoplasts by electrical stimulation[J]. Plant and Cell Physiology, 1979, 20(7): 1441-1443.

    [10] 謝廷棟, 汪和睦, 魯玉瓦, 等. 鼠B淋巴細(xì)胞和骨髓瘤細(xì)胞的電融合[J]. 中國(guó)生物醫(yī)學(xué)工程學(xué)報(bào), 1989, 8(1): 25-31.

    Xie Tingdong, Wang Hemu, Lu Yuwa, et al. Elec- trical fusion of mouse B lymphocytes and myeloma cells[J]. Chinese Journal of Biomedical Engineering, 1989, 8(1): 25-31.

    [11] 紀(jì)其雄, 魯玉瓦, 汪和睦. 運(yùn)用高效率電融合技術(shù)建立抗-hCG雜交瘤細(xì)胞株[J]. 生物物理學(xué)報(bào), 1993, 9(3): 507.

    Ji Qixiong, Lu Yuwa, Wang Hemu. A high- efficiency electrofusion technique was used to establish anti-hCG hybridoma cell line[J]. Acta Biophysica Sinica, 1993, 9(3): 507.

    [12] Usaj M, Flisar K, Miklavcic D, et al. Electrofusion of B16-F1 and CHO cells: the comparison of the pulse first and contact first protocols[J]. Bioelectro- chemistry, 2013, 89(2): 34-41.

    [13] Cardona-Costa J, García-Ximénez F, Espinós F J. Definition of fusion medium and electric parameters for efficient zygote electrofusion in the Pacific oyster (Crassostrea gigas)[J]. Theriogenology, 2010, 74(5): 828-834.

    [14] 李雁, 馮云, 孫貽娟, 等. 人卵母細(xì)胞與體細(xì)胞的電融合效率[J]. 上海交通大學(xué)學(xué)報(bào)(醫(yī)學(xué)版), 2007, 27(7): 837-840.

    Li Yan, Feng Yun, Sun Yijuan, et al. Efficiency of electrofusion of human somatic cells with oocytes[J]. Journal of Shanghai Jiaotong University (Medical Science), 2007, 27(7): 837-840.

    [15] Hu Ning, Yang Jun, Joo S W, et al. Cell electrofusion in microfluidic devices: a review[J]. Sensors and Actuators B: Chemical, 2013, 178: 63-85.

    [16] Sugahara K, Morimoto Y, Takamoei S, et al. A dynamic microarray device for pairing and electro- fusion of giant unilamellar vesicles[J]. Sensors and Actuators B: Chemical, 2020, 311: 127922.

    [17] Kanduser M, Imsirovic M K, Usaj M. The effect of lipid antioxidanta-tocopherol on cell viability and electrofusion yield of B16-F1 cells in vitro[J]. The Journal of Membrane Biology, 2019, 252(1): 105-114.

    [18] Qu Pengxiang, Shen Chong, Du Yue, et al. Melatonin protects rabbit somatic cell nuclear transfer (SCNT) embryos from electrofusion damage[J]. Scientific Reports, 2020, 10(1): 1-9.

    [19] Rems L, Usaj M, Kanduser M, et al. Cell electro- fusion using nanosecond electric pulses[J]. Scientific Reports, 2013, 3(1): 1-10.

    [20] Stacey M W, Fox P, Buescher S, et al. Nanosecond pulsed electric field induced cytoskeleton, nuclear membrane and telomere damage adversely impact cell survival[J]. Bioelectrochemistry, 2011, 82(2): 131-134.

    [21] 劉紅梅, 董守龍, 寧郡怡, 等. 納秒脈沖高頻透膜效應(yīng)優(yōu)先殺傷化療抗性腫瘤細(xì)胞的仿真與實(shí)驗(yàn)研究[J]. 電工技術(shù)學(xué)報(bào), 2019, 34(22): 4839-4848.

    Liu Hongmei, Dong Shoulong, Ning Junyi, et al. Simulation and experimental study on preferential killing of chemoresistance tumor cells induced by the high-frequency permeation effect of nanosecond pulse field[J]. Transactions of China Electrotechnical Society, 2019, 34(22): 4839-4848.

    [22] Hu Qingfang, Viswanadham S, Joshi R P, et al. Simulations of transient membrane behavior in cells subjected to a high-intensity ultrashort electric pulse[J]. Physical Review E, 2005, 71(3): 31914.

    [23] Silve A, Leray I, Mir L M. Demonstration of cell membrane permeabilization to medium-sized mole- cules caused by a single 10ns electric pulse[J]. Bioelectrochemistry, 2012, 87(10): 260-264.

    [24] Saulis G, Saul? R. Size of the pores created by an electric pulse: microsecond vs millisecond pulses[J]. Biochimica et Biophysica Acta (BBA)-Biomembranes, 2012, 1818(12): 3032-3039.

    [25] Usaj M, Kanduser M. The systematic study of the electroporation and electrofusion of B16-F1 and CHO cells in isotonic and hypotonic buffer[J]. The Journal of Membrane Biology, 2012, 245(9): 583-590.

    [26] Pucihar G, Miklavcic D, Kotnik T. A time-dependent numerical model of transmembrane voltage indu- cement and electroporation of irregularly shaped cells[J]. IEEE Transactions on Biomedical Engin- eering, 2009, 56(5): 1491-1501.

    [27] Krassowska W, Filev P D. Modeling electroporation in a single cell[J]. Biophysical Journal, 2007, 92(2): 404-417.

    [28] De Bruin K A, Krassowska W. Modeling electro- poration in a single cell. I. Effects of field strength and rest potential[J]. Biophysical Journal, 1999, 77(3): 1213-1224.

    [29] Liu Linying, Mao Zheng, Zhang Jianhua, et al. The influence of vesicle shape and medium conductivity on possible electrofusion under a pulsed electric field[J]. PloS One, 2016, 11(7): e0158739.

    [30] Neu J C, Smith K C, Krassowska W. Electrical energy required to form large conducting pores[J]. Bioelectrochemistry, 2003, 60(1): 107-114.

    [31] Yao Chenguo, Liu Hongmei, Zhao Yajun, et al. Analysis of dynamic processes in single-cell electro- poration and their effects on parameter selection based on the finite-element model[J]. IEEE Transa- ctions on Plasma Science, 2017, 45(5): 889-900.

    [32] Li Chengxiang, Ke Qiang, Yao Cheng, et al. Com- parison of bipolar and unipolar pulses in cell elec- trofusion: simulation and experimental research[J]. IEEE Transactions on Biomedical Engineering, 2018, 66(5): 1353-1360.

    [33] Polevaya Y, Ermolina I, Schlesinger M, et al. Time domain dielectric spectroscopy study of human cells: II. Normal and malignant white blood cells[J]. Biochimica et Biophysica Acta (BBA)-Biomembranes, 1999, 1419(2): 257-271.

    [34] Garner A L, Chen G, Chen Nianyong, et al. Ultrashort electric pulse induced changes in cellular dielectric properties[J]. Biochemical and Biophysical Research Communications, 2007, 362(1): 139-144.

    [35] 姚陳果, 鄭爽, 劉紅梅, 等. 面向臨床多針消融腫瘤的不可逆電穿孔量效關(guān)系模型[J]. 電工技術(shù)學(xué)報(bào), 2020, 35(11): 2491-2498.

    Yao Chenguo, Zheng Shuang, Liu Hongmei, et al. Model of irreversible electroporation dose-effect relationship for clinical multi-needle ablation tumors[J]. Transactions of China Electrotechnical Society, 2020, 35(11): 2491-2498.

    [36] Kotnik T, Miklav?i? D. Theoretical evaluation of voltage inducement on internal membranes of biological cells exposed to electric fields[J]. Bio- physical Journal, 2006, 90(2): 480-491.

    [37] 姚陳果, 劉紅梅, 趙亞軍, 等. 基于有限元分析的球形單細(xì)胞及不規(guī)則真實(shí)細(xì)胞電穿孔動(dòng)態(tài)過(guò)程仿真[J]. 高電壓技術(shù), 2016, 42(8): 2479-2486.

    Yao Chenguo, Liu Hongmei, Zhao Yajun, et al. Simulation on the process of electroporation in spherical single-cell and irregularly shaped real cells based on finite elements analysis[J]. High Voltage Engineering, 2016, 42(8): 2479-2486.

    [38] 姚陳果, 呂彥鵬, 趙亞軍, 等. 基于能量概率與微孔力模型的脈沖電場(chǎng)對(duì)細(xì)胞電穿孔動(dòng)態(tài)過(guò)程的仿真分析[J]. 電工技術(shù)學(xué)報(bào), 2016, 31(23): 141-149.

    Yao Chenguo, Lü Yanpeng, Zhao Yajun, et al. Simulation analysis on dynamic process of elec- troporation by the model based on energy probability and pore force in cell exposed to pulsed electric field[J]. Transactions of China Electrotechnical Society, 2016, 31(23): 141-149.

    [39] Pucihar G, Kotnik T, Teissié J, et al. Electro- permeabilization of dense cell suspensions[J]. European Biophysics Journal, 2007, 36(3): 173-185.

    [40] 米彥, 李盼, 劉權(quán), 等. 納秒脈沖電場(chǎng)聯(lián)合多壁碳納米管對(duì)皮膚癌細(xì)胞活性的劑量效應(yīng)研究[J]. 電工技術(shù)學(xué)報(bào), 2019, 34(22): 4849-4857.

    Mi Yan, Li Pan, Liu Quan, et al. Dose effect of the activity of skin cancer cells treated by nanosecond pulsed electric field combined with multi-walled carbon nanotubes[J]. Transactions of China Elec- trotechnical Society, 2019, 34(22): 4849-4857.

    [41] 米彥, 劉權(quán), 李盼, 等. 低強(qiáng)度納秒脈沖電場(chǎng)聯(lián)合靶向金納米棒對(duì)A375黑色素瘤細(xì)胞的殺傷效果研究[J]. 電工技術(shù)學(xué)報(bào), 2020, 35(12): 2534-2544.

    Mi Yan, Liu Quan, Li Pan, et al. Study of killing effects of A375 melanoma cells treated by low intensity nanosecond pulsed electric fields combined with targeted gold nanorods[J]. Transactions of China Electrotechnical Society, 2019, 35(12): 2534-2544.

    [42] 姚陳果, 寧郡怡, 劉紅梅, 等. 微/納秒脈沖電場(chǎng)靶向不同尺寸腫瘤細(xì)胞內(nèi)外膜電穿孔效應(yīng)研究[J]. 電工技術(shù)學(xué)報(bào), 2020, 35(1): 115-124.

    Yao Chenguo, Ning Junyi, Liu Hongmei, et al. Study of electroporation effect of different size tumor cells targeted by micro-nanosecond pulsed electric field[J]. Transactions of China Electrotechnical Society, 2020, 35(1): 115-124.

    Simulation and Experimental Research on Improving the Efficiency of Cell Electrofusion by Nano-Microsecond Pulsed Electric Fields

    (State Key Laboratory of Power Transmission Equipment & System Security and New Technology Chongqing University Chongqing 400030 China)

    Cell electrofusion technology is usually performed by short-term high-voltage electrical pulses with a duration of microseconds to induce cell electroporation. However, exposure to high amplitude microsecond pulses can easily lead to excessive cell electroporation while promoting cell electrofusion, resulting in the death of the fused cells, thereby limiting the efficiency of cell fusion. To increase the rate of cell electrofusion, this paper proposes a method of synergistic induction of cell electrofusion by combining nanosecond and microsecond pulsed electric fields, based on the phenomenon of “super-electroporation” induced by nanosecond pulsed electric field. In this paper, SP2/0 mouse myeloma cells are used as the experimental object. The COMSOL software is used to analyze the distribution of transmembrane voltage and electroporation density on the cell membrane under the exposure of the composite pulse, and preliminary electrofusion experiments are performed to explore the mortality and fusion rate of the fusion method. The simulation results show that under the traditional microsecond pulse, the transmembrane voltage across the two poles of the membrane and the contact area both exceed the perforation threshold, and the significant perforation area exceeds 10%. Besides, under the action of the same dose of nano/microsecond composite pulse, the transmembrane voltage only reaches a higher value in the cell contact area, and the electroporation is also concentrated in this area. Experimental results also prove that this method can guarantee a lower mortality rate (10%) while achieving a higher cell fusion efficiency (75%), which provides new ideas for the further development of cell electrofusion technology.

    Cell electrofusion, nano-microsecond pulsed electric fields, fusion efficiency, transmembrane voltage, electroporation density

    10.19595/j.cnki.1000-6753.tces.201315

    TM11; TM836; Q64

    李成祥 男,1979年生,研究員,博士生導(dǎo)師,研究方向?yàn)樯镝t(yī)學(xué)中的電工新技術(shù)、脈沖強(qiáng)磁場(chǎng)工業(yè)應(yīng)用技術(shù)及電力設(shè)備在線(xiàn)監(jiān)測(cè)及故障診斷。E-mail: lichengxiang@cqu.edu.cn(通信作者)

    李新皓 男,1997年生,碩士研究生,研究方向?yàn)樯镝t(yī)學(xué)中的電工新技術(shù)。E-mail: lixinhao@cqu.edu.cn

    2020-09-29

    2021-05-20

    國(guó)家自然科學(xué)基金專(zhuān)項(xiàng)基金(81727802)和國(guó)家自然科學(xué)基金面上(51677017)資助項(xiàng)目。

    (編輯 崔文靜)

    猜你喜歡
    融合實(shí)驗(yàn)
    記一次有趣的實(shí)驗(yàn)
    一次函數(shù)“四融合”
    微型實(shí)驗(yàn)里看“燃燒”
    村企黨建聯(lián)建融合共贏(yíng)
    融合菜
    從創(chuàng)新出發(fā),與高考數(shù)列相遇、融合
    寬窄融合便攜箱IPFS500
    《融合》
    做個(gè)怪怪長(zhǎng)實(shí)驗(yàn)
    NO與NO2相互轉(zhuǎn)化實(shí)驗(yàn)的改進(jìn)
    国产精品一区二区性色av| 亚洲在线观看片| av国产久精品久网站免费入址| 国产黄a三级三级三级人| 亚洲真实伦在线观看| 岛国毛片在线播放| 欧美日本视频| 尾随美女入室| 精品亚洲乱码少妇综合久久| 男人爽女人下面视频在线观看| 18禁在线播放成人免费| 我的女老师完整版在线观看| 99热国产这里只有精品6| 国产永久视频网站| 最近中文字幕2019免费版| 欧美激情国产日韩精品一区| 九九爱精品视频在线观看| 性色av一级| 日韩亚洲欧美综合| av在线播放精品| 色哟哟·www| 高清av免费在线| 丰满乱子伦码专区| 成人漫画全彩无遮挡| 啦啦啦中文免费视频观看日本| 91久久精品国产一区二区成人| 午夜福利视频精品| 国语对白做爰xxxⅹ性视频网站| 嫩草影院精品99| 狂野欧美激情性xxxx在线观看| 国产亚洲最大av| 久久ye,这里只有精品| 搡老乐熟女国产| 99九九线精品视频在线观看视频| 亚洲国产欧美人成| 自拍欧美九色日韩亚洲蝌蚪91 | 国产成人精品婷婷| 日韩伦理黄色片| 欧美bdsm另类| 久久精品久久久久久噜噜老黄| 少妇被粗大猛烈的视频| 汤姆久久久久久久影院中文字幕| 国产免费一级a男人的天堂| 天堂中文最新版在线下载 | 亚洲天堂av无毛| 婷婷色综合大香蕉| 九九久久精品国产亚洲av麻豆| 精品国产三级普通话版| 亚洲国产av新网站| 精品国产三级普通话版| 又大又黄又爽视频免费| av免费在线看不卡| 国产成人免费观看mmmm| 亚洲国产欧美在线一区| 一级a做视频免费观看| 日韩亚洲欧美综合| 日韩伦理黄色片| 久久国产乱子免费精品| 亚洲天堂av无毛| 少妇人妻久久综合中文| 亚洲国产精品专区欧美| 黑人高潮一二区| 欧美另类一区| 美女国产视频在线观看| 日本wwww免费看| 国产精品一及| 我的女老师完整版在线观看| 亚州av有码| 有码 亚洲区| av一本久久久久| 观看美女的网站| 六月丁香七月| 热99国产精品久久久久久7| 午夜福利视频1000在线观看| 国产精品熟女久久久久浪| 人妻少妇偷人精品九色| 狂野欧美白嫩少妇大欣赏| 青春草国产在线视频| 啦啦啦啦在线视频资源| 精品99又大又爽又粗少妇毛片| 亚洲av在线观看美女高潮| 日本一本二区三区精品| 成人欧美大片| 国产伦理片在线播放av一区| 99re6热这里在线精品视频| 国产成人精品福利久久| 黄色日韩在线| 国产有黄有色有爽视频| 日韩欧美 国产精品| 中文天堂在线官网| 免费观看a级毛片全部| 国产精品av视频在线免费观看| 国产精品一及| 国产精品一及| 国产精品爽爽va在线观看网站| 亚洲精品国产成人久久av| 久久97久久精品| 亚洲精品国产av成人精品| 国产精品久久久久久精品电影小说 | 99热这里只有是精品50| 久久久a久久爽久久v久久| 国产亚洲av片在线观看秒播厂| 一级爰片在线观看| 国产精品人妻久久久久久| 人妻制服诱惑在线中文字幕| 大香蕉久久网| 在线a可以看的网站| 成年女人看的毛片在线观看| 国产黄色免费在线视频| 国产午夜精品久久久久久一区二区三区| 免费看a级黄色片| 伦理电影大哥的女人| 插逼视频在线观看| 男女啪啪激烈高潮av片| 免费观看a级毛片全部| 国产成人91sexporn| 国产欧美另类精品又又久久亚洲欧美| 国产成人福利小说| 欧美日韩视频高清一区二区三区二| 男男h啪啪无遮挡| 久久久色成人| 亚洲国产精品成人久久小说| 国产乱来视频区| 99热这里只有精品一区| 久久精品熟女亚洲av麻豆精品| 亚洲性久久影院| 美女脱内裤让男人舔精品视频| 精品人妻偷拍中文字幕| 3wmmmm亚洲av在线观看| 亚洲国产欧美在线一区| 偷拍熟女少妇极品色| 亚洲精品国产av成人精品| 99热全是精品| 蜜桃久久精品国产亚洲av| 午夜精品一区二区三区免费看| 国产探花极品一区二区| 国产精品久久久久久久电影| 中文字幕av成人在线电影| 亚洲欧美成人综合另类久久久| 亚洲精品乱久久久久久| 亚洲av一区综合| 一本久久精品| 人妻系列 视频| 国产 一区 欧美 日韩| 国产午夜精品一二区理论片| 少妇高潮的动态图| 亚洲欧美精品自产自拍| 97人妻精品一区二区三区麻豆| 日韩大片免费观看网站| 亚洲精品久久午夜乱码| 婷婷色综合www| 成人二区视频| 亚洲在线观看片| 久久鲁丝午夜福利片| 天天躁夜夜躁狠狠久久av| 久久久久久伊人网av| 在线观看国产h片| 2021少妇久久久久久久久久久| 亚洲精品第二区| 777米奇影视久久| 国产成人福利小说| 嫩草影院入口| 99热全是精品| 国产色爽女视频免费观看| 一个人看视频在线观看www免费| 免费观看av网站的网址| 国产午夜精品久久久久久一区二区三区| 国产成人91sexporn| 激情 狠狠 欧美| 午夜福利视频1000在线观看| 香蕉精品网在线| 国产在线男女| 色网站视频免费| 只有这里有精品99| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 免费电影在线观看免费观看| 国产亚洲91精品色在线| 激情 狠狠 欧美| 18禁裸乳无遮挡免费网站照片| 欧美丝袜亚洲另类| 国产综合懂色| 久久久久久久久久成人| 亚洲欧美精品自产自拍| 最新中文字幕久久久久| 国产一区二区亚洲精品在线观看| 男男h啪啪无遮挡| 午夜爱爱视频在线播放| 高清毛片免费看| 亚洲无线观看免费| 欧美日韩一区二区视频在线观看视频在线 | 国产白丝娇喘喷水9色精品| 22中文网久久字幕| 黄片wwwwww| 精品久久久久久久久av| 在线观看美女被高潮喷水网站| 亚洲av免费高清在线观看| 欧美高清性xxxxhd video| 熟女av电影| 全区人妻精品视频| av国产免费在线观看| 国产一区二区亚洲精品在线观看| 免费观看无遮挡的男女| 乱系列少妇在线播放| av黄色大香蕉| 人人妻人人看人人澡| 我的老师免费观看完整版| 久久精品久久久久久噜噜老黄| 欧美国产精品一级二级三级 | 波野结衣二区三区在线| 少妇人妻一区二区三区视频| 深夜a级毛片| 97超视频在线观看视频| 国产精品女同一区二区软件| 99久国产av精品国产电影| 九九爱精品视频在线观看| 国产毛片在线视频| 午夜福利网站1000一区二区三区| 草草在线视频免费看| 日本爱情动作片www.在线观看| 久久久久网色| 亚洲不卡免费看| 午夜免费男女啪啪视频观看| 国产在线男女| av国产久精品久网站免费入址| 18禁在线播放成人免费| 69av精品久久久久久| 国产免费又黄又爽又色| 亚洲美女搞黄在线观看| 国产精品不卡视频一区二区| 91久久精品电影网| 亚洲内射少妇av| 黄色视频在线播放观看不卡| 欧美一级a爱片免费观看看| 精品国产露脸久久av麻豆| 欧美bdsm另类| av在线播放精品| 欧美成人精品欧美一级黄| 看黄色毛片网站| 永久网站在线| 九九久久精品国产亚洲av麻豆| 久久97久久精品| 久热久热在线精品观看| 在线播放无遮挡| 国产黄色免费在线视频| 在线观看国产h片| 狂野欧美白嫩少妇大欣赏| 成人综合一区亚洲| 人体艺术视频欧美日本| 日韩制服骚丝袜av| 日日摸夜夜添夜夜爱| 国产一级毛片在线| 国产伦理片在线播放av一区| 深夜a级毛片| 欧美少妇被猛烈插入视频| 一本一本综合久久| 国产久久久一区二区三区| 国产精品一区二区三区四区免费观看| 黄片wwwwww| 一区二区三区四区激情视频| 成人鲁丝片一二三区免费| 免费看不卡的av| 日日啪夜夜撸| 偷拍熟女少妇极品色| 成年免费大片在线观看| 中文资源天堂在线| 久久久亚洲精品成人影院| 国国产精品蜜臀av免费| 成年人午夜在线观看视频| 精品视频人人做人人爽| 国产爱豆传媒在线观看| 黄色配什么色好看| 婷婷色综合大香蕉| 久久6这里有精品| 特大巨黑吊av在线直播| 一个人看视频在线观看www免费| 久久国产乱子免费精品| 日本爱情动作片www.在线观看| 亚洲国产成人一精品久久久| 国产探花极品一区二区| 免费观看a级毛片全部| 亚洲精品久久午夜乱码| av线在线观看网站| 精品熟女少妇av免费看| 午夜激情久久久久久久| 最近的中文字幕免费完整| 免费av观看视频| 2022亚洲国产成人精品| 久久综合国产亚洲精品| 欧美极品一区二区三区四区| 精品酒店卫生间| 精品一区二区免费观看| 嫩草影院入口| 老师上课跳d突然被开到最大视频| 人体艺术视频欧美日本| 我的女老师完整版在线观看| 亚洲精品国产av成人精品| 日本三级黄在线观看| 老司机影院毛片| 日韩三级伦理在线观看| 国产在视频线精品| 亚洲在线观看片| 大话2 男鬼变身卡| 中文精品一卡2卡3卡4更新| 亚洲欧美一区二区三区国产| 97超碰精品成人国产| 久久99精品国语久久久| 听说在线观看完整版免费高清| 欧美丝袜亚洲另类| 欧美老熟妇乱子伦牲交| 亚洲自拍偷在线| av在线播放精品| 精品99又大又爽又粗少妇毛片| 99视频精品全部免费 在线| 亚洲av免费在线观看| 一区二区三区乱码不卡18| 熟妇人妻不卡中文字幕| h日本视频在线播放| 美女国产视频在线观看| 丰满少妇做爰视频| 又爽又黄a免费视频| av.在线天堂| 精品久久国产蜜桃| 久久精品国产亚洲网站| 国产成人一区二区在线| 日韩伦理黄色片| 久久久久九九精品影院| 天堂网av新在线| 永久免费av网站大全| 亚洲国产欧美人成| 韩国高清视频一区二区三区| 亚洲经典国产精华液单| 七月丁香在线播放| 色婷婷久久久亚洲欧美| 80岁老熟妇乱子伦牲交| 国产午夜精品久久久久久一区二区三区| 嫩草影院新地址| 亚洲人成网站在线播| 免费观看在线日韩| 青春草亚洲视频在线观看| 香蕉精品网在线| 99久久九九国产精品国产免费| 少妇人妻 视频| 最后的刺客免费高清国语| 成人一区二区视频在线观看| freevideosex欧美| 国产亚洲最大av| 欧美日韩亚洲高清精品| 永久网站在线| 国产一区二区亚洲精品在线观看| 超碰97精品在线观看| 一个人看的www免费观看视频| 最近最新中文字幕免费大全7| 成年人午夜在线观看视频| 黄色欧美视频在线观看| 免费高清在线观看视频在线观看| 91久久精品国产一区二区成人| tube8黄色片| 天天躁夜夜躁狠狠久久av| 在线观看美女被高潮喷水网站| 香蕉精品网在线| 欧美区成人在线视频| 国产亚洲午夜精品一区二区久久 | 大话2 男鬼变身卡| 亚洲精品国产色婷婷电影| 狠狠精品人妻久久久久久综合| 日本免费在线观看一区| 成人特级av手机在线观看| 女的被弄到高潮叫床怎么办| 高清午夜精品一区二区三区| 精品一区在线观看国产| 中文字幕制服av| 制服丝袜香蕉在线| 国产国拍精品亚洲av在线观看| 久久人人爽人人片av| 成年版毛片免费区| 麻豆国产97在线/欧美| 一级av片app| 天堂中文最新版在线下载 | 熟女人妻精品中文字幕| 亚洲最大成人中文| 国产av不卡久久| 又爽又黄无遮挡网站| 成年女人看的毛片在线观看| 久久人人爽av亚洲精品天堂 | 老女人水多毛片| 亚洲真实伦在线观看| 王馨瑶露胸无遮挡在线观看| 亚洲av中文字字幕乱码综合| 亚洲精品日韩av片在线观看| 国产精品熟女久久久久浪| 精品少妇黑人巨大在线播放| 亚洲综合色惰| 韩国av在线不卡| 大陆偷拍与自拍| 国产成人福利小说| 欧美 日韩 精品 国产| 精品熟女少妇av免费看| av在线天堂中文字幕| av黄色大香蕉| 免费看a级黄色片| 日韩欧美精品免费久久| 亚洲欧美成人精品一区二区| 国产成人免费观看mmmm| a级毛色黄片| 国产熟女欧美一区二区| 免费大片18禁| av在线app专区| 亚洲精品乱码久久久久久按摩| 六月丁香七月| 午夜免费男女啪啪视频观看| 亚洲av男天堂| 建设人人有责人人尽责人人享有的 | 国产高清国产精品国产三级 | 直男gayav资源| 亚洲av国产av综合av卡| 国产免费视频播放在线视频| 午夜免费观看性视频| 最新中文字幕久久久久| 水蜜桃什么品种好| 免费av不卡在线播放| 国产一区二区亚洲精品在线观看| 内地一区二区视频在线| 久久精品人妻少妇| 国产免费视频播放在线视频| 人妻 亚洲 视频| 亚洲av成人精品一二三区| 国产精品一二三区在线看| 亚洲欧美一区二区三区黑人 | 少妇 在线观看| 在线a可以看的网站| 欧美一区二区亚洲| 又黄又爽又刺激的免费视频.| 欧美日韩视频精品一区| 国产精品嫩草影院av在线观看| 不卡视频在线观看欧美| 免费不卡的大黄色大毛片视频在线观看| 色视频www国产| 成人毛片60女人毛片免费| 国产成人a∨麻豆精品| 18+在线观看网站| 18禁裸乳无遮挡动漫免费视频 | 黄片无遮挡物在线观看| 国产老妇伦熟女老妇高清| 国产精品一区二区三区四区免费观看| 国产黄a三级三级三级人| 国产精品久久久久久久久免| 亚洲av欧美aⅴ国产| 搡老乐熟女国产| 成人综合一区亚洲| 亚洲欧洲国产日韩| 免费看日本二区| 91午夜精品亚洲一区二区三区| 国模一区二区三区四区视频| 人人妻人人爽人人添夜夜欢视频 | 毛片女人毛片| 免费大片黄手机在线观看| 精品久久久久久久久av| 亚洲国产精品国产精品| 亚洲欧美精品专区久久| 亚洲成人久久爱视频| 国产高清不卡午夜福利| 偷拍熟女少妇极品色| 国产色婷婷99| 99热网站在线观看| 如何舔出高潮| 亚洲一区二区三区欧美精品 | 你懂的网址亚洲精品在线观看| 欧美激情国产日韩精品一区| 国产女主播在线喷水免费视频网站| 午夜福利视频精品| 91精品一卡2卡3卡4卡| 少妇的逼好多水| 国产综合懂色| 久久国产乱子免费精品| 久久久久网色| 成人高潮视频无遮挡免费网站| 综合色av麻豆| 免费少妇av软件| 亚洲怡红院男人天堂| 极品教师在线视频| 午夜福利在线在线| 国产精品久久久久久精品电影小说 | 国产精品一二三区在线看| 你懂的网址亚洲精品在线观看| 日日摸夜夜添夜夜爱| 搞女人的毛片| 97在线视频观看| 日本三级黄在线观看| av免费观看日本| 成年免费大片在线观看| 亚洲精品自拍成人| 国产片特级美女逼逼视频| 国产精品.久久久| 青春草亚洲视频在线观看| 成人一区二区视频在线观看| 丰满人妻一区二区三区视频av| 新久久久久国产一级毛片| 国产一区二区在线观看日韩| 天堂俺去俺来也www色官网| 日韩人妻高清精品专区| 啦啦啦中文免费视频观看日本| 少妇的逼水好多| 又粗又硬又长又爽又黄的视频| av在线app专区| 亚洲国产精品专区欧美| 国产成人aa在线观看| 麻豆久久精品国产亚洲av| 亚洲欧美清纯卡通| 亚洲经典国产精华液单| 国精品久久久久久国模美| 国产成人一区二区在线| 久久ye,这里只有精品| 欧美变态另类bdsm刘玥| 国产成人福利小说| 国产日韩欧美在线精品| av免费在线看不卡| 国产大屁股一区二区在线视频| 九九久久精品国产亚洲av麻豆| 欧美极品一区二区三区四区| 黄片无遮挡物在线观看| 亚洲精品成人av观看孕妇| 久久久亚洲精品成人影院| 精品一区二区三区视频在线| 大陆偷拍与自拍| 亚洲高清免费不卡视频| 黄色配什么色好看| 午夜免费鲁丝| 伦精品一区二区三区| 欧美一区二区亚洲| 伦精品一区二区三区| 在线观看美女被高潮喷水网站| 欧美高清成人免费视频www| 国产成人免费观看mmmm| 国产在线一区二区三区精| 97人妻精品一区二区三区麻豆| 成人二区视频| 国产一区二区三区av在线| 国产精品一二三区在线看| 国产午夜精品久久久久久一区二区三区| 免费黄色在线免费观看| 国产爱豆传媒在线观看| 亚洲国产精品成人久久小说| 国产永久视频网站| 最近2019中文字幕mv第一页| 亚洲欧美精品自产自拍| 日韩免费高清中文字幕av| 精品熟女少妇av免费看| 婷婷色综合www| 九九在线视频观看精品| 亚洲欧美清纯卡通| av.在线天堂| 97超碰精品成人国产| 欧美精品国产亚洲| 国产欧美日韩精品一区二区| 国产高清国产精品国产三级 | 国产成人freesex在线| 91在线精品国自产拍蜜月| 免费大片黄手机在线观看| 亚洲国产最新在线播放| 亚洲第一区二区三区不卡| 国产老妇伦熟女老妇高清| 亚洲精品自拍成人| 超碰97精品在线观看| 26uuu在线亚洲综合色| 久久精品国产a三级三级三级| 一本色道久久久久久精品综合| 欧美97在线视频| a级毛片免费高清观看在线播放| 午夜免费观看性视频| av在线天堂中文字幕| 人妻制服诱惑在线中文字幕| 嫩草影院入口| 日本黄大片高清| 免费黄色在线免费观看| 精品人妻视频免费看| 高清毛片免费看| 大码成人一级视频| 亚洲人成网站在线观看播放| 色视频在线一区二区三区| 国产白丝娇喘喷水9色精品| 亚洲欧美清纯卡通| 午夜福利高清视频| 国产69精品久久久久777片| 亚洲av男天堂| 内射极品少妇av片p| 丝袜脚勾引网站| 少妇人妻久久综合中文| 国产精品99久久久久久久久| 麻豆成人午夜福利视频| 九九久久精品国产亚洲av麻豆| 尾随美女入室| 五月天丁香电影| 一区二区av电影网| 国产成人一区二区在线| 中文乱码字字幕精品一区二区三区| 亚洲欧美成人综合另类久久久| 午夜免费观看性视频| 能在线免费看毛片的网站| 亚洲精品日韩av片在线观看| 日韩一区二区视频免费看| 一级毛片我不卡| 五月伊人婷婷丁香| 久久久欧美国产精品| 2021少妇久久久久久久久久久| 最新中文字幕久久久久| 熟女av电影| 国产淫片久久久久久久久| 菩萨蛮人人尽说江南好唐韦庄| 亚洲,欧美,日韩| 亚洲成人久久爱视频| 亚洲一级一片aⅴ在线观看| 岛国毛片在线播放| 亚洲精品第二区| 免费不卡的大黄色大毛片视频在线观看| 欧美激情久久久久久爽电影|