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

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