Multalin

前言

在寻找基因和致力于发现新蛋白的努力中，人们习惯于把新的序列同已知功能的蛋白序列作比对。由于这些比对通常都希望能够推测新蛋白的功能，不管它们是双重比对还是多序列比对，都可以回答大量的其它的生物学问题。举例来说，面对一堆搜集的比对序列，人们会研究隐含于蛋白之中的系统发生的关系，以便于更好地理解蛋白的进化。人们并不只是着眼于某一个蛋白，而是研究一个家族中的相关蛋白，看看进化压力和生物秩序如何结合起来创造出新的具有虽然不同但是功能相关的蛋白。研究完多序列比对中的高度保守区域，我们可以对蛋白质的整个结构进行预测，并且猜测这些保守区域对于维持三维结构的重要性。

显然，分析一群相关蛋白质时，很有必要了解比对的正确构成。发展用于多序列比对的程序是一个很有活力的研究领域，绝大多数方法都是基于渐进比对（progressive alignment）的概念。渐进比对的思想依赖于使用者用作比对的蛋白质序列之间确实存在的生物学上的或者更准确地说是系统发生学上的相互关联。不同算法从不同方面解决这一问题，但是当比对的序列大大地超过两个时（双重比对），对于计算的挑战就会很令人生畏。在实际操作中，算法会在计算速度和获得最佳比对之间寻求平衡，常常会接受足够相近的比对。不管最终使用的是什么方法，使用者都必须审视结果的比对，因为再次基础上作一些手工修改是十分必要的，尤其是对保守的区域。

方法原理

MultAlin方法也是基于用一系列双重比对开始的思想,然后基于双重比对的打分值进行一个分层次的聚类。当序列都分成类后，开始进行多序列比对，计算出多序列比对中的两个序列比对的新值，基于这些新值，重新构建一棵树。这个过程不断进行，直到分值不在上升，此时所序列比对也就结束了。

MultAlin可以在INRA Toulouse的一个环球网点上很容易地执行，要比对的序列按照FASTA的格式被粘贴到一个序列输入框内，然后从一系列下拉菜单中，用户定义适当的参数，比如输出格式，可选的输入格式，引用的分值矩阵以及空位开放和扩展罚分的分值。大多数用户只会根据输入序列的远近关系，选择不同的分值矩阵。然后，序列被提交到服务器上，当多序列比对返回时，会计算一个一致序列并显示在比对的下方。举例来说，如图8.1所示的用CLUSTAL W比对的同样的序列被提交给MultAlin服务器，接受缺省的比对参数，其结果如图8.2所示，在一致序列中，所有序列都匹配的残基相应的位置用此残基的大写字母表示，大多数都匹配的用小写字母表示，同样地，符号！、$、%和#分别表示保守取代，具体含义如图8.2上方的图例。

很明显，用两种方法分别得到的比对结果并不完全一样。主要区别在于CLUSTAL W在果蝇序列中开放了两个长度超过10的空位，而MultAlin只开放了一个长空位，而且，MultAlin可以得到比CLUSTAL W多20个完全相同残基的排队，当然，这并不以为这一种方法比另外一种方法优越，这有要重提本书的一个不便的话题，即从输入序列的性质出发，应用不同的方法会得到不同程度的成功。警慎的用户会选择若干个工具同时使用，并且对最终的比对结果作手工修正以期达到最佳效果。