В файле fullANDpart.fasta мной сохранены в формате FASTA небольшой фрагмент, а также полная последовательность заданного мне белка RECQ_Ecoli. Затем движением фрагмента относительно полной последовательности, с использованием программы GeneDoc я добилась полного совпадения букв. Заданный фрагмент соответствует позициям 91-115 в полной последовательности белка Recq_Ecoli.
Данное выравнивание сохранено в файле alignment1.msf, а картинка с данным выравниванием - aln1.gif
Скопировав обе короткие последовательности в файл shortseqs.fasta, я импортировала этот в файл в GeneDoc. Нужно выровнять вручную последовательности так, чтобы было сопоставлено максимальное число одинаковых букв при минимальном числе пропусков. С этой целью необходимо оценить вес выравнивания, для этого используем формулу: W = M - nG , где где M - число совпавших букв, G - штраф за пропуск, равен 2, n - общее число пропусков.
Итак, число а.о. в первом фрагменте равно 25, во втором - 26; длина выравнивания - 28; вес выравнивания: W= 12 - 2*3=6; процент идентичности двух выравненных последовательностей (отношение числа колонок, в которых стоят одинаковые буквы, к общему числу колонок, включая "гэповые", умноженное на 100): (12/28)*100%=42,86%.
Затем я получила картинку с этим выравниванием aln2.gif в рабочей директории).
Рассмотрим, как устроена наиболее популярная матрица весов замен аминокислотных остатков BLOSUM62. Будем считать близкородственными заменами те, для которых значение элемента матрицы положительно. Позиции выравнивания в данном фрагменте, в которых я обнаружила близкородственную замену аминокислотных остатков:
- № первой позиции выравнивания, в которой мы наблюдаем близкородственную замену а.о:9
- полные названия и однобуквенные обозначения поменявшихся а.о: Threonine, Thr, T заменяется на Serine, Ser, S.
- вес такой замены в соответствии с использованной матрицей:1
- комментарий: треонин - маленькая, гидрофобная, полярная аминокислота; серин - маленькая, полярная аминокислота; обе аминокислоты маленькие, полярные, поэтому вес замены положителен.
- Определить % сходства выровненных фрагментов.
- T --> S ( треонин заменяется на серин), из предыдущего задания мы знаем, что вес этой замены равен +1.
- E --> Q (глутаминовая кислота заменяется на глутамин), вес замены равен +2, что можно объяснить тем, что заменяется кислота на свой же амид.
- T --> S (+1)
- Построить иное выравнивание 2-х заданных фрагментов с весом не хуже, чем получилось в обязательном упр.2
- Придумать, как с помощью функций Excel из матрицы BLOSUM можно получить список замен вида "W-Y вес=...".
Дополнительные задания.
Процентом сходства будем считать отношение числа колонок со сходными буквами к общему числу, умноженное на 100%. Сходными буквами будем считать такие, для которых значение элемента матрицы BLOSUM62 положительно.
Итак, обратимся к полученному ранее выравниванию:
-
Из 28 колонок совпало 12 пар, а среди оставшихся 16 колонок сходными являются только 3, так как их замены, согласно таблице, положительны:
Так получается, что процент сходства можно вычислить по формуле: (12+3)/28*100%=53,57%.
Получилось несколько иное выравнивание с таким же весом, процентом идентичности и процентом сходства. Выравнивание представлено на рисунке aln2_2.gif:
Для начала экспортируем табллицу Blosum62 в Excel. В той же книге созданим лист, на котором в 1-м и во 2-м столбцах расположим аминокислоты из всех колонок моего выравнивания, включая гэповые. Затем в 3-м столбце воспользуемся формулой (ВПР и ПОИСКПО). Таким образом, мы получили вес для каждой замены. Полученную таблицу можно просмотреть здесь: Weight.xls