Образовательная программа «6D060722- молекулярная и клеточная биология»


Характеристика взаимодействия микроРНК с генами мишенями



бет8/11
Дата11.07.2016
өлшемі2.9 Mb.
#192326
түріОбразовательная программа
1   2   3   4   5   6   7   8   9   10   11
3.4 Характеристика взаимодействия микроРНК с генами мишенями

С помощью программы UNAFold.3.7 (http: // dinamelt.bioinfo.rpi.edu) были построены 2D структуры мРНК 13 генов мишеней микроРНК (ABCC2, ABCG2, ALCAM, APC, AXIN1, BAX, CDH1, FLCN, MLH3, MMP2, PTPN12, SRC, TP53). На вторичной структуре мишеней была изучена характеристика взаимодействия микроРНК с мРНК-мишенью. Как показано на рисунке 2, микроРНК взаимодействует с неспаренными нуклеотидами мРНК с 5' или 3'конца микроРНК. Это позволяет разрушить уотсон-криковское взаимодействие в 2D структуре мРНК для образования взаимодействия между микроРНК и мРНК. Неспаренные нуклеотиды мРНК служат затравкой для образования связи между микроРНК и мРНК. Из приведенных примеров видно, что затравкой может быть как 5' так и 3'конец микроРНК. В сайте для miR-21* в мРНК гена ALCAM затравкой служит взаимодействие 5’конца микроРНК с неспаренными нуклеотидами вторичной структуры мРНК (рисунок 9). А в сайте взаимодействия miR-1279 с мРНК гена ALCAM затравкой образования связи служит взаимодействие между 3‘концом микроРНК со свободными нуклеотидами мРНК (рисунок 9). Так как программа RNAhybrid предсказывает сайты на основе вторичной структуры мишени и микроРНК было обнаружено, что во всех изученных случаях микроРНК взаимодействует с неспаренными нуклеотидами вторичной структуры мишени.



Копмлементарные участки представлены непрерывной чертой, ондонуклеотидные мисматчи точкой, несколько некомплементарных нуклеотида в сайте связывания точкой-тире.


Рисунок 9 - мРНК с двумя сайтами микроРНК
Как видно из примеров генов CDH1 и ALCAM одна мРНК может связываться с несколькими микроРНК с высоким уровнем комплементарности (рисунки 9 и 10). Если микроРНК связываются с одной областью мРНК, то они конкурируют за область взаимодействия. miR-1587 и miR-4507 взаимодействуют в одном сайте с мРНК CDH1, то есть сайты перекрываются и конкурируют за мишень мРНК (рисунок 10). Тогда как в других мРНК есть сайты связывания для нескольких микроРНК в разных областях мРНК. Например, сайты связывания для miR-21* и miR-1279 в мРНК гена ALCAM (рисунок 9).

Комплементарные участки представлены непрерывной чертой, ондонуклеотидные мисматчи точкой, несколько некомплементарных нуклеотида в сайте связывания точкой-тире.


Рисунок 10 - мРНК с перекрывающими сайтами микроРНК
Из рассмотренных 63 сайтов в мРНК 13 генов 27% 5'-доминантных сайтов, 43% сайтов с центральным доминированием, 30% 3'-доминантных сайтов. На рисунках 9, 10, 11 можно найти сайты всех трех видов. Взаимодействия miR-21*:ALCAM, miR-1587:CDH1, miR-4508:FLCN1 являются сайтами с доминирующим 5'-концом микроРНК. Сайты miR-212:MMP2, miR-4472:ABCG2 являются сайтами с доминированием центрального участка микроРНК (рисунок 11). К 3'-доминирующим сайтам относятся сайты следующих пар: miR-1279:AlCAM, miR-4507:CDH1. Несмотря на то, какая часть микроРНК вносит больший вклад в энергию взаимодействия из рассмотренных примеров видно, что в основном затравкой для образования связи между микроРНК и мРНК служит 5' конец микроРНК. По итогам исследования была опубликована одна публикация [346].

Копмлементарные участки представлены непрерывной чертой, ондонуклеотидные мисматчи точкой, несколько некомплементарных нуклеотида в сайте связывания точкой-тире.


Рисунок 11 - Характеристика взаимодействия микроРНК с генами мишенями
Изучение распределения сайтов микроРНК в зависимости от длины микроРНК показало, что основная доля сайтов приходится на микроРНК длиной 22 нуклеотида. Пять микроРНК с длинной 16 нуклеотид не имеют сайтов связывания при установленных критериях взаимодействия. Из рисунка 12 видно, что микроРНК с короткой длиной имеют большое количество сайтов. Это связано с длиной микроРНК, потому что для отбора сайтов для микроРНК длиной 17 нуклеотид энергия гибридизации должна составлять 82,5% от максимальной энергии гибридизации, а для микроРНК длиной 27 нуклеотид 70%. Для длинных микроРНК критерий отбора были выше, так как основной задачей исследования было определение сайтов с высокой комплементарностью. Поэтому микроРНК с длиной выше 22 нуклеотид при установленных критериях имеют ограниченное количество сайтов. Только 10% сайтов составляют микроРНК с длиной от 23 до 27 нуклеотидов. Максимальное количество сайтов приходится на микроРНК длиной 22 нуклеотида, что соответствует отношению количества межгенных микроРНК.

Рисунок 12 - Распределение сайтов связывания в зависимости от длины микроРНК


3.5 Изучения различий между сайтами, предсказанных с помощью программ RNAhybrid и TargetScan

В следующем этапе исследования сайты, которые были найдены при помощи программы RNAhybrid, были сравнены с сайтами, которые предсказывает программа TargetScan 6.2 (www.targetscan.org). Отличительной особенностью сайтов отобранных нами по программе RNAhybrid является их высокая комплементарность во всей области взаимодействия, ранее ни в одной исследовательской работе не было показано наличие таких сайтов. Поэтому мы хотим проверить эффективность предсказывания сайтов с высокой комплементарностью с помощью программы TargetScan. Программа TargetScan разработана группой ученых под руководством Давида Бартела из Howard Hughes медицинского университета Америки и признана одной из широко используемых программ с низким уровнем предсказания ложно позитивных сайтов взаимодействия микроРНК [15]. Сравнения проведено для сайтов связывания микроРНК с мРНК генов CDH1, GSK3B, MLH3, MTHFR, PTPN12, SMAD4, SRC, TP53, ZEB1. В таблице 14 представлены сайты связывания микроРНК предсказанные программой RNAhybrid и TargetScan. Все сайты, которые имеют хоть одну Г:У пару в сайте связывания не предсказываются программой TargetScan. В таблице жирным шрифтов помечены сайты без Г:У пар в 'seed' районе, которые определены в одной области мРНК мишеней, но за счет специфических свойств алгоритма программы TargetScan не определена биологическая значимость этих сайтов. То есть комплементарность сайтов взаимодействия, которая в описанных сайтах близка к полной. Для сайтов между miR-1285:TP53, miR-4663:ZEB1, miR-4732-3p:GSK3B, miR-4710:GSK3B, miR-4711-3p:PTPN12 программа TargetScan не определила комплементарные пары в 3'конце микроРНК. Определение сайта микроРНК в мРНК гена TP53, во всех остальных сайтах, представленных в таблице 14, программа TargetScan не нашла сайты, которые ранее были предсказаны программой RNAhybrid и предлагает свои варианты сайтов в другой области мРНК. В базе TargetScan все сайты в зависимости от филогенетической консервативности подразделены на три группы: сайты консервативные у позвоночных, сайты консервативные у млекопитающих, сайты с низкой консервативностью. Все сайты, представленные в таблице 14, получены из он-лайн базы TargetScan, определены как сайты с низкой консервативностью. То есть экспериментально биологическая значимость этих сайтов в регуляции генов мишеней будет определяться в последнюю очередь, так как первый критерий для отбора сайтов для экспериментальных исследований основан на филогенетической консервативности сайтов.

Некоторые канонические сайты с высокой комплементарностью не предсказываются с помощью программы TargetScan. В таблице 15 представлены несколько сайтов, которые не были найдены при помощи программы TargetScan. Анализ структуры канонических сайтов выявил, что сайты с некомплементарными нуклеотидами мишени ко второй и третьей позиции 5'конца микроРНК (SMAD4:miR-1972; TP53:miR-2392) не определяются данной программой. Программа находит только канонические сайты: 5’-доминантный канонический сайт, сайт с 5’-доминантной seed областью и 3’-компенсаторным сайт (рисунок 6). Можно сделать вывод, что локализация seed является одним из основных критериев отбора сайтов в алгоритме программы TargetScan. Во вторых, сайты с очень высокой уровенью комплементарности, но с Г:У парами в seed области не определяются данной программой (SMAD4:miR-1268, SRC:miR-302f). Исходя из этого можно сделать вывод, что наличие Г:У пар в seed области также является критерием отбора.

Все эти данные свидетельствуют о преимуществе разработанной в нашей лаборатории методики отбора сайтов связывания. Методика позволяет отобрать сайты взаимодействия микроРНК с высоким уровнем комплементарности, то есть близких по свойствам к siRNA.


Таблица 14 - Сайты предсказанные RNAhybrid и TargetScan




RNAhybrid

TargetScan (Poorly conserved miRNA Families)

CCND1:hsa-miR-4487
target 5' G G A 3'

GCCCU CAGCCAGCUC

CGGGA GUCGGUCGAG

miRNA 3' GA A A 5'

5' ...CCAACGGCCCUGCAGCCAGCUCA...
            |||||    ||||||| 
3'        GACGGGAAGUCGGUCGAGA


3-utr,2022 ΔG = -37.6 kcal/mol

8mer, 3-utr,937

CCND1:hsa-miR-3180-5p

target5'U A A G3'

GGC GUGG GGUGGGGUGUUUGGGAG

CUG CACC CCGCCUCGCAGACCUUC

miRNA 3'G 5'


5'      ...AUGAUUGGAAUAGCUUCUGGAAU...
                          ||||||| 
3'       GCUGCACCCCGCCUCGCAGACCUUC 
5'     ...GUGAGAAAAAAACAAUCUGGAAG...
                         ||||||| 
3'      GCUGCACCCCGCCUCGCAGACCUUC

3-utr,2250 ΔG = -45.1 kcal/mol

3-utr,2578; 3-utr,1602

CDH1:hsa-miR-1285

target 5' G A A 3'



GGGUCUUGCU UGUUGCCCA

UCCAGAGUGA ACAACGGGU

miRNA 3' A CU 5'



5' ...AUGAGCCACUGCACCUGCCCAGC...
            ||||     ||||||| 
3'    UCCAGAGUGAAACA-ACGGGUCU
5' ...GGGCAUGAGCUGCUGUGCCCAGC...
                     ||||||| 
3'     UCCAGAGUGAAACAACGGGUCU

3-utr,3677 ΔG = -36.7 kcal/mol

750, 1009

MLH3:hsa-miR-378b

target 5' C C G 3'

UCUGCCUCCAGG CCAGU

AGACGGAGGUUC GGUCA

miRNA 3' A A 5'


5'  ...AUGGCAAAAAGGUGGGUCCAGAA...
                      ||||||  
3'        AAGACGGAGGUUCAGGUCA

3-utr,6094 ΔG = -37.6 kcal/mol

3-utr,3228

MTHFR:hsa-miR-1260

target 5' A U G 3'

GGUGGC GAGGUGGGA

CCACCG CUCCACCCU

miRNA 3' A U A 5'

5' ...UACUCAGGUGGCUGAGGUGGGAG...
            ||||||   ||||||| 
3'         ACCACCGUCUCCACCCUA




3-utr,6314 ΔG = -38.7 kcal/mol

7mer-m8, 3-utr,4130

MTHFR:hsa-miR-4665-5p

target 5' U C 3'

CUUGUGCUCACGC UCCCCCG

GAGCGCGAGUGCG AGGGGGU

miRNA 3' C C C 5


5'  ...AGAUUCCUGGGCCUGUCCCCCAA...
                      ||||||| 
3'     CGAGCGCGAGUGCGCAGGGGGUC

5' ...CACCCCGGCCUCCAC--UCCCCCAC...


               |||     ||||||| 
3'      CGAGCGCGAGUGCGCAGGGGGUC

3-utr,4252 ΔG = -43.1 kcal/mol

3-utr, 2188; 3-utr,116

TP53:hsa-miR-1285

target 5' U C 3'

GGGUCUCGCUUUGUUGCCCAGG

UCCAGAGUGAAACAACGGGUCU

miRNA 3' 5'

5' ...UGGGUCUCGCUUUGUUGCCCAGG...
               ||||||||||||| 
3'     UCCAGAGUGAAACAACGGGUCU



3-utr,2298 ΔG = -46.0 kcal/mol

7mer-m8,3-utr,1004

Продолжение таблицы 14




RNAhybrid

TargetScan (Poorly conserved miRNA Families)

ZEB1:hsa-miR-4663

target 5' U A 3'

UUGCAUGUUUAUGGAGCUCAGCU

GACGUGCAGGUACCUCGAGUCGA

miRNA 3' U 5'

5'  ...UUGCAUGUUUAUGGAGCUCAGCU...
                 |||||||||||| 
3'    UGACGUGCAGGUACCUCGAGUCGA


3-utr,6103 ΔG = -46.1 kcal/mol

7mer-m8, 3-utr, 2407

GSK3B:hsa-miR-4710

target 5' A U 3'

ACUUGCCCUCACCC

UGGACGGGAGUGGG

miRNA 3' UUGG 5'

5' ...UAAUAAUAAAAAGCUCCUCACCA...
                     ||||||| 
3'         UUGGUGGACGGGAGUGGG 
5' ...UUGUGUGUAACUUGCCCUCACCC...
                     ||||||| 
3'         UUGGUGGACGGGAGUGGG


3-utr,4046 ΔG = -34.6 kcal/mol

3-utr, 2573; 7mer-m8, 3-utr,1772

GSK3B:hsa-miR-4732-3p

target 5' A G A 3'

GGAA GGGACAGGUCAGGG

UCUU UCCUGUCCAGUCCC

miRNA 3' G G G 5'

5' ...AGAGGAAGGGGACAGGUCAGGGA...
               ||||||||||||| 
3'      GUCUUGUCCUGUCCAGUCCCG


3-utr,3170 ΔG =-41.0 kcal/mol

8mer,3-utr,901

PTPN12:hsa-miR-4711-3p

target 5' G U 3'

GAGCUAGAAGACAC

UUCGGUCUUCUGUG

miRNA 3' UAG C 5'

5' ...NNNUUCAGGGAGCUAGAAGACACU...
                      ||||||| 
3'          UAGUUCGGUCUUCUGUGC


3-utr,2438 ΔG =-28.9 kcal/mol

7mer-m8, 3-utr,14

Жирным шрифтов обозначены сайты обнаруженные в одной области мРНК.

Таблица 15 - Потенциальные сайта не найденные программой TargetScan




Сайты с высокой комплементарностью

Сайты без Г:У пар в 'seed' области

SMAD4 hsa-miR-1268 3-utr,4413 ΔG = -39.3 kcal/mol

target 5' G G 3'

CCCGCCACCAUGCCUG

GGGUGGUGGUGCGGGC

miRNA 3' GG 5'


SMAD4 hsa-miR-513b 3-utr,6663 ΔG = -25.7 kcal/mol

target 5' A C G 3'

GACACUUUC UGUGAA

CUGUGGAGG ACACUU

miRNA 3' UAUUUA A 5'


SRC hsa-miR-302f 3-utr,2950 ΔG = -24.4 kcal/mol

target 5' C A 3'

AGCAUGGAGGCAG

UUGUACCUUCGUU

miRNA 3' U AAU 5'


SMAD4 hsa-miR-1972 3-utr,4547 ΔG = -43.1 kcal/mol

target 5' G C C 3'

UGAGCCACUGU CCUGGCCU

ACUCGGUGACA GGACCGGA

miRNA 3' C CU 5'


VDR hsa-miR-4507 3-utr, 3066 ΔG = -41.6 kcal/mol

target 5' G G 3'

CCAGCCCAGCCCAGCU

GGUCGGGUCGGGUUGG

miRNA 3' G GUC 5'


TP53 hsa-miR-2392 3-utr,1540 ΔG = -36.5 kcal/mol

target 5' A C C 3'

GCCUC CACCCCCAUCU

UGGAG GUGGGGGUAGG

miRNA 3' G A AU 5'


GSK3B hsa-miR-4265 3-utr, 4020 ΔG =-38.2 kcal/mol 83.22

target 5' G U 3'

CAGAGCUGAGCCCAUGG

GUCUCGACUCGGGUGUC

miRNA 3' GG 5'


GSK3B hsa-let-7astar 3-utr,6823 ΔG = -26.6 kcal/mol 75.56

target 5' A U C 3'

AAGGAC GUGGGUUGUAUA

UUUCUG CAUCUAACAUAU

miRNA 3' C U C 5'


Используя выравнивание 3'UTR мРНК разных организмов, представленной в базе TargetScan, изучена филогенетическая консервативность сайтов связывания микроРНК локализованных в 3'UTR. Для семи мРНК были обнаружены сайты, предсказанные обеими программами в одной области мРНК мишени, но программа TargetScan во всех случаях не смогла показать все комплементарные пары в сайте взаимодействия. Для всех этих сайтов была изучена консервативность сайтов микроРНК среди различных видов. Если программой TargetScan определяется консервативность «seed» области (от 2 по 8 нуклеотид 5'конца микроРНК), мы изучили консервативность всей последовательности сайта, чтобы проверить уникальность таких сайтов для человека.

В таблице 16 представлены данные по консервативности сайтов в 3'UTR. Для трех изученных сайтов выявлена консервативность в одном геноме, во всех случаях для Pаn troglodytes (шимпанзе): GSK3B - miR-4732-3p, MTHFR - miR-1260, TP53 - miR-1285. Другие три сайта консервативны в двух геномах (у шимпанзе и кролика или у шимпанзе и галаго: CCND1 – miR-4487, PTPN12 – miR-4711-3p, ZEB1 – miR-4663. Сайт для miR-4710 в гене GSK3B был консервативен у шимпанзе, макаки-резуса и кролика. То есть, все эти сайты имеют низкую консервативность и не уникальны для генома человека.

Таблица 16 - Консервативность сайтов в 3'UTR




Локализация сайта

Нуклеотидная последовательность сайта взязывания

Сайт в 3'UTR CCND1 для miR-4487

......930........940.......950.......960..

hsa AACGGCCCUGCAGCCAGCUCACGUCCAGGUUCAACCCACAG

ptr AACGGCCCUGCAGCCAGCUCACGUCCAGGUUCAACCCACAG

oga AGUGGCCCUGCAGCCAGCUCACGUCCCGGUUCAACCCACAG



Сайт в 3'UTR ZEB1 для miR-4663

.....2380......2390......2400......2410..

hsa UUGAAUAAAAAUAUAAUUUGCAUGUUUAUGGAGCUCAGCUA

ptr UUGAAUAAAAAUAUAAUUUGCAUGUUUAUGGAGCUCAGCUA

ocu UAGAAUGACAUUACAAUUUGCAUGUUUAUGGAGCUCAGCUG



Сайт в 3'UTR PTPN12
для miR-4711-3p

1.......10.........20

hsa UUCAGGGAGCUAGAAGACACU

ptr UUCAGGGAGCUAAAAGACACU

ocu UUCAGGGAGCUAAAAGACACU




Сайт в 3'UTR GSK3B
для miR-4710

.....1770......1780......1790.

hsa AACUUGCCCUCACCCUCCUCAUUGCCACC

ptr AACUUGCCCUCACCCUCCUCAUUGCCACC

mml AACUUGCCCUCACCCUCCUCAUUGCCACC



ocu AACUUGCCCUCACCCUCCUCAUUGCCACC

Сайт в 3'UTR GSK3B для miR-4732-3p

.......890.......900.....

hsa CCCAGAGGAAGGGGACAGGUCAGGG

ptr CCCAGAGGAAGGGGACAGGUCAGGG


Сайт в 3'UTR TP53 для miR-1285

980.......990.......1000......1010.

hsa CUUUGAGACUGGGUCUCGCUUUGUUGCCCAGGCUG

ptr CUUUGAGACUGGGUCUCGCUUUGUUGCCCAGGCUG


Сайт в 3'UTR MTHFR
для miR-1260

4110......4120......4130......4140

hsa UCAGCUACUCAGGUGGCUGAGGUGGGAGGAUCGU

prt UCAGCUACUGAGGUGGCUGAGGUGGGAGGAUCGU



Достарыңызбен бөлісу:
1   2   3   4   5   6   7   8   9   10   11




©dereksiz.org 2024
әкімшілігінің қараңыз

    Басты бет