220 likes | 272 Views
Introduction to Bioinformatics - Tutorial no. 5. MEME – Discovering motifs in sequences MAST – Searching for motifs in databanks TRANSFAC – the Transcription Factor DB. WebLogo - Input. Aligned Sequences (e.g. output of ClulatlW). http://weblogo.berkeley.edu. RUN !. Genes:. Proteins:.
E N D
Introduction to Bioinformatics - Tutorial no. 5 MEME – Discovering motifs in sequences MAST – Searching for motifs in databanks TRANSFAC – the Transcription Factor DB
WebLogo - Input Aligned Sequences (e.g. output of ClulatlW) • http://weblogo.berkeley.edu RUN !
Genes: Proteins: WebLogo - Output
MEME • http://meme.sdsc.edu/ • Motif discovery from unaligned sequences • Genomic or protein sequences • Identifies profile motifs • Multiple motifs for any input • Flexible model of motif presence • Motif can be absent in some sequences • Can appear several times in one sequence
MEME Input Email address Multiple input sequences Range of motif lengths How many motifs? How many times in each sequence? How many times in total?
MEME Output (1) Like BLAST Motif length “Position-Specific Probability Matrix” = Motif Profile Number of times Most popular symbols Diversion of motif position from background
MEME Output (2) Position in sequence Strength of match Sequence names Reverse complement (genomic input only) Motif within sequence
MEME Output (3) Motif instance Original sequence lengths Overall strength of motif matches
MAST • Searches for motifs (one or more) in sequence databases: • Like BLAST but motifs for input • Similar to iterations of PSI-BLAST • Profile defines strength of match • Multiple motif matches per sequence • Combined E value for all motifs • MEME uses MAST to summarize results: • Each MEME result is accompanied by the MAST result for searching the discovered motifs on the given sequences.
MAST Input Email address Motif file (e.g. MEME output) Consider matched sequence length Database (like BLAST) E value threshold
MAST Output (1) Link to GenBank Matched accession Match E value Length of sequence
MAST Output (2) Motif diagram
MAST Output (3) Position of each instance Matched parts of sequence Motif and orientation P value of instance Motif ‘consensus’
TRANSFAC Database of eukaryotic DNA transcription regulation: • Individual regulatory sites (SITES table) • Genes to which they belong • Proteins which bind them • Proteins which bind sites (FACTORS table) • Cellular source of protein • Nucleotide motif profile for binding • Some grouping and classification • Classification of factors (CLASS table) • Position-specific matrices for select factors (MATRIX table) • Cell localization (CELL table)
Searching TRANSFAC • www.gene-regulation.com • Search a single table • By identifier, factor name, gene name • By species, author • Browse your way from table to table • Search within a sequence • MatInspector, TFScan (EMBOSS package)
TRANSFAC Factor (3) DT Date; author FA Factor name GE Encoding gene SF Structural features CP Cell specificity (positive) CN Cell specificity (negative) EX Expression pattern FF Functional features IN Interacting factors MX Matrix BS Binding SITE DR External databases References: RN Reference no. RX MEDLINE ID RA Reference authors RT Reference title RL Reference data
TRANSFAC Matrix Accession Position Specific Matrix Concensus (IUPAC subset symbols) Statistical basis
TRANSFAC complementary material at your free time
TRANSFAC Site (1) DNA or RNA Accession number Sequence of regulatory element Gene Gene region Position range of factor binding site
TRANSFAC Site (2) Binding factor accession Factor name Binding ‘quality’ Organism Cellular source External links Methods of identifying site
TRANSFAC Factor (1) AC: Accession number FA: Factor name HO: Homologs SX: Other names CL: Classification OS: Organism OC: Taxonomy SX: Amino acid sequence SZ: Size
TRANSFAC Factor (2) Protein sequence reference Features and positions Structural features Cell specificity