-
Notifications
You must be signed in to change notification settings - Fork 16
/
hybrid-only.html
125 lines (110 loc) · 5.54 KB
/
hybrid-only.html
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
<div class="row" id="introduction">
<div class="colW600">
The identification of an optimal RNA-RNA interaction <em>without</em> consideration
of intramolecular base pairs can be done efficiently in $O(nm)$ time and space
for two RNA sequences $S^{1}$ and $S^{2}$ of lengths $n$ and $m$, resp.,
<a href="https://dx.doi.org/10.1093%2Fnar%2Fgkl356">(Brian Tjaden et al., 2006)</a>.
For a <a href="index.jsp?toolName=Nussinov">Nussinov</a>-like
interaction prediction, a dynamic programming table $H$ is filled via the
<em>prefix-based</em> recursion
<div style="padding-left:45px;">
$H_{i,j}=\max\begin{cases} H_{i-1,j-1}+1 & \text{if }S^{1}_{i}, \overleftarrow{S^{2}_{j}}\text{ compl. base pair}\\ H_{i-1,j} \\ H_{i,j-1} \end{cases}$,
</div>
where an entry $H_{i,j}$ provides the maximal number of intermolecular
base pairs for the prefixes $S^{1}_{1..i}$ and $\overleftarrow{S^{2}_{1..j}}$, the reverse sequence of $S^2$.
Note, this recursion is a variant of a global
<a href="https://dx.doi.org/10.1016/0022-2836(70)90057-4">Needleman-Wunsch (1970)</a>
sequence alignment with special scoring, where only
Watson-Crick as well as GU base pairs are considered complementary and thus similiar.
A traceback from the entry $H_{i,j}$ with the highest value provides
an according optimal interaction.
<br>
<br>
Here, in order to prepare the computation of
<a href="index.jsp?toolName=accessibility">accessibility-based interaction prediction</a>,
the following <em>interval-based</em> recursion is used to identify the RNA-RNA interaction
with the maximal number of intermolecular base pairs. To this end,
the four-dimensional table $D$ is filled, where an entry $D^{i,k}_{j,l}$
represents the maximal number of base pairs for the sequence
intervals $S^{1}_{i..k}$ and $\overleftarrow{S^{2}_{j..l}}$ under the condition
that both interval ends are forming the intermolecular base
pairs $(S^{1}_{i},\overleftarrow{S^{2}_{j}})$ and $(S^{1}_{k},\overleftarrow{S^{2}_{l}})$.
<div id="rec_id" data-bind="text: latex()[0]"></div>
A list of up to 15 maximal base pair RNA-RNA interactions is provided
via exhaustive traceback for all table entries with maximal value,
i.e. $D^{i,k}_{j,l}=\max(D)$.
</div>
<div class="colW150">
<img alt="hybrid-only RNA-RNA interaction" src="hybrid-only-120x90.png" width=120 height=90 >
</div>
</div>
<div id="inputContainer">
<div class="row">
<div class="colW200" style="font-size: 120%; vertical-align: bottom; ">RNA sequence $S^{1}$:</div>
<div class="colW400">
<input data-bind="value: rawSeq" class="userInput"
placeholder="Example 'GCACG'" onkeydown="validate(event)"
style="text-transform:uppercase">
</div>
</div>
<div class="row">
<div class="colW200" style="font-size: 120%; vertical-align: bottom; ">RNA sequence $S^{2}$:<br>
<span style="font-size:60%;">(Computation uses reversed sequence $\overleftarrow{S^2}$)</span>
</div>
<div class="colW400">
<input data-bind="value: rawSeq2" class="userInput"
placeholder="Example 'GCACG'" onkeydown="validate(event)"
style="text-transform:uppercase">
</div>
</div>
<div class="row" style="display:none;">
<div class="colW200" style="font-size: 120%">Minimal loop length $l$:</div>
<div class="colW400">
<select data-bind="value: loopLength" id="ll" style="width:40px;">
<option value=0>0</option>
<option value=1>1</option>
<option value=2>2</option>
<option value=3>3</option>
</select>
<label for="ll" style="margin-left:10px;">i.e. minimal number of enclosed positions</label>
</div>
</div>
<div id="rec_select" style="display: none;">hybrid</div>
</div>
<div id="output">
<h2>Visualization of interacting base pairs (selected structure)</h2>
<div>
Due to the four-dimensionality of $D$, we only list the optimal
hybrid structures (up to 15). On selection, the intermolecular base pairs are
visualized.
</div>
<div class="row" style="margin:0;">
<div class="colW400" style="margin:0;overflow-x:auto;">
<table id="structures">
<thead>
<tr><th id="structTableHeading">Possible Structures</th></tr>
</thead>
<tbody data-bind="foreach: { data: tracebacks, as: 'tb' }">
<tr><td id="structTableCells" data-bind="text: tb.structure, event: {mousedown:$parent.clickStructure4d}"></td></tr>
</tbody>
</table>
</div>
<div class="colW400" style="margin:0;">
<div id="4dVisual" class="hybridVisual"></div>
</div>
</div>
<div style="font-weight:bold;">
Basepair number of selection:
<span id="output_value"></span>
</div>
<div>
The box provides an ASCII representation of the interacting
base pairs of the selected structure with $S^{1}$ on top and $\overleftarrow{S^{2}}$
on the bottom.
Note, sequence $S^{2}$ is reversed (running from right ($5'$-) to left
($3'$-end)) within this representation.
Note further, if no interacting
base pairs are present, no visualization is done.
</div>
</div>
<script src="js/visualize.js"></script>