Transcript Bioinformatics Tools

Introduction to Bioinformatics
236523
Lecturer: Dr. Yael Mandel-Gutfreund
Teaching Assistance:
Shula Shazman
Sivan Bercovici
Course web site :
http://webcourse.cs.technion.ac.il/236523
What is Bioinformatics?
2
Course Objectives
• To introduce the bioinfomatics discipline
• To make the students familiar with the major
biological questions which can be addressed
by bioinformatics tools
• To introduce the major tools used for
sequence and structure analysis and explain
in general how they work (limitation etc..)
3
Course Structure and Requirements
1.Class Structure
1.
2.
2 hours Lecture
1 hour tutorial
2. Home work
•
Homework projects will be given every third week
•
The homework will be done in pairs.
•
4/4 homework projects submitted
2. A final project will be conducted and submitted
in pairs
4
Grading
• 30 % Homework assignments
• 70% final project
5
Literature list
• Gibas, C., Jambeck, P. Developing Bioinformatics
Computer Skills. O'Reilly, 2001.
• Lesk, A. M. Introduction to Bioinformatics. Oxford
University Press, 2002.
• Mount, D.W. Bioinformatics: Sequence and Genome
Analysis. 2nd ed.,Cold Spring Harbor Laboratory
Press, 2004.
Advanced Reading
Jones N.C & Pevzner P.A. An introduction to
Bioinformatics algorithms MIT Press, 2004
6
What is Bioinformatics?
7
What is Bioinformatics?
“The field of science in which biology, computer
science, and information technology merge to
form a single discipline”
Ultimate goal: to enable the discovery of new
biological insights as well as to create a global
perspective from which unifying principles in
biology can be discerned.
8
Central Paradigm in Molecular Biology
Gene (DNA)
mRNA
Protein
21ST centaury
Genome
Transcriptome
Proteome
9
from purely lab-based science to an information science
Bioinformatics
Bio = Informatics
10
From DNA to Genome
Watson and Crick
DNA model
First protein
sequence
1955
1960
First protein
structure
1965
1970
1975
1980
1985
11
1990
First bacterial
genome
1995
Hemophilus Influenzae
Yeast genome
2000
First human
genome draft
12
Complete Genomes
Total
2008 2007
706 456
Eukaryotes
78
43
Bacteria
578 383
Archaea
50
29
13
What’s Next ?
The “post-genomics” era
Annotation
Comparative
genomics
Structural
genomics
Functional
genomics
Goal:
to understand the living cell
14
Annotation
CCTGACAAATTCGACGTGCGGCATTGCATGCAGACGTGCATG
CGTGCAAATAATCAATGTGGACTTTTCTGCGATTATGGAAGAA
CTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTC
AGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGA
AGAGCCAGTAAAAGACGCAGTGACGGAGATGTCTGATG CAA
TAT GGA CAA TTG GTT TCT TCT CTG AAT ......
.............. TGAAAAACGTA
15
Identify the genes within a
given sequence of DNA
Identify the sites
Which regulate the gene
Annotation
Predict the function
16
TF binding site
CCTGACAAATTCGACGTGCGGCATTGCATGCAGACGTGCATG
CGTGCAAATAATCAATGTGGACTTTTCTGCGATTATGGAAGAA
CTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTC
AGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGA
AGAGCCAGTAAAAGACGCAGTGACGGAGATGTCTGATG CAA
TAT GGA CAA TTG GTT TCT TCT CTG AAT .................................
Transcription
Start Site
promoter
.............. TGAAAAACGTA
ORF=Open Reading Frame Ribosome binding Site
CDS=Coding Sequence
17
Comparative
genomics
Human ATAGCGGGGGGATGCGGGCCCTATACCC
Chimp ATAGGGG - - GGATGCGGGCCCTATACCC
Mouse ATAGCG - - - GGATGCGGCGC -TATACCA
18
Perhaps not surprising!!!
How humans
are chimps?
Comparison between the full drafts of the human and chimp genomes
revealed that they differ only by 1.23%
19
Researchers have learned a great deal about the
function of human genes by examining their
counterparts in simpler model organisms such as
the mouse.
Conservation of the IGFALS (Insulin-like growth factor)
Between human and mouse.
20
Functional
genomics
21
Understanding the function of genes and other parts of the genome
22
A network of interactions can be built
for all proteins in an organism
A large network of 8184 interactions among 4140 S. Cerevisiae
proteins
23
Structural
Genomics
24
Assigning the structures of all proteins
protein complexes
Biologic processes
Shape and electrostatics
Active sites
fold
Evolutionary
relationship
Protein-ligand complexes
Functional sites
25
Resources and Databases
The different types of data are collected in
database
– Sequence databases
– Structural databases
– Databases of Experimental Results
All databases are connected
26
Sequence databases
•
•
•
•
Gene database
Genome database
SNPs database
Disease related mutation database
27
Gene database
• Give information into gene structure and
function
• Alternative splicing of genes
– Alternative pattern of exons included to create
gene product
28
Genome Databases
• Data organized by species
• Clones assembled into contigous pieces
‘contigs’ or whole chromosomes
• Information on non-coding regions
• Relativity
29
Genome Browsers
• Annotation adds value to sequence
• Easy “walk” through the genome
• Comparative genomics
30
Genome Browsers
• UCSC Genome Browser
http://genome.ucsc.edu/
• Ensembl Genome Browser (http://www.ensembl.org)
• WormBase: http://www.wormbase.org/
• AceDB: http://www.acedb.org/
• Comprehensive Microbial Resource:
http://www.tigr.org/tigr-scripts/CMR2/CMRHomePage.spl
• FlyBase: http://flybase.bio.indiana.edu/
31
SNP database
Single Nucleotide Polymorphisms (SNPs)
• Single base difference in a single position
among two different individuals of the same
species
• Play an important role in differentiation and
disease
32
Sickle Cell Anemia
• Due to 1 swapping an A for a T, causing inserted amino
acid to be valine instead of glutamine in hemoglobin
Image source: http://www.cc.nih.gov/ccc/ccnews/nov99/
33
Healthy Individual
>gi|28302128|ref|NM_000518.4| Homo sapiens hemoglobin, beta (HBB), mRNA
ACATTTGCTTCTGACACAACTGTGTTCACTAGCAACCTCAAACAGACACCATGGTGCATCTGACTCCTGA
GGAGAAGTCTGCCGTTACTGCCCTGTGGGGCAAGGTGAACGTGGATGAAGTTGGTGGTGAGGCCCTGGGC
AGGCTGCTGGTGGTCTACCCTTGGACCCAGAGGTTCTTTGAGTCCTTTGGGGATCTGTCCACTCCTGATG
CTGTTATGGGCAACCCTAAGGTGAAGGCTCATGGCAAGAAAGTGCTCGGTGCCTTTAGTGATGGCCTGGC
TCACCTGGACAACCTCAAGGGCACCTTTGCCACACTGAGTGAGCTGCACTGTGACAAGCTGCACGTGGAT
CCTGAGAACTTCAGGCTCCTGGGCAACGTGCTGGTCTGTGTGCTGGCCCATCACTTTGGCAAAGAATTCA
CCCCACCAGTGCAGGCTGCCTATCAGAAAGTGGTGGCTGGTGTGGCTAATGCCCTGGCCCACAAGTATCA
CTAAGCTCGCTTTCTTGCTGTCCAATTTCTATTAAAGGTTCCTTTGTTCCCTAAGTCCAACTACTAAACT
GGGGGATATTATGAAGGGCCTTGAGCATCTGGATTCTGCCTAATAAAAAACATTTATTTTCATTGC
>gi|4504349|ref|NP_000509.1| beta globin [Homo sapiens]
EEKSAVTALWGKVNVDEVGGEALGRLLVVYPWTQRFFESFGDLSTPDAVMGNPKVKAHGKKVLG
MVHLTP
AFSDGLAHLDNLKGTFATLSELHCDKLHVDPENFRLLGNVLVCVLAHHFGKEFTPPVQAAYQKVVAGVAN
ALAHKYH
34
Diseased Individual
>gi|28302128|ref|NM_000518.4| Homo sapiens hemoglobin, beta (HBB), mRNA
ACATTTGCTTCTGACACAACTGTGTTCACTAGCAACCTCAAACAGACACCATGGTGCATCTGACTCCTGA
GGTGAAGTCTGCCGTTACTGCCCTGTGGGGCAAGGTGAACGTGGATGAAGTTGGTGGTGAGGCCCTGGGC
AGGCTGCTGGTGGTCTACCCTTGGACCCAGAGGTTCTTTGAGTCCTTTGGGGATCTGTCCACTCCTGATG
CTGTTATGGGCAACCCTAAGGTGAAGGCTCATGGCAAGAAAGTGCTCGGTGCCTTTAGTGATGGCCTGGC
TCACCTGGACAACCTCAAGGGCACCTTTGCCACACTGAGTGAGCTGCACTGTGACAAGCTGCACGTGGAT
CCTGAGAACTTCAGGCTCCTGGGCAACGTGCTGGTCTGTGTGCTGGCCCATCACTTTGGCAAAGAATTCA
CCCCACCAGTGCAGGCTGCCTATCAGAAAGTGGTGGCTGGTGTGGCTAATGCCCTGGCCCACAAGTATCA
CTAAGCTCGCTTTCTTGCTGTCCAATTTCTATTAAAGGTTCCTTTGTTCCCTAAGTCCAACTACTAAACT
GGGGGATATTATGAAGGGCCTTGAGCATCTGGATTCTGCCTAATAAAAAACATTTATTTTCATTGC
>gi|4504349|ref|NP_000509.1| beta globin [Homo sapiens]
VEKSAVTALWGKVNVDEVGGEALGRLLVVYPWTQRFFESFGDLSTPDAVMGNPKVKAHGKKVLG
MVHLTP
AFSDGLAHLDNLKGTFATLSELHCDKLHVDPENFRLLGNVLVCVLAHHFGKEFTPPVQAAYQKVVAGVAN
ALAHKYH
35
Disease Databases
• Genes are involved in disease
• Many diseases are well studied
• Description of diseases and what is known
about them is stored
36
Structure Databases
• 3-dimensional structures of proteins, nucleic
acids, molecular complexes etc
• 3-d data is available due to techniques such
as NMR and X-Ray crystallography
37
38
Databases of Experimental Results
• Data such as experimental microarray
images- expression data
• Proteomic data
• Metabolic pathways, protein-protein
interaction data, regulatory networks
• ETC………….
39
Literature Databases
PubMed
http://www.ncbi.nlm.nih.giv/PubMed
Service of the National Library of Medicine
40
Putting it all Together
• Each Database contains specific
information
• Like other biological systems also these
databases are interrelated
41
PROTEIN
PIR
DISEASE
ASSEMBLED
GENOMES
LocusLink
SWISS-PROT
OMIM
GoldenPath
OMIA
WormBase
MOTIFS
TIGR
BLOCKS
Pfam
GENOMIC DATA
Prosite
GenBank
ESTs
dbEST
DDBJ
GENES
EMBL
RefSeq
unigene
AllGenes
SNPs
GENE
EXPRESSION
dbSNP
STRUCTURE
PDB
MMDB
SCOP
PATHWAY
Stanford MGDB
KEGG
NetAffx
COG
ArrayExpress
GDB
LITERATURE
PubMed
42