柯薩奇病毒A6型VP1蛋白的生物信息學分析①

劉洪波 陽廣菲 歐維琳 沈關心

(桂林醫學院第二附屬醫院檢驗科，桂林541199)

?

柯薩奇病毒A6型VP1蛋白的生物信息學分析①

劉洪波 陽廣菲②歐維琳③沈關心④

(桂林醫學院第二附屬醫院檢驗科，桂林541199)

目的：預測柯薩奇病毒A組6型(Coxsackievirus A6，CVA6)的衣殼蛋白VP1的基本理化性質、結構功能及線性B細胞表位。方法：應用Bioedit軟件、SubLoc、TargetP和生物信息學資源門戶ExPASy中的多種在線工具對CVA6 VP1的氨基酸序列進行分析。結果：CVA6 VP1為一親水性蛋白，其相對分子量為33.6 kD，等電點為7.92，含有24個可能的磷酸化位點，沒有信號肽、跨膜區和可能的脂酰化位點；其二級結構中以無規則卷曲居多，有48.52%的氨基酸殘基暴露于溶液界面；該分子內存在多個潛在的線性B細胞表位，其中的155～165位氨基酸殘基區域的抗原指數最高。結論：成功預測到CVA6 VP1的基本理化性質、結構功能特征及可能的線性B細胞表位，為該蛋白的進一步研究及疫苗和免疫診斷試劑的研制打下基礎。

柯薩奇病毒A組6型；VP1蛋白；生物學信息；B細胞表位；理化性質

手足口病(Hand,foot and mouth disease,HFMD)是兒童常見傳染病，其病原包括柯薩奇病毒(Coxsackievirus,CV)、埃柯病毒和新型人腸道病毒(Human enterovirus,HEV)的各部分成員，EV71型和CV A組(CVA)16型為最常見[1,2]。然而，近年其他HEV感染病例有增多趨勢[3]； CVA6型在多個國家和地區引起HFMD的暴發流行，也成為HFMD常見病原之一，并且可導致重癥發生[4]。目前，鮮有針對該病原的基礎性研究。

HEV的VP1具有與其血清型相對應的遺傳多樣性[5]，是其中和表位存在的最主要區域[6]，也是其毒力可能相關的區域之一[7]，具有重要研究價值。本研究運用生物信息學方法，對CVA6代表株VP1的基本理化性質、結構功能特征和含有的線性B細胞表位進行預測分析，為該病原的生物學特性研究以及疫苗和診斷試劑的研發提供理論基礎。

1 材料與方法

1.1 CVA6 VP1氨基酸序列和生物信息學工具 從NCBI數據庫檢索并下載CVA6的全長VP1氨基酸序列，選擇Gdula株(protein ID為AAD17701.1)的VP1序列為代表進行分析。除Bioedit軟件外，本研究主要利用SubLoc V1.0、TargetP 1.1和生物信息學資源門戶網站ExPASy(http://www.expa-sy.org/)提供的多種在線工具對VP1進行分析。這些在線工具的網址鏈接、參考文獻和選擇的參數見表1。

1.2 CVA6 VP1的理化性質預測 利用Bioedit對VP1的氨基酸組成進行分析；利用ProtParam預測相關VP1的原子組成、等電點、消光系數、不穩定系數以及在哺乳動物網織紅細胞、酵母和大腸桿菌中的半衰期等；利用ProtScale中Kyte和Doolittle方法預測共疏水性/親水性。

1.3 CVA6 VP1的結構功能和B細胞表位等的預測 利用SignalP預測VP1蛋白的信號肽；利用Tmpred和TMHMM進行跨膜螺旋預測；利用NetPhos預測磷酸化位點；利用Net-NGlyc和MotifScan預測糖基化、酯酰化等修飾位點，并用Big-PI Predictor預測糖基磷脂酰肌醇(Glycosyl phosphatidyl inositol，GPI)修飾位點；利用SubLoc V1.0和TargetP 1.1對VP1蛋白的亞細胞定位進行預測，并利用NetNES進行核定位信號分析；利用InterProScan、SMART和PROSITE對VP1一級結構中包含的結構和功能域特征序列進行分析；利用SOPMA和Predictprotein預測蛋白質VP1的二級結構和拓撲結構等；運用Bepipred預測VP1中可能的B細胞線性表位。

2 結果

2.1 理化性質

2.1.1 氨基酸的組成 應用Bioedit軟件對CVA6 VP1的氨基酸組成進行分析發現，由305個氨基酸組成的VP1中Ala(28個，9.2%)、Ser(28個，9.2%)、Thr(28個，9.2%)和Val(25個，8.2%)含量較多；所含堿性氨基酸和酸性氨基酸總數分別為26和25，分別占總氨基酸數的8.5%和8.2%(見圖1)。

2.1.2 基本理化性質 利用ProtParam工具預測了VP1蛋白的基本理化參數。該蛋白的原子組成為C1478H2269N415O461S12 ，相對理論分子量為33.6 kD，等電點為7.92。假設全部胱氨酸殘基以Cys形式形成二硫鍵，則其在水溶液中(A280nm)的摩爾消光系數為44 475 L/(mol·cm)；假設所有二硫鍵打開，其在水溶液中(A280nm)的摩爾消光系數則為44 350 L/(mol·cm)。若VP1 N端為Asn，其在哺乳動物體外網織紅細胞中的半衰期為1.4 h，在酵母體內為3 min，在大腸桿菌內大于10 h。其不穩定系數為43.50，屬于不穩定蛋白。

表1 本文使用工具

Tab.1 Tools used in this paper

ToolsnameWebsiteReferenceBig-PIPredictorhttp://mendel.imp.ac.at/sat/gpi/gpi_server.htmlTrendsBiochemSci,2000,25(7):340-341[J].Bepipredhttp://tools.immuneepitope.org/bcell/ImmunomeRes,2006,2:2[J].InterProScanhttp://www.ebi.ac.uk/interpro/scan.html/Bioinformatics,2014,30(9):1236-1240[J].MotifScanhttp://myhits.isb-sib.ch/cgi-bin/motif_scanUnpublishedNetPhoshttp://www.cbs.dtu.dk/services/NetPhos/JMolBiol,1999,294(5):1351-1362[J].Net-NGlychttp://www.cbs.dtu.dk/services/NetNGlyc/UnpublishedNetNEShttp://www.cbs.dtu.dk/services/NetNESProteinEngDesSel,2004,17(6):527-536[J].ProtParamhttp://www.expasy.ch/tools/protparam.htmlTheProteomicsProtocolsHandbook,2005,571-607[M].ProtScalehttp://www.Expasy.org/cgi-bin/protscale.plTheProteomicsProtocolsHandbook,2005,571-607[M].Predictproteinhttps://www.predictprotein.org/homeNucleicAcidsRes,2014,42:W337-43[J].PROSITEhttp://au.expasy.org/prosite/NucleicAcidsRes,2006,34:W362-365[J].SOPMAhttp://npsa-pbil.ibcp.fr/cgi-bin/npsa_automat.pl?page=npsa_sopma.htmlBMCBioinformatics,2006,7(1):255[J].SignalPV2.0.b2http://www.cbs.dtu.dk/services/SignalP-2.0/#submis-sionProcIntConfIntellSystMolBiol,1998,6:122-130[C].SubLochttp://www.bioinfo.tsinghua.edu.cn/SubLoc/Bioinformatics,2001,17(8):721-728[J].SMARThttp://smart.embl-heidelberg.de/NucleicAcidsRes,2015,43:D257-D260[J].Tmpredhttp://www.ch.embnet.org/software/TMPRED_form.htmlBiolChemHoppe-Seyler,1993,374:166[J].TMHMMhttp://www.cbs.dtu.dk/services/TMHMM-2.0/Bioinformatics,2001,17(7):646-653[J].TargetP1.1http://www.cbs.dtu.dk/services/TargetP/JMolBiol,2000,300(4):1005-1016[J].

Note：Tools are listed in alphabetical order;In this study,all tools were run in their default settings except MotifScan,which selected “mot_source:freq_pat” for application;The references for these tools were presented with some essential information for retrieval,but no publications are so far available for the two tools,MotifScan and Net-NGlyc.

2.1.3 疏水性分析 使用ProtParam預測VP1 氨基酸序列平均親水性，其脂肪指數61.08，GRAVY值為-0.459(范圍介于2和-2之間，正值表明為疏水性強，負值表明為親水性強)。運用ProtScale 中Kyte和Doolittle法預測VP1氨基酸序列的疏水性，從結果(圖2)可以看出，在整個肽鏈中，親水性氨基酸均勻分布，遠多于疏水性氨基酸；多肽鏈第161位的Arg分值-2.600最低(親水性最強)，第87位的Gly分值2.011最高(疏水性最強)。這些結果提示，VP1蛋白脂溶性差，親水性好，為一可溶的親水性蛋白。

2.2 功能相關序列/結構分析

圖1 CVA6 VP1的氨基酸組成分布圖Fig.1 Amino acid composition of CVA6 VP1

圖2 CVA6 VP1的疏水性分析Fig.2 Hydrophobicity analysis of CVA6 VP1

2.2.1 信號肽 應用Signal P對VP1蛋白進行信號肽分析，信預測值為0.002，遠低于cutoff值0.5(圖3)，提示VP1為非分泌性蛋白。

2.2.2 跨膜螺旋 使用Tmpred預測，其分值小于0，沒有發現可能的模型；用TMHMM預測，VP1整條肽鏈均處于膜外(圖4)。據此推測，VP1不存在跨膜結構域。

圖3 SignalP預測CVA6 VP1信號肽Fig.3 Prediction of CVA6 VP1 signal peptide by SignalP

圖4 CVA6 VP1的跨膜螺旋預測Fig.4 Prediction of transmembrane helices in CVA6 VP1

圖5 CVA6 VP1的磷酸化位點預測Fig.5 Prediction of phosphorylation sites in CVA6 VP1

圖6 CVA6 VP1的糖基化修飾位點預測Fig.6 Prediction of glycosylation sites in CVA6 VP1

圖7 CVA6 VP1的核定位信號分析Fig.7 Analysis of nuclear localization signal in CVA6VP1

2.2.3 磷酸化位點 應用NetPhos對VP1進行分析的結果顯示，VP1蛋白中共有24個可能的磷酸化位點，13個位于Ser，6個位于Thr，5個位于Tyr 殘基(見圖5)。

2.2.4 酯酰化和糖基化位點 用MotifScan預測，在VP1中沒有找到可能的脂酰化修飾位點；以Big-PI預測也證實，VP1中無可能的GPI錨定位點。用Net-NGlyc預測，在VP1中發現可能的糖基化位點有3個，分別位于第53、137、171位氨基酸殘基(見圖6)。然而，由于VP1不含信號肽序列，因此幾乎不可能產生糖基化修飾。

2.2.5 亞細胞定位 利用SubLoc預測，結果顯示VP1可能定位在細胞核上。而利用TargetP預測，結果顯示VP1 線粒體靶向肽(mTP)和分泌通路信號肽(SP)數值都不高，定位于其他任何位置均有可能。以NetNES進行核定位信號分析，各值均未超過閾值，提示該蛋白不含核定位信號(見圖7)。各工具預測結果不同，可能是不同工具利用的數據庫不同而產生的差異，說明針對病毒蛋白亞細胞定位的數據庫或預測工具尚不成熟，其亞細胞定位需以實驗驗證。

圖8 InterProScan預測CVA6 VP1蛋白一級結構中的保守功能域Fig.8 Prediction of conserved function domains in primary structure of CVA6 VP1 by InterProScan

圖9 CVA6 VP1的拓撲結構、二級結構和親水性特征Fig.9 Characteristics of topological structure,secondary structure and hydrophilicity of CVA6 VP1

圖10 Bepipred工具預測CVA6 VP1的線性B細胞表位Fig.10 Prediction of linear B-cell epitope of CVA6 VP1 by tool Bepipred

圖11 CVA6 VP1 潛在的線性B細胞表位位點Fig.11 Potential linear B-cell epitope sites on CVA6 VP1

2.2.6 一級結構中所含結構和功能域特征序列分析 對VP1序列使用InterProScan分析的結果顯示，VP1含有2個結構和功能域，一個是位于VP1中間與RhV具有同源性的小RNA病毒蛋白結構域，另一個是與DUF3724具有同源性的未知功能的RNA病毒蛋白結構域； 除VP1左端外的大部分序列屬于SSF88633超家族，含有非整合信號(見圖8)。以SMART預測發現，VP1序列第56～217位區域與Pfam數據庫中RhV結構域相似，第262～284位與DUF3724結構域相似。而運用PROSITE預測，未發現存在特異性的結構和功能域。

2.3 二級結構特征 以SOPMA預測，VP1中α-螺旋比例為27.87%，β-折疊為24.92%，β-轉角為8.52%，無規則卷曲為38.69%。以Predictprotein預測，α-螺旋比例為9.84%，β-折疊為25.90%，無規則卷曲為64.26%(見圖9)。結合兩種方法預測的結果可知，α-螺旋和β-折疊散布于整個VP1蛋白，成為蛋白的剛性支撐；而無規則卷曲為其中比例最高的結構元件，間雜少量的β-轉角，為蛋白的柔性區域，是抗原表位存在的可能部位。Predictprotein預測也顯示了VP1中每個氨基酸的溶劑可及性，其中48.52%的氨基酸殘基暴露于溶液界面，44.26%的氨基酸殘基包埋于蛋白質內部，其余7.21%的氨基酸處于中間態(見圖9)。

2.4 線性B細胞表位 運用Bepipred 法預測B細胞線性表位，結果顯示VP1 蛋白平均抗原性指數為0.442，分子內共有多個可能的抗原表位區域(圖10)。其中，分值最高的表位區域為第151-173位氨基酸殘基區段，分值達2.82(閾值0.350)；其次為第177-188，91-103，203-222和25-58位氨基酸區段(按分值高低排列)，其余位點見圖11。

3 討論

在基因組和蛋白質組學時代，運用生物信息學工具分析蛋白質序列已成為蛋白質研究的重要手段[8]。蛋白質理化性質預測可以用于蛋白質分離鑒定的實驗設計等，如分子量、凈電荷和親和力數據為蛋白質色譜分離提供先驗信息；蛋白質序列比對分析可以了解其進化和可能的結構和功能特征等。這些為蛋白質的進一步研究提供諸多信息，可大大提高蛋白質相關研究的工作效率。CVA6與EV71同屬于HEV A組成員，親緣關系很近，它們的VP1有著較高的序列一致性[9]，因此已有的諸多EV71 VP1相關研究可為CVA6 VP1的生物信息學分析提供預示信息[6,10-12]。比較俞海洋等[10]的研究結果，我們發現，CVA6和EV71的VP1蛋白具有很多相同或相似的理化和結構特征，例如，兩者均為可溶的親水性蛋白，無信號肽和跨膜螺旋結構，均存在與RhV相似的結構域，約50%氨基酸殘基暴露于溶液界面，二級結構以無規則卷曲為主等。這些相同或相似的理化和結構也預示它們VP1可能存在相似的功能。蛋白質的可逆磷酸化被認為是細胞調節和信號轉導的最頻繁和通用的機制[13]。病毒蛋白可以通過磷酸化/去磷酸化過程干擾宿主細胞的信號轉導，也使自身的生物學功能發生改變[14,15]。有研究證實，EV71可通過VP1激活細胞鈣調蛋白依賴的蛋白激酶II而對宿主細胞信號轉導和EV71自身的復制產生影響[11]。本研究預測發現，與EV71 VP1相似[10]，CVA6的VP1中也存在較多的磷酸化位點。因此我們推測，CVA6 VP1也是影響其自身復制循環和致病的重要分子。

運用生物信息學方法預測蛋白質的抗原優勢表位，并通過體外合成肽段進行實驗驗證，已廣泛應用于蛋白質抗原表位的鑒定[16]。通過該方法，已成功鑒定出包括EV71病毒在內的多種病原體抗原表位[6,16]。Bepipred方法是結合線性B細胞表位預測的一種經典算法和最佳預測模型-隱馬爾科夫模型發展的一種改良的傾向量表法；以HIV數據集驗證，其曲線下面積優于多種傳統方法[17]。本研究利用該方法預測發現，CVA6 VP1分子的抗原性很強(抗原性指數為0.442)，可能的線性B細胞表位分布于各個區段。與本文的其他預測結果非常相符，如較高的極性氨基酸比例、較強的親水性和較多的無規則卷曲結構等。另外，本研究預測到CVA6的可能B細胞表位25-58和203-222區段與已經實驗證實的EV71的SP12和SP70表位核心區域重疊且存在較多保守性氨基酸殘基[12]。結果提示，該兩肽段含B細胞表位的可能性較大，且EV71和CVA6之間可能存在交叉反應性；另一方面，這些表位中氨基酸殘基變異亦可能是導致這些病原抗原性轉變或其抗體親合力下降的結構基礎[18]，這些值得進一步研究。

總之，本研究對CVA6的主要衣殼蛋白VP1的生物信息學分析，成功預測到其基本的理化性質、結構和功能特征以及潛在的線性B細胞表位，這將有助于對其VP1蛋白的深入研究以及為疫苗和免疫診斷試劑的研制奠定理論基礎。

[1] 鄭 媛,周步宇,關路媛,等.2009-2013年陜西省手足口病流行病學及病原學特征分析[J].中華疾病控制雜志,2015,19(2):135-137.

[2] 海 巖,李忠延,躍 華,等.內蒙古2008—2013年手足口病疫情分析[J].中國公共衛生,2015,31(1):29-31.

[3] 田樹萍,何 菡,王文靜,等.2011-2013年手足口病病原譜及臨床特征分析[J].北京醫學,2014,36(6):447-449.

[4] Bian L,Wang Y,Yao X,etal.Coxsackievirus A6:a new emerging pathogen causing hand,foot and mouth disease outbreaks worldwide[J].Expert Review of Anti-infective Therapy,2015,13(9):1061-1071.

[5] Brown BA,Oberste MS,Alexander JP,etal.Molecular epidemiology and evolution of enterovirus 71 strains isolated from 1970 to 1998[J].J Virol,1999,73(12):9969-9975.

[6] Kirk K,Poh CL,Fecondo J,etal.Cross-reactive neutralizing antibody epitopes against Enterovirus 71 identified by an in silico approach[J].Vaccine,2012,30(49):7105-7110.

[7] Yang E,Zhao H,Zhang Y,etal.A comparative study of the characteristics of two Coxsackie A virus type 16 strains(genotype B)[J].Sci China Life Sci,2012,55(4):336-342.

[8] Gupta PS,Banerjee S,Islam RN,etal.PHYSICO:An UNIX based Standalone Procedure for Computation of Individual and Group Properties of Protein Sequences[J].Bioinformation,2014,10(2):105-107.

[9] Oberste MS,Penaranda S,Maher K,etal.Complete genome sequences of all members of the species human enterovirus A[J].J Gen Virol,2004,85(Pt 6):1597-1607.

[10] 俞海洋,陳 偉，倪德群，等.腸道病毒71型六安地區分離株VP1基因的生物信息學分析及B細胞表位預測[J].安徽醫科大學學報,2012,47(3):233-240.

[11] Haolong C,Du N,Hongchao T,etal.Enterovirus 71 VP1 activates calmodulin-dependent protein kinase II and results in the rearrangement of vimentin in human astrocyte cells[J].PLoS One,2013,8(9):e73900.

[12] Foo DG,Alonso S,Phoon MC,etal.Identification of neutralizing linear epitopes from the VP1 capsid protein of Enterovirus 71 using synthetic peptides[J].Virus Res,2007,125(1):61-68.

[13] Cesaro L,Pinna LA.The generation of phosphoserine stretches in phosphoproteins:mechanism and significance[J].Mol Biosyst,2015,11(10):2666-2679.

[14] Schulze zur Wiesch J,Schmitz H,Borowski E,etal.The proteins of the Hepatitis C virus:their features and interactions with intracellular protein phosphorylation[J].Archives Virol,2003,148(7):1247-1267.

[15] Ohman T,Soderholm S,Paidikondala M,etal.Phosphoproteome characterization reveals that Sendai virus infection activates mTOR signaling in human epithelial cells[J].Proteomics,2015,15(12):2087-2097.

[16] He L,Zhu J.Computational tools for epitope vaccine design and evaluation[J].Curr Opinion Virol,2015,11:103-112.

[17] Larsen JE,Lund O,Nielsen M.Improved method for predicting linear B-cell epitopes[J].Immunome Res,2006,2:2.

[18] Liu CC,Chou AH,Lien SP,etal.Identification and characterization of a cross-neutralization epitope of Enterovirus 71[J].Vaccine,2011,29(26):4362-4372.

[收稿2015-09-06 修回2015-09-15]

(編輯 許四平)

Bioinformatics analysis of VP1 protein of coxsackievirus A6

LIU Hong-Bo,YANG Guang-Fei,OU Wei-Lin,SHEN Guan-Xin.Department of Laboratory Medicine,the Second Affiliated Hospital of Guilin Medical University，Guilin 541199，China

Objective:To predict the basic physicochemical properties,structure and function,and linear B-cell epitopes of the capsid protein VP1 of coxsackievirus A6(CVA6).Methods: The amino acid sequence of the CVA6 VP1 was analyzed using Bioedit software and various online tools including SubLoc，TargetP and the others from ExPASy Bioinformatics Resource Portal.Results: The CVA6 VP1 protein was a hydrophilic protein with a relative molecular weight of 33.6 kD and an isoelectric point of 7.92.This protein containsed 24 phosphorylation sites,but no signal peptide,transmembrane domains and possible fatty acylation sites.Its secondary structure was characterized by the richest random coils,and 48.52 percent of its amino acid residues exposed at the solution interface.Epitope prediction by Bepipred showed a number of potential B cell epitopes in the protein,the highest antigenicity index among them located in the region of amino acids residue 155-165.Conclusion: The basic physicochemical properties,structure and function characteristics,and potential linear B-cell epitopes of CVA6 VP1 were successfully predicted,which laid foundations for the further study on the protein and the preparation of vaccines and immunological diagnostic reagents for CVA6 infection.

Coxsackievirus A6；VP1 protein；Bioinformatics；B-cell epitope；Physicochemical property

10.3969/j.issn.1000-484X.2016.04.019

①本文為國家自然科學基金項目(81460304；81360248)和廣西自然科學基金(2015GXNSFDA139020)和廣西壯族自治區衛生廳科研課題(Z2014298)。

劉洪波(1973年-)，男，博士，副主任技師，碩士生導師，主要從事新型疫苗、分子流行病學和分子診斷技術研究，E-mail:hbliu@glmc.edu.cn。

R373.2+3

A

1000-484X(2016)04-0536-06

②并列第一作者，桂林醫學院藥學院臨床藥學， 桂林541004。

③桂林醫學院附屬醫院兒科，桂林528403。

④華中科技大學同濟醫學院免疫學系，武漢430030。

指導教師：沈關心(1952年-)，男，碩士，教授，博士生導師，主要從事分子免疫與免疫相關疾病研究。