Difference between revisions of "Reference Sequence BCKDHA"
(→Sequence) |
(→Sequence Alignments) |
||
Line 31: | Line 31: | ||
[[Sequence_Alignments]] |
[[Sequence_Alignments]] |
||
− | |||
− | == Sequence Alignments == |
||
− | |||
− | === Sequence searches === |
||
− | * FASTA |
||
− | ../bin/fasta36 sequence.fasta database > FastaOutput.txt |
||
− | |||
− | * BLAST |
||
− | blastall -p blastp -d database -i sequence.fasta > BlastOutput.txt |
||
− | |||
− | * PSIBLAST |
||
− | blastpgp -i sequence.fasta -j iterations -h evalueCutoff -d database > PsiblastOutput.txt |
||
− | |||
− | * HHSearch |
||
− | hhsearch -i query -d database -o output.txt |
||
− | |||
− | |||
− | database = /data/blast/nr/nr |
||
− | |||
− | Sequences chosen for the multiple Alignment: |
||
− | {| border="1" style="text-align:center; border-spacing:0;" |
||
− | !SeqIdentifier |
||
− | !Seq Identity |
||
− | !source |
||
− | |- |
||
− | !colspan="3"| 99-90% Sequence Identity |
||
− | |- |
||
− | |gi|56967006|pdb|1X7Z||99%||PSI BLAST, 3 iterations, E-value cutoff 0.005 |
||
− | |- |
||
− | |gi|7546384|pdb|1DTW||95%||BLAST |
||
− | |- |
||
− | |gi|34810149|pdb|1OLU||99%||PSI BLAST, 3 iterations, E-value cutoff 10E-6 |
||
− | |- |
||
− | |gi|13277798|gb|AAH03787.1||95%||PSI BLAST, 3 iterations, E-value cutoff 10E-6 |
||
− | |- |
||
− | |gi|148727347|ref|NP_001092034.1||95%||BLAST |
||
− | |- |
||
− | !colspan="3"| 89-60% Sequence Identity |
||
− | |- |
||
− | |gi|196011048|ref|XP_002115388.1||66%||PSI BLAST, 3 iterations, E-value cutoff 0.005 |
||
− | |- |
||
− | |gi|149543950|ref|XP_001517857.1||67%||BLAST |
||
− | |- |
||
− | |gi|47227873|emb|CAG09036.1||82,5%||FASTA |
||
− | |- |
||
− | |gi|47196273|emb|CAF88112.1||81%||PSI BLAST, 5 iterations, E-value cutoff 0.005 |
||
− | |- |
||
− | |gi|12964598|dbj|BAB32665.1||88%||PSI BLAST, 5 iterations, E-value cutoff 10E-6 |
||
− | |- |
||
− | !colspan="3"| 59-40% Sequence Identity |
||
− | |- |
||
− | |gi|193290664|gb|ACF17640.1||47%||BLAST |
||
− | |- |
||
− | |gi|215431443|ref|ZP_03429362.1||40%||FASTA |
||
− | |- |
||
− | |gi|225557347|gb|EEH05633.1||51%||PSI BLAST, 3 iterations, E-value cutoff 10E-6 |
||
− | |- |
||
− | |gi|58267618|ref|XP_570965.1||50%||PSI BLAST, 5 iterations, E-value cutoff 0.005 |
||
− | |- |
||
− | |gi|162449842|ref|YP_001612209.1||41%||PSI BLAST, 5 iterations, E-value cutoff 10E-6 |
||
− | |- |
||
− | !colspan="3"| 39-20% Sequence Identity |
||
− | |- |
||
− | |gi|56966700|pdb|1W85||31%||PSI BLAST, 3 iterations, E-value cutoff 0.005 |
||
− | |- |
||
− | |gi|5822330|pdb|1QS0||38.1%||FASTA |
||
− | |- |
||
− | |gi|13516864|dbj|BAB40585.1||33%||PSI BLAST, 3 iterations, E-value cutoff 10E-6 |
||
− | |- |
||
− | |gi|284166853|ref|YP_003405132.1||35%||PSI BLAST, 5 iterations, E-value cutoff 0.005 |
||
− | |- |
||
− | |gi|76800932|ref|YP_325940.1||34%||PSI BLAST, 5 iterations, E-value cutoff 10E-6 |
||
− | |||
− | |} |
||
− | |||
− | Sequences for the Multiple Sequences Alignment were downloaded via NCBI, the sequence id can be changed in the link to retrieve the fasta format: |
||
− | http://www.ncbi.nlm.nih.gov/protein/76800932?report=fasta |
||
− | |||
− | === Multiple Alignments === |
||
− | * [[ClustalW]] |
||
− | clustalw sequences.fasta |
||
− | |||
− | * [[T-Coffee]] |
||
− | t_coffee -seq sequences.fasta |
||
− | |||
− | * [[T-Coffee(3D)]] |
||
− | t_coffee -seq sequences.fasta -mode expresso |
||
− | |||
− | * [[Muscle]] |
||
− | muscle -in sequences.fasta -out output.aln |
||
− | |||
− | * [[Cobalt]] |
||
− | |||
− | download [ftp://ftp.ncbi.nlm.nih.gov/pub/cobalt cobalt] |
||
− | |||
− | ./cobalt -i sequences.fasta -norps T > output.aln |
||
− | |||
− | === Conservation and Gaps === |
||
− | |||
− | {| border="1" style="text-align:center; border-spacing:0;" |
||
− | |'''Alignment methods''' |
||
− | |'''Gaps''' |
||
− | |'''Avg Gap Length''' |
||
− | |- |
||
− | |ClustalW |
||
− | |12 |
||
− | |3,75 |
||
− | |- |
||
− | |T-Coffee |
||
− | |25 |
||
− | |4,56 |
||
− | |- |
||
− | |T-Coffee (3D) |
||
− | |56 |
||
− | |4,75 |
||
− | |- |
||
− | |Cobalt |
||
− | |19 |
||
− | |3,26 |
||
− | |- |
||
− | |Muscle |
||
− | |17 |
||
− | |4,76 |
||
− | |} |
||
− | |||
− | back to [[Maple syrup urine disease]] main page |
Revision as of 00:02, 21 May 2011
Sequence
- Uniprot:
>sp|P12694|ODBA_HUMAN 2-oxoisovalerate dehydrogenase subunit alpha, mitochondrial OS=Homo sapiens GN=BCKDHA PE=1 SV=2
MAVAIAAARVWRLNRGLSQAALLLLRQPGARGLARSHPPRQQQQFSSLDDKPQFPGASAE FIDKLEFIQPNVISGIPIYRVMDRQGQIINPSEDPHLPKEKVLKLYKSMTLLNTMDRILY ESQRQGRISFYMTNYGEEGTHVGSAAALDNTDLVFGQYREAGVLMYRDYPLELFMAQCYG NISDLGKGRQMPVHYGCKERHFVTISSPLATQIPQAVGAAYAAKRANANRVVICYFGEGA ASEGDAHAGFNFAATLECPIIFFCRNNGYAISTPTSEQYRGDGIAARGPGYGIMSIRVDG NDVFAVYNATKEARRRAVAENQPFLIEAMTYRIGHHSTSDDSSAYRSVDEVNYWDKQDHP ISRLRHYLLSQGWWDEEQEKAWRKQSRRKVMEAFEQAERKPKPNPNLLFSDVYQEMPAQL RKQQESLARHLQTYGEHYPLDHFDK
Sequence info: [1] The Uniprot sequence is 445 aa long, as is contains the transit peptide sequence from position 1-45.
- PDB:
>1U5B:A|PDBID|CHAIN|SEQUENCE
SSLDDKPQFPGASAEFIDKLEFIQPNVISGIPIYRVMDRQGQIINPSEDPHLPKEKVLKLYKSMTLLNTMDRILYESQRQ GRISFYMTNYGEEGTHVGSAAALDNTDLVFGQYREAGVLMYRDYPLELFMAQCYGNISDLGKGRQMPVHYGCKERHFVTI SSPLATQIPQAVGAAYAAKRANANRVVICYFGEGAASEGDAHAGFNFAATLECPIIFFCRNNGYAISTPTSEQYRGDGIA ARGPGYGIMSIRVDGNDVFAVYNATKEARRRAVAENQPFLIEAMTYRIGHHSTSDDSSAYRSVDEVNYWDKQDHPISRLR HYLLSQGWWDEEQEKAWRKQSRRKVMEAFEQAERKPKPNPNLLFSDVYQEMPAQLRKQQESLARHLQTYGEHYPLDHFDK
Sequence info: [2]