• <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)定
    亚洲在线观看片| 国产日本99.免费观看| а√天堂www在线а√下载| 成人av在线播放网站| 免费观看人在逋| 别揉我奶头~嗯~啊~动态视频| 亚洲国产精品久久男人天堂| 日本熟妇午夜| 免费看a级黄色片| 精品乱码久久久久久99久播| 这个男人来自地球电影免费观看| 99久久精品国产亚洲精品| 国产69精品久久久久777片 | 18禁美女被吸乳视频| 久久中文字幕人妻熟女| 99精品久久久久人妻精品| 老司机福利观看| 法律面前人人平等表现在哪些方面| 看免费av毛片| 动漫黄色视频在线观看| 久久亚洲精品不卡| 天天添夜夜摸| 日韩av在线大香蕉| 中文字幕人妻丝袜一区二区| 成年女人看的毛片在线观看| www.精华液| 操出白浆在线播放| 成人无遮挡网站| 母亲3免费完整高清在线观看| 久久精品91无色码中文字幕| 日本a在线网址| 九九在线视频观看精品| a在线观看视频网站| 天堂动漫精品| 亚洲熟妇中文字幕五十中出| 日日干狠狠操夜夜爽| 亚洲欧美日韩高清专用| 亚洲天堂国产精品一区在线| 免费看光身美女| 大型黄色视频在线免费观看| 日韩欧美一区二区三区在线观看| 国产成人精品久久二区二区91| bbb黄色大片| 亚洲一区高清亚洲精品| 男插女下体视频免费在线播放| 中文字幕人妻丝袜一区二区| 男人舔奶头视频| 99久久99久久久精品蜜桃| 校园春色视频在线观看| 欧美xxxx黑人xx丫x性爽| 成年女人永久免费观看视频| 精品国产美女av久久久久小说| 国产成人系列免费观看| 一个人免费在线观看的高清视频| 亚洲最大成人中文| 窝窝影院91人妻| 亚洲国产色片| 午夜亚洲福利在线播放| 狠狠狠狠99中文字幕| 69av精品久久久久久| 亚洲五月天丁香| 精品久久久久久久末码| 嫩草影院入口| 国产成人精品无人区| 亚洲精华国产精华精| 丝袜人妻中文字幕| 国产真实乱freesex| 无人区码免费观看不卡| 国产探花在线观看一区二区| 久久国产精品影院| 久久精品aⅴ一区二区三区四区| 日本免费a在线| 狂野欧美白嫩少妇大欣赏| 久久久久久久久免费视频了| 亚洲人成网站在线播放欧美日韩| 在线永久观看黄色视频| 亚洲黑人精品在线| 亚洲九九香蕉| 久久精品国产清高在天天线| 日韩欧美 国产精品| 国产精品亚洲美女久久久| 又黄又粗又硬又大视频| www日本在线高清视频| 久久香蕉精品热| 每晚都被弄得嗷嗷叫到高潮| 好男人电影高清在线观看| 国产三级中文精品| 国产野战对白在线观看| www日本黄色视频网| av国产免费在线观看| 国产一区二区在线av高清观看| 99久久无色码亚洲精品果冻| 色综合欧美亚洲国产小说| 欧美日韩国产亚洲二区| 精品久久久久久,| 午夜视频精品福利| 一二三四社区在线视频社区8| 亚洲五月婷婷丁香| 制服人妻中文乱码| 老汉色∧v一级毛片| 国产伦精品一区二区三区视频9 | 日本五十路高清| 国产成年人精品一区二区| 欧美激情在线99| 非洲黑人性xxxx精品又粗又长| 九色国产91popny在线| 日韩精品中文字幕看吧| 国产精品综合久久久久久久免费| 国产综合懂色| 男插女下体视频免费在线播放| 天天添夜夜摸| 久久久久久国产a免费观看| 天天添夜夜摸| 99国产精品一区二区蜜桃av| 国产av不卡久久| 久久久久久久午夜电影| av在线天堂中文字幕| 男女之事视频高清在线观看| 欧美日韩乱码在线| 亚洲七黄色美女视频| 两个人视频免费观看高清| 亚洲五月婷婷丁香| xxx96com| 午夜两性在线视频| 亚洲无线在线观看| 老司机午夜十八禁免费视频| 国产午夜精品论理片| 亚洲专区中文字幕在线| 欧美日韩中文字幕国产精品一区二区三区| xxxwww97欧美| 国产免费男女视频| 色噜噜av男人的天堂激情| 国内毛片毛片毛片毛片毛片| 女警被强在线播放| 日韩欧美在线乱码| 国产伦精品一区二区三区视频9 | 欧美一区二区精品小视频在线| 国产私拍福利视频在线观看| 成年人黄色毛片网站| 国产久久久一区二区三区| 国内精品一区二区在线观看| 成人三级黄色视频| 日本免费a在线| 欧美中文日本在线观看视频| 亚洲va日本ⅴa欧美va伊人久久| 国产男靠女视频免费网站| 九九热线精品视视频播放| 欧美乱妇无乱码| 中文字幕最新亚洲高清| 不卡一级毛片| 午夜福利免费观看在线| 99国产精品一区二区蜜桃av| 听说在线观看完整版免费高清| 亚洲国产精品sss在线观看| 国产野战对白在线观看| 丰满人妻一区二区三区视频av | 国产精品一区二区免费欧美| 欧美极品一区二区三区四区| 欧美成人免费av一区二区三区| 丁香六月欧美| 一区二区三区高清视频在线| 久久久久久久午夜电影| 老司机午夜十八禁免费视频| 又黄又爽又免费观看的视频| 九九久久精品国产亚洲av麻豆 | 在线播放国产精品三级| 精品福利观看| 国产主播在线观看一区二区| 欧美xxxx黑人xx丫x性爽| av黄色大香蕉| 亚洲国产高清在线一区二区三| 黄片小视频在线播放| 一边摸一边抽搐一进一小说| 国产在线精品亚洲第一网站| 国产乱人视频| 99久久精品热视频| 欧美成人一区二区免费高清观看 | 波多野结衣高清无吗| 哪里可以看免费的av片| 精品无人区乱码1区二区| 久久精品国产综合久久久| 国产精品影院久久| av视频在线观看入口| 亚洲一区高清亚洲精品| 免费av不卡在线播放| 久久婷婷人人爽人人干人人爱| 欧美黑人巨大hd| 法律面前人人平等表现在哪些方面| 国产毛片a区久久久久| 热99在线观看视频| 久久九九热精品免费| 日韩大尺度精品在线看网址| 女人高潮潮喷娇喘18禁视频| 久久亚洲真实| 亚洲在线自拍视频| 一进一出抽搐gif免费好疼| 免费高清视频大片| 欧美色欧美亚洲另类二区| 亚洲午夜理论影院| 老汉色∧v一级毛片| 国产成人福利小说| 欧美乱妇无乱码| 亚洲在线观看片| 91字幕亚洲| 偷拍熟女少妇极品色| 精品无人区乱码1区二区| 国产伦在线观看视频一区| 久久国产乱子伦精品免费另类| 两个人视频免费观看高清| 法律面前人人平等表现在哪些方面| 日本成人三级电影网站| 亚洲精品色激情综合| 亚洲成av人片免费观看| 可以在线观看毛片的网站| 九九热线精品视视频播放| x7x7x7水蜜桃| 日本熟妇午夜| www日本黄色视频网| 亚洲精品中文字幕一二三四区| 亚洲男人的天堂狠狠| 九九在线视频观看精品| 国产欧美日韩精品亚洲av| 一进一出抽搐动态| 亚洲午夜理论影院| 女人高潮潮喷娇喘18禁视频| 免费看十八禁软件| 国产成人av激情在线播放| 国产高潮美女av| 国产成人福利小说| 操出白浆在线播放| 久久久久久久久中文| 97碰自拍视频| 哪里可以看免费的av片| 国产单亲对白刺激| 99视频精品全部免费 在线 | av女优亚洲男人天堂 | 亚洲专区中文字幕在线| 精品免费久久久久久久清纯| av片东京热男人的天堂| 老司机在亚洲福利影院| 在线免费观看的www视频| 美女cb高潮喷水在线观看 | 啦啦啦免费观看视频1| 性色av乱码一区二区三区2| 亚洲av电影不卡..在线观看| 身体一侧抽搐| 我要搜黄色片| 久久国产精品人妻蜜桃| 亚洲精品在线观看二区| 舔av片在线| 91麻豆av在线| 黄色日韩在线| 亚洲av美国av| av片东京热男人的天堂| 亚洲最大成人中文| 老汉色av国产亚洲站长工具| 国产成人av教育| 少妇熟女aⅴ在线视频| avwww免费| 十八禁人妻一区二区| 一区二区三区国产精品乱码| 一级a爱片免费观看的视频| 夜夜爽天天搞| 国产精品野战在线观看| 女人高潮潮喷娇喘18禁视频| 免费搜索国产男女视频| 小蜜桃在线观看免费完整版高清| 国内少妇人妻偷人精品xxx网站 | 色播亚洲综合网| 国产综合懂色| 久久久国产欧美日韩av| 麻豆av在线久日| 成人一区二区视频在线观看| 婷婷精品国产亚洲av在线| 成人午夜高清在线视频| 亚洲在线观看片| 国产三级中文精品| a级毛片在线看网站| 婷婷丁香在线五月| 国产精品国产高清国产av| 深夜精品福利| 欧美成人性av电影在线观看| 日韩欧美国产在线观看| 日本五十路高清| 中亚洲国语对白在线视频| 免费观看精品视频网站| av欧美777| 午夜精品在线福利| 国产免费男女视频| 久久久国产精品麻豆| 国内精品久久久久久久电影| 免费无遮挡裸体视频| 热99在线观看视频| xxxwww97欧美| 亚洲国产高清在线一区二区三| 女生性感内裤真人,穿戴方法视频| 色av中文字幕| 精品人妻1区二区| 噜噜噜噜噜久久久久久91| 特大巨黑吊av在线直播| 99国产精品一区二区三区| 久久久成人免费电影| 久99久视频精品免费| 免费看光身美女| 久久久久久久久久黄片| 国产精品亚洲av一区麻豆| 国产精品久久久久久久电影 | 国产精品影院久久| 色播亚洲综合网| 看免费av毛片| 国产伦精品一区二区三区视频9 | 村上凉子中文字幕在线| 亚洲中文av在线| 亚洲av日韩精品久久久久久密| 最近最新免费中文字幕在线| 国产欧美日韩一区二区精品| 美女 人体艺术 gogo| 好看av亚洲va欧美ⅴa在| 99re在线观看精品视频| 男女视频在线观看网站免费| 韩国av一区二区三区四区| 久久久成人免费电影| 精品无人区乱码1区二区| 久久国产精品影院| 搡老妇女老女人老熟妇| 啦啦啦观看免费观看视频高清| 亚洲真实伦在线观看| 中亚洲国语对白在线视频| 国产av麻豆久久久久久久| 亚洲国产看品久久| 国产真实乱freesex| 老熟妇乱子伦视频在线观看| 国产精品日韩av在线免费观看| 成人鲁丝片一二三区免费| 成人特级黄色片久久久久久久| 一个人免费在线观看的高清视频| 亚洲第一电影网av| 男人舔女人的私密视频| 午夜福利视频1000在线观看| 在线观看舔阴道视频| 婷婷六月久久综合丁香| 国产日本99.免费观看| 一级黄色大片毛片| 51午夜福利影视在线观看| e午夜精品久久久久久久| 国产精品综合久久久久久久免费| 一夜夜www| 亚洲在线观看片| 精品久久久久久,| 最近视频中文字幕2019在线8| 欧美大码av| 欧美在线黄色| 国产美女午夜福利| 亚洲自偷自拍图片 自拍| 亚洲熟女毛片儿| 麻豆国产97在线/欧美| 欧美成人免费av一区二区三区| 人妻夜夜爽99麻豆av| 精品国产乱码久久久久久男人| 女警被强在线播放| 国产99白浆流出| 亚洲九九香蕉| 国产精品野战在线观看| 国产成年人精品一区二区| 久久久国产成人精品二区| 三级男女做爰猛烈吃奶摸视频| 亚洲狠狠婷婷综合久久图片| 国产探花在线观看一区二区| 日韩欧美三级三区| 午夜福利成人在线免费观看| 精品电影一区二区在线| 精品久久久久久久末码| 在线十欧美十亚洲十日本专区| 久久99热这里只有精品18| 中文字幕精品亚洲无线码一区| 午夜两性在线视频| 色综合站精品国产| 国产aⅴ精品一区二区三区波| 美女免费视频网站| 国产精品美女特级片免费视频播放器 | 精品无人区乱码1区二区| 在线a可以看的网站| 亚洲av电影在线进入| 国产综合懂色| 丁香欧美五月| xxxwww97欧美| www.自偷自拍.com| 国产成人aa在线观看| 国产精品久久久久久久电影 | 露出奶头的视频| 观看免费一级毛片| 一a级毛片在线观看| 在线免费观看不下载黄p国产 | 免费无遮挡裸体视频| 最新中文字幕久久久久 | 国产精品影院久久| 这个男人来自地球电影免费观看| 又黄又粗又硬又大视频| 国产黄a三级三级三级人| 国产成人啪精品午夜网站| 村上凉子中文字幕在线| 亚洲国产色片| 欧美日韩黄片免| 69av精品久久久久久| 中文资源天堂在线| 亚洲熟妇中文字幕五十中出| 精品一区二区三区av网在线观看| 日韩欧美在线乱码| 国产蜜桃级精品一区二区三区| 欧美丝袜亚洲另类 | 亚洲avbb在线观看| 精品国产亚洲在线| 精品欧美国产一区二区三| 久久精品亚洲精品国产色婷小说| 长腿黑丝高跟| 九色国产91popny在线| 久久人人精品亚洲av| 在线观看日韩欧美| 身体一侧抽搐| 黄色片一级片一级黄色片| 国内精品美女久久久久久| 久久久久国内视频| 精品久久久久久久人妻蜜臀av| 亚洲美女黄片视频| 国产一区二区在线av高清观看| 国产精品久久久久久人妻精品电影| 99久久成人亚洲精品观看| 97超视频在线观看视频| av女优亚洲男人天堂 | 亚洲欧美精品综合一区二区三区| 欧美一区二区精品小视频在线| 久久中文看片网| 久久久国产成人免费| 波多野结衣巨乳人妻| 国产一级毛片七仙女欲春2| 国产主播在线观看一区二区| 一本一本综合久久| 国产亚洲精品一区二区www| 草草在线视频免费看| 可以在线观看的亚洲视频| 特级一级黄色大片| 久久久久久久久久黄片| 国产精品亚洲av一区麻豆| 日本成人三级电影网站| 日韩欧美三级三区| 国产亚洲欧美98| 欧美日韩精品网址| 亚洲精品粉嫩美女一区| 校园春色视频在线观看| 搡老岳熟女国产| 亚洲成a人片在线一区二区| 免费在线观看影片大全网站| 欧美xxxx黑人xx丫x性爽| 国产精品一区二区精品视频观看| cao死你这个sao货| 亚洲专区国产一区二区| 亚洲精品中文字幕一二三四区| 一本精品99久久精品77| 免费av毛片视频| 在线观看免费午夜福利视频| 亚洲美女黄片视频| 色视频www国产| 日韩精品青青久久久久久| 热99re8久久精品国产| 老熟妇乱子伦视频在线观看| 免费人成视频x8x8入口观看| 成人性生交大片免费视频hd| 午夜精品在线福利| 国产av麻豆久久久久久久| 国产真人三级小视频在线观看| 久久久久九九精品影院| 亚洲在线自拍视频| 日韩欧美精品v在线| 91av网站免费观看| 亚洲精品456在线播放app | 白带黄色成豆腐渣| 国产黄a三级三级三级人| 亚洲国产中文字幕在线视频| 2021天堂中文幕一二区在线观| 亚洲人成伊人成综合网2020| 悠悠久久av| 又紧又爽又黄一区二区| 九九热线精品视视频播放| 亚洲avbb在线观看| 成人永久免费在线观看视频| 亚洲精品在线美女| 亚洲成人久久性| 欧美在线一区亚洲| 国内精品久久久久精免费| 97碰自拍视频| 亚洲色图 男人天堂 中文字幕| 美女 人体艺术 gogo| 亚洲午夜理论影院| 九九久久精品国产亚洲av麻豆 | 99久久久亚洲精品蜜臀av| 国产精品精品国产色婷婷| 可以在线观看毛片的网站| 91字幕亚洲| a级毛片a级免费在线| 中文在线观看免费www的网站| 99热精品在线国产| 在线永久观看黄色视频| 男人和女人高潮做爰伦理| 国产成人精品久久二区二区免费| 99热只有精品国产| 久久这里只有精品19| 波多野结衣巨乳人妻| 波多野结衣高清作品| 亚洲国产中文字幕在线视频| 91av网一区二区| 狂野欧美白嫩少妇大欣赏| 国产三级中文精品| 久99久视频精品免费| 久久久久国产精品人妻aⅴ院| 久久精品影院6| tocl精华| 婷婷亚洲欧美| 久久中文字幕一级| 天堂av国产一区二区熟女人妻| 国产一区二区在线观看日韩 | 亚洲七黄色美女视频| 两个人看的免费小视频| av欧美777| 非洲黑人性xxxx精品又粗又长| 亚洲中文字幕日韩| 综合色av麻豆| 一级毛片精品| 久久99热这里只有精品18| 黄色片一级片一级黄色片| 在线十欧美十亚洲十日本专区| 在线视频色国产色| 性色av乱码一区二区三区2| 国产激情偷乱视频一区二区| 日本 av在线| 日本黄色片子视频| 黄色 视频免费看| 国内毛片毛片毛片毛片毛片| 国产激情欧美一区二区| 国产毛片a区久久久久| 久久婷婷人人爽人人干人人爱| 在线看三级毛片| 亚洲自偷自拍图片 自拍| АⅤ资源中文在线天堂| 熟女人妻精品中文字幕| 白带黄色成豆腐渣| av在线天堂中文字幕| 欧美丝袜亚洲另类 | 午夜福利成人在线免费观看| 黄色成人免费大全| 精品国产乱码久久久久久男人| 精品国产乱子伦一区二区三区| 国产男靠女视频免费网站| 亚洲av第一区精品v没综合| 麻豆一二三区av精品| 一进一出好大好爽视频| 在线观看日韩欧美| 在线观看66精品国产| 国产高清视频在线观看网站| 国产成人一区二区三区免费视频网站| 一二三四在线观看免费中文在| АⅤ资源中文在线天堂| 校园春色视频在线观看| 在线十欧美十亚洲十日本专区| 少妇丰满av| 日本黄色片子视频| 曰老女人黄片| 国产伦精品一区二区三区视频9 | 成人18禁在线播放| 欧美日韩精品网址| 久久中文字幕一级| 看片在线看免费视频| 丰满人妻一区二区三区视频av | 亚洲精品色激情综合| 欧美日韩亚洲国产一区二区在线观看| 亚洲va日本ⅴa欧美va伊人久久| 99国产精品一区二区三区| 在线免费观看的www视频| 午夜成年电影在线免费观看| 国产野战对白在线观看| 日韩有码中文字幕| 黄色 视频免费看| 日本 欧美在线| 欧美日本视频| 国产欧美日韩精品亚洲av| 亚洲avbb在线观看| 香蕉av资源在线| 欧美日韩亚洲国产一区二区在线观看| 国产精品,欧美在线| 国产激情偷乱视频一区二区| 噜噜噜噜噜久久久久久91| 一级毛片高清免费大全| 日韩 欧美 亚洲 中文字幕| 久久伊人香网站| 日韩中文字幕欧美一区二区| 国产美女午夜福利| 一区二区三区激情视频| 狠狠狠狠99中文字幕| 国产午夜精品论理片| 女生性感内裤真人,穿戴方法视频| 久久久久国产一级毛片高清牌| 欧美日韩国产亚洲二区| 久久久久久久久久黄片| 少妇的丰满在线观看| 国产视频一区二区在线看| 久久亚洲真实| 无限看片的www在线观看| 国产精品久久久人人做人人爽| 国产成人一区二区三区免费视频网站| 日本黄大片高清| 91久久精品国产一区二区成人 | 亚洲第一电影网av| 国产一区在线观看成人免费| 免费在线观看成人毛片| 欧美日韩国产亚洲二区| 国产av不卡久久|