С помощью команды программы einverted из пакета EMBOSS (einverted -sequence 1EIY.seq -gap 12 -threshold 10 -match 3 -mismatch -3 -outfile outfile -outseq seqout) ищу инвертированные участки нуклеотидных последовательностей. При заданных параметрах находился один стебель, поэтому я уменьшала значение threshold, но даже при значении 0 оставался всего один стебель:
Score 21: 7/7 (100%) matches, 0 gaps
1 gccgagg 7
72 cggctcc 66
С помощью алгоритма Зукера построила возможную вторичную сткрутру тРНК.
Последовательность команд для демонстрации разных множеств в PyMOL (после использования скрипта выше):
hide all
show lines, chain D+E
show spheres, set1
hide spheres, set1
show spheres, set2
hide spheres, set2
show spheres, set3
Контакты атомов белка с | Полярные | Неполярные | Всего |
---|---|---|---|
Остатками 2'-дезоксирибозы | 7 | 30 | 37 |
Остатками фосфорной кислоты | 11 | 22 | 33 |
Остатками азотистых оснований со стороны большой бороздки | 10 | 15 | 25 |
Остатками азотистых оснований со стороны малой бороздки | 5 | 6 | 11 |
Табл. 2. Контакты разного типа в комплексе 1hw2.pdb.
На основе данных таблицы видно, что наибольшее количество контактов наблюдается у остатков дезоксирибозы и фосфорной кислоты. Они находятся снаружи относительно азотистых оснований, поэтому белок взаимодействует с ними больше, чем с основанями. Также видно, что больше контактов азотистых оснований со стороны большой бороздки, чем со стороны малой, что ожидаемо при увеличении длины бороздки.
получение схемы ДНК-белковых контактов с помощью программы nucplot
Применила следующие команды:
1) wget "https://files.rcsb.org/download/1hw2.pdb" -O "1hw2.pdb"
2) remediator --old ''1hw2.pdb'' > ''1hw2_old.pdb
3) nucplot 1hw2_old.pdb
Среди взаимодействующих аминокислот наибольшее число контактов с ДНК имеет 49 аргинин цепи B белка, что отображено на Рис. 5.
Наиболее важным для рапознавания ДНК является 65 His цепи А белка, взаимодействующий с 16 аденином цепи D и остатком дезоксирибозы 17 гуанина цепи D. Аминокислотный остаток полностью проникает в ДНК.