计算蛋白质分子量大小（MW）、等电点（pI）和消光系数

以下代码由PHP站长网 52php.cn收集自互联网

现在PHP站长网小编把它分享给大家，仅供参考

# -*- coding: UTF-8 -*-

AA="ACDEFGHIKLMNPQRSTVWY"
AA_residue={'A': 71.0788,'C': 103.1388,'D': 115.0886,'E': 129.1155,'F': 147.1766,'G': 57.0519,'H': 137.1411,'I': 113.1594,'K': 128.1741,'L': 113.1594,'M': 131.1926,'N': 114.1038,'P': 97.1167,'Q': 128.1307,'R': 156.1875,'S': 87.0782,'T': 101.1051,'V': 99.1326,'W': 186.2132,'Y': 163.176}
pI_e={'C': 9.0,'D': 4.0,'E': 4.5,'H': 6.4,'K': 10.4,'R': 12.0,'Y': 10.0}
COOH="CDEY"
NH2="HKR"


P=raw_input(u"请输入蛋白序列:n",)

number={k:P.count(k) for k in AA}

a=0
for k in AA:
    a += number[k]*AA_residue[k]

E=(number['Y']*1490+number['W']*5500+number['C']/2*125)/a

def f(x) :
    b=0
    c=0
    for m in COOH:
        b += (number[m]*(10**x))/(10**x+10**pI_e[m])
    for n in NH2:
        c += (number[n]*10**pI_e[n])/(10**x+10**pI_e[n])
    return b+10**x/(10**x+10**3.2)-10**8.2/(10**x+10**8.2)-c
    
r=3.2
s=12.0
x=(r+s)/2
for i in range(11) :
    if f(x)>0:
        s=x
        x=(r+s)/2
    elif f(x)<0:
        r=x
        x=(r+s)/2


print(u"分子量大小为: "+str(round((a+18.01524)/1000,2))+u" kD    pI值为: "+str(round(x,2)))
print(u"消光系数："+str(round(E,2)))

以上内容由PHP站长网【52php.cn】收集整理供大家参考研究

如果以上内容对您有帮助，欢迎收藏、点赞、推荐、分享。

（编辑：李大同）

【声明】本站内容均来自网络，其相关言论仅代表作者个人观点，不代表本站立场。若无意侵犯到您的权利，请及时与联系站长删除相关内容!