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
|
Достарыңызбен бөлісу: |