蛋白质氨基酸的组合问题=.pptVIP

蛋白质氨基酸的组合问题 一.问题的提出 生命蛋白质是由若干种氨基酸的不同组合构成的. 18种氨基酸的已知分子量ai (i=1,2,3,4,5……..)分别如下: 57,71,87,97,99,101,103,113,114,115,128,129,131,137,147,156,163,186. 给定某一蛋白质的分子量下X(X≤1000且为正整数)设计出数学模型给出该蛋白质的所有可能的组成。即确定该蛋白质是哪几种氨基酸组成及每种氨基酸的数目. 二.问题的分析 根据给定的分子量X及ai 测定蛋白质的组成,实际是求多元线性方程: ∑ai xi =X 的所有整数解的问题. 一般采用枚举法求解, 组合共有 ∏([X/ ai ]+1) 种. 计算量大,耗费时间长 给出的解的个数过多反而失去了解的意义. 模型的设计和改进：围绕着减少运算和缩小解的范围的思路展开 三.模型假设. 1.给定的蛋白质分子量 ---和氨基酸已知分子量---是准确的,没有测试误差; 2.假设所有被测定的蛋白质均由给定分子量的这几种氨基酸构成,而不含有其他种类的氨基酸.实际中,构成生命蛋白质的主要氨基酸有20种----,其中两对氨基酸的分子量相等; 3.假设蛋白质分子是构成过程中,各个氨基酸分子之间相互结合的方式不影响蛋白质的分子量. 4.假设被测定的蛋白质所含氨基酸的个数=2,即x=114. 5.假设氨基酸分子结合构成过程中是任意排列组合的,不存在互斥和互补现象,即任何两种氨基酸都可以同时纯在同一个蛋白质中,没有任何一种氨基酸的存在是以其他氨基酸的纯在为前提的. 6.假设在蛋白质中,每种氨基酸存在的概率是相等的,不存在某种必须存在的氨基酸. 7.假设该试验拥有测定化学性质的仪器. 四.最一般的模型 在没有任何其他补充信息和约束条件的情况下,最一般的模型可以表示为 ∑ ai x i =X ； x i为非负整数，i =1,2,3……n； 采用了深度优先算法,遍历了整个解空间,由于采用了分支限定界法,其实际最坏的时间效率也是远小于∏([X/ai]+1)的. 在化学中,我们知道,生命蛋白质氮的含量约占总量的16%左右(其波动范围15%--17%).蛋白质测定的凯式定氮法---就是利用了这个性质.我们考虑了含氮量的模型 模型A 对该模型有两点说明: 1.常见20的种氨基酸中:有两对的分子量相等.其中亮氨酸与异亮氨酸为同分异构,分质量与分子式均相同,因而不会影响该模型的计算.而另一对谷酰氨酸与赖氨酸仅仅是分子量相同而分子式不同.因此在模型中,把含硫的两种氨基酸作特殊处理后,还剩下16种分子量不同,然后加入一个变量,用以区分谷酰氨酸与赖氨酸.最后将结果合并. 模型B 已知蛋白质中某些氨基酸是存在的. 在实际的蛋白质的一级结构测定--中,通常可以对蛋白质经过充分水解后所得到的氨基酸混合液作离子交换层析.纸层析或薄层层析,定性研究的结果可以确定该蛋白质所含的全部或部分氨基酸种类. 假设b 已知被测定的蛋白质中肯定含有其中的k种氨基酸,其分子量为a j(j=1,2,3…k). 可假设 X‵=X-∑a j 已经存在的氨基酸,其对应的个数 x’ j= xj-1 则模型表为 ∑aix’ i=X’ ； x’i为非负正整数（i=1,2,..18）. 模型C： 已知蛋白质中只含有某几种氨基酸. 在比较成功的氨基酸定性分析中,可以得到被测定的蛋白质完全水解生成的氨基酸的全部种类 假设c 假设某蛋白质由且仅由k种已知的氨基酸构成. 对应的xi=1,这样我们可以令X‵=X-∑ai . 模型可表示为 ∑ki=1 aix’i=X’ ； x’i为非负正整数（i =1,2,..k）. 模型D 18种已知氨基酸分子量的平均值为118.5,因而平均来看对于X1000的蛋白质来说其所含氨基酸的分子数在8-9之间. 不妨设:每种氨基酸分之的数目仅为0或1,因而模型表示为 ∑aixi=X xi= 0或 1； 模型E 如实验室不拥有微机, 但可能拥有较为先进的化学分析设备,设实验室对完全水解后的氨基酸混合液作定性的分析,并可以通过质谱仪侧得蛋白质的分子式. 设构成被测氨基酸的构成分别为第==种,则模型可以