# #——环境,回声= FALSE,消息= FALSE,警告= FALSE - - - - - - - - - - - - - - - - - - - - - -库(“刀”)库(UniProt.ws)库(“大脑”)# # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -库(“刀”)# # - - - - - eval = FALSE - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - #帮助(“打通”)# # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - # #裂开,裂开(“LAAGKVEDSD”, enzym =“胰蛋白酶”)# #把乳沟范围cleavageRanges (“LAAGKVEDSD”, enzym =“胰蛋白酶”)# #只有乳沟网站cleavageSites (“LAAGKVEDSD”, enzym =“胰蛋白酶”)# # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - # #小姐一个乳沟位置裂开(“LAAGKVEDSD enzym =“胰蛋白酶missedCleavages = 1) cleavageRanges (“LAAGKVEDSD enzym =“胰蛋白酶missedCleavages = 1) # # 0或1小姐乳沟位置裂开(“LAAGKVEDSD enzym =“胰蛋白酶missedCleavages = 0:1) cleavageRanges (“LAAGKVEDSD enzym =“胰蛋白酶missedCleavages = 0:1) # # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - # #创建AAStringSet对象p < AAStringSet (c (gaju =“LAAGKVEDSD pnm =“AGEPKLDAGV”)) # #裂开,裂开(p, enzym =“胰蛋白酶”)cleavageRanges (p, enzym =“胰蛋白酶”)cleavageSites (p, enzym =“胰蛋白酶”)# # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - # # UniProt负载。ws库库(UniProt.ws) # #选择物种智人< - UniProt.ws (taxId = 9606) # #下载序列的胰岛素/生长激素抑制素s < -选择(= c,键(“P01308”、“P61278”),列= c(“序列”),keytype =“UniProtKB”) # # < - setNames只获取序列序列(序列,新元新元条目)# # < - gsub去除序列(模式=“[[:空间:]]”,替换= " x =序列)# # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -裂开(序列,enzym =“胃蛋白酶”)# # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - # #加载大脑库库(“大脑”)# #裂开胰岛素cleavedInsulin < -裂开(序列[1],enzym =“胰蛋白酶”)[[1]]# #建立空图形区情节(NA xlim = c (150、4300), ylim = c (0, 1), xlab =“质量”,ylab =“相对强度”,主要=“胰蛋白酶的消化胰岛素-同位素分布”)# #循环通过肽(我在seq (= cleavedInsulin)){# #计算c、H、N, O,在当前肽s原子原子< -大脑::getAtomsFromSeq (cleavedInsulin[[我]])# #计算同位素分布d <——useBRAIN(原子)# #画峰线(质量,美元d isoDistr美元,类型=“H”,上校= 2)}# #——sessioninfo,回声= FALSE - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - sessioninfo ()