С помощью инструментов пакета 3DNA были построены модели структур A-, B- и Z-формы дуплекса ДНК, последовательность одного из цепей которого — пять повторов GATC.
Файлы с полученными моделями:
Для построения была использована команда ‹fiber -a-seq=GATCGATCGATCGATCGATC gatc-a.pdb› в случае с А-формой, опция -а заменялась на -b и -z для B- и Z-формы соответственно.
В цепи B экспериментально полученной структуры А-формы ДНК (3V9D) было выбрано 7 основание — тимин, который выделяется скриптом с участком сахаро-фосфатного остова в апплете ниже.
Далее были определены атомы, явно обращённые в сторону большой бороздки (C4, C5, C6, O4), и атомы, явно обращённые в малую (N1, C2, O2). На Изображении 1 представленная нарисованная в MarvinSketch молекула тимина, атомы которой окрашен в зависимости от обращённости к определённой бороздке. Синим покрашены направленные в сторону малой, а красной в сторону большой.
Также были измерены параметры различных форм ДНК из экспериментально полученных структур (1tne для Z-формы, 1bna для B-формы, 3v9d для A-формы). Результаты измерений приведены в таблице ниже.
A-форма | B-форма | Z-форма | |
---|---|---|---|
Тип спирали | правая | правая | левая |
Шаг спирали (Å) | 28.09 | 33.75 | 43.5 |
Число оснований на виток | 11 | 10 | 12 |
Ширина большой бороздки | 7.98(3:A.P 32:B.P) | 17.21 (29:B:P-9:A:P) | 7.2 (33:B.P-11:A.P) |
Ширина малой бороздки | 16.81 (10:A:P-33:B:P) | 13.1 (30:B:P-14:A:P) | 16.08 (28:B.P-12:A.P) |
В этой части практикума работа ведётся с т-РНК 1G59. К сожалению, пакет 3DNA позволяет работать только со старым форматом pdb, поэтому сначала необходимо для конвертации воспользоваться командой remediator ‹--old 1g59.pdb > 1g59_old.pdb›.
Далее с помощью команды ‹find_pair -t 1g59_old.pdb stdout | analyze› были получены файлы, в одном из которых (1g59_old.out) можно найти информацию о торсионных углах и водородных связях.
Ниже представлена таблица созначениями торсионных углов для первой цепи.
Strand I base alpha beta gamma delta epsilon zeta chi 1 G --- -157.0 59.6 98.9 -153.8 -79.8 -171.3 2 G -84.9 -162.0 51.6 96.1 -155.3 -68.9 -154.4 3 C -69.6 -179.4 51.3 92.5 -152.5 -74.0 -148.7 4 C -60.1 176.5 47.0 90.4 -161.9 -83.7 -149.4 5 C -60.4 168.9 51.9 89.1 -157.7 -69.5 -161.6 6 C -66.0 -176.8 55.6 96.2 -161.6 -76.4 -165.7 7 A -73.5 174.6 73.2 128.9 -65.1 -113.2 -144.7 8 G 59.7 151.9 31.3 97.1 -148.8 -79.0 -170.4 9 G -75.8 176.5 63.6 91.3 -160.6 -64.4 -165.4 10 G -63.8 -170.7 41.7 90.7 -162.9 -97.4 -156.2 11 G 130.2 -154.5 -168.9 104.6 -118.0 -26.4 -178.1 12 G -109.6 56.1 178.5 110.9 -115.0 -87.5 167.2 13 U 0.9 155.6 -12.8 97.3 -136.3 -77.9 -152.5 14 U -69.9 156.0 60.4 88.9 -110.2 -75.4 -153.1 15 A -69.8 173.1 54.9 94.0 -146.6 -61.6 -167.3 16 G -69.5 177.9 49.6 93.0 -157.1 -77.9 -156.7 17 G -61.8 178.8 44.4 93.8 -157.3 -74.7 -154.8 18 C -81.4 175.2 64.5 87.3 -159.0 -69.1 -157.2 19 C -70.5 -179.2 55.4 87.5 -157.7 -73.2 -154.6 20 G -66.2 172.9 62.0 87.9 -158.6 -79.5 -164.1 21 A -60.2 -176.6 42.0 88.8 -141.2 -66.9 -145.2 22 G -154.9 159.7 45.6 93.8 -144.5 -73.1 -175.6 23 U -68.1 -171.7 46.2 93.1 -153.8 -79.1 -166.8 24 C -59.7 -176.0 41.7 91.1 -151.4 -71.1 -147.9 25 U 120.0 -154.7 -146.4 92.3 -161.4 -79.7 -170.3 26 U -37.9 -148.7 56.4 96.3 -148.2 -73.8 -140.9 27 A -68.2 -169.7 52.4 91.3 -134.5 -68.0 -162.7 28 G -78.8 162.3 72.4 99.5 24.2 165.6 -147.4 29 G -147.7 -111.4 52.7 121.2 58.6 139.0 -62.9 30 G --- -161.2 55.0 98.3 -155.1 -75.7 -163.9 31 G -78.8 -165.2 44.9 96.0 -151.7 -71.0 -155.8 32 C -65.5 174.8 48.7 92.0 -152.0 -75.1 -146.4 33 C -57.5 175.0 44.6 91.2 -159.7 -84.2 -151.5 34 C -64.4 170.1 55.5 88.1 -158.5 -69.8 -162.5 35 C -64.6 -177.7 55.1 96.5 -162.4 -76.1 -163.7 36 A -75.6 175.3 74.8 130.5 -66.4 -112.1 -145.1 37 G 60.1 151.0 30.0 93.2 -147.9 -78.5 -173.7 38 G -69.3 177.2 61.9 92.2 -158.5 -66.6 -166.9 39 G -61.2 -173.6 40.7 90.2 -159.8 -99.2 -157.4 40 G 126.3 -151.0 -166.9 104.6 -119.8 -27.5 179.1 41 G -110.5 56.4 178.7 112.4 -118.9 -89.0 166.5 42 U 9.3 148.6 -19.7 92.5 -131.1 -75.9 -152.2 43 U -71.9 158.4 68.3 89.9 -114.2 -74.7 -146.5 44 A -67.6 172.6 51.8 92.5 -147.5 -58.4 -167.8 45 G -79.3 -179.5 57.9 90.6 -156.5 -77.8 -157.2 46 G -60.6 177.9 43.0 92.2 -156.5 -75.3 -157.2 47 C -76.6 174.7 60.8 88.5 -159.3 -69.1 -153.4 48 C -68.9 -178.7 54.0 87.1 -156.5 -73.1 -153.3 49 G -64.2 171.0 62.2 89.1 -155.9 -81.0 -163.1 50 A -59.9 -178.7 41.1 87.0 -141.1 -67.4 -146.0 51 G -163.4 166.4 48.1 97.2 -147.7 -72.6 -174.0 52 U -66.5 -172.1 46.0 93.2 -153.5 -78.6 -159.6
α | β | γ | δ | ε | ζ | χ | |
---|---|---|---|---|---|---|---|
A-форма (теор.) | 62 | 173 | 52 | 88/3 | 178 | -50 | -160 |
B-форма (теор.) | 63 | 171 | 54 | 123/131 | 155 | -90 | -117 |
Из отрицательных значений углов α и β цепей РНК очевидно, что сравнение по средним значениям этих углов с формами ДНК не представляется возможным. Однако можно заметить, что модули этих углов зачастую достаточно близки к таковым у A- и B-форм ДНК.
Углы дельта и хи ближе к A-форме, а эпсилон по модулю к B-форме.
Ниже представлена таблица с информацией о водородных связях.
RMSD of the bases (----- for WC bp, + for isolated bp, x for helix change) Strand I Strand II Helix 1 (0.011) ....>B:.501_:[..G]G-----C[..C]:.572_:B<.... (0.010) | 2 (0.008) ....>B:.502_:[..G]G-*---U[..U]:.571_:B<.... (0.015) | 3 (0.011) ....>B:.503_:[..C]C-----G[..G]:.570_:B<.... (0.007) | 4 (0.010) ....>B:.504_:[..C]C-----G[..G]:.569_:B<.... (0.008) | 5 (0.010) ....>B:.505_:[..C]C-----G[..G]:.568_:B<.... (0.007) | 6 (0.011) ....>B:.506_:[..C]C-----G[..G]:.567_:B<.... (0.007) | 7 (0.012) ....>B:.507_:[..A]A-----U[..U]:.566_:B<.... (0.014) | 8 (0.007) ....>B:.549_:[..G]G-----C[..C]:.565_:B<.... (0.010) | 9 (0.006) ....>B:.550_:[..G]G-----C[..C]:.564_:B<.... (0.011) | 10 (0.007) ....>B:.551_:[..G]G-----C[..C]:.563_:B<.... (0.010) | 11 (0.007) ....>B:.552_:[..G]G-----C[..C]:.562_:B<.... (0.010) | 12 (0.009) ....>B:.553_:[..G]G-----C[..C]:.561_:B<.... (0.010) | 13 (0.016) ....>B:.554_:[..U]U-**--A[..A]:.558_:B<.... (0.008) | 14 (0.017) ....>B:.555_:[..U]U-**+-G[..G]:.517_:B<.... (0.015) x 15 (0.008) ....>B:.538_:[..A]A-**--C[..C]:.532_:B<.... (0.012) | 16 (0.009) ....>B:.539_:[..G]G-----C[..C]:.531_:B<.... (0.010) | 17 (0.012) ....>B:.540_:[..G]G-----C[..C]:.530_:B<.... (0.010) | 18 (0.012) ....>B:.541_:[..C]C-----G[..G]:.529_:B<.... (0.011) | 19 (0.011) ....>B:.542_:[..C]C-----G[..G]:.528_:B<.... (0.006) | 20 (0.008) ....>B:.543_:[..G]G-----C[..C]:.527_:B<.... (0.010) | 21 (0.007) ....>B:.544_:[..A]A-**--G[..G]:.526_:B<.... (0.006) | 22 (0.010) ....>B:.510_:[..G]G-----C[..C]:.525_:B<.... (0.010) | 23 (0.015) ....>B:.511_:[..U]U-----A[..A]:.524_:B<.... (0.013) | 24 (0.012) ....>B:.512_:[..C]C-----G[..G]:.523_:B<.... (0.007) | 25 (0.017) ....>B:.513_:[..U]U-*---G[..G]:.522_:B<.... (0.010) | 26 (0.014) ....>B:.508_:[..U]U-**--A[..A]:.546_:B<.... (0.008) | 27 (0.009) ....>B:.514_:[..A]A-**--A[..A]:.521_:B<.... (0.009) | 28 (0.008) ....>B:.515_:[..G]G-**+-C[..C]:.548_:B<.... (0.009) x 29 (0.008) ....>B:.518_:[..G]G-----C[..C]:.556_:B<.... (0.015) + 30 (0.009) ....>D:.501_:[..G]G-----C[..C]:.572_:D<.... (0.010) | 31 (0.010) ....>D:.502_:[..G]G-**--U[..U]:.571_:D<.... (0.016) | 32 (0.010) ....>D:.503_:[..C]C-----G[..G]:.570_:D<.... (0.008) | 33 (0.009) ....>D:.504_:[..C]C-----G[..G]:.569_:D<.... (0.006) | 34 (0.009) ....>D:.505_:[..C]C-----G[..G]:.568_:D<.... (0.009) | 35 (0.008) ....>D:.506_:[..C]C-----G[..G]:.567_:D<.... (0.008) | 36 (0.009) ....>D:.507_:[..A]A-----U[..U]:.566_:D<.... (0.015) | 37 (0.011) ....>D:.549_:[..G]G-----C[..C]:.565_:D<.... (0.012) | 38 (0.008) ....>D:.550_:[..G]G-----C[..C]:.564_:D<.... (0.010) | 39 (0.009) ....>D:.551_:[..G]G-----C[..C]:.563_:D<.... (0.014) | 40 (0.007) ....>D:.552_:[..G]G-----C[..C]:.562_:D<.... (0.011) | 41 (0.011) ....>D:.553_:[..G]G-----C[..C]:.561_:D<.... (0.009) | 42 (0.017) ....>D:.554_:[..U]U-**--A[..A]:.558_:D<.... (0.008) | 43 (0.019) ....>D:.555_:[..U]U-**+-G[..G]:.517_:D<.... (0.022) x 44 (0.011) ....>D:.538_:[..A]A-**--C[..C]:.532_:D<.... (0.013) | 45 (0.013) ....>D:.539_:[..G]G-----C[..C]:.531_:D<.... (0.008) | 46 (0.008) ....>D:.540_:[..G]G-----C[..C]:.530_:D<.... (0.009) | 47 (0.010) ....>D:.541_:[..C]C-----G[..G]:.529_:D<.... (0.008) | 48 (0.010) ....>D:.542_:[..C]C-----G[..G]:.528_:D<.... (0.006) | 49 (0.009) ....>D:.543_:[..G]G-----C[..C]:.527_:D<.... (0.012) | 50 (0.009) ....>D:.544_:[..A]A-**--G[..G]:.526_:D<.... (0.009) | 51 (0.010) ....>D:.510_:[..G]G-----C[..C]:.525_:D<.... (0.012) | 52 (0.018) ....>D:.511_:[..U]U-----A[..A]:.524_:D<.... (0.009) | 53 (0.011) ....>D:.512_:[..C]C-----G[..G]:.523_:D<.... (0.008) | 54 (0.016) ....>D:.513_:[..U]U-*---G[..G]:.522_:D<.... (0.007) | 55 (0.015) ....>D:.508_:[..U]U-**--A[..A]:.546_:D<.... (0.008) | 56 (0.010) ....>D:.514_:[..A]A-**--A[..A]:.521_:D<.... (0.010) | 57 (0.010) ....>D:.515_:[..G]G-**+-C[..C]:.548_:D<.... (0.008) x 58 (0.012) ....>D:.518_:[..G]G-----C[..C]:.556_:D<.... (0.032) + Note: This structure contains 18[15] non-Watson-Crick base-pairs.
По данным в таблице выше можно определить номера нуклетидов, образующих стебли, и сопоставить их с существующими названиями (с чем помогла эта статья):
Также можно найти неканонические пары оснований:
В файле 1g59_old.out в разделе Overlap area были найдены пары азотисных оснований с наибольшим перекрыванием. Они находились на втором шаге.
Далее с помощью команды ‹ex_str -2 stacking.pdb step2.pdb› была вырезана соответствующая 2 шагу структура.
На следующем шаге было построено изображение этой структуры с помощью команды ‹stack2img -cdolt step2.pdb step2.ps›. Результат после конвертации представлен на изображении ниже: