680 likes | 883 Views
Exercises. Pairwise alignment Homology search (BLAST) Multiple alignment (CLUSTAL W) Iterative Profile Search: Profile Search Pfam Prosite PSI-BLAST SAM. Overview. Query Sequence Unknown. Blast. Blast Sequence to search for close homologs. ClustalW. PFAM PROSITE.
E N D
Exercises • Pairwise alignment • Homology search (BLAST) • Multiple alignment (CLUSTAL W) • Iterative Profile Search: • Profile Search • Pfam • Prosite • PSI-BLAST • SAM
Overview Query Sequence Unknown Blast Blast Sequence to search for close homologs ClustalW PFAM PROSITE Make a multiple sequence alignment Search pFAM, Prosite for conserved motifs HMMer, PSSM PSI-blast Generate profile or HMM You detected homology with an annotated protein family Profile Search Search database for remote homologs
Terminal Oxidases Cytc e- Fe Fe Fe OUT e- e- CuB Fe Fe IN • Unknown protein is a heme cupper oxidase • Enzyme that reduces O2 to H2O in respiratory chain O2 I H2O • Subunit contains 2 hemes and a Cu prosthetic group • The residues that are ligands of these groups have been conserved in all types of terminal oxidase complexes
Terminal Oxidases O2 cytbc1 ? Cytc H2O NADH Nadh.dh e- cytochrome c oxidase e- e- quinol O2 succinate succ.dh ? H2O e- quinol oxidase e- e- O2 H2O H+ H+
Multiple Alignment Multiple alignment: standard gap cost Ligands Cu center Prosite pattern Ligands hemes
Multiple Alignment Ligands Cu center Multiple alignment: large gap cost Prosite pattern Ligands hemes
Phylogenetic Tree Tree based on subselection
Prosite Prosite domain
Pfam COX family
Overview Query Sequence Unknown Blast Blast Sequence to search for close homologs PFAM PROSITE Make a multiple sequence alignment Search pFAM, Prosite for conserved motifs HMMer, PSSM PSI-blast Generate profile or HMM You detected homology with an annotated protein family Profile Search Search database for remote homologs
PSI-BLAST • PSI BLAST • Start from a single sequence • Blast it against NCBI • Select high scoring hits • Perform multiple alignment • Construct profile • Iterate and find remote homologs • Usually cut the sequence in pieces • Avoid to give as input multi domain proteins
SAM http://www.cse.ucsc.edu/research/compbio/HMM-apps/HMM-applications.html