Improved RNA
Download
1 / 22

Improved RNA Secondary Structure Prediction by Maximizing Expected Pair Accuracy - PowerPoint PPT Presentation


  • 97 Views
  • Uploaded on
  • Presentation posted in: General

Improved RNA Secondary Structure Prediction by Maximizing Expected Pair Accuracy. Zhi John Lu, Jason Gloor, and David H. Mathews University of Rochester Medical Center, Rochester, New York. RNA Secondary and Tertiary Structure:. AAUUGCGGGAAAGGGGUCAA CAGCCGUUCAGUACCAAGUC

loader
I am the owner, or an agent authorized to act on behalf of the owner, of the copyrighted work described.
capcha

Download Presentation

Improved RNA Secondary Structure Prediction by Maximizing Expected Pair Accuracy

An Image/Link below is provided (as is) to download presentation

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.


- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -

Presentation Transcript


Improved RNA

Secondary Structure Prediction by Maximizing

Expected Pair Accuracy

Zhi John Lu, Jason Gloor, and David H. MathewsUniversity of Rochester Medical Center, Rochester, New York


RNA Secondary and Tertiary Structure:

AAUUGCGGGAAAGGGGUCAA

CAGCCGUUCAGUACCAAGUC

UCAGGGGAAACUUUGAGAUG

GCCUUGCAAAGGGUAUGGUA

AUAAGCUGACGGACAUGGUC

CUAACCACGCAGCCAAGUCC

UAAGUCAACAGAUCUUCUGU

UGAUAUGGAUGCAGUUCA

Cate, et al. (Cech & Doudna).

(1996) Science 273:1678.

Waring & Davies. (1984) Gene 28: 277.


Gibbs Free Energy Change:

Ki =

=

= Ki/Kj =

The structure with the lowest DG° is the

most favored at a given temperature.


Nearest Neighbor Model for Free Energy Change of a Sample Hairpin Loop:

Mathews et al., J. Mol. Biol., 1999, 288: 911.

Mathews et al., PNAS, 2004, 101: 7287.


RNA Secondary Structure Prediction Accuracy:

Percentage of Known Base Pairs Correctly Predicted:

Mathews, Disney, Childs, Schroeder, Zuker, & Turner. 2004. PNAS 101: 7287.


Limitations to Prediction of the Minimum Free Energy Structure:

  • A minimum free energy structure provides the single best guess for the secondary structure.

  • Assumes that:

    • RNA is at equilibrium

    • RNA has a single conformation

    • RNA thermodynamic parameters are without error

      • Non-nearest neighbor effects

      • Some sequence-specific stabilities are averaged


A Method that Looks at the Probability of a Structure could be more Informative:

  • A partition function can be used to determine the probability of a structure at equilibrium.


The Partition Function, Q:


So, what is Q good for?

where k is the sum over all structures with the i-j base pair.


Accuracy:

  • Sensitivity – what percentage of known pairs occur in the predicted structure.

  • Positive Predictive Value (PPV) – what percentage of predicted pairs occur in the known structure.

  • PPV ≤ Sensitivity because the structures determined by comparative sequence analysis do not have all pairs and there is a tendency to over-predict base pairs by free energy minimization.


Applying Pi,j to Structure Prediction:

PPV

PBP≥ 90%

PPV

PBP≥ 70%

PPV

PBP> 50%

Positive

Predictive

Value (PPV)

PPV

PBP≥ 99%

PPV

PBP≥ 95%

Sensitivity

Mathews. RNA. 10: 1178. (2004).


Percent of Predicted BP above Threshold:

PPV

PBP≥ 99%

PPV

PBP≥ 95%

PPV

PBP≥ 90%

PPV

PBP≥ 70%

PPV

PBP> 50%

Mathews. RNA. 10: 1178. (2004).


Color Annotation:

E. coli 5S rRNA


Structures Constructed from Highly Probable Pairs:

PBP≥ 99%

PBP≥ 90%

PBP≥ 70%

PBP> 50%


“Maximizing Expected Accuracy:”


CONTRAfold:

  • “Statistical learning method” to predict Pi,j

  • Generate structures:

Where:

Bioinformatics. 22: e90-e98. (2006).


Implement Maximum Expectation:

  • Zhi John Lu, Jason Gloor, David Mathews

  • Implement dynamic programming algorithm

    using partition function prediction of Pi,j.

  • Also implement suboptimal structure prediction.

    • Alternative hypotheses.


Sensitivity and PPV vs. g:


Comparison:


Summary:

  • Maximizing expected accuracy can predict structures with greater sensitivity and positive predictive value than free energy minimization.

  • Maximizing expected accuracy using an underlying thermodynamic model is more accurate than an underlying statistical model.


Methanococcus thermolithotrophicus 5S rRNA (Szymanski et al., 1998):

MaxExpect Predicted Structure:


Minimum Free Energy Structure:

CONTRAfold

Predicted Structure:


ad
  • Login