Модели эволюции нуклеотидных и аминокислотных последовательностей BLOSUM 62

Модели эволюции нуклеотидных и аминокислотных последовательностей. BLOSUM 62 Matrix AAA AAG AAA AGA GGA AGA

Допущения при построении филогенетического дерева ветвь внутренний узел лист • Топология T : бифуркационная – Листья - 1…N – Внутренние узлы N+1 … 2 N-2 • Длины t = { ti } для каждой ветви • Филогенетическое дерево = (Топология, Длины) = (T, t )

Метод максимального правдоподобия Пусть филогенетическое дерево представляет собой стохастический процесс AAA AAG AAA AGA GGA AGA Ненаблюдаемые Наблюдаемые Вероятность перехода от символа a к символу b равна b|a. Вероятность найти символ a в корне равна qa. Эти параметры определяются через скорости изменения за единицу времени, умноженное на время. Имея полное дерево, вероятность данных определяются через значения b|a и qa’.

Вероятностные методы • Филогенетическое дерево представляет собой обобщенную вероятностную модель для наблюдаемых последовательностей. • Корневые вероятности: q( a ) • Вероятности мутаций: P( a | b, t ) • Модели мутаций – Модель Джукса-Кантора (Jukes Cantor, 1969) JC 69 – Двупараpsqметрическая модель Кимуры (Kimura, 1980) K 2 P – Модель Хасегава-Кишино-Яно (Hasegawa, Kishino and Yano, 1985) HKY 85

Накопление замен • Две последовательности, происшедшие от общего предка. D - процент сайтов, которые различаются между собой Предполагается, что замены происходят случайным образом с постоянной скоростью Число видимых замен меньше или равно числу реальных замен

Отношение D и t • На малых временах D пропорционально t • На больших временах D растет медленнее, чем t • Очень большое время -> очень много замен -> получились случайные последовательности • Выравненные случайные последовательности с равновероятным распределением нуклеотидов -> 3/4 сайтов будут различными • D=3/4 при t -> бесконечности • Почему D плоха для использования? üнелинейная зависимость от t üнеаддитивна (D 12 не равно D 01+D 02)

Что необходимо? • Определить эволюционное расстояние, которое было бы аддитивно и линейно зависило от времени • Пусть d - среднее число замен в каждом сайте между двумя последовательностями • Если замены происходят случайным образом с постоянной скоростью, то d пропорционально t • d аддитивно, так как d 12=d 01+d 02 по определению • Но d не подлежит наблюдению, в отличие от D • Чтобы оценить d нужна эволюционная модель

Оценка эволюционного расстояния d • Пусть скорость замены одного нуклеотида. • 3 скорость замены одного нуклеотида на любой другой • t - время, прошедшее после расхождения двух последовательностей • d=3 t=6 t

Модель Джукса-Кантора (1969) • Модель для скоростей мутации • Мутации происходят с постоянной скоростью • Каждый нуклеотид РАВНОВЕРОЯТНО может перейти в любой другой со скоростью .

Модель Джукса-Кантора Нужно получить формулу эволюции ДНК через вероятности Prob(y | x, t) , где x и y взяты из {A, C, G, T} и t временной интервал. Модель Джукса-Кантора подразумевает равные скорости переходов нуклеотидов: A C -3 G T

Модель Джукса-Кантора Пусть S(t) - матрица вероятностей переходов (транзиций): Мы предполагаем, что матрица мультипликативна, то есть S(t+s)=S(t)S(s) для любых интревалов времени s или t.

Модель Джукса-Кантора Для короткого интервала времени : Из свойства мультипликативности следует: S(t+ ) = S(t) S( ) S(t)(I+R ) Отсюда: [S(t+ ) - S(t)] / S(t) R Приходим к линейному дифференциальному уравнению: S’ (t) S(t)R С дополнительным условием, что в пределе при t, стремящемся к бесконечности:

Модель Джукса-Кантора Подставляя S(t) в дифференциальное уравнение, получаем систему: Система имеет единственное решение, известное под названием модели Джукса-Кантора:

Модель Джукса-Кантора вероятность мутации 1 -3 вероятность отсутствия мутации Процент идентичных сайтов - i(t); процент неидентичных сайтов - n(t): i(t)+n(t)=1 За период времени (t+1), каждый идентичный сайт остается без изменения с вероятностью 1 -3 Вероятность, что обe последовательности не мутировали (1 -3 )^2=1 -6 ^2~1 -6 <<1 Вероятность, что неидентичные сайты мутируют обратно к идентичному состоянию =2 n(t) i(t+1)=i(t)(1 -6 2 n(t)=i(t)(1 -6 2 (1 -i(t)) i(t+1)-i(t)= 2 i(t) =1/4(1+3 exp(-8 t)) n(t)=1 -i(t) d=6 t

Зависимость D и d от t

Двупараметрическая модель Кимуры • Учитывает разные скорости для транзиций и трансверсий.

Модель Кимуры K 2 P (1980) Модель Джукса-Кантора не принимает во внимание, что скорости транзиций (между пуринами) A G и (между пиримидинами) C T отличаются от скоростей трансверсий A C, A T, C G, G T. Кимура предложил следующую матрицу скоростей:

Модель Кимуры K 2 P (1980) Матрицу скоростей S(t) можно проинтегрировать, так же как и для модели Джукса-Кантора, и получится общая, зависящая от времени форма где И эта модель еще далека от реальности, так как частоты нуклеотидов равны между собой: q. A=q. C=q. G=q. T=1/4 У многих организмов AT и GC сильно сдвинуто.

Модель Кимуры K 2 P • S - процент неидентичных сайтов, получившихся в результате транзиции • V - процент неидентичных сайтов, получившихся в результате трансверсии • D=S+V • Аналогично подходу в модели Джукса-Кантора D=0. 3 d(JC 69)=0. 383 S=0. 2 d(K 2 P)=0. 402 V=0. 1

Модель Хасегава-Кишино-Яно HKY 85 (Hasegawa-Kishino-Yano, 85) • Частоты встречаемости нуклеотидов неодинаковы -

Модели эволюции аминокислотных последовательностей. • Аминокислотные последовательности (20 аминокислот) • Протеин-кодирующие последовательности ДНК

Модели эволюции аа последовательностей 1. Эмпирические матрицы 20 x 20 – – – Dayhoff et al. (1978) матрица основана на наблюдениях 1572 фиксированных мутаций в 34 суперсемействах достаточно близких последовательностях JTT матрица (Jones et al. 1992; Gonnett et al. 1992): методология такая же как и у Dayhoff, но данные взяты из современных баз данных (позднее были сделаны модификации для трансмембранных белков Jones et al. 1994) mt. REV (Adachi and Hasegawa 1995, 1996) matrix derived from maximum likelihood-inferred replacements in mitochondrial proteins of 20 vertebrate species WAG (Whelan and Goldman 2001) matrix derived from maximum likelihood improvement of JTT Poisson assumes equal stationary state frequencies and equal substitution rates (equivalent to JC model for DNA). Not really empirical, but it is fixed

Матрица счета PAM 250

Матрицы счета S = [sij] дает число выравнивания символа i с символом j для каждой пары i, j. C 12 S 0 T -2 1 3 P -3 1 0 6 A -2 1 1 1 2 C T P A STPP CTCA 2 S 0 + 3 + (-3) + 1 =1

Типы матриц • • • PAM BLOSUM Gonnet JTT DNA matrices матрицы PAM, Gonnet, JTT, и DNA PAM основаны на явных эволюционных моделях; • Матрицы BLOSUM основаны на неявных эволюционных моделях

Матрицы PAM основаны на простой эволюционной модели GAATC GAGTT Последовательность предка? Два изменения GA(A/G)T(C/T) • Разрешены только мутации • Сайты эволюционируют независимо

Матрицы замен PAM • Point Accepted Mutation: точечные мутации, которые распространяются и фиксируются. • Единица эволюционного изменения для белковых последовательностей [Dayhoff 78]. • Единица PAM - это такая единица эволюции, при которой в среднем изменяется 1% аминокислот внутри белковых последовательностей.

Матрица замен PAM • Первая PAM Dayhoff et al. 1978 • 71 семейство • 15 % - расхождение аа последовательностей • 1572 замены • Дерево построено методом максимальной бережливости. Пример для 6 аа

Число фиксированных точечных мутаций 1572 мутации A ij

Определение PAM • Пусть P(i, j) - вероятность находиться в состоянии j во время t, при условии, что в при t=0 сайт находился в состоянии i. • Определим матрицу PAM 1 как где t - время одной единицы PAM Ni - общее количество данной аминокислоты

Итак • Пусть частота аминокислот =N /N i i Atot - общее число элементов в матрице A Вероятность, что данная аминокислота осталась неизмененной за единицу времени PAM 1. total

PAM 1 Значения умножены на 100 000, для удобства. Каждый ряд суммируется в 100 000. m - способность к мутациям (вероятность мутировать >1% или <1%)

Экстраполяция модели на большие растояния -- PAM 1 матрица вероятностей мутаций. -- PAM 2 матрица вероятностей мутаций? -- Мутации, происшедшие за эволюционное расстояние в две единицы PAM 1

В два интервала PAM 1: • {A R} = {A A and A R} or {A N and N R} or {A D and D R} or … or {A V and V R}

Матрица вероятностей замен PAM 2

Матрица замен PAM-k

Матрица весов PAM-k

Матрица весов PAM 250

Соответствия PAM-k проценту схожести белков • PAM 60— 60%, PAM 80— 50%, • PAM 120— 40% • матрица PAM-250 дает лучший вес выравнивания по сравнению с более низкими версиями PAM для белков со схожестью 14 -27%

Построение матриц BLOSUM Blocks Substitution Matrices

Матрицы BLOSUM заданной схожести • отбираются последовательности выше порогового значения. • Если схожесть равна 62%, окончательная матрица называется BLOSUM 62

Так считали PAM Так считают BLOSUM Считаем пары выравненных аминокислот, а не число замен Пример для 6 аа

Пример построения матрицы BLOSUM на основе 4 последовательностей

Построение матрицы BLOSUM. 1. Подсчет пар

2. Cуммирование частот мутаций

3. Матрица мутационных вероятностей.

4. Посчитать избыточность каждого символа (маргинальные вероятности)

5. Получаем матрицу BLOSUM

Построение реальной матрицы BLOSUM 62

1. 2. 3. Таблица частот мутаций

4. Расчет избыточности аминокислот

5. Получение матрицы BLOSUM 62

Оригинальная публикация матрицы BLOSUM • S. Henikoff and J. Henikoff (1992). “Amino acid substitution matrices from protein blocks”. PNAS 89: 10915 -10919 • Данные для обучения: ~2000 консервативных блоков из базы данных BLOCKS. Выравненные белковые участки без пробелов. Каждый блок представляет консервативный участок семейства белка.

Выбор матрицы счета

Сравнение матриц счета BLOSUM • Основана на последовательностях с разными эволюционными периодами. • Консервативные блоки • Поиск консервативных доменов PAM • Основана на экстраполяции последовательностей с малым эволюционным периодом • Отслеживает эволюционное происхождение • Гомологичные последовательности