250 likes | 401 Views
Setting Up a Replica Exchange Approach to Motif Discovery in DNA . Jeffrey Goett Advisor: Professor Sengupta. RNA polymerase. Binding Proteins. Regulation. Transcription. Binding sites. gene. Translation to Proteins. Protein Synthesis from DNA. Binding Site. Binding protein “A”.
E N D
Setting Up a Replica Exchange Approach to Motif Discovery in DNA Jeffrey Goett Advisor: Professor Sengupta
RNA polymerase Binding Proteins Regulation Transcription Binding sites gene Translation to Proteins Protein Synthesis from DNA
Binding Site Binding protein “A” Binding protein “A” T - T - C - A - A - C - C - A - A - A - C - G - A - C - Sequence A: T - T - G - C - T - G - A - A - G - T - T - G - G - T - code for protein C - G - T - T - G - C - T - C - A - A - G - G - A - C - Sequence B: T - T - C - C - T - G - G - C - A - A - C - G - A - G - code for protein Binding Sites
Discovering New Binding Motifs Motif: GCTCAG …ATCG GCTCAG CTAG… …CACT GATCAG AGTA… …TTCC GCTCTG TAAC… …GCTA GCTCAA ATCG… Motif Probability Model
Modeling Motifs in Sequences Assume: Break into N sequences Each sequence has one instance of motif embedded in random background Variations of motif by point mutation, but not insertion or deletion ATATCCGTA AATCGAGAC TCGATGTGT CCACCTGCA
Modeling Motifs in Sequences The “Alignment:”Starting position of motif in each sequence AT ATC CGTA A ATC GAGAC TCG ATG TGT CC ACC TGCA The “Motif Probability Distribution:”Probability of each letter occurring at each motif position
0 0 0 0 0 0 0 0 0 0 pC pC pG pT pA pA pA pT pC pG p1,T p2,A p3,T 0 ATATCCGTA AATCGAGAC TCGATGTGT CCACCTGCA pA p1,A p2,G p3,A 0 pC 0 0 0 pT pC pG 0 0 0 pT pG pT p1,A p2,T p3,G 0 0 0 0 0 0 p1,C p2,C p3,A pC pC pT pG pC pA Scoring a Model “Log-likelihood” score:
{3, 2, 4, 3} {2, 4, 7, 3} Example Models A TAT CCGTA AAT CGA GAC TCGATG TGT CC ACC TGCA AT ATC CGTA A ATC GAGAC TCG ATG TGT CC ACC TGCA
The Gibbs Sampler that maximizes We want to find
The Gibbs Sampler Times visited Over time, the frequency distribution approaches
Optimization Technique If we assume areas of local maximization contribute the most during “integration” to the local maximizations of Biasing our search to these areas may discover the pj,ro values which maximize faster.
Multiple Gibbs Samplers By combining results from Gibbs Samplers begun at random positions, find maximizing sooner
Replica Exchange/Parallel Tempering “Low-sensitivity” samplers which “scout out area” periodically swap with “high-sensitivity” samplers good at focused searches if swap appears promising.
Controlling Sensitivity Adjust the relative probability of sampling an xi by adjusting a new parameter in distribution: Large Small Search breadth of space Focused search of region
Betas 2 1 .9 .1 Testing the Sensitivity Running on randomly generated sequences to see motifs found, different sensitivity samplers converge to different scores.
Predicting Convergence Score Measure of Similarity: magnetization Ex: m=.5 “Configuration Score:” energy m=0 m=0 m=1 m=0 m=.5 E=2J E=2J E=-6J E=2J E=0
Alignment Analogue m=1 E=-9J A: B: m=.77 E=-5J C: m=.77 m=.77 E=-5J E=-5J
Test Results L<|alphabet|w
Test Results L>|alphabet|w
Hidden Motifs: Gibbs Sampler Beta = .1 Beta = .5 Beta = .9 Beta = 1.3 Beta = 1.7 Beta = 2 W=5, l=500
Betas .8 .9 .93 .96 1 1.5 Hidden Motifs: Replica Exchange