Bioinformatyka

Wykład Bioinformatyka

Wykład 4, 2008

Bioinformatyka

Wykład 4.

E. Banachowicz

Zakład Biofizyki Molekularnej

IF UAM

http://www.amu.edu.pl/~ewas

Formaty danych - GenBank

Poco wprowadza się dane do komputerów?
1. żeby je pobrać 
2. żeby coś odkryć

Jeśli baza danych nie pozwala na wyszukanie potrzebnej
informacji, jest by bezużyteczna. (Nawet największa baza!)

Wykład Bioinformatyka

Wykład 4, 2008

Formaty danych - GenBank

1. dane muszą mieć jednoznaczną strukturę

i zdefiniowane powiązania

2. dane muszą być stabilne

Model danych w NCBI oparty jest na sekwencji DNA
(stabilność), i daje możliwość śledzenia informacji od literatury
do sekwencji.

Stabilność danych

sekwencja

DNA

literatura

PubMed

sekwencja

białka

struktura

4 podstawowe dane

Full Text

el. Journal

OMIM

Mapy i

genomy

Taksonomia

Wykład Bioinformatyka

Wykład 4, 2008

pliki ASCII

• większość programów do analizy sekwencji

nie akceptuje znaków spoza zestawu ASCII

(różna interpretacja, problemy z transferem)

• Poza sekwencją DNA lub białka (raw sequence) –

odpowiedni format

• Kod DNA i białka ujednolicony został przez

NC-IUB (

Nomenclature Committee of the International Union of

Biochemistry and Molecular Biology -

http://www.chem.qmul.ac.uk/iubmb/

Zasady/kwasy nukleinowe

- ujednolicony kod

 adenozyna

M  A C (amino)

 cytozyna

S  G C (strong)

 guanina

W  A T (weak)

 tymidyna

B  G T C

 urydyna

D  G A T

R  G A (puryna)

H  A C T

Y  T C (pyrymidyna)

V  G C A

K  G T (keto)

N  A G C T (dowolna)
- 

gap of indeterminate length

NC – IUP (1984)

Wykład Bioinformatyka

Wykład 4, 2008

Standardowy kod aminokwasów

A Ala alanina

Pro prolina

Asx kw. asparaginowy/asparagina Q

Gln glutamina

C Cys cysteina

Arg arginina

D Asp kw. asparaginowy

Ser seryna

E Glu kw. glutaminowy

Thr treonina

F Phe fenyloanina

selenocysteina

G Gly glicyna

V Val

walina

H His histydyna

W Trp tryptofan

I Ile izoleucyna

Y Tyr tyrozyna

K Lys lizyna

Glx kw.glutaminowy/glutamina

L Leu leucyna

Xxx dowolny

M Met metionina

* stop translacji

N Asn asparagina

- gap of indeterminate length

Abstract Syntax Notation Sequence Format ASN.1

ASN.1 (skrót od Abstract Syntax Notation One
abstrakcyjna notacja składniowa numer jeden)

język opisu danych przejęty i rozwijany przez NCBI

Wykład Bioinformatyka

Wykład 4, 2008

Integracja danych z wielu różnych źródeł

• np. PubMed

(np. wyszukiwanie według autorów)

Tag

Name

Abstract

Affiliation

AID

Article Identifier

Author

CIN

Comment In

Corporate Author

CON

Comment On

CRF

Corrected and republished from

CRI

Corrected and republished in

Date Created

DCOM

Date Completed

DEP

Date of Electronic Publication

Publication Date

EDAT

Entrez Date

EFR

Erratum For

EIN

Erratum In

FAU

Full Author Name

FIR

Full Investigator

FPS

Full Personal Name as Subject

General Note

Grant Number

Gene Symbol

Issue

Investigator

IRAD

Investigator Affiliation

ISSN

JID

NLM Unique ID

Full Journal Title

Language

LID

Location ID

MEDLINE Display

Last Revision Date

MeSH Terms

MHDA

MeSH Date

OAB

Other Abstract

OCI

Other Copyright Information

OID

Other ID

ORI

Original Report In

Other Term

OTO

Other Term Owner

OWN

Owner

Pagination

PHST

Publication History Status Date

Place of Publication

PMID

PubMed Unique Identifier

PRIN

Partial Retraction In

PROF

Partial Retraction Of

Personal Name as Subject

PST

Publication Status

Publication Type

PUBM

Publishing Model

Number of References

RIN

Retraction In

EC/RN Number

ROF

Retraction Of

RPF

Republished From

RPI

Republished In

Subset

SFM

Space Flight Mission

Secondary Source Identifier

Source

SPIN

Summary For Patients In

STAT

Status Tag

Journal Title Abbreviation

Title

Transliterated Title

UIN

Update In

UOF

Update Of

Volume

Wykład Bioinformatyka

Wykład 4, 2008

streszczenie

brak streszczenia

dostępny w PMC

dostępny pełen teskt

cytowanie

nr stron

autorzy

tytuł

identyfikator

czasopismo

data publikacji

DDBJ/GenBank/EMBL

Podobna struktura i identyfikatory: A12345=A12345

PIR/ Swiss-Prot

Różne identyfikatory: A12345≠A12345

Seq-id klasa obiektów

Wykład Bioinformatyka

Wykład 4, 2008

GenBank:

http://www.ncbi.nlm.nih.gov/

nazwa lokusa
(locus)

długość i typ
sekwencji

klasyfikacja
organizmu

data
wprowadzenia

nazwa lokusa
(locus)

długość i typ
sekwencji

klasyfikacja
organizmu

data
wprowadzenia

Wykład Bioinformatyka

Wykład 4, 2008

GenBank:

http://www.ncbi.nlm.nih.gov/

opis objektu

ACCESSION

numer dostępu do oryginalnego

źródła
VERSION numer kolejnej wersji
KEYWORDS

słowa kluczowe (cross reference)

SOURCE

organizm, z którego pochodziło DNA

ORGANISM opis organizmu
REFERENCE bibliografia

GenBank:

http://www.ncbi.nlm.nih.gov/

COMMENT np.funkcja biologiczna
FEATURES informacje o sekwencji przez

podanie położenia zasad lub przedziału położeń

sourece, misc_signal, mRNA, CDS,
intron, mutation

ORIGIN

początek sekwencji

koniec sekwencji

Wykład Bioinformatyka

Wykład 4, 2008

EMBL:

http://www.ebi.ac.uk/embl/index.html/

European Molecular Biology Laboratory

Wygląd strony w 2006

Wykład Bioinformatyka

Wykład 4, 2008

EMBL:

http://www.ebi.ac.uk/embl/index.html/

European Molecular Biology Laboratory

numer identyfikacyjny w bazie danych

numer dostępowy do pierwotnej sekwencji

wersja

data wprowadzenia lub modyfikacji

opis

OS,OC

– organizm pochodzenia DNA

(

,…) bibliografia

– informacje o sekwencji (FEATUREs)

- początek i koniec sekwencji

Wykład Bioinformatyka

Wykład 4, 2008

Format sekwencji FASTA

>embl|DQ423612|DQ423612 Influenza A virus (A/Cygnus olor/Astrakhan/Ast05-2-
10/2005(H5N1)) polymerase basic protein 1 (PB1) gene, complete cds. ...
caaaccatttgaatggatgtcaatccgactttacttttcttgaaagtaccagtgcaaaat
gctataagtaccacattcccttatactggagaccctccatacagccatgggacagggaca
ggatacaccatggacacagtcaacagaacacaccaatattcagaaaaggggaagtggaca
acaaacacagagactggagcaccccaactcaacccgattgatggaccactacctgaggat
aatgagcccagtggttatgcacaaacagattgtgtattggaagcaatggctttccttgaa
gaatcccacccagggatctttgaaaactcgtgtcttgaaacgatggaaattgttcaacaa
acaagagtggataaactgacccaaggtcgtcagacctatgactggacattgaatagaaac
caaccggctgcaaccgctttggccaacactatagaaatcttcagatcgaacggtctaaca
gccaatgaatcgggacggctaatagatttcctcaaggatgtgatggaatcaatggataag
gaagaaatggagataacaacacacttccagagaaagagaagagtgagagacaacatgacc
aaaaagatggtcacacaaagaacaatagggaagaaaaagcaaaggctgaacaaaaagagc
tacctgataagagcactgacactgaatacaatgacaaaagatgcagaaagaggcaaattg
aagaggcgagcaattgcaacacccggaatgcaaatcagaggattcgtgtactttgttgaa
acattagcgaggagtatctgtgagaaacttgagcaatctggactcccagttggagggaat
gaaaagaaggctaaattggcaaacgtcgtgaggaagatgatgactaactcacaagatact
gaactctcctttacaattactggagacaatactaaatggaatgagaatcagaatcctagg

> jednoliniowy opis
wszystkie linie tekstu nie powinny być dłuższe niż 80 znaków

>gi|89213215|gb|ABD64049.1| polymerase basic protein 1 [Influenza A virus (A/Cygnus olor/Astrakhan/Ast05-2-
10/2005(H5N1))]
MDVNPTLLFLKVPVQNAISTTFPYTGDPPYSHGTGTGYTMDTVNRTHQYSEKGKWTTNTETGAPQLNPID
GPLPEDNEPSGYAQTDCVLEAMAFLEESHPGIFENSCLETMEIVQQTRVDKLTQGRQTYDWTLNRNQPAA
TALANTIEIFRSNGLTANESGRLIDFLKDVMESMDKEEMEITTHFQRKRRVRDNMTKKMVTQRTIGKKKQ
RLNKKSYLIRALTLNTMTKDAERGKLKRRAIATPGMQIRGFVYFVETLARSICEKLEQSGLPVGGNEKKA
KLANVVRKMMTNSQDTELSFTITGDNTKWNENQNPRMFLAMITYITRNQPEWFRNVLSIAPIMFSNKMAR
LGRGYMFESKSMKLRTQIPAEMLANIDLKYFNELTKKKIEKIRPLLIDGTASLSPGMMMGMFNMLSTVLG
VSILNLGQKRYTKTTYWWDGLQSSDDFALIVNAPNHEGIQAGVDRFYRTCKLVGINMSKKKSYINRTGTF
EFTSFFYRYGFVANFSMELPSFGVSGINESADMSIGVTVIKNNMINNDLGPATAQMALQLFIKDYRYTYR
CHRGDTQIQTRRSFELKKLWEQTRSKAGLLVSDGGPNLYNIRNLHIPEVCLKWELMDEDYQGRLCNPLNP
FVSHKEIESVNNAVVMPAHGPAKGMEYDAVATTHSWIPKRNRSILNTSQRGILEDEQMYQKCCNLFEKFF
PSSSYRRPVGISSMVEAMVSRARIDARIDFESGRIKKEEFAEIMKICSTIEELRRPK

Wykład Bioinformatyka

Wykład 4, 2008

NBRF/PIR (National Biomedical Research

Foundation/Protein Information Resource

http://www-nbrf.georgetown.edu/pirwww/dbinfo/

>P1;gi|89213215|gb|ABD64049_1|
gi|89213215|gb|ABD64049_1| 757 bases

MDVNPTLLFL KVPVQNAIST TFPYTGDPPY SHGTGTGYTM DTVNRTHQYS
EKGKWTTNTE TGAPQLNPID GPLPEDNEPS GYAQTDCVLE AMAFLEESHP
GIFENSCLET MEIVQQTRVD KLTQGRQTYD WTLNRNQPAA TALANTIEIF
RSNGLTANES GRLIDFLKDV MESMDKEEME ITTHFQRKRR VRDNMTKKMV
TQRTIGKKKQ RLNKKSYLIR ALTLNTMTKD AERGKLKRRA IATPGMQIRG

READSEQ – konwersja formatów

http://www.ebi.ac.uk/cgi-bin/readseq.cgi

Wykład Bioinformatyka

Wykład 4, 2008

Name

Read

Write

Int'leaf

Features

Sequence

Content-type

Suffix

GenBank|gb

yes

biosequence/genbank

.gb

EMBL|em

yes

biosequence/embl

.embl

Pearson|Fasta|fa

yes

biosequence/fasta

.fasta

GCG

yes

biosequence/gcg

.gcg

MSF

yes

biosequence/msf

.msf

Clustal

yes

biosequence/clustal

.aln

NBRF

yes

biosequence/nbrf

.nbrf

PIR|CODATA

yes

biosequence/codata

.pir

ACEDB

yes

biosequence/acedb

.ace

Phylip3.2

yes

biosequence/phylip2

.phylip2

Phylip|Phylip4

yes

biosequence/phylip

.phylip

Plain|Raw

yes

biosequence/plain

.seq

PAUP|NEXUS

yes

biosequence/nexus

.nexus

XML

yes

biosequence/xml

.xml

FlatFeat|FFF

yes

biosequence/fff

.fff

GFF

yes

biosequence/gff

.gff

BLAST

yes

biosequence/blast

.blast

Pretty

yes

biosequence/pretty

.pretty

SCF

yes

biosequence/scf

.scf

DNAStrider

yes

biosequence/strider

.strider

IG|Stanford

yes

biosequence/ig

.ig

Fitch

yes

biosequence/fitch

.fitch

ASN.1

yes

biosequence/asn1

.asn

znane formaty sekwencji

Anatomia danych SwissProt/TrEMBL

http://www.expasy.ch/

Wykład Bioinformatyka

Wykład 4, 2008

Wykład Bioinformatyka

Wykład 4, 2008

MeCP2

NCBI
•

http://www.ncbi.nlm.nih.gov/entrez/viewer.fcgi?d
b=protein&val=1708973

EMBL-EBI
•

http://www.ebi.ac.uk/

PIR
•

http://pir.georgetown.edu/cgi-bin/textsearch.pl

ReadSeq
•

http://www.ebi.ac.uk/cgi-bin/readseq.cgi

PDB

Wykład Bioinformatyka

Wykład 4, 2008

plik PDB

Wykład Bioinformatyka

Wykład 4, 2008

plik PDB

Ser

Lys

Val

Wykład Bioinformatyka

Wykład 4, 2008

plik PDB

Identyfikacja sekwencji w BD

• Identyfikacja przez porównanie z innymi

sekwencjami

Zestawienia sekwencji = uliniowienie =
=porównanie = alignment

Wykład Bioinformatyka

Wykład 4, 2008

Porównywanie sekwencji

• Pierwsze pytanie biologa molekularnego,

kiedy odkryje nową sekwencję:

Czy w bazie sekwencji są już sekwencje podobne

do mojej?



sekwencje są identyczne –

nic nowego….

 sekwencja jest podobna (ma „krewnych”)

– nowy członek

znanej rodziny

 sekwencja ma kilka podobnych regionów, motywów lub

domen

– można zaproponować funkję

 Nie ma znaczącego podobieństwa

– dużo pracy…..

Porównywanie sekwencji

• Celem porównania białek jest między

innymi przypisanie informacji znanej dla
jednej cząsteczki drugiej cząsteczce

Wykład Bioinformatyka

Wykład 4, 2008

Pokrycie sekwencji

• dopasowanie globalne – dopasowanie

wzdłuż całej sekwencji

(zastosowanie: do białek

składających się z pojedynczej domeny lub
homologicznych słabo zróżnicowanych)

• dopasowanie lokalne – uwzględnia

domenową naturę białek, szuka
subsekwencji (

zastosowanie: do białek

wielodomenowych, mRNA z sekwencją genomową

)

BLAST

Wykład Bioinformatyka

Wykład 4, 2008

Wykład Bioinformatyka

Wykład 4, 2008

CDN …..

…..na ćwiczeniach