html5-img
1 / 82

Proteins Determine Function

Proteins Determine Function. Proteins make us tick Problems occur when Proteins are missing Proteins are malfunctioning Proteins are present that should not be there Controlling disease Understanding protein function Stopping protein function Supplying missing/desired protein function.

jimbo
Télécharger la présentation

Proteins Determine Function

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. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Proteins Determine Function • Proteins make us tick • Problems occur when • Proteins are missing • Proteins are malfunctioning • Proteins are present that should not be there • Controlling disease • Understanding protein function • Stopping protein function • Supplying missing/desired protein function Sequence Alignment

  2. The Central Dogma Sequence Alignment

  3. The Human Genome • The entire collection of our DNA • Consists of about 3.5x109 base pairs • Our genome is split across chromosomes • Makes packaging easier/more efficient Sequence Alignment

  4. Human Genome Project • Began in 1990, a 13-year effort coordinated by the DOE and the NIH. Project goals included: • identify all the approximately 30,000 genes in human DNA, • determine the sequences of the 3 billion chemical base pairs that make up human DNA, • store this information in databases, • improve tools for data analysis, • transfer related technologies to the private sector, and. • address the ethical, legal, and social issues. • The Human Genome Project ended in 2003 with the completion of the human genetic sequence. Sequence Alignment

  5. DNA Sequencing • Determining the sequence of nucleotides in a strand of DNA • 1978 Sanger • Determined the DNA sequence for phi-x 174 • 3.5x103 base pairs • Took about 2 years to sequence • 2001 Venter/Celera • Sequenced the human genome • 3x109 base pairs • Took about 9 months to sequence Sequence Alignment

  6. Sequence Data TCTAGGAGGGAAGCACCCACCTCCCCTAAGCTCCATCTCCCTGAGCACTCATTTCCCAATGACCATACCAGGTTTTGGCCCTAGAGAGTTTATTACAAAATAAGAAAGAGAAGTCTGGGGAAGGTTCACTCATCATAGAATTTTGGCAGTTCATTGCCCAAGATGACTCGATGGTCCACACCGGCAGCTGTAATAGTGACCAGGTAGATGACACCCCCGCTTGAGCCATCCCGGCTCATGGCCAGAGCAATAGCTGCAGAGGGTTTCAAGTTGGAGAGAGGGAGAGAGAGGATGGCTTAGCTTCAAAAATCTTTTTACTCCCCCTCCATCCATATGCCTACTACCACTTTCACCTCAAAACTCATCTTCCAGGAAGGCATATTTAGTGGTGTGCTGGTAAATCAGTTTTTTTACAAAAAGGCTTCCATATGTGGCATCTGCTGATGTCCGTGGTGTAAATGCTCCCGCTATGATGAATTGCAAGTTACAAATAGCTAAGCAGTTCACAAATCCTTGACTATTTAACAGTCCGCTCTCATGAGTGGTCCCAAGCCAGCCTCAGCACACCTCAGCACACCACTGGTTCTTTTTTTTTTTTTTTTTCTCCAGACAGGGTCTCTCTCTGTCACCTAGGCTGCAGCGCAGTGGTGCAATCACCGCTCACTACAGCCTTGATCTCCCCGGCTCAGATGATCTTTCCACCTCAGCCTCCTGAGTAGCTGGGACTACAGGTGTGCACCACTATGCCCAGTTCATTTTTTTTTTTACTTTTTTTTATTGTTTTTTGTGGAGACAGGGTTTCACCATCTTGCCTAGGCTGGCCTCAAACTCCTGGGCTCAAGTAATCCTCCTGCCTCAGCCTCCCAAATTGTTGGCATTACAGGTGTGAGCCACTGTGCTTAGCACACCACTGGTTCTCACAGTGACTGTGTATCCTCATTTGATTTACTCAGAACAGCCCTGGTTTATCCGTATTGCCCAAGAACCCCATTGAGCTTTGCATTTGTCCTGCCCCTTTTCACTCTTAAAAGTGTACCAGGCCCGGCATTAACTTAAATGGCCACCCCTGTATTTCTCTTCCTGTTCCTCATAATCTACTTCCTTCCCATGTTTCAAAGCCCTCCCCAGGTACCCTTCCACTTGGCTGGTTACCGTCTGTGGTGAAGCGCCTGCACTCCTCGGGAGACATGCCTGGCTTATATGCTGCATCCACATAACCATAGATAAAGGTGCTGCCGGAGCCACCAATGGCAAAAGGCTGTCGAGTCAGCATTCCTCCCAGGGTTCCATATACCTGGGAAAGGGATCCTCAGGTTAAAGAATCATCAAGCCCTTCCTTCCCACTGAGACATTAAGTGGTCTCTGCACCCTGCAATGAAGCCCTGGTATCTCATATCCCCAAAGTACTATGCTTTCAGAGGTAGTGTCCTTGGAACTCATTGCTAGAATGACATAGGACTTCCATCTTCCTCTGCAGGAGAGTGGGGAAGCCCAGAGGAGAGAGTGCTTTGGGAGAAACTCACCTGACCTCCTTCACGTTGGTCCCAGCCAGCTACCATGAGATGTGCAGACAAGTCCTCTCGATATTTATAGCTGATATTTCTCACCACATTTGCAGCAGCCAAAACAAGTGGAGGTTCCTCCAGTTCTATCCTGAGGGAAATATTAGGAATAAAGGTTGATAGAATTTTAAGTCTCATTCTCCTATACTGTTACCATCATCCCTGCTAAACGACCCCTGAAAACTGTAACTGCAATAGCTCAAACTGCAGCCTCCCTCCCACATGTACAGGGGAACCAGAGTCCCACACCACCAACTGGTAAGAAGCTTTCAATTGCTCACTCTTTTGCTCAGCCCCACCCACATAACTTTCTTTTGGCTGCAAGGACCCTGCTCTTATGGGGAAAAGCAGATAAGGTTCACTCGGTTCACCACCGCCTCGCTGTCAGGAGGGAGTCAACAGTCACCAAGTTAAAACTCAGGTTTTTTTTTTTTTTTTTTTTTTTTGAGACAGTCTCACTCTGTCACCCAGGCTGGAGTGCAGTGGATCAATCTTGGGCTCACTGCAAACTTCGCCTCCCTGGTTCAAGTGATTCTCCTGCCTCAGCCTCCCGAATAGCTGGGATTACAGGCACCCACCACCAAGCCCAGCTAATGTTTGTATTTTCAGTAGAGACAAGGTCCCAACATGTTGGCCAGGCTGGTCTCAAACTCCTGACCTCAAATATCTGCCCACCTCGGCATCCCAAAGTGCTGAGATTATAGATGTGAGCCACTGCACCCAACCAGAACTCAGGAATTTTTGAGGGTGATCATTCAATGTCTCTCAAATTTCTTTGACAAGAGAATAGCATGAAGTTTAATGCTTGGATTAAAGCAGGAGGCAAATAATCATCTCAGATATTATTAATCACTGCAGATGTTAATCAAAATTAGGCTTATTTTTCAGGCTTAGATTTTATAACAAAGCAAAAAATGCTAAGGTAAGAAAAATATGCCTCATCAATTTTCTTTGCTATTAACAATCTTGAGAGAGTTATGTTCTATGGAACATAATGTCAGTAATATTGACCTAACCCCATATACTCATTTTGCATGTGAGGAAATTGGTTAGGAGTGGGAGAAGAGACAAAATAGTTCAATATATGGTAAATGAGAAACCAGGTATCTGCTTGACAGAATCATCTTTTTGATCCCTAAGCACAGATGGAAAGAAGACCCTCAAAAATCTATCTCCTGTCCCCCTCTCAGACCCTATTCCTTTACTCATCCCTGTACACTACTGGGACAGGTCACATACACATTCAGACCCCAGATCCTCCTCCACAAATTCAGAGACCCAAGCACCCACCAAATAGCTTATCATAGTGGCTTTTGGGGAAGGTCAACTCCATTCCTCCAAGGCTCCAGTTTGCCAGTCTTTTCATGAATGGGTAAGGAAAGTGTGTATTTGAGGCCATTAGCTTCTTTCCAAATGCATACATCTTCACTTTTACTCACCCTGCAGACACTCGGGAATCAGAACCCATCACAACGCCCCCGTCAAACTCCACTGCCATGATGGTGGTCTGCAGAGACACAGAATATGGAATGTCAGGGCAAGAACAGCCTTGATGCCCTCATGTTAGAGAAGAAGAAACATTCCCAGAGAGGCGAAGTGACTGGCTCAAAGATTACACAGTAACAGGCCAGAGCTGACTGTCAGTACAGGCTTTTTTTCCCTTCATCTTTCCACTTTCTCTATTGCTTCATCCGGCTGCAGGGGAATGCCACAGCCCAGCTGTGATACAACACAGAAAGAACTGTGTCCCTAAGTTCCAACTTGCCTAGTGGAATCCTCTCCACTGTAGAGAGGTGGAG …….. 19,000 base pairs omitted!!! Sequence Alignment

  7. Genes • Definition varies • To a geneticist a gene is the region of a chromosome that confers a particular trait • To a molecular biologist a gene is a sequence of DNA which encode a protein or RNA and includes all of the relevant regulatory sequences • Genes can be turned on/off, up/down • Not all parts of the DNA sequence are involved with the production of proteins Sequence Alignment

  8. Gene Structure Sequence Alignment

  9. Locate Genes Gn.Ex Type S .Begin ...End .Len Fr Ph I/Ac Do/T CodRg P.... Tscr.. ----- ---- - ------ ------ ---- -- -- ---- ---- ----- ----- ------ 2.00 Prom + 5833 5872 40 -14.22 2.01 Init + 6023 6620 598 1 1 57 87 371 0.621 27.34 2.02 Intr + 7157 7271 115 0 1 122 41 76 0.997 5.81 2.03 Intr + 7420 7550 131 1 2 83 9 151 0.979 7.04 2.04 Intr + 8510 8715 206 0 2 123 91 98 0.826 12.52 2.05 Intr + 9142 9339 198 0 0 69 39 276 0.998 20.35 2.06 Intr + 10541 10669 129 1 0 89 78 131 0.992 13.09 2.07 Intr + 10819 11007 189 0 0 137 87 125 0.999 17.08 2.08 Intr + 11567 11740 174 1 0 101 80 233 0.966 23.94 2.09 Intr + 11984 12146 163 1 1 103 78 108 0.999 10.85 2.10 Intr + 12455 12591 137 0 2 101 94 127 0.999 14.89 2.11 Intr + 13874 14050 177 1 0 113 16 188 0.043 14.22 2.12 Intr + 16570 16717 148 0 1 98 64 116 0.984 10.01 2.13 Intr + 16876 16987 112 2 1 103 115 134 0.999 17.04 2.14 Intr + 17396 17525 130 2 1 100 98 174 0.999 20.30 2.15 Intr + 17924 18128 205 1 1 81 55 295 0.999 24.27 2.16 Term + 18612 18700 89 1 2 96 39 148 0.999 8.52 2.17 PlyA + 18919 18924 6 1.05 Sequence Alignment

  10. Codons Alanine, Arginine, Aspartic Acid, Asparagine, Cystinine, Glutamic Acid, Glutamine, Glycine, Histidine, Isoleucine, Leucine, Lysine, Methionine, Phenylalanine, Proline, Serine, Threonine, Tryptophan, Tyrosine, Valine Sequence Alignment

  11. Conversion to Protein Sequence Alignment

  12. Predict Protein Sequence >18:21:37|GENSCAN_predicted_peptide_2|966_aa MASSRCPAPRGCRCLPGASLAWLGTVLLLLADWVLLRTALPRIFSLLVPTALPLLRVWAV GLSRWAVLWLGACGVLRATVGSKSENAGAQGWLAALKPLAAALGLALPGLALFRELISWG APGSADSTRLLHWGSHPTAFVVSYAAALPAAALWHKLGSLWVPGGQGGSGNPVRRLLGCL GSETRRLSLFLVLVVLSSLGEMAIPFFTGRLTDWILQDGSADTFTRNLTLMSILTIASAV LEFVGDGIYNNTMGHVHSHLQGEVFGAVLRQETEFFQQNQTGNIMSRVTEDTSTLSDSLS ENLSLFLWYLVRGLCLLGIMLWGSVSLTMVTLITLPLLFLLPKKVGKWYQLLEVQVRESL AKSSQVAIEALSAMPTVRSFANEEGEAQKFREKLQEIKTLNQKEAVAYAVNSWTTSISGM LLKVGILYIGGQLVTSGAVSSGNLVTFVLYQMQFTQAVEVLLSIYPRVQKAVGSSEKIFE YLDRTPRCPPSGLLTPLHLEGLVQFQDVSFAYPNRPDVLVLQGLTFTLRPGEVTALVGPN GSGKSTVAALLQNLYQPTGGQLLLDGKPLPQYEHRYLHRQVAAVGQEPQVFGRSLQENIA YGLTQKPTMEEITAAAVKSGAHSFISGLPQGYDTEVDEAGSQLSGGQRQAVALARALIRK PCVLILDDATSALDANSQLQVEQLLYESPERYSRSVLLITQHLSLVEQADHILFLEGGAI REGGTHQQLMEKKGCYWAMPTEFFQSLGGDGERNVQIEMAHGTTTLAFKFQHGVIAAVDS RASAGSYISALRVNKVIEINPYLLGTMSGCAADCQYWERLLAKECRLYYLRNGERISVSA ASKLLSNMMCQYRGMGLSMGSMICGWDKKGPGLYYVDEHGTRLSGNMFSTGSGNTYAYGV MDSGYRPNLSPEEAYDLGRRAIAYATHRDSYSGGVVNMYHMKEDGWVKVESTDVSDLLHQ YREANQ Sequence Alignment

  13. hello.exe 011110100100101111111110000001111110010101000011110100010001101010100011101110101101110100000010101100011010011101001110001010011010000001000100011101 100010001001011110010100111000101101101000010101010111001101000110000000101101001011111100000111011111001011100110011010101111111010101011111101111010 001010101000010101011100011110110110001010010110101011001111100010010110100110110110100011001000111101001111011110110101101111111110101010110011010011 001011011010101001010000011011101101100010110100001110110011001001000101101111001010111001101000011000110001110111011111111001101110000100110100011001 111010000011111011000010011010011110010010001100101110111001000100100011000001000111101100110100110101001111000111001100000010000010000101110110010110 111011100110110100010000010100001001001010001100010010000010101000000110111011111111111010110110000001010001011110000101011011101101000011000110010101 011001110001001100100101111100110001010000011000010010110011010000111011110100000011000110111000000010000101010011101101010010111000011110101100101101 111010001101011111011100000001011101010101000101110101011110000011001010000011001000000101011000101100011101101010100011111100010000010001001100011100 101111000101000101100110010001100100110010011010100100111001010001010000101000011001101110100001010001010010010111000101110010100000111100100011111100 001100110010001010001100000100100111111100001110001010111000110000011100011100001101100010000111111001111011111101001011000100100100011001011011111111 100011000000111101010000110010110110001011101110101111000001111101111110111001011010110111001010000101001010101001100100000111010001101100100101101001 101111111011000101011111010010011001010110111100111111011100100010000101110100000100111110001001111010110101000001000111101110001100100011010011100111 011101100101100001111001110111010010101001111111001011111101100101110110111000000100101100110110100010001001011011101011001100001000000110100010011001 000010110110111111101000011000010001010000111000010001111001000111111011110001011100100011010111001011101000111011101110000100001011000011010110001110 111001100101111010110111001010110010010010001010100010001100010000001001110101100011011010100111000010001010000011101011111110111010010101010110001101 010111110010011101010111010011110110000101011110001111000101010011101110101001101100011010001000000101001010011001000110101001010011010110010111100101 001110100001010101010111000100000001101110111000010001101001111111100010100010011111001000001101001110101000110010110010001010100100011111100010111010 101011101100111010101010101001001000110111110011001100100101001011000101100111000110110011100101100111010001010101010100111101011100011010011001101100 000001100100001000010110001100111001101010001001001110101100110010011101001101001011110110100111000110011001100001011011000111011111000000011010010100 111011000001101001001101010101000110110110100011101011001000100111100111111111010000000011100010100110010110010111110011011010100000111000001010010010 000110000101101110111010110010011010001111100110001000100001011011001000101011101111010111111100001111110100111001001000101100110100000100001100110101 101011010110100010111011111010101100101001110100101100011010001111110110001000110100111010110010100010010001100100110000101111001100100011000100010010 Sequence Alignment

  14. Find Similar Proteins gi|549042|sp|Q03518|TAP1_HUMAN ANTIGEN PEPTIDE TRANSPORTER ... 1221 0.0 gi|2506117|sp|P36370|TAP1_RAT ANTIGEN PEPTIDE TRANSPORTER 1... 831 0.0 gi|2506116|sp|P21958|TAP1_MOUSE ANTIGEN PEPTIDE TRANSPORTER... 827 0.0 gi|1172602|sp|P28062|PRCY_HUMAN PROTEASOME COMPONENT C13 PR... 455 e-127 Sequence Alignment

  15. Sequence Alignment • Sequences of genes and proteins are compared to infer • Structural, functional, and evolutionary relationships between the sequences • Over time evolution causes • Substitutions which change residues in a sequence • Insertions/deletions add or remove residues (gaps) Sequence Alignment

  16. Sequence Alignment • Aligning two sequences is the cornerstone of Bioinformatics • What are we looking for when aligning sequences? • Identity: Two sequence that have a certain number of positions in common at aligned positions • Similarity: Often a number of positions will be replaced by ones of similar chemical properties • Homology: Two sequences that are evolutionarily related and stem from a common ancestor Sequence Alignment

  17. Possible Alignments HEAGAWGHE-E -PA--W-HEAE HEAGAWGHE-E ---PAW-HEAE HEAGAWGHE-E --P-AW-HEAE substitution insertion deletion Sequence Alignment

  18. How Do We Choose? • Obviously there are many alignments, how do we choose? • The Biologists provide a scoring mechanism which can be used to determine which alignment is best • A simple scoring mechanism: • 0 for a match • 1 mismatch or gap • Lowest score is the best Sequence Alignment

  19. Alignments With Score HEAGAWGHE-E -PA--W-HEAE Score: 6 Sequence Alignment

  20. Computing The Result • We know how to generate different alignments • We know how to score alignments • One algorithm might be to generate all possible alignments and choose the one with the best score • Not feasible!! Sequence Alignment

  21. Alignment • At every step in the process of aligning two sequences, S1 and S2, you have to make one of three decisions • Match/mismatch • Add gap to S1 • Add gap to S2 • For example, when aligning ABC with BC ABC-ABC ABC BCBC-BC Sequence Alignment

  22. Using An Array • An array could be used to keep track of the possible moves • Diagonal: match/mismatch • Across: gap in sequence 2 • Down: gap in sequence 1 Sequence Alignment

  23. Example Sequence Alignment

  24. Example Sequence Alignment

  25. Example A- bc -B c Sequence Alignment

  26. Example A bc B c Sequence Alignment

  27. Example -A bc B- c Sequence Alignment

  28. Example A- bc (2) -B c Abc (1) Bc -A bc (2) B- c Sequence Alignment

  29. Example AB- c --B c Sequence Alignment

  30. Example AB c -B c Sequence Alignment

  31. Example A-Bc -B-c ABc B-c -ABc -B-c Sequence Alignment

  32. Example AB- c (3) --B c A-Bc (3) -B-c ABc (2) B-c -ABc (2) -B-c AB c (1) -B c Sequence Alignment

  33. Dynamic Programming • The word Programming in the name has nothing to do with writing computer programs. • Mathematicians use the word to describe a set of rules which anyone can follow to solve a problem. • They do not have to be written in a computer language. • Dynamic programming was the brainchild of an American Mathematician Richard Bellman • Store the results for small sub-problems and looks them up, rather than recomputing them, when they are needed later to solve larger sub-problems Sequence Alignment

  34. Fib(6) 6 4 5 2 3 3 4 0 1 1 2 1 2 2 3 0 1 0 1 0 1 1 2 0 1 Sequence Alignment

  35. Keep Only the Best Match • There is no need to remember every possible extension • We only need to keep the best one • There may be ties which is okay • How do you determine which one is the best? Sequence Alignment

  36. Keep Only the Best Match • There is no need to remember every possible extension • We only need to keep the best one • There may be ties which is okay • How do you determine which one is the best? • Ask a Biologist!!!! Sequence Alignment

  37. Example A- Score==2 -B A Score==1 B -A Score==2 B- Sequence Alignment

  38. Example AB- Score==3 --B AB Score==1 -B AB Score==2 B- Sequence Alignment

  39. Example ABC- Score==4 ---B ABC Score==3 --B ABC Score==2 -B- Sequence Alignment

  40. Example Note: Down or across always adds one to the score Diagonal will add either one or a zero depending on whether or not the bases match ABC Score==1 -BC Sequence Alignment

  41. Dynamic Programming • Three steps • Initialization • Matrix Fill • Traceback • These steps apply in general whether you are doing sequence alignment or folding predictions Sequence Alignment

  42. Recurrence • A mathematical relationship that defines fn as some combination of fi with i<n • For our next alignment we will use the recurrence • Mi,j = Maximum of • Mi-1,j-1 + Si,j • Mi,j-1 + w (gap in sequence 1) • Mi-1,j + w (gap in sequence 2) • Where • Si,j (1 if match, 0 if mismatch) • w = 0 (gap penalty) Sequence Alignment

  43. Initialization Sequence Alignment

  44. Matrix Fill Sequence Alignment

  45. Matrix Fill Sequence Alignment

  46. Matrix Fill Sequence Alignment

  47. Matrix Fill Sequence Alignment

  48. Matrix Fill Sequence Alignment

  49. Matrix Fill Sequence Alignment

  50. Matrix Fill Sequence Alignment

More Related