蛋白质分子量
蛋白质是由许多氨基酸残基通过肽键（一个氨基酸的α-羧基与另一个氨基酸的α-氨基脱水缩合形成的化学键）连接而成。蛋白质的分子量（molecular weight）为各个氨基酸的分子量之和，是蛋白质的重要理化参数。单位为Dalton（道尔顿，缩写Da，或D），定义为碳12原子质量的1/12，1D=1/N g，N为阿弗加德罗常数。通常，我们可以使用氨基酸残基数*110大概计算蛋白质的分子量。
目录
一、四个蛋白质分子量计算网站
1）UniProt数据库
2）Expasy
3) EMBOSS
4) PIR
二、python 代码版
1）Biopython包
2）python代码从头计算
附1：Average vs monoisoform.. 8
附2：氨基酸分子量表格... 8
一、四个蛋白质分子量计算网站
1）UniProt数据库
打开https://www.uniprot.org，在输入框输入“P05130”，然后点击“Search”按钮。
点击结果页面左侧的“Sequence & Isoform”，链接到序列处。此处有长度679，质量77,695 Da。粗略等于679*110。
点击“Download”下载序列备用。
如果仅有几条蛋白质，并且有对应的UniProtKB id，建议直接检索获得结果。
如果有N条蛋白质，并且有对应的UniProtKB id，可以使用UniProt的idmapping工具检索，批量下载结果，然后进行解析以获得分子量。
2）Expasy
打开https://web.expasy.org/protparam/，在上面的输入框中输入蛋白质id或者序列id，例如P05130或者KPC1_DROME；或者在下面的输入框中输入蛋白质单字母序列。
点击“compute parameters”提交，弹出来参数，全部默认，点击“Submit”按钮。
分子量就计算出来了。
也可以输入我们在UniProt上获得的序列（这里仅序列）计算蛋白质分子量。略，请大家自行尝试。
3）EMBOSS
打开https://www.ebi.ac.uk/Tools/seqstats/emboss_pepstats/，在输入框中输入在UniProtKB上获得的序列（fasta格式），点击“Submit”按钮提交。
经过约10秒后，会返回结果。
也可以一次性提交N条序列（fasta格式）。略，请大家自行尝试。
4）PIR
打开https://proteininformationresource.org/pirwww/search/comp_mw.shtml，在上面的输入框输入UniProtid或者在下面的输入框中输入在UniProt上获得的序列，点击“Submit”按钮。
约1秒后返回结果。红框中为计算所用的公式和每个氨基酸的分子量（带水分子的分子量）。
也可以一次性提交N条序列（换行分割）。略，请大家自行尝试。
总结：
二、python代码版
1）Bipython包
2）python代码从头计算
虽然Biopython 3行代码就可以搞定蛋白质分子量计算，但是我们需要知道计算原理。最重要的就是获得氨基酸对应的分子量表格。
代码解释：
输入1：单字母的氨基酸序列
输入2：每个氨基酸对应的分子量表格
原理：根据每个氨基酸对应的分子量表格，将所有字母对应的分子量数值加起来，然后再加上水分子的分子量，就是最终的分子量。
注意：这里的分子量是不带水分子的，如果你用的表格是带水分子的，需要用PIR结果中的计算公式。
总结：1）网站上使用的基本都是平均质量。
2）使用上述6种方法计算出来的P05130的分子量基本一样，小数点细微差别可能是由于精度不同或者使用的分子量表格（见附2）不同。
附1：Average vs monoisoform
蛋白质/化合物的平均质量(average mass，红线处)，是由其组成的元素质量的加和，每种元素的质量选取所有同位素的平均质量。
蛋白质/化合物的单同位素质量（monoisoform，最高峰处），是由其组成元素质量的加和，每种元素质量选择其最大丰度同位素的质量。
绝大多数情况下，平均质量>单同位素质量（见附2）。
附2：氨基酸分子量表格
网友评论