分子量分解问题的优化设计.doc

分子量分解法研究 杨艳林 周旭斌 （河海大学理学院） 摘要：没有蛋白质就没有生命。蛋白质是由C、H、O、N、P、S等元素组成的一类高分子化合物。氨基酸为其主要的组成物质，蛋白质水解后亦生成氨基酸，可以说氨基酸是蛋白质的基石。研究蛋白质的组成，最重要的就是研究其是由那些氨基酸组成的。对于给定的18中氨基酸，以及分子量X，我们利用循环的思想，对于所有可能的情况加以考虑。当有计算机时，使用计算机穷尽所有的可能组合计算得出结果；当没有计算机时，将问题转换成求多元一次方程解的问题，通过一系列的矩阵变化及限定条件便可将解一一得出。经过计算，我们发现当X=500时，解的个数N=158个；当X=1000时，N=28268；而当X=2000时，N从画出的X与N的关系图中，可以很简单的看出X与N呈现指数次幂的增长关系；而在X与程序运行时间T的关系图中，亦显示两者呈指数关系（具体表达式由下文给出）。 关键词：蛋白质 循环 n元一次多项式 问题重述 众多周知，蛋白质分子的基本组成单位是氨基酸，各种氨基酸分子经过不同的结合方式形成了各种各样的生命蛋白质分子。测量某种物质的分子量，是我们在研究化学物质的一种基本手段。在实验室中，我们通常质谱实验测定某蛋白质分子量x (正整数)，然后将分子量x分解为n个已知分子量a[i](i=1,.......,n)氨基酸的和的形式。现在某实验室所研究的问题中： n=18, x1000 a[i](i=1,.......,18)分别为57, 71, 87, 97, 99, 101, 103, 113, 114, 115, 128, 129, 131, 137, 147, 156, 163, 186 要求针对该实验室拥有或不拥有计算机的情况作出解答。 问题分析 本问题已知蛋白质分子以及18种基本氨基酸分子的分子量，需要计算得到蛋白质分子由这18种分子的组成安排。组成中，可以有多种组合，可以有某些氨基酸分子大量重复使用，亦可以有一些氨基酸分子不用到。将每个氨基酸分子使用的次数乘上其对应分子量若是恰好等于给出值X，那便得到了蛋白质分子的一种组合，穷尽所有可能组合，便得到了最终的结果。 模型假设 1）蛋白质分子在由氨基酸分子的结合过程中，由于产生肽腱所失去的水分子分子量忽略不计。 2）蛋白质分子由氨基酸组合中产生的一些特定腱不单独考虑，而计算在单独氨基酸中。 3）忽略蛋白质分子在产生过程中参加反应的其他分子。 模型建立 在假设的情况下，对于给定个分子量x，将按题给顺序的氨基酸分子分别编号1--18，假设每种氨基酸分子分子量分别为（=1、2、3.....18）（最大公因子是1）分别有（i=1、2、3.....18）种，根据我们的基本设计模型： 解出所有的bi即可得到所有组合。 以下即为在实验室有计算机和无计算机两种情况下的算法解法。在有计算机时，主要运用“循环”的思想，考虑所有可能情况以得出解；在无计算机时，通过矩阵变化和求出一个特解及其参数，即可表示得到所有的解。 在有计算机的情况下求解 在有计算机时，我们认为使用VB语言来编程可以形成一个应用程序非常方便运行，故此我们采用VB来编程。我们通过计算机强大的计算能力，通过“循环”的思想：先假象确定第i种氨基酸分子的个数，然后对第i+1种氨基酸的所有可能情况分别进行计算，得到多有可能性，当其中满足基本模型的组合即为所要求的解，最后输出解的情况、个数及运行时间即可。 以下是我们根据不同的x所求出的不同解的个数及运行时间情况所制成的一张表格： x的取值 解的个数N 运行时间T 50 0 0 100 0 0 150 0 0 200 4 0 250 3 0 300 14 0 350 7 0 400 45 0 450 44 0 500 158 0 550 202 0 600 522 0 650 694 0 700 1508 0.01 750 2197 0.03 800 4291 0.05 850 6337 0.07 900 11249 0.09 950 16885 0.11 1000 28268 0.16 1100 67339 0.38 1200 154143 0.82 1300 338158 1.85 1400 716481 4.02 1500 1467221 8.39 1600 2915738 17.19 1700 5633990 34.08 180065.69 1900123.76 2000228.3 当X不断变大时，X与解的个数N所呈现的大体态势如图（1）： 图（1） X与N的变化关系图 由上图可以看出，随着X的增大