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

    基于免疫信息學(xué)的SARS-CoV-2表位篩選及疫苗分子設(shè)計(jì)分析①

    2022-07-21 07:30:44鄧輝雄王革非李雁嫘宋鑫利谷李銘
    中國(guó)免疫學(xué)雜志 2022年7期
    關(guān)鍵詞:表位抗原氨基酸

    鄧輝雄 王革非 李雁嫘 宋鑫利 谷李銘 李 蕊

    (汕頭大學(xué)醫(yī)學(xué)院微生物學(xué)與免疫學(xué)教研室,病原與免疫學(xué)研究中心,汕頭 515041)

    新型冠狀病毒(SARS-CoV-2)感染引起的新冠肺炎(Corona Virus Disease 2019,COVID-19)疫情在全球200多個(gè)國(guó)家和地區(qū)流行,對(duì)全球公共衛(wèi)生安全造成巨大威脅,給社會(huì)生活及經(jīng)濟(jì)發(fā)展造成巨大損失。我國(guó)與其他國(guó)家的科研人員針對(duì)SARS-CoV-2 的表征、候選抗原和表位的鑒定、動(dòng)物模型的建立、免疫應(yīng)答的表征以及疫苗設(shè)計(jì)研究取得了實(shí)質(zhì)性進(jìn)展[1-2]。截至2020年11月,至少已有5種SARS-CoV-2疫苗進(jìn)入Ⅲ期臨床試驗(yàn)。但SARS-CoV-2 的抗原特性、抗原變異對(duì)宿主的保護(hù)性免疫反應(yīng)的影響有待進(jìn)一步驗(yàn)證[2]。

    利用免疫信息學(xué)篩選、分析鑒定具有免疫保護(hù)性的表位,可有效提高保護(hù)性抗體親和力和細(xì)胞免疫平衡。病毒抗原,包括Spike 蛋白及RBD 區(qū)域存在多種表位。其中部分表位具有免疫保護(hù)作用,能夠誘導(dǎo)中和抗體和抑制病毒復(fù)制的抗體,以及抑制或殺傷感染細(xì)胞及輔助免疫平衡應(yīng)答的效應(yīng)T 細(xì)胞。大多數(shù)表位所誘導(dǎo)的免疫應(yīng)答不具有中和或抑制病毒復(fù)制的作用,反而可能會(huì)誘導(dǎo)細(xì)胞免疫應(yīng)答失衡并加重炎癥反應(yīng),存在免疫病理?yè)p傷、抗體依賴(lài)的感染增強(qiáng)(ADE)等負(fù)面風(fēng)險(xiǎn),甚至可能減弱具有免疫保護(hù)力表位的免疫保護(hù)效應(yīng)[3-4]。

    Spike 蛋白是SARS-COV-2 的關(guān)鍵抗體靶點(diǎn),其同源三聚體在受體結(jié)合和病毒入侵的過(guò)程中發(fā)揮重要作用。全長(zhǎng)三聚體Spike 蛋白通常具有很高的免疫原性(約600 kD),其不僅包含RBD 及強(qiáng)效中和抗體的主要靶點(diǎn),還包含能誘導(dǎo)中和抗體或保護(hù)性抗體的非RDB 區(qū)域,如N 端結(jié)構(gòu)域。目前已報(bào)道了一種基于結(jié)構(gòu)設(shè)計(jì)的疫苗,確保MERS-CoV Spike蛋白穩(wěn)定維持在具有最佳抗病毒效果的融合前構(gòu)象,進(jìn)而提高基于全長(zhǎng)Spike蛋白的疫苗保護(hù)效果,該策略能維持宿主體液免疫和細(xì)胞免疫的有效平衡[5-6]。

    為確保合成肽疫苗效果,本研究基于免疫信息學(xué)進(jìn)行SARS-CoV-2 優(yōu)勢(shì)保護(hù)性抗原表位的篩選、鑒定及對(duì)接分析,計(jì)算機(jī)克隆及免疫模擬預(yù)測(cè)分析抗原表位,構(gòu)建候選多價(jià)表位疫苗,現(xiàn)將結(jié)果報(bào)道如下。

    1 資料與方法

    1.1 SARS-CoV-2候選疫苗毒株的確定 選擇NCBI數(shù)據(jù)庫(kù)發(fā)表的surface glycoprotein[Severe acute respiratory syndrome coronavirus 2](NC_045512.2)為細(xì)胞表位預(yù)測(cè)模板。通過(guò)NCBI(https://www.ncbi.nlm.nih.gov/protein)數(shù)據(jù)庫(kù)收集了579 條序列不同的SARS-COV-2 Spike 蛋白氨基酸序列作為保守性分析模板,所有蛋白數(shù)據(jù)以FASTA文件格式收集。

    1.2 細(xì)胞表位預(yù)測(cè) 通過(guò)ABCpred(http://crdd.osdd. net/raghava/abcpred/)和BCPreds(http://ailab.ist.psu.edu/bcpred/)預(yù)測(cè)候選新冠病毒毒株的線性B 細(xì)胞表位。ABCpred 是基于遞歸神經(jīng)網(wǎng)絡(luò)的學(xué)習(xí)算法,默認(rèn)閾值為0.51,表位長(zhǎng)度為16,其準(zhǔn)確率為65.93%[7]。而基于內(nèi)核方法實(shí)現(xiàn)的BCPreds 準(zhǔn)確率為75.8%[8]。IEDB 服務(wù)器(http://www.iedb.org/)和RANKPEP 服務(wù)器(http://imed. med. ucm. es/Tools/rankpep. html)用來(lái)預(yù)測(cè)MHCⅡ型表位[9-10]。選擇IEDB 默認(rèn)參數(shù)設(shè)置,并選擇一個(gè)MHC HLA 等位基因中人群覆蓋度可達(dá)最小等位基因99%覆蓋率的等位基因參考集[11]。RANKPEP 使用定位特異性評(píng)分矩陣(PSSMs)或已知與給定MHC 分子結(jié)合的一組對(duì)齊肽作為MHC肽結(jié)合的預(yù)測(cè)因子,可重點(diǎn)預(yù)測(cè)保守表位,避免突變引起的免疫逃逸[10],表位與MHCⅡ的結(jié)合潛力通過(guò)EpiTOP(http://www.pharm‐fac.net/EpiTOP/index.php)評(píng)估[12]。IFNepitope(http://crdd. osdd. net/raghava/ifnepitope/)和IL4pred(http://crdd.osdd.net/raghava/il4pred/)預(yù)測(cè)IFN-γ 誘導(dǎo)性抗原表位和非IL-4誘導(dǎo)性抗原表位[13-14]。MHCⅠ結(jié)合CTL表位預(yù)測(cè)分別用NetMHCpan 4.1服務(wù)器(http://www. cbs. dtu. dk/services/NetMHCpan/)、IEDB 服務(wù)器(http://tools. iedb. org/mhci/)、Rankpep 服務(wù)器(http://imed.med.ucm.es/Tools/rankpep.html)、CTL‐pred(http://crdd. osdd. net/raghava/ctlpred/)、NetCTL server 1.2(http://www.cbs.dtu.dk/services/NetCTL/)進(jìn)行分析。選擇的所有服務(wù)器HLA 等位基因的參數(shù)均可達(dá)到97%的群體覆蓋率[15]。NetMHCpan4.1 server 利用人工神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)表位與任何已知序列的MHC分子的結(jié)合[16],本研究只選擇了強(qiáng)結(jié)合能力的抗原表位進(jìn)行進(jìn)一步分析。IEDB 服務(wù)器使用ANN、SMM 和CombLib 組成的共識(shí)推薦方法[17]。Rankpep 使用位置特異性評(píng)分矩陣(PSSM)[10]。CTLpred 利用T細(xì)胞表位模式代替MHC 結(jié)合物預(yù)測(cè)CTL 表位[18]。NetCTL 預(yù)測(cè)方法綜合預(yù)測(cè)MHCⅠ類(lèi)結(jié)合肽、蛋白酶體C端裂解、TAP 轉(zhuǎn)運(yùn)效率。所有預(yù)測(cè)的表位均由IEDB(http://tools.iedb.org/immunoge‐nicity/)進(jìn)行免疫原性檢測(cè),并通過(guò)ToxinPred(https://webs. iiitd. edu. in/raghava/toxinpred/design.php)確定其無(wú)毒性[19-20]。在本研究的表位預(yù)測(cè)策略中,使用多種工具可比單獨(dú)使用一種工具獲得更可靠的結(jié)果。

    1.3 保護(hù)性表位的鑒定 為了選擇合適的表位進(jìn)行進(jìn)一步分析,通過(guò)比較服務(wù)器生成的B細(xì)胞表位、Th 表位和CTL 表位,找出重疊度高的最佳片段,使表位同時(shí)具有誘導(dǎo)體液免疫反應(yīng)和細(xì)胞免疫反應(yīng)的能力。選擇MHCⅠ結(jié)合物,CTL 和B 細(xì)胞表位重疊區(qū)域作為潛在的表位,使用IEDB 服務(wù)器(http://tools. immuneepitope. org/mhcii/)確定最佳結(jié)合物的IEDB評(píng)分[9]。

    1.4 候選抗原表位保守性評(píng)估及T 細(xì)胞(HTL&CTL)表位等位基因人群覆蓋分析 通過(guò)IEDB(http://tools.iedb.org/conservancy/)上的表位保守性分析工具評(píng)估病毒血清型的保守水平[21]。根據(jù)已確定的候選表位與579條不同的SARS-COV-2 Spike蛋白質(zhì)序列比對(duì),計(jì)算抗原表位的保守性情況。分析類(lèi)型和序列識(shí)別閾值分別設(shè)置為“線性”和“100%”,選擇去除相同模板序列。IEDB 人群覆蓋度工具(http://tools. iedb. org/population/)用于評(píng)估預(yù)測(cè)的T細(xì)胞表位(MHCⅠ和MHCⅡ)的人群覆蓋率[22]。

    1.5 疫苗構(gòu)建及抗原性、過(guò)敏原性、溶解度的預(yù)測(cè) 基于上述免疫信息學(xué)分析得到的表位,使用KK、AAY 和GPGPG linker 將B 細(xì)胞表位、HTL 和CTL 表位按順序連接在一起,50S 核糖體蛋白L7/L1分子佐劑通過(guò)EAAAK 在N 端連接[23]。采用Vexigen(http://www. ddg-pharmfac. net/vaxijen/VaxiJen/Vaxi‐Jen.html)和ANTIGENpro(http://scratch.proteomics.ics.uci.edu/)預(yù)測(cè)構(gòu)建的多表位疫苗的抗原性[24-25]。AllerTop server(https://www. ddg-pharmfac. net/Aller‐TOP/)用于預(yù)測(cè)最終疫苗構(gòu)建的過(guò)敏原性,與其他服務(wù)器進(jìn)行的過(guò)敏原預(yù)測(cè)相比,AllerTOPV2.0 的敏感度高達(dá)94%[26]。SOLpro(http://scratch.proteomics.ics.uci.edu/)用于預(yù)測(cè)構(gòu)建的疫苗蛋白的溶解度[27]。利用ExPASy 工具(https://web.expasy.org/protparam/)對(duì)疫苗的分子量、等電點(diǎn)(PI)、半衰期、穩(wěn)定性等理化性質(zhì)進(jìn)行預(yù)測(cè)[28]。

    1.6 疫苗二級(jí)、三級(jí)結(jié)構(gòu)的預(yù)測(cè)及二硫鍵設(shè)計(jì) 以Raptor X(http://raptorx. uchicago. edu/StructPredV2/predict/)和PSIPRED 4.0(http://bioinf.cs.ucl.ac.uk/psipredtest)預(yù)測(cè)蛋白的二級(jí)結(jié)構(gòu)[29-30]。以trRosetta(https://yanglab. nankai. edu. cn/trRosetta/)完成對(duì)最終蛋白疫苗的三級(jí)結(jié)構(gòu)的預(yù)測(cè)[31],trRosetta 是一種快速、準(zhǔn)確地從頭預(yù)測(cè)蛋白質(zhì)結(jié)構(gòu)的算法,其基于一個(gè)限制性的Rosetta 模塊直接能量最小化構(gòu)建蛋白質(zhì)穩(wěn)定構(gòu)象,約束了包括殘基間距離和方向分布,由深度殘基神經(jīng)網(wǎng)絡(luò)進(jìn)行預(yù)測(cè)。在CASP1 和CAMEO 派生集的基準(zhǔn)測(cè)試中,trRosetta 的性能優(yōu)于所有之前描述的方法?;贠mega、Phi、Distance、Theta、Contact 參數(shù)以及TM-score。PyMOL 程序用于三維模型的可視化。在進(jìn)行下一步之前,有必要提高蛋白模型的穩(wěn)定性。利用Disulfide by Design2.0(http://cptweb. cpt. wayne. edu/DbD2/index. php)對(duì)候選蛋白疫苗進(jìn)行二硫鍵設(shè)計(jì)[32]。

    1.7 三級(jí)結(jié)構(gòu)的重修與驗(yàn)證 GalaxyRefine(http://galaxy. seoklab. org/cgi-bin/submit. cgi?type=REFINE)用于重修三級(jí)結(jié)構(gòu)[33]。為了驗(yàn)證三級(jí)結(jié)構(gòu)模型的質(zhì)量,使用RAMPAGE 在線服務(wù)器(http://mordred.bioc. cam. ac. uk/~rapper/rampage. php)和ProSA-web(https://prosa. services. came. sbg. ac. at/prosa. php)對(duì)模型進(jìn)行檢測(cè)[34-35]。

    1.8 B 細(xì)胞構(gòu)象表位的預(yù)測(cè)及蛋白分子對(duì)接分析 利用DiscoTope2.0(http://www. cbs. dtu. dk/ser‐vices/DiscoTope/)預(yù)測(cè)疫苗蛋白三級(jí)結(jié)構(gòu)的構(gòu)象B細(xì)胞表位[36]。蛋白質(zhì)間的相互作用對(duì)于理解細(xì)胞功能和組織結(jié)構(gòu)非常重要。分子對(duì)接是預(yù)測(cè)受體與配體在穩(wěn)定構(gòu)象下相互作用的一種方法。配體和受體間的相互作用取決于他們的親和力,親和力越高,配體和受體分子間的相互作用越穩(wěn)定。使用ClusPro 2.0(https://cluspro. bu. edu/login. php)預(yù)測(cè)作為受體的TLR-4(PDB ID:4G8A)和作為配體的蛋白疫苗間的相互作用[37]。TLR-4(PDB ID:4G8A)受體使用pyMol 去掉水分子及其他雜質(zhì)。該服務(wù)器產(chǎn)生的對(duì)接配合物具有良好的脫溶自由能和靜電相互作用。LigPLot+用來(lái)確定疫苗配體和TLR-4 受體是否具有氫鍵和疏水相互作用的氨基酸殘基。

    1.9 計(jì)算機(jī)克隆及免疫模擬分析 計(jì)算機(jī)模擬克隆通過(guò)SnapGene 軟件完成,而多表位疫苗的免疫原性和免疫應(yīng)答模式通過(guò)C-ImmSim 服務(wù)器(http://150.146.2.1/C-IMMSIM/index.php)完成,該服務(wù)器基于位置特異性評(píng)分矩陣(PSSM)原理,預(yù)測(cè)候選疫苗模型抗原表位,并基于經(jīng)典免疫學(xué)原理使用機(jī)器學(xué)習(xí)模型預(yù)測(cè)免疫應(yīng)答相互作用[38]。在本實(shí)驗(yàn)中,基于群體覆蓋率使用了以下MHC 等位基因:HLA A*11:01、HLA A*31:01、HLA B*35:01、HLA B*53:01、HLA DRB5*01:01、HLA DRB1*09:01。按典型免疫過(guò)程間隔1 周注射,共免疫3 次。所有模擬參數(shù)均設(shè)置為默認(rèn)值,時(shí)間步長(zhǎng)設(shè)置為1、20 和40(每個(gè)時(shí)間步長(zhǎng)為8 h,時(shí)間步長(zhǎng)1 為第1 次免疫時(shí)間點(diǎn)),抗原分子的氨基酸序列作為疫苗。

    2 結(jié)果

    2.1 SARS-CoV-2 表位的初步篩選 表1 為經(jīng)鑒定的5 個(gè)B 細(xì)胞表位,表2 為6 個(gè)Th 表位,表3 為5 個(gè)CTL 表位。所有候選抗原表位在Spike 蛋白的構(gòu)象位置如圖1A,獲得能夠中和Spike蛋白與ACE2結(jié)合域的中和抗體具有重要意義,本研究選擇的B-1 和B-2 候選抗原表位如圖1B,在Spike 蛋白R(shí)BD 和ACE2 結(jié)合域上,具有能夠誘導(dǎo)主要中和抗體、破壞其與ACE2結(jié)合的潛力。候選Th抗原表位和CTL抗原表位結(jié)合的HLA 等位基因?qū)?yīng)的不同地區(qū)人群覆蓋度分析結(jié)果如表4 所示,候選表位可覆蓋全世界95.4%的人群,盡管對(duì)應(yīng)不同地區(qū)人群有差異,但依然具有較高的人群覆蓋度,如目前新冠疫情的高發(fā)地區(qū)南亞有98.39%的人群覆蓋度,歐洲有99.67%的人群覆蓋度,南美有99.95%的人群覆蓋度。對(duì)于全球性暴發(fā)流行的新冠病毒,開(kāi)發(fā)能夠適應(yīng)高人群覆蓋度的疫苗具有重要意義。

    表4 候選Th和CTL表位等位基因中不同地區(qū)人群覆蓋率Tab.4 Population coverage calculation of Th and CTL epitopes specific alleles

    圖1 SARS-CoV-2 Spike蛋白候選抗原表位空間展示Fig.1 Epitopes of SARS-CoV-2 Spike protein candidate antigen

    表1 候選B細(xì)胞表位Tab.1 Candidate B cell epitopes

    表2 候選Th細(xì)胞表位Tab.2 Candidate Th cell epitopes

    表3 候選CTL細(xì)胞表位Tab.3 Candidate CTL cell epitopes

    2.2 SARS-CoV-2表位疫苗構(gòu)建 基于上述免疫信息學(xué)分析得到的表位,使用KK、AAY 和GPGPG linker將B 細(xì)胞表位、HTL 和CTL 表位按順序連接在一起,B 細(xì)胞表位之間使用雙賴(lài)氨酸(KK)linker 保持其獨(dú)立的免疫原性活性,MHC 分子在MHCⅡ分子介導(dǎo)的抗原提呈到Th0之前,作用于抗原的KK序列,使用AAY 和GPGPG linker 可增強(qiáng)疫苗亞基的識(shí)別能力。Linker 可確保良好的穩(wěn)定性、結(jié)構(gòu)域取向和折疊率,使用EAAAK linker 連接表位和50S L7/L12分子佐劑氨基酸序列,示意圖如圖2所示。

    圖2 SARS-CoV-2表位疫苗氨基酸序列Fig.2 Amino acid residues sequence of SARS-CoV-2 epitopes-based vaccine

    2.3 SARS-CoV-2表位疫苗理化性質(zhì)分析 疫苗全長(zhǎng)392 個(gè)氨基酸,分子量大小為41.649 92 kD,pI:9.35,在體外哺乳動(dòng)物網(wǎng)織紅細(xì)胞中半衰期為30 h,體內(nèi)酵母細(xì)胞中大于20 h,大腸桿菌中大于10 h,抗原穩(wěn)定指數(shù)為21.82,歸類(lèi)為穩(wěn)定蛋白質(zhì)。Vaxijen v2.0 預(yù)測(cè)保護(hù)性抗原在設(shè)定閾值為0.4 時(shí)分?jǐn)?shù)為0.416 2,可作為抗原;ANTIGENpro 預(yù)測(cè)抗原免疫原性為0.964 894;預(yù)測(cè)過(guò)表達(dá)溶解度為0.942 995,提示可溶。脂肪族指數(shù)為84.52,Grand average of hydropathicity(GRAVY):?0.117,球狀蛋白結(jié)構(gòu)中的脂肪族鏈代表了該蛋白的熱穩(wěn)定性。蛋白質(zhì)脂肪族指數(shù)的升高表明該蛋白質(zhì)在性質(zhì)上是熱穩(wěn)定的,GRAVY代表蛋白質(zhì)序列疏水性的平均值。

    2.4 SARS-CoV-2表位疫苗二級(jí)結(jié)構(gòu)預(yù)測(cè) PSIPRED和Raptox 對(duì)二級(jí)結(jié)構(gòu)的預(yù)測(cè)結(jié)果表明,蛋白疫苗由11.1%的螺旋、35.3%的折疊和53.6%的卷曲組成。溶劑可及性(ACC)由26%的埋藏、25%的居中和48%的外側(cè)組成。二級(jí)結(jié)構(gòu)預(yù)測(cè)如圖3所示。

    圖3 SARS-CoV-2表位疫苗二級(jí)結(jié)構(gòu)預(yù)測(cè)分析Fig.3 Secondary structure prediction of SARS-CoV-2 epitopes-based vaccine

    2.5 疫苗三級(jí)結(jié)構(gòu)建模 最佳疫苗三級(jí)結(jié)構(gòu)建模結(jié)果如圖4F 所示。進(jìn)一步使用Disulfide by De‐sign2.0進(jìn)行二硫鍵設(shè)計(jì),在默認(rèn)參數(shù)下,分析得到的B-factor 均為0,認(rèn)為通過(guò)trRosetta 能量最小化的蛋白質(zhì)建模已足夠穩(wěn)定,不需要進(jìn)一步的二硫鍵設(shè)計(jì)。預(yù)測(cè)模型的Ramachandran 質(zhì)量分析表明,預(yù)測(cè)模型中分別有379 個(gè)(97.2%)、8 個(gè)(2.1%)和3 個(gè)氨基酸(0.8%)的氨基酸殘基分別位于有利、允許和異常區(qū)域(圖5)。

    圖4 SARS-CoV-2疫苗三級(jí)結(jié)構(gòu)建模Fig.4 SARS-CoV-2 epitopes-based vaccine tertiary struc?ture modeling

    圖5 疫苗三級(jí)結(jié)構(gòu)模型refine 前Ramachandran 分析模型質(zhì)量Fig.5 Validation of vaccine tertiary structure model before refine using Ramachandran plot

    2.6 疫苗三級(jí)結(jié)構(gòu)的refine 和驗(yàn)證 應(yīng)用Galaxy‐Refine 程序?qū)ν葱阅P徒Y(jié)構(gòu)進(jìn)行了改進(jìn)重修。GalaxyRefine 為模型構(gòu)建了5 個(gè)重修模型。重修前蛋白質(zhì)模型各種參數(shù)為GDT-HA(1.000 0)、RMSD(0.000)、Molprobity(2.509)、Clash score(48.7)、poor rotamers(2.0)以及Ramachandran plot(97.2),重修后確認(rèn)模型1 為最佳候選模型,其具體參數(shù)為GDT-HA(0.9458)、RMSD(0.440)、Molprobity(1.995)、Clash score(19.4)、poor rotamers(0.7)以及Ramachandran plot(97.9)。使用Pymol 軟件繪制的卡通圖比對(duì)結(jié)果如圖6A 所示。Ramachandran 對(duì)重修三級(jí)結(jié)構(gòu)模型的驗(yàn)證表明,重修后的模型質(zhì)量獲得提升(圖6B)。Ramachandran 圖表明,在優(yōu)化后,有利區(qū)域含382個(gè)氨基酸(97.9%),允許區(qū)域含6個(gè)氨基酸(1.5%),不允許區(qū)域含2 個(gè)氨基酸(0.5%)。而ProSA 網(wǎng)絡(luò)(蛋白質(zhì)結(jié)構(gòu)分析)顯示,蛋白質(zhì)疫苗結(jié)構(gòu)準(zhǔn)確度分析結(jié)果在類(lèi)似大小的天然蛋白質(zhì)常見(jiàn)的分?jǐn)?shù)范圍內(nèi),z值為?7.07(圖6C)。該圖通過(guò)繪制氨基酸序列位置函數(shù)的能量分布來(lái)顯示局部模型質(zhì)量(圖6D),一般來(lái)說(shuō),正值對(duì)應(yīng)輸入結(jié)構(gòu)存在問(wèn)題部分。結(jié)果提示,基于同源性建模的多價(jià)抗原表位腸道病毒蛋白質(zhì)結(jié)構(gòu)模型精度可滿(mǎn)足進(jìn)一步分析。

    圖6 候選疫苗三級(jí)結(jié)構(gòu)的重修和驗(yàn)證Fig.6 Refinement and validation of vaccine tertiary struc?ture model

    2.7 計(jì)算機(jī)模擬克隆候選多價(jià)表位疫苗及構(gòu)象B細(xì)胞表位分析 計(jì)算機(jī)模擬克隆選擇pGEX-KG 原核表達(dá)載體,使用Snapgene 軟件對(duì)疫苗的氨基酸序列進(jìn)行反向翻譯得到疫苗表達(dá)基因序列,通過(guò)密碼子的優(yōu)化獲得更好的表達(dá)。通過(guò)選擇BamHⅠ和HindⅢ限制性?xún)?nèi)切酶酶切位點(diǎn)構(gòu)建pGEX-S 新冠疫苗原核表達(dá)載體,示意圖如圖7A所示。候選疫苗構(gòu)象B 細(xì)胞表位分析surface 結(jié)構(gòu)如圖7B 所示,共有44個(gè)構(gòu)象B細(xì)胞表位。

    圖7 原核表達(dá)載體構(gòu)建和構(gòu)象B細(xì)胞表位分析Fig.7 Diagram of silico cloning and conformational B-cell epitope prediction

    2.8 疫苗與TLR4-MD2 復(fù)合體對(duì)接 在平衡模式下選擇model.000.05 作為最佳對(duì)接模型,對(duì)接結(jié)果如圖8A 所示,在該對(duì)接模式下可以獲得穩(wěn)定的結(jié)合。LigPLot+確定疫苗配體和TLR-4/MD2 受體具有氫鍵的氨基酸殘基為Asn365、Asn181、Gly363、Arg184、Asn210、Lys170、Tyr371、His350、Gly391、Tyr383、Thr389、Glu33、Arg368。具有疏水相互作用的氨基酸殘基為Ile366、Ile362、Ile364、Gly213、Trp374、Val353、Ala357、Phe375、Ala392、Phe390、Trp388、His379、Tyr371、Thr35、Phe375、His376。

    圖8 候選疫苗蛋白與人TLR4/MD2受體對(duì)接復(fù)合體結(jié)果Fig.8 Complexes of candidated vaccine protein docking with human TLR4/MD2 receptor

    2.9 計(jì)算機(jī)模擬免疫應(yīng)答刺激 多表位疫苗的免疫原性和免疫應(yīng)答由C-ImmSim 服務(wù)器進(jìn)行,按典型免疫過(guò)程間隔1 周注射,共免疫3 次。所有模擬參數(shù)均設(shè)置為默認(rèn)值,時(shí)間步長(zhǎng)設(shè)置為1、20 和40(每個(gè)時(shí)間步長(zhǎng)為8 h,時(shí)間步長(zhǎng)1為第1次免疫時(shí)間點(diǎn)),抗原分子的氨基酸序列作為疫苗。通過(guò)計(jì)算機(jī)模擬免疫應(yīng)答刺激可快速評(píng)估疫苗的免疫原性,快速篩選候選疫苗,節(jié)約時(shí)間和成本。計(jì)算機(jī)模擬免疫應(yīng)答刺激可以確定疫苗符合誘導(dǎo)免疫應(yīng)答規(guī)律,具有良好的免疫原性,可有效激活體液免疫和細(xì)胞免疫,結(jié)果如圖9~12所示。

    圖9 模擬免疫刺激免疫B細(xì)胞和Th數(shù)量變化展示Fig.9 B cell and Th counts shown in silico immune simu?lation

    圖10 模擬免疫刺激免疫CTL、NK、MA、DC、EP 細(xì)胞數(shù)量變化展示Fig.10 CTL,NK,MA,DC,EP cell counts shown in silico immune simulation

    圖11 免疫球蛋白和免疫復(fù)合物變化趨勢(shì)Fig.11 Trends of immunoglobulins and immunocomplexes

    圖12 細(xì)胞因子和白細(xì)胞介素的濃度變化趨勢(shì)Fig.12 Concentration trends of cytokines and interleukins

    3 討論

    借助免疫信息學(xué)和結(jié)構(gòu)生物學(xué)開(kāi)發(fā)新冠病毒疫苗是一種有效的策略,目前已有一些借助免疫信息學(xué)和結(jié)構(gòu)生物學(xué)設(shè)計(jì)的新冠疫苗研究和本文類(lèi)似[39-40],不同的是,本研究策略側(cè)重于體液免疫和細(xì)胞免疫的平衡。疫苗的保護(hù)作用除了能夠誘導(dǎo)中和抗體外,還需依賴(lài)于細(xì)胞毒性CD8 T 細(xì)胞和輔助性CD4 T細(xì)胞來(lái)完全清除病毒。免疫應(yīng)答的重要特征之一就是輔助性T細(xì)胞亞型的激活及相對(duì)應(yīng)的特異性類(lèi)型細(xì)胞因子的釋放。加載在MHCⅡ上的抗原表位可激活Th1 和Th2 型輔助細(xì)胞。Th2 主要引發(fā)針對(duì)細(xì)胞外病原體的體液反應(yīng),而Th1 激活針對(duì)細(xì)胞內(nèi)病原體(病毒、癌癥)的細(xì)胞免疫。然而,Th1和Th2 分別與介導(dǎo)細(xì)胞免疫和體液免疫不嚴(yán)格相等[41],例如Th1 途徑也可刺激適度水平的基于抗體的反應(yīng)。Th1細(xì)胞因子傾向于產(chǎn)生促炎反應(yīng),而Th2與抗炎反應(yīng)相關(guān)。不平衡的Th1/Th2 反應(yīng)可能引起免疫病理學(xué)并發(fā)癥,例如廣泛的組織損傷炎癥或強(qiáng)烈的過(guò)敏反應(yīng)。因此,在疫苗開(kāi)發(fā)過(guò)程中應(yīng)考慮適當(dāng)平衡Th1和Th2反應(yīng)。Th與APC相互作用可依靠IL-4 這樣的細(xì)胞因子形成局部正反饋,特異性的激活遞呈了有效表位的APC,同時(shí)也擴(kuò)增了能夠識(shí)別有效表位的Th,避免泥沙俱下的胡亂遞呈[14]。針對(duì)每個(gè)不同表位的Th,可以確立反應(yīng)方式,是奮力殺死敵人,還是選擇免疫耐受,與之和平相處。理想抗原應(yīng)由一般人群中的多個(gè)MHCⅠ和MHCⅡ等位基因呈遞,并包含與中和抗體相關(guān)的線性B 細(xì)胞表位。很多抗原并沒(méi)有足夠有效的抗原表位,無(wú)法有效激活細(xì)胞毒性T細(xì)胞,所以效率低,準(zhǔn)確找到有效的抗原表位具有重要意義。本研究通過(guò)免疫信息學(xué)篩選出具有平衡體液免疫和細(xì)胞免疫特征的新冠病毒候選B 細(xì)胞、Th 和CTL 細(xì)胞表位,設(shè)計(jì)具有多表位的高免疫原性候選疫苗。通過(guò)計(jì)算機(jī)免疫模擬刺激實(shí)驗(yàn)分析發(fā)現(xiàn)候選疫苗可誘導(dǎo)CD8+T 細(xì)胞,同時(shí)產(chǎn)生CD4+T 輔助性1型細(xì)胞免疫應(yīng)答,其具有良好的抗病毒特性。佐劑不僅可用于提高蛋白質(zhì)疫苗引起的抗體反應(yīng)的強(qiáng)度和持久性,還可影響T 細(xì)胞來(lái)源的細(xì)胞因子模式,從而調(diào)節(jié)免疫反應(yīng),50S L7/L12 分子佐劑已被證明可通過(guò)激活TLR4 信號(hào)通路增強(qiáng)免疫應(yīng)答,實(shí)驗(yàn)通過(guò)分子對(duì)接確定候選疫苗與TLR4 的結(jié)合模式和相互作用的氨基酸殘基。在設(shè)計(jì)疫苗時(shí)需要在激發(fā)最大的免疫原性的同時(shí)盡可能地減少潛在的副作用。冠狀病毒株之間的多樣性為我們提出了另一個(gè)挑戰(zhàn),因?yàn)樵诠跔畈《緦僦兄辽儆? 個(gè)不同的亞組。盡管SARS-CoV和SARS-CoV-2 處在相同的亞組,但二者的基因組只有77%相同[42]。了解SARS-CoV-2 的共有和獨(dú)特表位有助于為患者設(shè)計(jì)更好的疫苗或開(kāi)展診斷性臨床試驗(yàn)。另一方面,新冠病毒疫苗可能存在的抗體依賴(lài)增強(qiáng)(antibody-dependent enhancement,ADE)效應(yīng)也是一個(gè)巨大的挑戰(zhàn),傳統(tǒng)的ADE 效應(yīng)主要認(rèn)為抗體不能中和病毒,反而充當(dāng)了一個(gè)特洛伊木馬,讓病毒更容易感染單核細(xì)胞。決定抗體是否引起ADE 的因素主要包括抗體的特異性、滴度、親和力以及抗體的亞型[43]。在猴子模型中,針對(duì)SARSCoV 的刺突蛋白的不同表位的抗體,其誘導(dǎo)的反應(yīng)也各不相同,有些可起到很好的保護(hù)作用,有些則容易引起ADE 效應(yīng)[44]。使用有效的表位設(shè)計(jì)疫苗可認(rèn)為是有效避免ADE效應(yīng)的策略。

    4 結(jié)論

    通過(guò)免疫信息學(xué)和結(jié)構(gòu)生物學(xué)分析,進(jìn)行基于抗原表位的新冠病毒分子疫苗設(shè)計(jì),信息學(xué)分析結(jié)果初步顯示,具備高保守性、免疫原性的B 細(xì)胞表位、MHC 等位基因高人群覆蓋度的IL-4、IFN-γ 誘導(dǎo)型Th 表位及CTL 候選表位構(gòu)建的疫苗具有良好的平衡體液免疫和細(xì)胞免疫應(yīng)答能力。

    猜你喜歡
    表位抗原氨基酸
    H3N2流感病毒HA保守Th表位對(duì)CD4+T細(xì)胞活化及分化的影響
    月桂酰丙氨基酸鈉的抑菌性能研究
    UFLC-QTRAP-MS/MS法同時(shí)測(cè)定絞股藍(lán)中11種氨基酸
    中成藥(2018年1期)2018-02-02 07:20:05
    聯(lián)合T、B細(xì)胞表位設(shè)計(jì)多肽疫苗的研究進(jìn)展①
    一株Nsp2蛋白自然缺失123個(gè)氨基酸的PRRSV分離和鑒定
    梅毒螺旋體TpN17抗原的表達(dá)及純化
    小反芻獸疫病毒化學(xué)合成表位多肽對(duì)小鼠的免疫效果研究
    結(jié)核分枝桿菌抗原Lppx和MT0322人T細(xì)胞抗原表位的多態(tài)性研究
    APOBEC-3F和APOBEC-3G與乙肝核心抗原的相互作用研究
    氨基酸分析儀測(cè)定玉米漿中17種游離氨基酸的不確定度評(píng)定
    tube8黄色片| 免费在线观看完整版高清| 亚洲熟女精品中文字幕| 欧美成人精品欧美一级黄| svipshipincom国产片| 久久99精品国语久久久| 18禁观看日本| 国产成人欧美在线观看 | 爱豆传媒免费全集在线观看| 尾随美女入室| 少妇猛男粗大的猛烈进出视频| 中文字幕人妻丝袜一区二区 | 国产av精品麻豆| 久久精品熟女亚洲av麻豆精品| 久久精品国产综合久久久| 成人黄色视频免费在线看| 久久久欧美国产精品| 少妇人妻精品综合一区二区| 激情五月婷婷亚洲| 国产精品av久久久久免费| 国产日韩一区二区三区精品不卡| 亚洲国产av影院在线观看| 熟女少妇亚洲综合色aaa.| 一级黄片播放器| 国产av国产精品国产| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美黑人欧美精品刺激| 夜夜骑夜夜射夜夜干| 亚洲欧洲精品一区二区精品久久久 | 亚洲 欧美一区二区三区| 一级片'在线观看视频| 母亲3免费完整高清在线观看| 国产欧美亚洲国产| 999精品在线视频| 在线观看免费日韩欧美大片| 黄频高清免费视频| 国产精品.久久久| 国产成人午夜福利电影在线观看| 日本91视频免费播放| 热99国产精品久久久久久7| 亚洲欧洲精品一区二区精品久久久 | 日韩电影二区| 亚洲精品在线美女| 欧美日韩国产mv在线观看视频| 黑人欧美特级aaaaaa片| 黄片无遮挡物在线观看| 国产又色又爽无遮挡免| 女人被躁到高潮嗷嗷叫费观| 9191精品国产免费久久| 日韩不卡一区二区三区视频在线| 久久久久视频综合| 中文字幕人妻熟女乱码| 99精国产麻豆久久婷婷| av视频免费观看在线观看| 十八禁高潮呻吟视频| 国产欧美日韩一区二区三区在线| 免费观看性生交大片5| 久久久久视频综合| 婷婷色综合大香蕉| 久久免费观看电影| 男女之事视频高清在线观看 | 久久这里只有精品19| 国产女主播在线喷水免费视频网站| 欧美人与善性xxx| 午夜福利视频在线观看免费| 国产精品国产三级国产专区5o| 99久久综合免费| 欧美中文综合在线视频| 国产精品秋霞免费鲁丝片| 交换朋友夫妻互换小说| 国产有黄有色有爽视频| 欧美黑人欧美精品刺激| 黑人巨大精品欧美一区二区蜜桃| 日本欧美国产在线视频| 丝袜美足系列| 丰满迷人的少妇在线观看| 美女主播在线视频| 成人午夜精彩视频在线观看| 国产 精品1| 成人影院久久| 成人黄色视频免费在线看| 少妇被粗大的猛进出69影院| 丝袜在线中文字幕| 久久精品国产综合久久久| av女优亚洲男人天堂| 国产精品欧美亚洲77777| 日韩中文字幕视频在线看片| 女人精品久久久久毛片| 亚洲av成人精品一二三区| 国产欧美日韩综合在线一区二区| 女人久久www免费人成看片| 亚洲成人免费av在线播放| 中文字幕制服av| 女性被躁到高潮视频| 我的亚洲天堂| 国精品久久久久久国模美| 欧美日韩成人在线一区二区| 亚洲一区二区三区欧美精品| 日韩av不卡免费在线播放| 国产女主播在线喷水免费视频网站| 中文字幕人妻丝袜制服| 男女之事视频高清在线观看 | 国产精品秋霞免费鲁丝片| 日韩大码丰满熟妇| 波野结衣二区三区在线| 亚洲,欧美精品.| 欧美在线黄色| 黑人巨大精品欧美一区二区蜜桃| 国产精品久久久久久久久免| 咕卡用的链子| 久久青草综合色| 最近最新中文字幕免费大全7| 日本色播在线视频| 午夜日本视频在线| 国产精品.久久久| 性高湖久久久久久久久免费观看| 欧美变态另类bdsm刘玥| av视频免费观看在线观看| 亚洲少妇的诱惑av| 老司机影院成人| 亚洲 欧美一区二区三区| 女人精品久久久久毛片| 日本一区二区免费在线视频| 精品一区二区三区四区五区乱码 | 丁香六月天网| 国产精品亚洲av一区麻豆 | 精品酒店卫生间| 美女午夜性视频免费| 久久久久久久大尺度免费视频| 九色亚洲精品在线播放| 黄色 视频免费看| 人人妻人人添人人爽欧美一区卜| 丰满迷人的少妇在线观看| 久久久久国产一级毛片高清牌| 成人国语在线视频| www.精华液| 美国免费a级毛片| 悠悠久久av| 午夜福利视频在线观看免费| 精品人妻一区二区三区麻豆| 国产日韩欧美在线精品| 一区二区日韩欧美中文字幕| 19禁男女啪啪无遮挡网站| 中文字幕另类日韩欧美亚洲嫩草| h视频一区二区三区| 少妇 在线观看| 黄色毛片三级朝国网站| 欧美国产精品va在线观看不卡| 99久久99久久久精品蜜桃| 一级毛片 在线播放| 国产精品偷伦视频观看了| 日本色播在线视频| 大香蕉久久成人网| 99热国产这里只有精品6| 最新的欧美精品一区二区| 1024视频免费在线观看| 国产黄色免费在线视频| 在线观看免费日韩欧美大片| 最黄视频免费看| 亚洲国产日韩一区二区| 亚洲四区av| e午夜精品久久久久久久| 国产精品av久久久久免费| 亚洲精品国产av蜜桃| 国产精品女同一区二区软件| 国产精品久久久久久人妻精品电影 | svipshipincom国产片| 亚洲伊人久久精品综合| 黑丝袜美女国产一区| www.自偷自拍.com| 啦啦啦中文免费视频观看日本| 日韩av不卡免费在线播放| 亚洲成av片中文字幕在线观看| 乱人伦中国视频| 男女午夜视频在线观看| 久久精品国产a三级三级三级| 熟妇人妻不卡中文字幕| videos熟女内射| 一本一本久久a久久精品综合妖精| 九草在线视频观看| 伦理电影大哥的女人| 捣出白浆h1v1| www.av在线官网国产| 性高湖久久久久久久久免费观看| 五月开心婷婷网| 久热这里只有精品99| 久久久久久免费高清国产稀缺| 国产乱来视频区| 成年av动漫网址| 高清在线视频一区二区三区| 国产一级毛片在线| 天天躁狠狠躁夜夜躁狠狠躁| 免费黄色在线免费观看| 亚洲美女搞黄在线观看| 在线观看人妻少妇| 国产亚洲av片在线观看秒播厂| 一区福利在线观看| 日本黄色日本黄色录像| 亚洲国产看品久久| 五月开心婷婷网| 亚洲国产欧美一区二区综合| 欧美日韩一级在线毛片| 少妇人妻 视频| 制服诱惑二区| 精品国产乱码久久久久久小说| 亚洲一级一片aⅴ在线观看| 国产精品三级大全| 91aial.com中文字幕在线观看| 亚洲一区二区三区欧美精品| 一区在线观看完整版| 免费观看a级毛片全部| 国产麻豆69| 亚洲av福利一区| 亚洲国产精品一区三区| 日韩欧美精品免费久久| 国产精品99久久99久久久不卡 | 91精品三级在线观看| 亚洲精品美女久久久久99蜜臀 | 免费黄色在线免费观看| 日本欧美国产在线视频| 一二三四中文在线观看免费高清| 99九九在线精品视频| 国产日韩一区二区三区精品不卡| 亚洲情色 制服丝袜| 亚洲 欧美一区二区三区| 亚洲av电影在线观看一区二区三区| 精品一区二区三区av网在线观看 | 国产免费现黄频在线看| 你懂的网址亚洲精品在线观看| av.在线天堂| 伊人亚洲综合成人网| 欧美日韩视频精品一区| 极品人妻少妇av视频| 亚洲欧洲精品一区二区精品久久久 | 别揉我奶头~嗯~啊~动态视频 | 国产成人欧美| 极品少妇高潮喷水抽搐| 国产精品香港三级国产av潘金莲 | 久久狼人影院| 19禁男女啪啪无遮挡网站| 肉色欧美久久久久久久蜜桃| 亚洲av欧美aⅴ国产| 69精品国产乱码久久久| 一本色道久久久久久精品综合| 丰满乱子伦码专区| 啦啦啦在线免费观看视频4| 99九九在线精品视频| www日本在线高清视频| 欧美日韩精品网址| 成年av动漫网址| 男人添女人高潮全过程视频| 国产又色又爽无遮挡免| 日韩熟女老妇一区二区性免费视频| 51午夜福利影视在线观看| 丝袜人妻中文字幕| 别揉我奶头~嗯~啊~动态视频 | 一二三四在线观看免费中文在| 亚洲情色 制服丝袜| 亚洲综合精品二区| 老熟女久久久| 国产午夜精品一二区理论片| 看免费成人av毛片| av.在线天堂| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美精品av麻豆av| 欧美日韩视频高清一区二区三区二| 日日爽夜夜爽网站| 午夜日韩欧美国产| 自线自在国产av| 天天影视国产精品| 人人妻人人澡人人爽人人夜夜| 天天躁夜夜躁狠狠久久av| 色婷婷av一区二区三区视频| www.熟女人妻精品国产| 咕卡用的链子| 丝袜在线中文字幕| 日韩一区二区视频免费看| 欧美日本中文国产一区发布| 国产精品av久久久久免费| 又大又爽又粗| 国产麻豆69| 精品少妇久久久久久888优播| 美国免费a级毛片| 女性被躁到高潮视频| 亚洲欧美激情在线| 久久精品国产a三级三级三级| 中文字幕亚洲精品专区| 午夜福利乱码中文字幕| 久久女婷五月综合色啪小说| 亚洲激情五月婷婷啪啪| 久久久久视频综合| 91精品国产国语对白视频| 国产精品蜜桃在线观看| 亚洲欧美色中文字幕在线| 十八禁网站网址无遮挡| 国产伦人伦偷精品视频| 又大又爽又粗| 精品一区在线观看国产| 国产精品久久久久久精品电影小说| 国产精品久久久av美女十八| 亚洲人成网站在线观看播放| 久久精品人人爽人人爽视色| 久久免费观看电影| 日韩不卡一区二区三区视频在线| 国产成人欧美在线观看 | 无限看片的www在线观看| 日韩中文字幕欧美一区二区 | 飞空精品影院首页| 国产av一区二区精品久久| 国产高清不卡午夜福利| 午夜日本视频在线| 人人妻人人澡人人看| 制服人妻中文乱码| 亚洲精品国产av成人精品| 亚洲精品国产区一区二| 欧美av亚洲av综合av国产av | 女人精品久久久久毛片| 久久国产亚洲av麻豆专区| 超碰97精品在线观看| 80岁老熟妇乱子伦牲交| 精品人妻熟女毛片av久久网站| 中文天堂在线官网| 亚洲人成电影观看| 久久久久视频综合| 国产成人免费无遮挡视频| 97人妻天天添夜夜摸| 国产亚洲av片在线观看秒播厂| 亚洲,欧美,日韩| 人成视频在线观看免费观看| 亚洲国产av影院在线观看| 欧美精品av麻豆av| 亚洲自偷自拍图片 自拍| 999精品在线视频| avwww免费| 大陆偷拍与自拍| 美女主播在线视频| svipshipincom国产片| 欧美日韩国产mv在线观看视频| 性色av一级| 国产一区有黄有色的免费视频| 欧美日本中文国产一区发布| 交换朋友夫妻互换小说| 少妇人妻久久综合中文| 国产99久久九九免费精品| 一区二区三区激情视频| 国产xxxxx性猛交| 国产精品二区激情视频| 亚洲色图 男人天堂 中文字幕| 欧美精品av麻豆av| 亚洲精品日本国产第一区| 午夜精品国产一区二区电影| a级片在线免费高清观看视频| 1024视频免费在线观看| 亚洲国产欧美在线一区| 国产成人免费无遮挡视频| 亚洲精品av麻豆狂野| 美女福利国产在线| 丁香六月欧美| 精品少妇一区二区三区视频日本电影 | 日韩人妻精品一区2区三区| 国产精品国产三级专区第一集| 成人三级做爰电影| 老司机靠b影院| 宅男免费午夜| 丝袜人妻中文字幕| 自拍欧美九色日韩亚洲蝌蚪91| 精品久久久精品久久久| 中文字幕制服av| 欧美日本中文国产一区发布| 国产熟女午夜一区二区三区| 久久影院123| 最近2019中文字幕mv第一页| 青春草视频在线免费观看| 精品免费久久久久久久清纯 | 免费日韩欧美在线观看| 男女午夜视频在线观看| 色94色欧美一区二区| 男女边摸边吃奶| 欧美少妇被猛烈插入视频| 欧美日韩av久久| videos熟女内射| 欧美在线一区亚洲| 日韩制服骚丝袜av| 校园人妻丝袜中文字幕| 亚洲精品成人av观看孕妇| 熟女av电影| 国产福利在线免费观看视频| 免费人妻精品一区二区三区视频| 成人影院久久| 少妇 在线观看| 少妇被粗大的猛进出69影院| 1024视频免费在线观看| 伦理电影免费视频| 亚洲欧美激情在线| 亚洲人成网站在线观看播放| 街头女战士在线观看网站| 欧美日韩亚洲国产一区二区在线观看 | 18在线观看网站| 成人漫画全彩无遮挡| 黄网站色视频无遮挡免费观看| 在线精品无人区一区二区三| 少妇被粗大猛烈的视频| 午夜91福利影院| 日本91视频免费播放| 你懂的网址亚洲精品在线观看| 一本—道久久a久久精品蜜桃钙片| 午夜福利影视在线免费观看| 中文字幕高清在线视频| 男女高潮啪啪啪动态图| 老司机亚洲免费影院| 十八禁高潮呻吟视频| 午夜福利在线免费观看网站| 久久久久精品久久久久真实原创| 亚洲色图综合在线观看| 晚上一个人看的免费电影| 日本黄色日本黄色录像| 操出白浆在线播放| 毛片一级片免费看久久久久| 国产亚洲精品第一综合不卡| 国产免费福利视频在线观看| 成人国语在线视频| 欧美 亚洲 国产 日韩一| 亚洲第一av免费看| 亚洲精品aⅴ在线观看| 中文精品一卡2卡3卡4更新| 丝袜脚勾引网站| 亚洲成av片中文字幕在线观看| 亚洲成国产人片在线观看| 久久久久精品性色| 精品一区二区三卡| 在线看a的网站| 别揉我奶头~嗯~啊~动态视频 | 老司机靠b影院| 国产免费现黄频在线看| 亚洲精品乱久久久久久| 国产片特级美女逼逼视频| 亚洲久久久国产精品| 国产欧美日韩综合在线一区二区| 亚洲欧美一区二区三区久久| 日本欧美视频一区| 亚洲国产精品一区二区三区在线| 亚洲精品国产色婷婷电影| 999精品在线视频| 性高湖久久久久久久久免费观看| 国产av一区二区精品久久| 青春草国产在线视频| 久久毛片免费看一区二区三区| 99热国产这里只有精品6| 国产 精品1| 亚洲成色77777| 国产精品秋霞免费鲁丝片| 一本久久精品| 久久久国产欧美日韩av| 亚洲第一区二区三区不卡| xxx大片免费视频| 好男人视频免费观看在线| 另类精品久久| 国产不卡av网站在线观看| 国产熟女欧美一区二区| 日本vs欧美在线观看视频| 色婷婷久久久亚洲欧美| 国产欧美亚洲国产| e午夜精品久久久久久久| 亚洲欧美精品综合一区二区三区| 欧美黄色片欧美黄色片| 亚洲av成人精品一二三区| 欧美 亚洲 国产 日韩一| 午夜福利,免费看| 国产精品亚洲av一区麻豆 | 女人被躁到高潮嗷嗷叫费观| 叶爱在线成人免费视频播放| www.精华液| 免费看av在线观看网站| 女性生殖器流出的白浆| 可以免费在线观看a视频的电影网站 | 9热在线视频观看99| 日韩,欧美,国产一区二区三区| 人人妻人人澡人人看| 成年人午夜在线观看视频| 男女床上黄色一级片免费看| 亚洲伊人久久精品综合| 一边摸一边抽搐一进一出视频| 欧美老熟妇乱子伦牲交| 99热网站在线观看| av不卡在线播放| 在线精品无人区一区二区三| 日韩制服骚丝袜av| 久久久久精品性色| 亚洲中文av在线| 久久久久精品国产欧美久久久 | 三上悠亚av全集在线观看| videosex国产| 久久久久国产一级毛片高清牌| 少妇 在线观看| 欧美激情极品国产一区二区三区| 久久这里只有精品19| 在线免费观看不下载黄p国产| 女人久久www免费人成看片| 自拍欧美九色日韩亚洲蝌蚪91| 久久久久久人人人人人| 叶爱在线成人免费视频播放| 欧美少妇被猛烈插入视频| 久久久久久人人人人人| 亚洲七黄色美女视频| 黄色怎么调成土黄色| 极品少妇高潮喷水抽搐| 老鸭窝网址在线观看| 黄色毛片三级朝国网站| 免费少妇av软件| 亚洲美女黄色视频免费看| 一级毛片我不卡| 亚洲欧美成人精品一区二区| 女人精品久久久久毛片| 亚洲精品在线美女| 国产av码专区亚洲av| 成人影院久久| 久久精品国产a三级三级三级| 国产av码专区亚洲av| 观看美女的网站| 一级爰片在线观看| 亚洲av成人精品一二三区| 国产精品免费视频内射| 亚洲中文av在线| 国产精品女同一区二区软件| 国产亚洲av片在线观看秒播厂| 一区二区三区精品91| 亚洲欧美一区二区三区久久| 国产成人午夜福利电影在线观看| 最近最新中文字幕免费大全7| 激情视频va一区二区三区| 狠狠婷婷综合久久久久久88av| 黄片播放在线免费| 国产免费又黄又爽又色| 亚洲熟女精品中文字幕| 一二三四在线观看免费中文在| 美女脱内裤让男人舔精品视频| 日日撸夜夜添| 99久久99久久久精品蜜桃| 欧美黄色片欧美黄色片| 黄片无遮挡物在线观看| 精品少妇内射三级| 亚洲精品一二三| e午夜精品久久久久久久| 亚洲av综合色区一区| a级毛片在线看网站| 婷婷色综合www| 人体艺术视频欧美日本| 日日摸夜夜添夜夜爱| 国产av一区二区精品久久| 亚洲免费av在线视频| 热re99久久精品国产66热6| 久久精品国产亚洲av高清一级| 国产成人免费无遮挡视频| 赤兔流量卡办理| 桃花免费在线播放| 欧美日韩亚洲高清精品| 久久天躁狠狠躁夜夜2o2o | 十分钟在线观看高清视频www| 在线观看www视频免费| 男的添女的下面高潮视频| 国产免费一区二区三区四区乱码| a级片在线免费高清观看视频| 国产成人欧美在线观看 | 亚洲av中文av极速乱| 又黄又粗又硬又大视频| 亚洲成人av在线免费| 国产高清国产精品国产三级| 久久国产精品男人的天堂亚洲| 欧美日韩一级在线毛片| 精品国产超薄肉色丝袜足j| 国产极品天堂在线| 一级a爱视频在线免费观看| 男的添女的下面高潮视频| www日本在线高清视频| 亚洲av电影在线进入| 人人妻人人添人人爽欧美一区卜| 80岁老熟妇乱子伦牲交| 久久久久国产一级毛片高清牌| 97在线人人人人妻| 精品久久久久久电影网| 另类精品久久| 成人手机av| 午夜福利视频在线观看免费| 别揉我奶头~嗯~啊~动态视频 | 久久精品人人爽人人爽视色| 91老司机精品| 日韩欧美一区视频在线观看| 在线观看三级黄色| 日韩电影二区| 中文精品一卡2卡3卡4更新| 丝袜喷水一区| e午夜精品久久久久久久| www.自偷自拍.com| 精品久久久久久电影网| 狠狠精品人妻久久久久久综合| 午夜av观看不卡| 国产成人欧美| 日韩大码丰满熟妇| 亚洲成人手机| h视频一区二区三区| 69精品国产乱码久久久| 国产片特级美女逼逼视频| 色视频在线一区二区三区| 男人舔女人的私密视频| 欧美精品一区二区大全| av国产精品久久久久影院| 免费少妇av软件| 成人手机av| 不卡av一区二区三区| 亚洲久久久国产精品| 精品久久蜜臀av无| 啦啦啦啦在线视频资源| 国产精品av久久久久免费| 精品第一国产精品| 国产精品麻豆人妻色哟哟久久|