WikiDer > ProbCons

ProbCons представляет собой вероятностную согласованность с открытым исходным кодом, основанную на множественном согласовании аминокислота последовательности. Это один из самых эффективных белков. множественное выравнивание последовательностей программ, поскольку он неоднократно демонстрировал статистически значимое преимущество в точности перед аналогичными инструментами, в том числе Clustal и MAFFT.^[1][2]

Алгоритм

Ниже описана основная схема алгоритма ProbCons.^[3]

Шаг 1. Надежность выравнивающей кромки

Для каждой пары последовательностей вычислите вероятность того, что буквы ${displaystyle x_ {i}}$ и ${displaystyle y_ {i}}$ в паре в ${displaystyle a ^ {*}}$ выравнивание, которое создается моделью.

${displaystyle {egin {align} P (x_ {i} sim y_ {i} | x, y) & {stackrel {def} {=}} Pr [x_ {i} sim y_ {i} {ext {в некоторых }} | x, y] & = sum _ {{ext {выравнивание a с}} x_ {i} -y_ {i}} Pr [a | x, y] & = sum _ {ext {alignment a} } mathbf {1} {x_ {i} -y_ {i} in a} Pr [a | x, y] конец {выровнено}}}$

(Где ${displaystyle mathbf {1} {x_ {i} sim y_ {i} in a}}$ равно 1, если ${displaystyle x_ {i}}$ и ${displaystyle y_ {i}}$ находятся в выравнивании и 0 в противном случае.)

Шаг 2: максимальная ожидаемая точность

Точность совмещения ${displaystyle a ^ {*}}$ по отношению к другому раскладу ${displaystyle a}$ определяется как количество общих выровненных пар, деленное на длину более короткой последовательности.

Рассчитайте ожидаемую точность каждой последовательности:

${displaystyle {egin {align} E_ {Pr [a | x, y]} (acc (a ^ {*}, a)) & = sum _ {a} Pr [a | x, y] acc (a ^ { *}, a) & = {frac {1} {min (| x |, | y |)}} cdot sum _ {a} mathbf {1} {x_ {i} sim y_ {i} in a} Pr [a | x, y] & = {frac {1} {min (| x |, | y |)}} сумма cdot _ {x_ {i} -y_ {i}} P (x_ {i} sim y_ {j} | x, y) конец {выровнено}}}$

Это дает выравнивание с максимальной ожидаемой точностью (MEA):

${displaystyle E (x, y) = arg max _ {a ^ {*}}; E_ {Pr [a | x, y]} (acc (a ^ {*}, a))}$

Шаг 3: преобразование вероятностной согласованности

Все пары последовательностей x, y из множества всех последовательностей ${displaystyle {mathcal {S}}}$ теперь переоцениваются с использованием всех промежуточных последовательностей z:

${displaystyle P '(x_ {i} -y_ {i} | x, y) = {frac {1} {| {mathcal {S}} |}} sum _ {z} sum _ {1leq kleq | z |} P (x_ {i} sim z_ {i} | x, z) cdot P (z_ {i} sim y_ {i} | z, y)}$

Этот шаг можно повторить.

Шаг 4: Расчет дерева направляющих

Постройте направляющее дерево путем иерархической кластеризации, используя оценку MEA в качестве оценки сходства последовательностей. Сходство кластера определяется с использованием средневзвешенного значения попарного сходства последовательностей.

Шаг 5: вычислить MSA

Наконец, вычислите MSA, используя прогрессивное или итеративное выравнивание.

Смотрите также

• Программное обеспечение для выравнивания последовательностей
• Clustal
• МЫШЦЫ
• КАРТА
• Т-кофе
• Probalign

Рекомендации

внешняя ссылка

• Официальный веб-сайт