NP-completeness of the direct energy barrier problem without - - PowerPoint PPT Presentation
NP-completeness of the direct energy barrier problem without - - PowerPoint PPT Presentation
NP-completeness of the direct energy barrier problem without pseudoknots Jn Mauch, Anne Condon, Ladislav Stacho, Chris Thachuk Motivation: structure, pathway, barrier Motivation: structure , pathway, barrier Motivation: structure, pathway ,
Motivation: structure, pathway, barrier
Motivation: structure, pathway, barrier
Motivation: structure, pathway, barrier
Yin et al., 2008
(multiple copies
- f “fuel”)
Motivation: structure, pathway, barrier
Yin et al., 2008
Motivation: structure, pathway, barrier
Yin et al., 2008
Motivation: structure, pathway, barrier
Yin et al., 2008
Motivation: structure, pathway, barrier
Yin et al., 2008
Motivation: structure, pathway, barrier
Yin et al., 2008
Motivation: structure, pathway, barrier
unfolded
true
Motivation: structure, pathway, barrier
energy
unfolded
true
MFE barrier from true to MFE
Motivation: literature on energy barriers
Motivation: literature on energy barriers
Chen, S.J., Dill, K.A.: RNA folding energy landscapes. Proc. Nat.
- Acad. Sci. 97(2) (January 2000) 646–651
Russell, R., Zhuang, X., Babcock, H., Millett, I., Doniach, S., Chu, S., Herschlag, D.: Exploring the folding landscape of a structured
- RNA. Proc. Nat. Acad. Sci. 99 (2002) 155–160
Shcherbakova, I., Mitra, S., Laederach, A., Brenowitz, M.: Energy barriers, pathways, and dynamics during folding of large, multidomain RNAs. Curr. Opin. Chem. Biol. (2008) 655–666 Treiber, D.K., Williamson, J.R.: Beyond kinetic traps in RNA
- folding. Curr. Opin. Struc. Biol. 11 (2001) 309–314
Motivation: literature on energy barriers
Flamm, C., Fontana, W., Hofacker, I.L., Schuster, P.: RNA folding at elementary step resolution. RNA (2000) 325–338 Tang, X., Thomas, S., Tapia, L., Giedroc, D.P., Amato, N.M.: Simulating RNA folding kinetics on approximated energy
- landscapes. J. Mol. Biol. 381 (2008) 1055–1067
Wolfinger, M.T.: The energy landscape of RNA folding. Master’s thesis, University Vienna (2001) Flamm, C., Hofacker, I.L., Stadler, P.F., Wolfinger, M.T.: Barrier trees of degenerate landscapes. Zeitschrift f¨ur Physikalische
Our contribution
Q: is there a polynomial-time algorithm to calculate the min-barrier folding pathway between two structures? A: probably not
Overview
- clear statement of our result
- some details about our proof
- open questions
Notation: secondary structures
Notation: secondary structures
I F
Notation: secondary structures
I F
8 3 8 10 7
Notation: direct folding pathway
8 3 8 10 7
I F
+
- +
Notation: energy barrier
+
- +
- 3
- 8
- 18
- 13
b a r r i e r 1 5 8 3 8 10 7
I F
Notation: energy barrier
+ +
b a r r i e r 1 5 8 3 8 10 7
I F
- 3
- 8
- 18
- 13
The EB-DPKF problem
(Energy Barrier for Direct PseudoKnot Free pathways)
two (pkfree) structures I and F, and integer k is there a direct folding pathway from I to F whose energy barrier is at most k?
Given: Q:
Example
what is the min-barrier pathway?
4 8 7 5 10 8 3 10
Example
what is the min-barrier pathway?
4 8 7 5 10 8 3 10
+
- +
- +
Example
what is the min-barrier pathway?
4 8 7 6 10 8 3 10
+
- +
- +
15 10 5 b a r r i e r 1 2
Main result
The EB-DPKF problem is NP-complete
Main result
The EB-DPKF problem is NP-complete Proof: we show a polynomial time reduction from a known NP-complete problem, namely 3-Partition, to EB-DPKF
3-Partition problem: example
integers 10, 9, 8, 7, 7, 7 can we partition the integers into triples, with the integers in each triple summing to 24? Given: Q:
3-Partition problem: example
integers 10, 9, 8, 7, 7, 7 can we partition the integers into triples, with the integers in each triple summing to 24? yes! {9,8,7} and {10,7,7} Given: Q: A:
3-Partition problem
integers a1, . . . ,a3m and A (in unary), where A/4 < ai < A/2 and the sum of the ai’s is mA can we partition the integers into m triples, with the integers in each triple summing to A? Given: Q:
Reduction: 3-partition → EB-DPKF
Reduction: 3-partition → EB-DPKF
➓➒➑➐➐➐ ➉➈➇➆➆➆
triple1 triple2
→
10, 9, 8, 7, 7, 7
k
Reduction: 3-partition → EB-DPKF
➓➒➑➐➐➐ ➉➈➇➆➆➆
triple2
→
10, 9, 8, 7, 7, 7
k
Correctness: Suppose that the reduction maps instance x of 3-Partition to instance y of EB-DPKF. Then x is a “yes”-instance → y is a “yes”-instance, and x is a “no”-instance → y is a “no”-instance. triple1
Reduction: 3-partition → EB-DPKF
➓ ➈ ➇ ➆ ➐ ➐ triple1
triple2 ➒➑➐➉ ➆➆
➓➒➑➐➐➐ ➉➈➇➆➆➆
triple2
→
10, 9, 8, 7, 7, 7
Correctness (yes → yes): consider folding pathway:
k
triple1
Reduction: 3-partition → EB-DPKF
➓ ➈ ➇ ➆ ➐ ➐ triple1
triple2 ➒➑➐➉ ➆➆
{
triple-choosing
{
triple-validating
{
clean-up
➓➒➑➐➐➐ ➉➈➇➆➆➆
triple2
→
10, 9, 8, 7, 7, 7
Correctness (yes → yes): consider folding pathway:
k
triple1
Reduction: 3-partition → EB-DPKF
➓ ➈ ➇ ➆ ➐ ➐ triple1
triple2 ➒➑➐➉ ➆➆
➓➒➑➐➐➐ ➉➈➇➆➆➆
triple2
→
10, 9, 8, 7, 7, 7
Correctness (yes → yes): consider folding pathway:
k k
triple1
Reduction: 3-partition → EB-DPKF
➈ ➆ ➐ ➐ triple1
triple2 ➒
➐ ➆➆ ➒ ➐➐➐ ➈ ➆➆➆
triple2
→
9, 9, 9, 7, 7, 7
Correctness (no → no): all folding pathways exceed barrier:
➒ ➒ ➈ ➈ ➒ ➈ ➒ ➈
k k
triple1
Observations
- the EB-DPKF instance obtained from a
“yes”-instance of 3-Partition has an energy landscape such that
- there are exponentially many “initial”
paths from I that stay within barrier
- few of these continue to F within
barrier
Open problems
- is the problem of determining the min-