



全文預(yù)覽已結(jié)束
下載本文檔
版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)
文檔簡(jiǎn)介
VAMA:aversatileweb-basedtoolforvariabilityanalysisinmultiply-alignedaminoacidsequencesVAMA:VariabilityAnalysisofMultipleAlignmentsAditiGupta,AridamanPanditandSomdattaSinhaMathematicalModelingandComputationalBiologyGroup,CentreforCellularandMolecularBiology(CSIR)Hyderabad500007,IndiaEmail:sinhaccmb.res.inAbstractQuantifyingresiduevariabilityateachcolumninamultiplesequencealignmentofaminoacidshelpsinindicatingtheirsimilarities,andisusefultohighlightinformationaboutthesignificancesofeachpositionfromtheperspectiveoftheirstructure,function,andevolution.Itisbecomingincreasinglyclearthatthegroupsofaminoacidsthatallowconservedreplacementvarywiththepositionoftheresidueintheprotein.Mostmultiplealignmentalgorithmscatertogeneralusersandhencedonotaddressthisspecificfeature.Atoolforscoringvariabilityinmultiply-alignedaminoacidsequences,thatallowsdifferentconservationgroups,ishighlydesirable.VAMA(VariabilityAnalysisofMultipleAlignments)isasimpleyetversatileprogramthatcalculatesandplotsresiduevariabilityinagivensetofalignedsequencesbasedonknownconservationgroupsspecificfordifferentfunctionallyimportantregionsofaprotein,andalsoallowsuser-definedgroupsfornewdiscoveries.VAMAisavailableat84/VAMA/Overview.htmlKeywords-MultipleSequenceAlignment,VariabilityAnalysis,ResidueConservationGroupsI.INTRODUCTIONAlignmentofaminoacidsequencesiswidelyappliedtoidentifyconservationofresidues1.Variabilityisameasureoftheextentofvariationofaminoacidsatapositioninmultiplesequencealignment.Functionallyimportantresiduesareknowntoexhibithigherconservation,orlowervariability.Multiplesequencealignment(MSA)toolsalignsequencesbasedonsomepredeterminedclassificationofaminoacids,dependingontheirphysicochemicalproperties2.Recentstudieshaveshownthatthesamegroupofaminoacidsmaynotalwaysbeuseful,assequenceconservationclassificationsvarywiththestructureandfunctionaloftheprotein.Forexample,differentclassificationshavebeenshowntoexistforresiduesinvolvedinligandbinding3,inprotein-proteininteractioninterfaces4,inmaintainingstructure5,andindeterminingproteinfunctionalspecificity6-8.ClassificationbyMirnyandShakhnovich(MS)wasintendedforproteinstructurecores6;Williamsons(W)wastailoredtodealwithtransporterproteins8,whileGuharoyandChakrabarti(GC)classificationwasmeantforthestudyofinterfacialaminoacids4.ThissuggeststhatdifferentsubstitutionclassificationsshouldbeappliedinMSAbaseduponstructure/functionoftheproteins.Programmesavailableforstudyingconservationorvariabilityinmultiplealignmentsbasedondifferentscoringmethodsarerathercomplex9,andnoneofthemconsideralltheabove-mentionedposition-specificfeaturesinproteins10-15.Atoolforscoringvariabilityinmultiply-alignedaminoacidsequences,thatissimplebutaddressthisspecificfeatureofallowingdifferentconservationgroups,ishighlydesirable.Wehavedevelopedaweb-basedprogram(VAMA)thatquantifiesthevariabilityoftheresiduesateachalignedpositioninMSAusingasimplesymboldiversityscore16.Here,alldifferentaminoacidresiduespresentinaparticularcolumnofthealignedsequencesareconsidered,andthescore,v,iscalculatedusingthesumofthefrequencyofeachresidueas()1211iinNvN=(1)where,niisthefrequencyofeachresidue,andNthetotalnumberofresiduesatthisposition.Thesumisoverallthedifferenttypesofaminoacidspresentinthecolumn.ForcompleteconservationofaminoacidsataparticularcolumnintheMSA(i.e.,n=N),thescore,v=0,indicatingnovariability;whereas,fornoconservationortotaldiversity,thescoreisv=1.Since,vvariesbetween0and1,anormalizedscoreisobtained,whichisusefulforcomparisonofdifferentMSAs.BecauseofthegenericnatureofscoringinVAMA,whichisnotbasedonanystereochemicalpropertyoftheaminoacids,vcanbeusedtocalculatetherelativefrequenciesofaminoacidsforanygivenclassification.Here,alongwiththebasicscoringmethodologyadoptedbycommonlyusedMSAprogrammeCLUSTAL17,severaloptionsareprovidedforfunction-specificclassificationsasmentionedearlier.Thus,VAMAprovidestheuserwithflexibilitytoquantifyvariabilityusingthesedifferentclassificationsasrequired.II.FEATURESOFVAMAFig.1showstheVAMAinterface.TherearethreewaysinwhichvariabilityanalysisisdoneinVAMA:Basic,GroupBasedandReferenceSequenceBased.TheinputforVAMAcanbemultiplealignmentfilesinCLUSTALandFASTA978-1-4244-4713-8/10/$25.002010IEEEFigure1.VAMAinterface.formats.ThesefilescanbepastedontotheVAMAworkwindow,ormaybeuploadedusingthe“Browse”button.VAMAalsoaddressesthecommonproblemofexistenceof“gaps”inthealignment.TheusercandefineaGapcutoffforincludingresiduepositionshavinggapsinthealignment.Forcolumns,havingmoregapsthanGapcutoff,variabilityisnotcalculatedandablankspaceisdisplayedintheoutput.ForcolumnshavinggapslessthanorequaltothedefinedGapcutoff,thevalueofNisadjustedbysubtractingfromN,thenumberofgapsatthatposition.Thesepositionsareindicatedbya#intheoutput.VAMAcalculatesthevariabilityscoreforeachcolumninthealignmentbasedontheconservationgroupschosenbytheuser,andtheoutputfileconsistsofthefollowingparts-(i)thevariabilityscoreateachalignedposition;(ii)statisticsofthevariabilitydatadisplayingthemean,standarddeviationandrangeforthesame;and(iii)theplotofvariabilityvaluesversusalignmentpositions.ThedatacanbesavedbothintextandEXCELformatforfurtheranalysis.ThevariabilityanalysisinVAMAcanbedoneinthefollowingthreewaysA.BasicVariabilityAnalysisHereallaminoacidsareconsideredtobeinseparateclasses.Hence,itdoesnottakeintoaccountanysubstitutions,andanynon-identitycontributestothevariability.B.GroupBasedVariabilityAnalysisConservativesubstitutionscanbeaccountedforbyclassifyingaminoacidsaccordingtotheirphysicochemicalpropertiesandpositionalattributes.Variabilityisassigned0foraparticularpositioniftheaminoacidsbelongtoagroupinthespecificclassification.VAMAprovidesthefollowinggroupbasedanalysisoptions:i.DefaultClassification:ThisclassificationissameastheoneusedinCLUSTAL17.AminoacidsareclassifiedintostrongandweakgroupsbasedonthephysicochemicalpropertiesandtheGonnetPam250matrix18.ii.MS,GC,andWClassifications:Severaldifferentclassificationsareproposeddependinguponthefunctionalconstraintsapplicable.MSclassificationisapplicabletoproteinstructurecores6,Wclassificationisapplicabletotransporterproteins8,andGCclassificationisapplicabletotheinterfaceaminoacids4.TheclassificationsaredescribedintheUserGuide.iii.UserDefinedClassification:Variousotherclassificationschemeshavebeenproposed7.VAMAoffersuserstheoptiontouseanyortheirownclassification.C.ReferenceSequenceBasedVariabilityAnalysisReferencesequenceisthesequencewithrespecttowhichtheresiduesarenumbered.Herex=0inthevariabilityplotcorrespondstofirstresidueofthereferencesequence.Thisisusefulwhenthe3-dimensionalstructureofthereferencesequenceisavailabletohelpanalyzetheresultsforposition-specificchangesinotheraminoacidsequencesinthealignment.Thecalculationisfirstdonebasedupondifferentclassifications,andthentheresiduesarenumberedaccordingtotheReferenceSequencegivenbytheuser.BythismethodtheusercanaccessthealignmentagainsttheReferenceSequence.AnExampleofanalysis,usingVAMA,isshowninFig.2.Asetof25aminoacidsequencesrepresentingtheRosmannfold19wereextractedfromProteinDataBank20.CLUSTAL-alignedsequenceswerepastedasinputtoVAMA.WeusedMSclassificationtocalculatethevariability,andcompareditwiththeDefault(CLUSTAL)classification.Fig.2AisthevariabilityplotforMSandDefaultclassifications,showinglowerscoreforMSclassificationthanthatoftheDefaultclassification.InFig.2A,Gapcutoffof0isusedtocalculatethescores.Thus,forpositionswith1ormoregaps,variabilityisnotcalculated.Rosmannfoldbeinganexampleofproteinstructurecore,theresultswithMSaremorereliable.Fig.2BalsoshowsthattheMSclassificationgivestheloweststatisticswhencomparedtoothers.Hereweshowonlythefirstsubsetofpositions(92-100)forwhichthevariabilityscoreswerecalculated.Clearly,theminimumvariabilityscoregivenbyMSclassificationadvocatestheapplicabilityofthistool.Importantly,incaseofaproteinlackingwell-definedfunction,VAMAallowscalculationofthevariabilityscoreusingthegivenclassificationstoidentifyfunctionallyandstructurallyimportantresiduesbasedontheircomparativescore.Thisfeature,whereseveraldifferentclassificationscanbeusedtocalculatethevariabilityinMSA,isuniquetoVAMA.Figure2.AnalysisofVAMAoutputfor25aminoacidsequencesoftheRossmanFold.(A)VariabilityplotcomparingMSandDefaultclassificationscores,(B)Variabilityvaluesofresidues92to100,usingdifferentclassifications,alongwithmeanandstandarddeviation(SD).VAMAincludesseveralusefulfeaturesfortheuser.The“UserGuide”explainsallfeaturesclearlywithexample.The“RelatedLinks”provideslinkstootherMSAtools(e.g.CLUSTALW,T-Coffee,CINEMA,etc.),andusefulwebsitessuchas,NCBI,PDB,Swiss-Prot,andKEGG.A“Search”buttonallowssearchwithinVAMAandtheWorldWideWeb.III.DISCUSSIONVAMAisasimple,user-friendly,yetversatile,toolforcalculatingthevariabilityinmultiplyalignedproteinsequences.VAMAsupportsbothbasicandgroupbasedanalysis,byconsideringthephysicochemicalpropertiesofaminoacids,aswellastheirdifferentialusagefordifferenttopologicaldeterminantsintheprotein.Itquantifiesvariabilityinthesequencesbasedonaminoacidclassificationgroupsdependingontheabovefactors.Gap-cutofffeatureactsasanadditionaltooltoscorethesequences.Statisticalanalysisperformedonthevariabilitydata,likemean,standarddeviationandrange,canalsobehelpfulinfurtheranalysis.VAMAis,thus,ausefulfunction-specificvariabilityanalysistoolthatallowsacomparativeanalysisafeaturelackinginothersimilartools.ACKNOWLEDGMENTSSthanksDepartmentofBiotechnology,Indiaforfinancialsupport.AGthankstheIndianAcademyofSciencesforasummerfellowshiptoworkattheCCMB.APthankstheCouncilofScientificandIndustrialResearch(CSIR)forfellowship.REFERENCES1T.F.Smith,andM.S.Waterman,“IdentificationofCommonMolecularSubsequences.”J.Mol.Biol.,vol.147,pp.195-197,1981.2D.J.Lipman,S.F.Altschul,andJ.D.Kececioglu,“Atoolformultiplesequencealignment.”Proc.Natl.Acad.Sci.USA,vol.86,pp.4412-4415,1989.3T.J.Magliery,andL.Regan,“Sequencevariationinligandbindingsitesinproteins.”BMCBioinformatics,vol.6,p.240,2005.4M.Guharoy,andP.Chakrabarti,“Conservationandrelativeimportanceofresiduesacrossprotein-proteininterfaces.”Proc.Natl.Acad.Sci.USA,vol.102,pp.15447-15452,2005.5O.Schueler-Furman,andD.Baker,“Conservedresidueclusteringandproteinstructureprediction.”Proteins:StructureFunctionandBioinformatics,vol.52,pp.225-235,2003.6L.A.Mirny,andE.I.Shakhnovich,“Evolutionaryconservationofthefoldingnucleus.”J.Mol.Biol.,vol.308,pp.123-129,2001.7W.R.Taylor,“Theclassificationofaminoacidconservation.”J.Theor.Biol.,vol.119,pp.205-218,1986.8R.M.Williamson,“Informationtheoryanalysisoftherelationshipbetweenprimarysequencestructureandligandrecognitionamongaclassoffacilitatedtransporters.”J.Theor.Biol.,vol.174,pp.179-188,1995.9W.S.J.Valdar,“Scoringresidueconservation.”Proteins:Structure,FunctionandBioinformatics,vol.48,pp.227-241,2002.10J.A.Capra,andM.Singh,“Predictingfunctionallyimportantresiduesfromsequenceconservation.”Bioinformatics,vol.23,pp.1875-1882,2007.11M.Clamp,J.Cuff,S.M.Searle,andG.J.Barton,“TheJalviewJavaAlig
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫(kù)網(wǎng)僅提供信息存儲(chǔ)空間,僅對(duì)用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。
最新文檔
- 車輛無償借用及緊急救援服務(wù)協(xié)議
- 彩票站行業(yè)規(guī)范與自律管理合作協(xié)議
- 右心衰竭臨床診療規(guī)范
- 嘔血窒息護(hù)理要點(diǎn)與實(shí)施規(guī)范
- 天然藥物學(xué):附子專題研究
- 2025年逆回購(gòu)協(xié)議
- 小兒骨折護(hù)理要點(diǎn)
- 積水治療與護(hù)理
- 低鉀血癥護(hù)理
- 細(xì)胞環(huán)境與互作
- 2025工程建設(shè)項(xiàng)目多測(cè)合一成果報(bào)告書范本
- 麻醉科麻精藥品PDCA管理
- 語言習(xí)得神經(jīng)機(jī)制探究-深度研究
- 兒童發(fā)展問題的咨詢與輔導(dǎo)-案例1-5-國(guó)開-參考資料
- 2025年河北石家莊市市屬國(guó)有企業(yè)招聘筆試參考題庫(kù)含答案解析
- 2025年國(guó)航股份地面服務(wù)部校園招聘筆試參考題庫(kù)含答案解析
- 宣傳物料制作合同范本
- 2025年度安徽白帝集團(tuán)限公司社會(huì)招聘高頻重點(diǎn)提升(共500題)附帶答案詳解
- 公益招貼設(shè)計(jì)課件
- 電廠灰?guī)旖Y(jié)構(gòu)工程施工方案
- 靜脈治療小組開展工作匯報(bào)
評(píng)論
0/150
提交評(píng)論