Предскажите минимальную свободную энергию вторичная структура последовательности RNA
rnafold(
Seq
)
RNAbracket
= rnafold(Seq
)
[RNAbracket, Energy
] =
rnafold(Seq
)
[RNAbracket, Energy, RNAmatrix
]
= rnafold(Seq
)
...
= rnafold(Seq
,
...'MinLoopSize', MinLoopSizeValue
, ...)
...
= rnafold(Seq
,
...'NoGU', NoGUValue
, ...)
...
= rnafold(Seq
,
...'Progress', ProgressValue
, ...)
Seq | Любое из следующего:
|
MinLoopSizeValue | Целое число, задающее минимальный размер циклов (в базисах), чтобы быть рассмотренным при вычислении свободной энергии. Значением по умолчанию является 3 . |
NoGUValue | Средствам управления или GU или парам UG запрещают сформироваться. Выбором является true или false (значение по умолчанию). |
ProgressValue | Управляет отображением индикатора выполнения во время расчета минимальной свободной энергии вторичная структура. Выбором является true или false (значение по умолчанию). |
RNAbracket | Вектор символов точек и скобок, указывающих на обозначение скобки для энергии без минимума вторичная структура последовательности RNA. В обозначении скобки каждая точка представляет непарную основу, в то время как пара одинаково вложенных, открывающих и закрывающих скобок представляет пару оснований. |
Energy | Значение, задающее энергию (в kcal/mol) минимальной свободной энергии вторичная структура последовательности RNA. |
RNAmatrix | Матрица смежности, представляющая минимальную свободную энергию вторичная структура последовательности RNA. Двоичный файл, верхняя треугольная матрица, где если и только если i остаток th в последовательности RNA Seq соединяется с j остаток th Seq . |
rnafold(
предсказывает и отображает вторичную структуру (в обозначении скобки) сопоставленный с минимальной свободной энергией для последовательности RNA, Seq
)Seq
, использование термодинамического подхода ближайшего соседа.
Примечание
Для длинных последовательностей это предсказание может быть трудоемким. Например, последовательность с 600 нуклеотидами может занять несколько минут и последовательностей, больше, чем 1 000 нуклеотидов могут принять 1 час, в зависимости от вашей системы.
предсказывает и возвращает вторичную структуру, сопоставленную с минимальной свободной энергией для последовательности RNA, RNAbracket
= rnafold(Seq
)Seq
, использование термодинамического подхода ближайшего соседа. Возвращенная структура, RNAbracket
, находится в обозначении скобки, которое является вектором из точек и скобок, где каждая точка представляет непарную основу, в то время как пара одинаково вложенных, открывающих и закрывающих скобок представляет пару оснований.
[
также возвращает RNAbracket, Energy
] =
rnafold(Seq
)Energy
, энергетическая ценность (в kcal/mol) минимальной свободной энергии вторичная структура последовательности RNA.
[
также возвращает RNAbracket, Energy, RNAmatrix
]
= rnafold(Seq
)RNAmatrix
, матрица смежности, представляющая вторичную структуру, сопоставлена с минимальной свободной энергией. RNAmatrix
верхняя треугольная матрица где
если и только если RNAmatrix
(i, j) = 1i
остаток th в последовательности RNA Seq
соединяется с j
остаток th Seq
.
... = rnafold (
вызовы Seq
PropertyName
', PropertyValue
, ...)rnafold
с дополнительными свойствами, которые используют имя свойства / пары значения свойства. Можно задать одно или несколько свойств в любом порядке. Каждый PropertyName
должен быть заключен в одинарные кавычки и нечувствительный к регистру. Это имя свойства / пары значения свойства следующие:
задает минимальный размер циклов (в базисах), чтобы быть рассмотренным при вычислении свободной энергии. Значением по умолчанию является ...
= rnafold(Seq
,
...'MinLoopSize', MinLoopSizeValue
, ...)3
.
средствам управления или GU или парам UG запрещают сформироваться. Выбором является ...
= rnafold(Seq
,
...'NoGU', NoGUValue
, ...)true
или false
(значение по умолчанию).
управляет отображением индикатора выполнения во время расчета минимальной свободной энергии вторичная структура. Выбором является ...
= rnafold(Seq
,
...'Progress', ProgressValue
, ...)true
или false
(значение по умолчанию).
Определите минимальную свободную энергию вторичная структура (и в скобке и в матричном обозначении) и энергетическая ценность следующей последовательности RNA:
seq = 'ACCCCCUCCUUCCUUGGAUCAAGGGGCUCAA';
[bracket, energy, matrix] = rnafold(seq);bracket
bracket =
..(((((...((....))...))))).....
[1] Wuchty, S., Фонтана, W., Hofacker, я., и Шустер, P. (1999). Завершите субоптимальное сворачивание RNA и устойчивость вторичных структур. Биополимеры 49, 145–165.
[2] Мэтьюс, D., Сабина, J., Zuker, M. и Токарь, D. (1999). Расширенная зависимость последовательности термодинамических параметров улучшает предсказание RNA вторичная структура. J. Молекулярная масса Biol. 288, 911–940.